FITS形式のイメージを読み込む

image.fits というイメージファイルを読ませてみましょう。

  IDL> im = readfits('image.fits')
                :

im という変数は、宣言をしなくても2次元配列になります。 もし、FITS のヘッダも

  IDL> im = readfits('W49B_Fe_image256_deco30_Lucy.fits',hdr)

とすれば、hdrは文字列の変数になります。

  IDL> print,hdr
          :

とすると、hdr の内容をみることができます。