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()