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

Grigoris Maravelias gr.maravelias at gmail.com
Tue Feb 25 09:34:20 EST 2014


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





More information about the AstroPy mailing list