[AstroPy] pywcs problem on Linux
OKUMURA, Akira
oxon@ceres.phys.s.u-tokyo.ac...
Wed May 6 01:30:23 CDT 2009
Hello all,
I have installed pywcs in two environments.
(1) Scientific Linux 4.7 (64 bit)
+ Python 2.5.4 (r254:67916, May 6 2009, 07:09:17)
+ GCC 3.4.6 20060404 (Red Hat 3.4.6-10)] on linux2
+ pyfits-1.3
+ pywcs-1.3a1-4.3
(2) Red Hat Enterprise Linux 4 (32 bit).
+ Python 2.5.1 (r251:54863, Aug 24 2007, 16:22:50)
+ GCC 3.4.3 (TWW)] on linux2
+ pywcs-1.3a1-4.3
+ pyfits-2.1.1
But I got following errors even though it runs properly on Mac OS X
10.5.
=== Error on (1) ===
>>> hdulist = pyfits.open("cmap.fits")
>>> wcs = pywcs.WCS(hdulist[0].header)
zsh: segmentation fault python
=== Error on (2) ===
>>> hdulist = pyfits.open("cmap.fits")
>>> wcs = pywcs.WCS(hdulist[0].header)
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "/u/gl/oxon/tmp/lib/python2.5/site-packages/pywcs/pywcs.py",
line 155, in __init__
self.footprint = self._calcFootprint(header)
File "/u/gl/oxon/tmp/lib/python2.5/site-packages/pywcs/pywcs.py",
line 181, in _calcFootprint
return self.all_pix2sky(corners)
File "/u/gl/oxon/tmp/lib/python2.5/site-packages/pywcs/pywcs.py",
line 307, in all_pix2sky
return self._array_converter(self._all_pix2sky, *args, **kwargs)
File "/u/gl/oxon/tmp/lib/python2.5/site-packages/pywcs/pywcs.py",
line 294, in _array_converter
return func(args[0], **kwargs)
SystemError: Python/getargs.c:1246: bad argument to internal function
Is there any tips on running pywcs with Linux ?
The FITS is as follows.
>>> print hdulist[0].header
SIMPLE = T / File conforms to NOST standard
BITPIX = 32 / Bits per pixel
NAXIS = 2 / No data is associated with this header
NAXIS1 = 160 / Length of data axis 1
NAXIS2 = 160 / Length of data axis 2
EXTEND = T / Extensions may be present
COMMENT FITS (Flexible Image Transport System) format is defined in
'Astronomy
COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...
376..359H
CTYPE1 = 'GLON-AIT' / RA---%%%, %%% represents the
projection method
CRPIX1 = 80.5 / Reference pixel
CRVAL1 = 211. / RA at the reference pixel
CDELT1 = -0.25 / X-axis incr per pixel of physical
coord at posi
CUNIT1 = 'deg ' / Physical unit of X-axis
CTYPE2 = 'GLAT-AIT' / DEC---%%%, %%% represents the
projection method
CRPIX2 = 80.5 / Reference pixel
CRVAL2 = -17. / DEC at the reference pixel
CDELT2 = 0.25 / Y-axis incr per pixel of physical
coord at posi
CUNIT2 = 'deg ' / Physical unit of Y-axis
CROTA2 = 0. / Image rotation (deg)
DATE = '2009-04-01T15:34:22' / file creation date (YYYY-MM-
DDThh:mm:ss U
FILENAME= 'cmap.fits' /
TELESCOP= 'GLAST ' / name of telescope generating data
INSTRUME= 'LAT ' / name of instrument generating data
DATE-OBS= '2009-03-11T21:33:18.000' / start date and time of the
observation (U
DATE-END= '2009-03-17T16:26:38.000' / end date and time of the
observation (UTC
NDSKEYS = 5
EQUINOX = 2000. / Equinox of RA & DEC specifications
OBSERVER= 'Peter Michelson' / GLAST/LAT PI
CREATOR = 'gtbin ' / Software and version creating file
HISTORY LatCountMapTemplate,v 1.3 2005/04/05
21:06:39 peachey
HISTORY Exp
CHECKSUM= 'c8JZf8JXc8JXc8JX' / HDU checksum updated
2009-04-01T06:34:22
DATASUM = '182648 ' / data unit checksum updated
2009-04-01T06:34:22
DSTYP1 = 'TIME '
DSUNI1 = 's '
DSVAL1 = 'TABLE '
DSREF1 = ':GTI '
DSTYP2 = 'POS(RA,DEC)'
DSUNI2 = 'deg '
DSVAL2 = 'CIRCLE(86.81257852,-6.005536809,20)'
DSTYP3 = 'CTBCLASSLEVEL'
DSUNI3 = 'dimensionless'
DSVAL3 = '3: '
DSTYP4 = 'ENERGY '
DSUNI4 = 'MeV '
DSVAL4 = '200:20000'
DSTYP5 = 'ZENITH_ANGLE'
DSUNI5 = 'deg '
DSVAL5 = '0:100 '
Regards,
OKUMURA, Akira oxon@ceres.phys.s.u-tokyo.ac.jp
Department of Physics, The University of Tokyo
7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033
TEL/FAX +81 3-5841-4173/4059
Skype : okumura.akira
More information about the AstroPy
mailing list