これまでボスのサンプルプログラムを参考にして、34天体のfitsファイルの平均とメジアンをとるプログラムを作ったわけですが、同じことをIRAFでやるにはimcombineというコマンドを用いればよいようです。
目的：複数のfitsファイル（画像リストinput.list）の平均を取る
im > imcombine @input.list output.fits combine="average"
目的：複数のfitsファイル（画像リストinput.list）のメジアンを取る
im > imcombine @input.list output.fits combine="median"
