[AstroPy] convert equinox issue with nicmos/radio data?

[Thomas Robitaille]

Hi Grigoris,

Could you send me the problematic FITS file and a short script to
reproduce the problem? (off-list) Once I have the files I can look
into it.

Thanks,

Tom (APLpy developer)


On 25 February 2014 15:34, Grigoris Maravelias <gr.maravelias at gmail.com> wrote:
> Hello to all!
>
> A colleague of mine experiences a problem with astropy when she tries to
> plot a nicmos image combined with radio contours.
> Although the scripts runs flawlessly with most objects, when it fails
> the error output is this:
>
> ERROR: Exception: Cannot convert between equatorial coordinates for
> equinoxes j2000 and j2000 [aplpy.wcs_util]
> Traceback (most recent call last):
>    File "./nicmos.py", line 98, in <module>
>      f.show_contour(a149, colors="black",
> levels=np.logspace(math.log10(sa[hrs]),math.log10(maxa[hrs]),
> num=10),linestyles=['--'],slices=[0,0],linewidths=1)
>    File "<string>", line 2, in show_contour
>    File "/Library/Python/2.7/site-packages/aplpy/decorators.py", line
> 25, in _auto_refresh
>      f(*args, **kwargs)
>    File "/Library/Python/2.7/site-packages/aplpy/aplpy.py", line 916, in
> show_contour
>      contour_util.transform(c, wcs_contour, self._wcs, filled=filled,
> overlap=overlap)
>    File "/Library/Python/2.7/site-packages/aplpy/contour_util.py", line
> 26, in transform
>      output=(system_out, equinox_out))
>    File "/Library/Python/2.7/site-packages/aplpy/wcs_util.py", line 213,
> in convert_coords
>      raise Exception("Cannot convert between equatorial coordinates for
> equinoxes %s and %s" % (equinox_in, equinox_out))
> Exception: Cannot convert between equatorial coordinates for equinoxes
> j2000 and j2000
>
> Is there any idea why it fails ?
>
> Best
> Grigoris
>
>
> _______________________________________________
> AstroPy mailing list
> AstroPy at scipy.org
> http://mail.scipy.org/mailman/listinfo/astropy