radio-astro-tools / pvextractor

Position-Velocity diagram extract (but it can slice out other things too)
BSD 3-Clause "New" or "Revised" License
16 stars 13 forks source link

Header choke causes a ds9 pvextractor failure. #20

Closed low-sky closed 10 years ago

low-sky commented 10 years ago

The following header barfs with a

ERROR: ValueError: No valid coordinate format found. [pvextractor.pvregions]

though the file ingests happily into SpectralCube.

SIMPLE  =                    T
BITPIX  =                  -32
NAXIS   =                    3
NAXIS1  =                  329
NAXIS2  =                  260
NAXIS3  =                  680
CTYPE1  = 'RA      '
CRVAL1  =  1.8000000000000E+02
CRPIX1  = -9.2495000000100E+03
CDELT1  = -1.6666700607100E-02
CTYPE2  = 'DEC     '
CRVAL2  =  0.0000000000000E+00
CRPIX2  = -1.6755000000100E+03
CDELT2  =  1.6666700607100E-02
CTYPE3  = 'VELO-LSR'
CRVAL3  = -4.0008280000300E+05
CRPIX3  =  1.0000000000000E+00
CDELT3  =  7.3612501160100E+02
BUNIT   = 'K (Tb)  '
TELESCOP= 'Arecibo 305m'
DATE-OBS= '2000-01-01T12:00:00.0'
RESTFREQ=  1.4204058647200E+09
CELLSCAL= 'CONSTANT'
DATAMAX =  4.2133846282959E+01
DATAMIN = -6.0426634550095E-01
OBSERVER= 'Turn On Galfa Su'
OBJECT  = 'GALFA-HI RA+DEC Tile 22.7500 +30'
EPOCH   =  2.0000000000000E+03
ORIGIN  = 'Miriad Fits: version 1.1 09-Feb-07'
HISTORY kvis: Module version 1.0.0  Karma v1.7.18  compiled with v1.7.17
astrofrog commented 10 years ago

@low-sky - are you still seeing this issue? pvextractor should be using SpectralCube to read in the data, so it should work now.

low-sky commented 10 years ago

Looks resolved to me. Nice work!