カタログ作成
目的:整約済みのFITSデータからカタログを作成する
狭帯域バンド
広帯域バンド
1pixel=0.117"
GAIN=3.15 e-/ADU
- 位置合わせ
すべてのマップKバンドのマップに位置を合わせる
GT2のみで行う
KバンドのマップからGT2の領域[(1,1969)(4017,4438)]だけ切り取る
imcopy K.fits[1:4017,1969:4438] K_cut.fits
KバンドとNB119バンドの共通の星の座標を得る
imexam K_cut.fits
ds9上でrを押すと輝度プロファイル、aを押すと端末上に掃き出す
qで終了
このときシーイングも測定しておく
シーイング
Kバンド:5.26[pixel]
Jバンド:4.82[pixel]
NB119:6.5
10~20個程度の星の座標のテキストを作成
Xk Yk X209 Y209
x1 y1 x2 y2
epar geomap
input:作成した天体の座標テキスト
database:出力テキスト名
xmin,xmax,ymin,ymax:基本となるマップの座標を入力
今回の場合KバンドのGT2のマップに合わせるので
1,4017,1,2470
:goで実行
変換するためのテキストが作成できた
epar geotran
input:座標変換するマップのfitsを入力
output:
database:geomapで作成したテキストを入力
transfer:星の座標のテキストを入力
- 検出・測光
位置を合わせた画像の不必要な領域をimcopyを使って切り取る
imcopy NB119_trans.fits[355:3905,254:2079] NB119_set.fits
シーイングの差をなくすためにガウシアンでスムージングを行う
epar gauss
sex [検出する画像], [測光する画像] -c [sex]
今回は全てKバンドで検出した天体を測光
測光パラメータについて
DETECT_MINAREA は25か16ピクセルで
DETECT_THRESH は2σか2.5σで
NBデータ解析(1次整約)
NB209解析手順
0. ファイルリストの準備
chip1とchip2それぞれのファイルリストを作成する
cl> listprep /(ディレクトリ) chip 120789 121976
chip_1.lst : chip1のファイルリスト
chip_2.lst : chip2のファイルリスト
ファイルリストからほしいフレームのみを取り出す
cl> hselect @chip_2.lst $I 'FILTER03==NB2095' > object.list
ファイルリストからヘッダを見る
cl> hselect @object.list FRAMEID,DATE-OBS,FILTER03,OBJECT,EXPTIME yes
1. ダーク引き
準備すべきリスト
- 積分時間の等しいダークフレームのリスト
- 積分時間の等しいオブジェクトフレームのリスト
積分時間の等しいダークフレームからメジアンをとったダークフレームを作成する
cl> imcomb @dark_900.list std_dark_900.fits combine=median
オブジェクトフレームから積分時間の等しいダークフレームを引く
cl> !sed 's/MCSA/d_MCSA/g' object_900.list > object_900_d.list
cl> imarith @object_900.list - std_dark_900.fits @object_900_d.list
2. 天体マスク・フラット作成
ヘッダの書き込み禁止になっている場合、chmodで可能にしておく
# chmod -R a+rwx /ディレクトリ
観測日ごとにオブジェクトフレームを分けてリストを作成
タスクmcs_mksflatでセルフフラットを作成
PACKAGE = mcsred
TASK = mcs_mksflat
infile = @object_0419.list The list of input images
nflat = flat_0419.fits The name of output flatframe
(bis = no) Subtract dark frame?
(fbis = ) The name of dark frame
(keepbis= yes) Keep dark-subtracted frames after the process?
(flt = no) Execute flatfielding for each frame?
(ownflat= no) Use user-supplied flat for making mask? (not recommended)
(fownfla= ) The name of the user-supplied flat for mask
(crrej = yes) Execute cosmicray cleaning? (for NB data)
(skip = no) Skip the mask making process (but add maskfile name to the BPM keyword).
(second = no) 2nd-path process?
(sthresh= 1.5) Threshold for object mask.
(msksize= 4.) Object expansion size (pix).
(rail = no) Do the rail remain on chip 1?
(cmask1 = ) Optional bad region mask for chip 1
(cmask2 = ) Optional bad region mask for chip 2
(disp = yes) display the skysub process?
(mode = q)
天体マスクがヘッダに書き込まれたfitsが生成されるので
このファイルを今後使わないと天体マスクに関する情報がないのでエラーがでる
3. フラット引き
cl> !sed 's/d_MCSA/f_MCSA/g' object_0419.list > object_0419_f.list
cl> imarith @bisobject_0419.list / flat_0419.fits @object_0419_f.list
4. 1次スカイ引き
リストにあるファイルが異なるディレクトリにあると出力ファイルを作成するときにエラーがでる
PACKAGE = mcsred
TASK = sbselfsky
infile = object_0419_f.list The list of input images.
nsf = 2 2 x nsf frames are used for making sky
(skytype= median) Combine algorithm for sky (median/average)?
(zero = yes) zero matching the sky when subtraction?
(sc = no) scale the sky when subtraction?
(sreject= avsig) Rejection type for sky (none/minmax/avsig)?
(lrej = 5000) minmax: Lower rejection number.
(hrej = 14000) minmax: Higher rejection number.
(nkeep = 3) minmax: Minimum to keep.
(lsig = 3.) avsig: Lower sigma rejection value.
(hsig = 3.) avsig: Higher sigma rejection value.
(disp = no) Display and check result by display?
(tm = 0.) sleep time during image display [sec].
(ars = 0.8) relative size for statistics region
(mode = q)
5. 2次スカイ引き
PACKAGE = mcsred
TASK = qmsepskysb
infile = @object_0419_f.list An image or @list format
(ng = 5) parameter ngrow in imsurfit
(xo = 3) parameter xorder in imsurfit
(yo = 3) parameter yorder in imsurfit
(upper = 3) Upper limit for sigma rejection in imsurfit
(lower = 3) Lower limit for sigma rejection in imsurfit
(niter = 3) Maximum number of rejection cycles in imsurfit
(ff = spline3) fitting function in imsurfit
(ao = yes) Replace AO frange region to average value for sky?
(rail = no) Do the rail remain on chip 1?
(cmask1 = ) Optional chip-1 common bad region mask.
(cmask2 = ) Optional chip-2 common bad region mask.
(disp = yes) Display the anomaly subtraction process?
(hkeep = no) Keep object mask in the header? (usually no!)
(mode = q)
6. 溝埋め
7. 歪み補正
PACKAGE = mcsred
TASK = mcsgeocorr
infile = @qcor_object.list Input image name (allow @list format)
(config = dir_mcsred$DATABASE/ana_feb08.cfg) the name of mcsred config file
(crrej = yes) Execute cosmicray cleaning? (time consuming...)
(mode = q)
ここまでが1ショット画像
画像合成
1.各ショットのPSFを測定
2.PSFをそろえる
gauss MCSA00120924.fits g_0924.fits 1.6
3.画像端などの使えない領域は極端に小さな値を入れておく
imreplace g_0924.fits[600:2048,1:2048] 0.001
4.gsextcat,gmkgtrimagesを用いて位置合わせを行う
出来上がった仮の画像を用いてMODSのデータに位置合わせを行う
5.3で作成したオブジェクトフレームとマスクを両方geotranを用いてMODSに合わせる
最終更新:2009年10月29日 20:45