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

[Thomas Robitaille]

Hi Grigoris,

(just replying to the list in case anyone has the same issue)

The issue was due to the fact that one image had the first axis be the
declination and the second the right ascension (as opposed to RA then
Dec). I have fixed it so that this is now properly taken into account
in APLpy. You can update with

git clone https://github.com/aplpy/aplpy.git
cd aplpy
python setup.py install

Let me know if this does not fix the issue. I did see another error,
but I think this is unrelated.

Cheers,
Tom


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