[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