[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