Hi, <br>I have a problem with pywcs. I tried to get a feeling for the library by playing around with it. <br>At the end of the email, there is my modified &quot;test&quot;, which is taken from here &quot;<a href="http://stsdas.stsci.edu/astrolib/pywcs-1.0a1-4.3_api/">http://stsdas.stsci.edu/astrolib/pywcs-1.0a1-4.3_api/</a>&quot;<br>
<br>Now when I let this test run, I get the following output:<br>&gt;-------------------USING TEST HEADER-------------------<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;[[ 99.322&nbsp;&nbsp; 36.3123]]<br>
<br>To me this seems rather strange, since in the testheader, there is this part here:<br>CRPIX1&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 151.0 / Reference pixel on axis 1<br>
CRPIX2&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 151.0 / Reference pixel on axis 2<br>
CRVAL1&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 279.322 / Value at ref. pixel on axis 1<br>
CRVAL2&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -7.3123 / Value at ref. pixel on axis 2<br><br>and my test uses [150,150] as test pixel position.<br>Should I not get something like:~ [[279, -7]] ?<br><br>I hope someone can help me here.<br><br>Thanks, <br>
Stefan<br><br><br><br>#!/usr/bin/env python<br>import numpy<br>import pywcs<br>import pyfits<br>import sys<br><br>print &quot;-------------------USING TEST HEADER-------------------&quot;<br>h = &quot;&quot;&quot;<br>SIMPLE&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; T / file does conform to FITS standard<br>
BITPIX&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -32 / number of bits per data pixel<br>NAXIS&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 / number of data axes<br>NAXIS1&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 301 / length of data axis 1<br>NAXIS2&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 301 / length of data axis 2<br>
EXTEND&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; T / FITS dataset may contain extensions<br>COMMENT&nbsp;&nbsp; FITS (Flexible Image Transport System) format is defined in &#39;Astronomy<br>COMMENT&nbsp;&nbsp; and Astrophysics&#39;, volume 376, page 359; bibcode: 2001A&amp;A...376..359H<br>
CTYPE1&nbsp; = &#39;RA---CAR&#39;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; / Type of co-ordinate on axis 1<br>CTYPE2&nbsp; = &#39;DEC--CAR&#39;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; / Type of co-ordinate on axis 2<br>CRPIX1&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 151.0 / Reference pixel on axis 1<br>CRPIX2&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 151.0 / Reference pixel on axis 2<br>
CRVAL1&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 279.322 / Value at ref. pixel on axis 1<br>CRVAL2&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -7.3123 / Value at ref. pixel on axis 2<br>CDELT1&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -0.02 / Pixel size on axis 1<br>CDELT2&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.02 / Pixel size on axis 2<br>
MJD-OBS =&nbsp;&nbsp;&nbsp;&nbsp; 51544.4992571308 / Modified Julian Date of observation<br>DATE-OBS= &#39;2000-01-01T11:58:55.816&#39;/ Date of observation<br>RADESYS = &#39;FK5&nbsp;&nbsp;&nbsp;&nbsp; &#39;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; / Reference frame for RA/DEC values<br>EQUINOX =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2000.0 / [yr] Epoch of reference equinox<br>
OBJECT&nbsp; = &#39;SOMETHING&#39;<br>BUNIT&nbsp;&nbsp; = &#39;Count&nbsp;&nbsp; &#39;<br>BSCALE&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.0<br>BZERO&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0<br>&quot;&quot;&quot;<br><br>wcs = pywcs.WCS(h)<br><br># Print out the &quot;name&quot; of the WCS, as defined in the FITS header<br>
print <a href="http://wcs.name">wcs.name</a><br><br># Some interesting pixel coordinates<br>pixcrd = numpy.array([[150,150]], numpy.double)<br><br># Convert pixel coordinates to world coordinates<br>world = wcs.pixel2world(pixcrd)<br>
print world<br><br>