Reading IMAPS Data Using the IDL
Astronomy User's Library Routines

The following examples show how to use three different FITS readers (available from the IDL Astronomy Users library) to read IMAPS images.

1) READFITS and TBGET

written by J. Woffard & Wayne Landsman from the UIT project:

Read the corrected image and wavelengths from a coadded, corrected file:

image = readfits('imaps1_198-199.fits',h)      ; read image
tab = readfits('imaps1_198-199.fits',e,/exten) ; read binary table
w0 = tbget(e,tab,'wavelength',0)               ; starting wavelengths
dw = tbget(e,tab,'deltaw',0)                   ; wavelength increment
m = tbget(e,tab,'order',0)                     ; order numbers

2) FXREAD and FXBREAD

written by Bill Thompson of the SOHO project:

Read a coadded corrected image and wavelength vectors:

fxread,'imaps1_306-309.fits',image,h           ; read image
fxbopen,un,'imaps1_306-309.fits',1,h           ; open table 
fxbread,un,w0,'wavelength',1                   ; read specified fields
fxbread,un,dw,'deltaw',1
fxbread,un,m,'order',1

3) MRDFITS

written by Tom McGlynn from HEASARC:

Extract the coadded image and wavelength information:

image=mrdfits('imaps1_337-339.fits',0,h) ; read image
struc=mrdfits('imaps1_337-339.fits',1,e) ; read binary table fields
help,st,/struc                           ; display table entries
m = st.order       ; order numbers
w0 = st.wavelength ; starting wavelengths
dw = st.deltaw     ; wavelength increments