ERROR: MWCS on TKRS data
jkennefick wrote on Nov 13, 2009
Hello,
When trying to SPLOT some 1d TKRS spectra, I get the following error:
ERROR: MWCS: dimension mismatch (mw_open)
Frank Valdes previously mentioned that you could try a wcsreset to wcs = world, but that gives the same error.
The complete header is:
spec1d.KTRS16.070.10574.fits[0][short]:
No bad pixels, min=0., max=0. (old)
Line storage mode, physdim [0], length of user area 41 s.u.
Created Fri 13:44:01 13-Nov-2009, Last modified Fri 13:44:01 13-Nov-2009
Pixel file "spec1d.KTRS16.070.10574.fits" [NO PIXEL FILE]
EXTEND = T /Extensions may (will!) be present
and that's it, obviously not the whole header.
Any suggestions?
Thanks!
When trying to SPLOT some 1d TKRS spectra, I get the following error:
ERROR: MWCS: dimension mismatch (mw_open)
Frank Valdes previously mentioned that you could try a wcsreset to wcs = world, but that gives the same error.
The complete header is:
spec1d.KTRS16.070.10574.fits[0][short]:
No bad pixels, min=0., max=0. (old)
Line storage mode, physdim [0], length of user area 41 s.u.
Created Fri 13:44:01 13-Nov-2009, Last modified Fri 13:44:01 13-Nov-2009
Pixel file "spec1d.KTRS16.070.10574.fits" [NO PIXEL FILE]
EXTEND = T /Extensions may (will!) be present
and that's it, obviously not the whole header.
Any suggestions?
Thanks!
Mike Fitzpatrick wrote on Nov 13, 2009
What you sent is the primary header of a multi-extension FITS file, presumaby the data and the WCS information is in one of the extension headers. Do you specify an extension number to splot/imheader/wcsreset or just the filename (or is it the '[0]' primary header)? If you still have problems, could you upload the image to ftp://iraf.noao.edu/pub so we can have a look?
-Mike
-Mike
Francisco Valdes wrote on Nov 13, 2009
Hi Julia,
You sent me a separate email with a header listing of the extension. I am going to reply through the forum.
The information you sent indicates to me that the spectrum is in a table. Currently the IRAF spectral formats supported by tasks such as SPLOT consist only of images. So you would have to convert the data into a spectral image or use tools intended for tables.
I found the following link with some format information
http://tkserver.keck.hawaii.edu/tksurvey/data_products/type.html
So you can use IGI or SGRAPH to plot things from the table format. You can dump the wavelength and flux data with tdump (or tprint) and import it into and IRAF spectrum image with RSPECTEXT.
Good luck,
Frank Valdes
P.S. The error message does sound like things which have been noted earlier in the forums having to do with problems with the spectral WCS. However, in this case the error is really due to the data not being an image.
You sent me a separate email with a header listing of the extension. I am going to reply through the forum.
The information you sent indicates to me that the spectrum is in a table. Currently the IRAF spectral formats supported by tasks such as SPLOT consist only of images. So you would have to convert the data into a spectral image or use tools intended for tables.
I found the following link with some format information
http://tkserver.keck.hawaii.edu/tksurvey/data_products/type.html
So you can use IGI or SGRAPH to plot things from the table format. You can dump the wavelength and flux data with tdump (or tprint) and import it into and IRAF spectrum image with RSPECTEXT.
Good luck,
Frank Valdes
P.S. The error message does sound like things which have been noted earlier in the forums having to do with problems with the spectral WCS. However, in this case the error is really due to the data not being an image.
Hello Dr. Valdes,
When trying to SPLOT some 1d TKRS spectra, I get the following error:
ERROR: MWCS: dimension mismatch (mw_open)
You previously mentioned on iraf.net that you could try a wcsreset to wcs = world, but that gives the same error.
The complete header is:
spec1d.KTRS16.070.10574.fits[0][short]:
No bad pixels, min=0., max=0. (old)
Line storage mode, physdim [0], length of user area 41 s.u.
Created Fri 13:44:01 13-Nov-2009, Last modified Fri 13:44:01 13-Nov-2009
Pixel file "spec1d.KTRS16.070.10574.fits" [NO PIXEL FILE]
EXTEND = T /Extensions may (will!) be present
and that's it, obviously not the whole header.
Any suggestions?
---
XTENSION= 'BINTABLE' /Binary table written by MWRFITS v1.2
BITPIX = 8 /Required value
NAXIS = 2 /Required value
NAXIS1 = 106520 /Number of bytes per row
NAXIS2 = 1 /Number of rows
PCOUNT = 0 /Normally 0 (no varying arrays)
GCOUNT = 1 /Required value
TFIELDS = 15 /Number of columns in table
COMMENT
COMMENT *** End of mandatory fields ***
COMMENT
COMMENT
COMMENT *** Column names ***
COMMENT
TTYPE1 = 'SPEC ' /
TTYPE2 = 'LAMBDA ' /
TTYPE3 = 'IVAR ' /
TTYPE4 = 'CRMASK ' /
TTYPE5 = 'BITMASK ' /
TTYPE6 = 'ORMASK ' /
TTYPE7 = 'NBADPIX ' /
TTYPE8 = 'INFOMASK ' /
TTYPE9 = 'OBJPOS ' /
TTYPE10 = 'FWHM ' /
TTYPE11 = 'NSIGMA ' /
TTYPE12 = 'R1 ' /
TTYPE13 = 'R2 ' /
TTYPE14 = 'SKYSPEC ' /
TTYPE15 = 'IVARFUDGE ' /
COMMENT
COMMENT *** Column formats ***
COMMENT
TFORM1 = '4096E ' /
TFORM2 = '4096E ' /
TFORM3 = '4096E ' /
TFORM4 = '4096I ' /
TFORM5 = '4096I ' /
TFORM6 = '4096I ' /
TFORM7 = '4096I ' /
TFORM8 = '4096I ' /
TFORM9 = 'E ' /
TFORM10 = 'E ' /
TFORM11 = 'E ' /
TFORM12 = 'J ' /
TFORM13 = 'J ' /
TFORM14 = '4096E ' /
TFORM15 = 'E ' /
COMMENT
COMMENT *** End of mandatory fields ***
COMMENT
EXTNAME = 'Bxspf-R ' /Extension name
DATE = '2003-09-05' /Creation date
OBJNO = '10574 ' /Object number
OBJPOS = 27.0040 /Object Position from Spatial Profile
FWHM = 9.35104 /FWHM from Spatial Profile
CAT_OBJP= 24.1813 /Design Object Position
CAT_FWHM= 0.00000 /FWHM from PCAT
COR_FWHM= 9.35104 /FWHM from PCAT corrected for seeing diff
EXT_FWHM= 9.52629 /FWHM employed in object extraction
RA_OBJ = '12:37:21.25' /RA of object
DEC_OBJ = '+62:12:47.0' /DEC of object
XMM = 50.4740 /X position of slit
YMM = 84.2920 /Y position of slit
MAGB = 999.000 /B magnitude
MAGR = 999.000 /R magnitude
MAGI = 999.000 /I magnitude
FRAMENO = 80 /
OUTFILE = 'd0528_ ' /
EXPTIME = 3600.14941406 /Total exposure time
DARKTIME= 1200.54882812 /
OBSERVER= 'faber, davis, conselice, wirth' /
OBJECT = 'KTRS16 spectral' /
OBSTYPE = 'Object ' /
ROTATVAL= 14.3313827500 /
DATE-OBS= '2003-05-28' /
UT = 0 /
AIRMASS = 1.36917934000 /
TARGNAME= 'KTRS15 ' /
EQUINOX = 2000.00000000 /
DEC = '+62:11:46.3' /
RA = '12:37:19.13' /
AZ = -8.60179084000 /
EL = 46.8671886300 /
HA = '+00:50:54.82' /
ST = '13:28:24.23' /
UTC = '07:28:43.42' /
MJD-OBS = '52787.311518' /
PARANG = 162.392763840 /
SYNOPSIS= '1200 s, Object Spectral KTRS16/600ZD/GG455 Off' /
DWFILNAM= 'GG455 ' /
SLMSKNAM= 'KTRS16 ' /
GRATEPOS= 4 /
HPLOGTIM= 'Tue May 27 21:46:35 HST 2003' /
SP2DVERS= 'v1_1_4 ' /Version of spec2d
AUTHOR = 'Finkbeiner & Davis' /
CHIPNO = 3 /chip number for this slitlet
SLITNO = 70 /slit number
SLITX0 = 489.747 /center of trace
SLITX1 = 552.451 /center of trace
DX = 0.0328321 /offset from calibslit trace
WAVESIG = 0.0105905000000 /RMS error in lambda solution, in AA (1&10 are fl
SHIFT0 = -0.616827 /shift to fix low-x slit end
SHIFT1 = 0.230574 /shift to fix high-x slit end
SKYSIGMA= 2.18104 /Gaussian sigma of skylines (AA)
COMMENT
COMMENT *** Column names ***
COMMENT
COMMENT
COMMENT *** Column formats ***
COMMENT
COMMENT
COMMENT *** Column dimensions (2 D or greater) ***
COMMENT
COMMENT FLUX array in units of counts/hour
COMMENT spslit frames combined, sky subtracted
COMMENT CR removal by means of chisqr analysis
OBJPA = 999.000 /Object PA on sky
SLITPA = 74.0000 /Slit PA on sky
MASKPA = 44.0000 /Mask PA on sky
WAVETYPE= 'POLYFLAG' /wavelength solution method
CTYPE1 = 'LAMBDA ' / Coordinate Type
CTYPE2 = 'LAMBDA ' / Coordinate Type
CD1_1 = 0.631918549538 / Degrees / Pixel
CD2_1 = 0.00000000000 / Degrees / Pixel
CD1_2 = 0.00000000000 / Degrees / Pixel
CD2_2 = 0.00000000000 / Degrees / Pixel
CRPIX1 = 0.00000 / Reference Pixel in X
CRPIX2 = 0.00000 / Reference Pixel in Y
CRVAL1 = 4627.66268329 / R.A. (Degrees)
CRVAL2 = 0.00000000000 / Dec (Degrees)
HISTORY PUTAST: Sep 5 02:26:00 2003 CD, CRPIX and CRVAL parameters written
END
jkennefick wrote on Nov 13, 2009
That seems right. We have been able to see them with IGI. We'll give that a try. Thanks a bunch.
Julia
Julia
Last post on Nov 13, 2009