Reading binary table images
Philip Massey wrote on Jul 06, 2018
I have some images produced by the UKIRT WFCAM pipeline. An example can be found at:
http://www.lowell.edu/users/massey/w20171208_00206_sf_st.fit
I can load this image into ds9 standalone, but no matter what I've tried, I can't access any of it using IRAF. If I do a stsdas "catfits" on it, I see:
EXT# FITSNAME FILENAME EXTVE DIMENS BITPI OBJECT
0 w20171228_002 8 M31
1 BINTABLE 1Fx4135R
2 BINTABLE 1Fx4135R
3 BINTABLE 1Fx4135R
4 BINTABLE 1Fx4135R
I can understand why there are four extensions (there are four chips involved) but nothing I would normally do in these circumstances (such as incopy image.fit[1] image1.fits
works. I'd like to be able to implot it, imexamine it, run apphot, etc.. Any suggestions on how to transform this into something that IRAF can read?
I can do an imhead on it; see below.
Any help would be appreciated!
Best regards,
Phil
vocl> imhead *206*sf_st.fit l+
w20171208_00206_sf_st.fit[0][short]: M31
No bad pixels, min=0., max=0. (old)
Line storage mode, physdim [0], length of user area 8627 s.u.
Created Thu 15:28:11 05-Jul-2018, Last modified Thu 15:28:11 05-Jul-2018
Pixel file "w20171208_00206_sf_st.fit" [NO PIXEL FILE]
EXTEND = T / FITS dataset may contain extensions
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H
DATE = '2017-12-08T17:46:54' / file creation date (YYYY-MM-DDThh:mm:ss UT)
ORIGIN = 'Starlink Software' / Origin of this FITS file
TELESCOP= 'UKIRT ' / Telescope name
INSTRUME= 'WFCAM ' / Instrument
DHSVER = 'UKDHS 2002 Oct 31 ' / Data handling version
HDTFILE = 'wfcam.hdt ' / Name of global hdt file
OBSERVER= ' ' / Observers names
USERID = ' ' / Userid logged in as
OBSREF = ' ' / PATT or other reference
PROJECT = 'U/17B/UA03' / Time-allocation code
SURVEY = ' ' / Survey name
SURVEY_I= ' ' / Pointing ID within survey
MSBID = '3442a7fc55865fff66ad6401aa700342' / Id min.-schedulable block
MSBTID = 'UKIRT_1512707887_911449' / MSB Transaction ID
MSBTITLE= 'M31 ' / MSB Title
QUEUE = 'UA ' / Queue the MSB belongs to
RQ_MINSB= ' ' / [mag/arcsec**2] requested min J sky brightness
RQ_MAXSB= ' ' / [mag/arcsec**2] requested max J sky brightness
RQ_MNSEE= '0.00 ' / [arcsec] requested min seeing
RQ_MXSEE= '1.0 ' / [arcsec] requested max seeing
RQ_MINCL= '0 ' / [percentage] requested min cloud coverage
RQ_MAXCL= '20 ' / [percentage] requested max cloud coverage
RQ_MNTAU= '0 ' / requested min tau constraint
RQ_MXTAU= ' ' / requested max tau constraint
RQ_MINMN= '0 ' / [percentage] req. min illumination moon
RQ_MAXMN= '100 ' / [percentage] req. max illumination moon
RMTAGENT= ' ' / Name of remote agent
AGENTID = ' ' / Unique identifier for remote agent
OBJECT = 'M31 ' / Object name from telescope
RECIPE = 'JITTER ' / Data reduction recipe to be used
OBSTYPE = 'OBJECT ' / BIAS|DARK|SKYFLAT|DOMEFLAT|OBJECT|SKY|FOCUS
OBSNUM = 206 / Observation number
GRPNUM = 206 / Group number applied to all members
GRPMEM = T / Group membership
TILENUM = 0 / Tile number applied to all members
STANDARD= F / Is the target a standard star observation?
NJITTER = 5 / Number of positions in tel jitter pattern
JITTER_I= 1 / Serial number in this tel jitter pattern
JITTER_X= 0.00 / [arcsec] X (RA) offset in tel jitter pattern
JITTER_Y= 0.00 / [arcsec] Y (Dec) offset in tel jitter pattern
NUSTEP = 4 / Number of positions in microstep pattern
USTEP_I = 1 / Serial number in this microstep pattern
USTEP_X = 0.00 / [arcsec] X (RA) offset in microstep pattern
USTEP_Y = 0.00 / [arcsec] Y (Dec) offset in microstep pattern
NFOC = 0 / Number of positions in focus scan
NFOCSCAN= 0 / Number of focus scans in focus test
UTDATE = '20171208' / UT date as integer in yyyymmdd format
DATE-OBS= '2017-12-08T04:38:19.510878' / Date time (UTC) of start of observation
DATE-END= '2017-12-08T04:38:30.819151' / Date time (UTC) of end of observation
GPS-OBS = '2017-12-08T04:38:17' / Date time (UTC) of start of observation from G
GPS-END = '2017-12-08T04:38:28' / Date time (UTC) of end of observation from G
MJD-OBS = 58095.19327500 / DATE-OBS as Modified Julian Date
WCSAXES = 2 / Number of axes in world co-ordinate system
RADESYS = 'FK5 ' / Mean IAU 1984 equatorial co-ordinates
EQUINOX = 2000.000 / [yr] Equinox of object position
TRACKSYS= 'J2000 ' / Telescope tracking coordinate system
RABASE = 0.6603389 / [h] Right ascension of base position
DECBASE = 40.7929639 / [deg] Declination of base position
TELRA = 0.6603389 / [h] Current telescope right ascension
TELDEC = 40.7929639 / [deg] Current telescope declination
GSRA = 0.6591111 / [h] Right ascension of guide star
GSDEC = 40.7711083 / [deg] Declination of guide star
LST_OBS = 0.975576278 / [days] / LST at start of observation [days]
LST_END = 0.975709745 / [days] / LST at end of observation [days]
HABASE = -18.9280939 / [deg] Hour angle
TRAOFF = 0.000 / [arcsec] Telescope offset in RA
TDECOFF = 0.000 / [arcsec] Telescope offset in Dec
AMSTART = 1.118 / Airmass at start of observation
AMEND = 1.118 / Airmass at end of observation
RARATE = 0.000000000E+00 / Non-sidereal tracking rate in RA [rad]
DECRATE = 0.000000000E+00 / Non-sidereal tracking rate in DEC [rad]
APER_X = 0.000 / [arcsec] Aperture X axis coordinate
APER_Y = 0.000 / [arcsec] Aperture Y axis coordinate
ISU2PORT= 'Lost ' / ISU2 port name
ISU2POFF= 0 / ISU2 port offset
ISU2TOFF= 0 / ISU2 tilt offset
GUIDING = 'ON ' / Guider status: OFF|ON|PAUSED
GUIDFFOC= ' ' / Fast guider CCD fine focus
AGFREQ = 16.000 / Autoguider: Frequency
FGMODE = 'Fast Guider' / Fast guider mode
IRPOLARM= ' ' / Position of IRPOL arm
NINT = 1 / Number of integrations in observation
READMODE= 'CDS ' / Name of camera readmode
EXP_TIME= 10.000000 / Integration time per exposure
NEXP = 1 / Number of exposures in integration
READINT = 0.000000 / Interval between reads
NREADS = 2 / Number of reads per exposure
AIRTEMP = 3.100 / [degC] Air temperature
BARPRESS= 615.214 / [mbar] Ambient pressure
DEWPOINT= -18.100 / [degC] Dewpoint
DOMETEMP= 3.800 / [degC] Dome temperature
HUMIDITY= 16.300 / Relative Humidity
MIRR_NE = 896.475 / [degC] Mirror temperature NE
MIRR_NW = 3.374 / [degC] Mirror temperature NW
MIRR_SE = 883.019 / [degC] Mirror temperature SE
MIRR_SW = 841.256 / [degC] Mirror temperature SW
MIRRBTNW= 1.740 / [degC] Mirror bottom temp. SW
MIRRTPNW= 1.628 / [degC] Mirror top temp. SW
SECONDAR= 7.371 / [degC] Temperature of secondary
TOPAIRNW= 5.239 / [degC] Top air NW
TRUSSENE= 6.737 / [degC] Truss leg ENE
TRUSSWSW= 5.032 / [degC] Truss leg WSW
TRUSSTOP= 7.264 / [degC] Temperature at the top of the truss
TRUSSBOT= 8.010 / [degC] Temperature at the bottom of the truss
WIND_DIR= 343.679 / [deg] Wind direction, azimuth
WIND_SPD= 0.000 / [km/h] Wind speed
CSOTAU = 0.109 / Tau at 225 GHz from CSO
TAUDATE = '2017-12-08T04:21:00' / Time and date of Tau reading
TAUSRC = 'CSO225GHZ' / Source of opacity data
IZ5 = 987.000 / Astigmatism: IZ5
HZ5SH = -763.000 / Astigmatism: HZ5SH
HZ5CH = 9.000 / Astigmatism: HZ5CH
HZ5SD = 1492.000 / Astigmatism: HZ5SD
HZ5CD = -362.000 / Astigmatism: HZ5CD
IZ6 = 44.000 / Astigmatism: IZ6
HZ6SH = 2273.000 / Astigmatism: HZ6SH
HZ6CH = 140.000 / Astigmatism: HZ6CH
HZ6SD = -805.000 / Astigmatism: HZ6SD
HZ6CD = -160.000 / Astigmatism: HZ6CD
Z7 = -405.000 / Coma: Z7
Z8 = 1450.000 / Coma: Z8
WEBEAM = -405.000 / Topend: W-E beam
NSBEAM = 45.000 / Topend: N-S beam
IZ9 = -630.000 / Trefoil: IZ9
HZ9SH = 60.000 / Trefoil: HZ9SH
HZ9CH = -30.000 / Trefoil: HZ9CH
HZ9SD = 124.000 / Trefoil: HZ9SD
HZ9CD = 206.000 / Trefoil: HZ9CD
IZ10 = 404.000 / Trefoil: IZ10
HZ10SH = -155.000 / Trefoil: HZ10SH
HZ10CH = -100.000 / Trefoil: HZ10CH
HZ10SD = 7.000 / Trefoil: HZ10SD
HZ10CD = 36.000 / Trefoil: HZ10CD
Z11 = 250.000 / Spherical: Z11
AAMP = 2295.437 / Astigmatism: amplitude
APHI = -15.988 / Astigmatism: angle
TAMP = 590.139 / Trefoil: Amplitude
TPHI = 46.121 / Trefoil: angle
CASSFOC = ' ' / Tel_inst focus
M2_X = -0.968 / [mm] M2 X (E-W) raw position
M2_XD = -0.968 / [mm] M2 X (E-W) demand position
M2_Y = 0.493 / [mm] M2 Y (N-S) raw position
M2_YD = 0.493 / [mm] M2 Y (N-S) demand position
M2_Z = 0.698 / [mm] M2 Z (focus) raw position
M2_ZD = 0.698 / [mm] M2 Z (focus) demand position
M2_U = 5.613 / [mrad] M2 U (N-S tilt) raw position
M2_UD = 5.613 / [mrad] M2 U (N-S tilt) demand position
M2_V = 3.550 / [mrad] M2 V (E-W tilt) raw position
M2_VD = 3.545 / [mrad] M2 V (E-W tilt) demand position
M2_W = 0.000 / [mrad] M2 W (axial rotation) raw position
TCS_FOC = 0.360 / telescope (TCS) focus
TEL_FOFF= -3.400 / focus offset for inst. filter
FOC_TOFF= 0.000 / [mm] Offset from nominal telescope focus
FOC_POSN= 0.060 / [mm] Internal focus position
FOC_ZERO= 3.500 / [mm] Focus zero-point position
FOC_OFFS= 0.000 / [mm] Focus offset
FOC_FOFF= 0.060 / [mm] Focus filter offset
FILTER = 'J ' / Filter name
DETTEMPB= 74.7 / [K] DET_TEMP alternate sensor
DETTEMPS= 75.0 / [K] Det temperature setpoint
DET_HTP = 75.0 / [percent] Det heater level
CCD_TEMP= 160.1 / [K] AG CCD temperature
CCDTEMPS= 160.0 / [K] AG CCD temperature setpoint
CCD_HTP = 14.0 / [percent] AG CCD heater level
FS_TEMP = 72.0 / [K] CCC 1st stage temperature
LNC_TEMP= 96.5 / [K] LN can temperature
SS_TEMP = 16.9 / [K] CCC 2nd stage temperature
M3_TEMP = 113.5 / [K] M3 temperature
TR_TEMP = 121.2 / [K] Top ring temperature
RS_TEMP = 136.0 / [K] Rad shield temperature
CP_TEMP = 192.1 / [K] Cor plate temperature
SC_TEMP = 1.4 / [K] SCSU cab temperature
CC_PRES = 3.00000e-06 / [mbar] Cryostat Pressure
DATE = '2017-12-08T17:46:54' / file creation date (YYYY-MM-DDThh:mm:ss UT)
ORIGIN = 'Starlink Software' / Origin of this FITS file
UTSTART = 4.638611 / [h] Start time of integration
UTEND = 4.641667 / [h] End time of integration
FOC_I = 0 / Serial number in focus scan
FOC_OFF = 0.0 / [mm] Offset from nominal inst. focus pos
HISTORY 20171209 15:49:41
HISTORY $Id: cir_compress_fits.c,v 1.14 2010/10/01 09:26:57 jim Exp $
CIR_VERS= '$Revision: 1.202 $' / CIRDR Version
SOFTVERS= '$Revision: 1.202 $' / CIRDR Version
SOFTDATE= '$Date: 2014/08/08 11:42:28 $' / CIRDR release date
SOFTAUTH= 'Jim Lewis ' / Contact for bug reports
SOFTINST= 'CASU: http://casu.ast.cam.ac.uk' / CASU URL
PREVIEWV= '$Id: preview,v 1.77 2014/05/06 15:22:57 jim Exp $' / Version of previ
WSA_TIME= ' ' / WSA time stamp
WSA_MFID= 0 / WSA Multiframe ID
Last post on Jul 06, 2018