[AstroPy] Bug in wcs_world2pix?

Maik Riechert maik.riechert at arcor.de
Fri Sep 19 01:53:54 EDT 2014


Hi Dominik,

Nope, same results.

Cheers
Maik

Am 19.09.2014 00:50, schrieb Dominik Klaes:
> Hi Maik,
> 
> try instead of this
> 
> x,y = wcs.wcs_world2pix(211, -26, 0)
> 
> this
> 
> x,y = wcs.wcs_world2pix([[211, -26]], 0), also in the transformation back. Does this help? I remember that I had a similar problem.
> 
> Cheers,
> 
> Dominik
> 
> On 18. September 2014 23:04:17 MESZ, Maik Riechert <maik.riechert at arcor.de> wrote:
>> Hi,
> 
>> I'm having some great trouble currently and am a bit confused too:
> 
>> wcs = WCS(header)
>> x,y = wcs.wcs_world2pix(211, -26, 0)
>> print x,y # 836.316942718 658.26364248
>> ra,dec = wcs.wcs_pix2world(x, y, 0)
>> print ra, dec # 31.0 26.0
> 
>> It seems like world2pix gives me bogus results. I was expecting values
>> way outside my image area (image size is 4256x2832). Am I doing
>> something obvious wrong?
> 
>> The WCS headers can be seen at http://pastebin.com/JrdiL949
> 
>> Thanks for the help,
>> Maik
>> _______________________________________________
>> AstroPy mailing list
>> AstroPy at scipy.org
>> http://mail.scipy.org/mailman/listinfo/astropy
> 




More information about the AstroPy mailing list