FOCAS画像からプレイメージを作る・座標から天体を選ぶ
2005年8月1日 とうこく
---------------------------------------------------------------------
FOCAS合成画像を作る
---------------------------------------------------------------------
idlたちあげ
>idlde &
下のコマンドラインから合成ソフト立ち上げ
>bigimage
「overscan subtraction」「distortion correction」を選択。
フラット画像があるときは「flat fielding」も選択。
歪みを補正した合成画像ができる。
暗い天体のときはフラットして画像を重ね合わせたほうがよい。
明るい領域なら不要。
---------------------------------------------------------------------
ディザーした画像を重ね合わせる
---------------------------------------------------------------------
手順
1.fits=ref
2.fits=input
sex 天体検出
xyxymatch in ref 1.fitsに対して2.fitsの対応天体を求める
geomap in ref 1.fitsに対して2.fitsのずれ量を求める
geotran 2.fits 2.fitsずれを直す
imcombine 重ね合わせる
>sex FCSA_Q1623_1.fits -c default.sex
>mv sex.cat FCSA_Q1623_1.cat
>awk '{if ($5>1.00 && $5<8.0) print $1, $2}' FCSA_Q1623_1.cat > FCSA_Q1623_1a.cat
>sex FCSA_Q1623_2.fits -c default.sex
>mv sex.cat FCSA_Q1623_2.cat
>awk '{if ($5>1.00 && $5<8.0) print $1, $2}' FCSA_Q1623_2.cat > FCSA_Q1623_2a.cat
ずれ量をだいたい求めておく(xin・yinはref原点に対してin原点のずれ量)
>xyxymatch FCSA_Q1623_2a.cat FCSA_Q1623_1a.cat xyxy.list toler=5.0 xin=50 yin=-1
>geomap xyxy.list xyxygmp.dbs 1 4029 1 4095
>geotran FCSA_Q1623_2.fits FCSA_Q1623_2tran.fits xyxygmp.dbs xyxy.list
重ね合わせる画像をリスト(imcomb_in.lst)に入れる
FCSA_Q1623_1.fits
FCSA_Q1623_2tran.fits
imcombineのためにヘッダ修正(WCSDIM=2)
>hedit FCSA_Q1623_1.fits WCSDIM 2 add+
重ね合わせ
>imcomb @imcomb_in.lst FCSA_Q1623_12.fits combine="average"
---------------------------------------------------------------------
座標から天体を選びたいとき
---------------------------------------------------------------------
WCS座標付けをするため、FOCAS画像とカタログ天体の対応リストを作る。
--- Example ------------------------------------------------------------------
画像:
http://archive.stsci.edu/cgi-bin/dss_form
POSS2/UKSTUのデータがよい (Palmar Obs. Sky Survey)
座標:
http://irsa.ipac.caltech.edu/
2MASSデータ
持ってきた座標データから、RA、DEC、等級などを抜き出して、「*.gaia」という
ファイルにする。columnはタブで区切る。タイトルとリストの間は「---」で区切る。
ID RA DEC J H K
--- --- --- --- --- ---
.
.
.
skycatで画像を表示してから、「data-servers → Local catalog → Load from file」
で *.gaia を読み込む。「Option → Set Plot symboles」でプロット種類を設定すれば
カタログの位置をプロットできる。
------------------------------------------------------------------------------
FOCASの画像でさちっていない天体とカタログにある天体を
並べてリストファイル *.coo にする(ここではFCSA_Q1623.coo)。
X[pix] Y[pix] RA DEC
.
.
.
パラメータファイルを作る
>ccmap FCSA_Q1623.coo CSA_Q1623.param subpix以下でフィットさせる
WCSを張りつける(fitsヘッダが書き変わるので注意)
>ccsetwcs FCSA_Q1623_12f.fits CSA_Q1623.param FCSA_Q1623.coo
skycatでWCSを付けたFOCAS画像を表示させ、観測したい天体のリストを
先の*.gaiaと同じように作って表示させれば完成。