PyFITS

PyFITSのメモ
fitsファイルをオープン

fitsdata = pyfits.open(files)

ヘッダーの読みだし。fitsのサイズであるNAXISの値を引っ張る。

X1 = fitsdata[0].header['NAXIS1']
Y1 = fitsdata[0].header['NAXIS2']

データをリストとして読み込み。その表示。

scidata = fitsdata[0].data
linedata=scidata[Y1/2-2:Y1/2+2,0:X1]

ファイルクローズ

fitsdata.close()