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

[Grigoris Maravelias]

Hi Tom! Indeed this fixed the issue.

Thanks for the quick response.

Best
Grigoris

On 02/26/2014 09:24 PM, Thomas Robitaille wrote:
> 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