spacer link to MAST page spacer logo image spacer
 
link to STScI page


Planned outage -- MAST Services will be unavailable Tuesday 16 January 2018, 09:00 to 17:00 EST (14:00 to 22:00 GMT) for system upgrades. We apologize for the inconvenience.

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