[AstroPy] pywcs problem on Linux

OKUMURA, Akira oxon@ceres.phys.s.u-tokyo.ac...
Wed May 6 09:07:27 CDT 2009


Hello Mike and Thomas,

I modified test.py to use the second or third argument, but still I  
get a seg fault after calling wcs_pix2sky. And wcs.wcs.fix() always  
cause a seg fault. I have uploaded my sample FITS and modified test.py.

http://www-utheal.phys.s.u-tokyo.ac.jp/~oxon/model_total_ICHII_bin00.fits
http://www-utheal.phys.s.u-tokyo.ac.jp/~oxon/test.py

Could you investigate it?

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

On 2009/05/06, at 22:30, Michael Droettboom wrote:

> It looks as if your new traceback is related to an API change -- all  
> of the transformation functions now require an "offset" argument,  
> which should be 0 for C/Numpy indexing or 1 for Fortran/FITS indexing.
>
> If that change doesn't work for you, can you send me a FITS file and  
> script to reproduce this off list?  I'll see if I can reproduce what  
> you're seeing here.
>
> Mike
>
>
>
> OKUMURA, Akira wrote:
>> Hello Nadia,
>>
>> Thank you for your reply. I didn't know that new version has been   
>> released recently. I installed pywcs-1.4-4.3.tar.gz and tried  
>> test.py  in the package. But I got errors again.
>>
>> 1)
>> $ ./test.py .
>> = = =  
>> = 
>> = 
>> = 
>> =====================================================================
>> ./cmap.fits
>> zsh: segmentation fault  ./test.py .
>>
>> the line wcs.wcs.fix() causes the seg fault.
>>
>> 2) comment out wcs.wcs.fix()
>> $ ./test.py .
>> = = =  
>> = 
>> = 
>> = 
>> =====================================================================
>> ./cmap.fits
>> Traceback (most recent call last):
>>   File "./test.py", line 82, in <module>
>>     run_directory(directory)
>>   File "./test.py", line 74, in run_directory
>>     test_file(filepath)
>>   File "./test.py", line 38, in test_file
>>     world = wcs.wcs_pix2sky(data3)
>>   File "/usr/local/lib/python2.5/site-packages/pywcs/pywcs.py",  
>> line  419, in wcs_pix2sky
>>     *args, **kwargs)
>>   File "/usr/local/lib/python2.5/site-packages/pywcs/pywcs.py",  
>> line  378, in _array_converter
>>     raise TypeError("Expected 2 or 3 arguments, %d given" %  
>> len(args))
>> TypeError: Expected 2 or 3 arguments, 1 given
>>
>> I tried several FITS files but I got the same errors. Do you have  
>> any  idea?
>>
>> 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
>>
>> On 2009/05/06, at 20:47, Nadezhda Dencheva wrote:
>>
>>
>>> Hi,
>>>
>>> You are running an older version of pywcs.
>>> I believe an upgrade to the latest release will fix this.
>>>
>>>
>>> Nadia Dencheva
>>>
>>> ---- Original message ----
>>>
>>>> Date: Wed, 6 May 2009 15:30:23 +0900
>>>> From: astropy-bounces@scipy.org (on behalf of "OKUMURA, Akira" <oxon@ceres.phys.s.u-tokyo.ac.jp
>>>>> )
>>>>>
>>>> Subject: [AstroPy] pywcs problem on Linux
>>>> To: astropy@scipy.org
>>>>
>>>> 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
>>>>
>>>> _______________________________________________
>>>> AstroPy mailing list
>>>> AstroPy@scipy.org
>>>> http://mail.scipy.org/mailman/listinfo/astropy
>>>>
>>
>> _______________________________________________
>> AstroPy mailing list
>> AstroPy@scipy.org
>> http://mail.scipy.org/mailman/listinfo/astropy
>>
>
> -- 
> Michael Droettboom
> Science Software Branch
> Operations and Engineering Division
> Space Telescope Science Institute
> Operated by AURA for NASA
>



More information about the AstroPy mailing list