座標の分かってる天体を、IRAFのタスクphotを使って測光しよう

一年前：http://d.hatena.ne.jp/yoshiyuki_kitne/20070530/p1

から少しだけ上達w

#! /bin/csh -f

rm -f tmp_phot NB816spec_ch1mag.data

cl <<EOF 
noao
digi
ap
phot \
hogehoge_ch1.fits \          # 測光したい画像
coords=NB816spec_X_Y.coo \   # 測光したい天体の座標ファイル（col1:X(Image),col2:Y(Image)）
wcsin=tv \
verify- \
inter- \
apertur=5.0 \                # 測光するアパーチャの半径[pixs]
zmag=33.001 \                # 等級原点
calgorithm="none" \          # 中心位置を求めるアルゴリズム
output="tmp_phot"
txdump tmp_phot MAG yes > NB816spec_ch1mag.data
logout
EOF

このスクリプトにより、NB816spec_ch1mag.dataに等級がずらっと出力されます。