How do I read a RHESSI image cube FITS file?
To read an image (single or multi) image FITS file into the image object, use the im_input_fits parameter as follows: o = hsi_image() o -> set, im_input_fits = ‘hsi_imagecube_19tx2e_20020527_180350.fits’ The control and info parameters, and images will be set in the object. You can retrieve images or info for single images, or for the entire array. You may not change any control parameters until you set im_input_fits back to a blank string, or set the im_allow_reprocess parameter. See How to work with the Image Object for more details. To read the FITS file without creating an image object, you can use hsi_image_fitsread , hsi_fits2map, or mrdfits. If you simply want to retrieve the array of images from the FITS file, the mrdfits routine suffices: images = mrdfits(filename,0,header) However, hsi_image_fitsread and hsi_fits2map are tailored to the RHESSI single image and image cube files. Look at the header documentation for all of the options. hsi_image_fitsread examples: images = hsi_
