手順3: bias引き

biasフレームを作ります。 biasフレームリストを作って、imcombineで足し上げます。

#biasフレームのリスト作成。
ls kls0004096[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 kls00040979-cut.fits > obs-cut.list
ls kls0004098[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