手順3: bias引き
biasフレームを作ります。 biasフレームリストを作って、imcombineで足し上げます。 ここで作成するbiasの代わりにMaster Biasを使用することも出来ますが、 biasの時間変動は未検証なのでデータを確認しながら使用してください。
#biasフレームのリスト作成。
ls kif0004096[5-9]-cut.fits > bias-cut.list
…
#imcombineのオプションはcombine=medianかcombine=average reject=sigclipなどが良いだろう。
#出力ファイルがあれば、事前に削除しておく (以下同様)。
rm -f bias-comb-cut.fits
imcombine @bias-cut.list bias-comb-cut.fits combine=median
天体+較正フレームリストを作って、biasフレームを引きます。
ls kif00040979-cut.fits > obs-cut.list
ls kif0004098[0-1]-cut.fits >> obs-cut.list
…
sed -e s/-cut.fits/-bias.fits/ obs-cut.list > obs-bias.list
awk '{print "rm -f " $1}' obs-bias.list | sh
imarith @obs-cut.list - bias-comb-cut.fits @obs-bias.list