iraf.net

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\$this->_split2($m[0]) 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

These images are FITS compressed which DS9 natively understands, but IRAF doesn't. You'll need to use the 'FUNPACK' tool from CFITSIO to unpack the image, or the FUNPACK task if you have the FITSUTIL package installed.

Mike, thanks! FUNPACK, huh! I'll try it!