[AstroPy] pywcs problem on Linux

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


Hello again,

ver 1.4-4.3 runs properly with the sample FITS and the modified script  
on my Mac (10.5). So the problem seems to happen only with Linux.


On 2009/05/06, at 23:07, OKUMURA, Akira wrote:

> 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
>>
>
> _______________________________________________
> AstroPy mailing list
> AstroPy@scipy.org
> http://mail.scipy.org/mailman/listinfo/astropy




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