2009.12/11
今フォーマットは2048*2048になっているが、標準星画像(flFSMCSA38???.fits)は見れば分かるようにもっと小さい。なのでその部分だけ切り出してからmedian sky引きをする。標準星画像は枚数少ない上にchip1とchip2のどちらかにしか星が写ってないわけだからflat-fielding,sky引きだけで一次処理終わりにしてしまう。
その後、次週本格的に行うzero点決めの練習をする。
☆其の一 flat fieldingした標準光源(flFSMCSA???.fits @flFS〜リスト)のmedian sky引き
◎標準星をデータ別に分類してリスト化
$ more standardodd
MCSA00038343 "FS27 Standard" KS 8.000 1
MCSA00038345 "FS27 Standard" KS 8.000 1
MCSA00038347 "FS27 Standard" KS 8.000 1
MCSA00038349 "FS27 Standard" KS 8.000 1
MCSA00038351 "FS27 Standard" KS 8.000 1
MCSA00038353 "FS27 Standard" KS 8.000 1
MCSA00038355 "FS27 Standard" KS 8.000 1
MCSA00038357 "FS103 Standard" KS 8.000 1
MCSA00038359 "FS103 Standard" KS 2.000 1
MCSA00038361 "FS103 Standard" KS 3.000 1
MCSA00038363 "FS103 Standard" KS 3.000 1
MCSA00038365 "FS103 Standard" KS 3.000 1
MCSA00038367 "FS103 Standard" KS 3.000 1
MCSA00038369 "FS103 Standard" KS 3.000 1
MCSA00038371 "FS103 Standard" KS 3.000 1
MCSA00038373 "FS103 Standard" KS 3.000 1
MCSA00038375 "FS103 Standard" KS 3.000 1
MCSA00038377 "FS103 Standard" KS 3.000 1
MCSA00038595 "FS103 Standard" KS 2.000 1
MCSA00038597 "FS103 Standard" KS 3.000 1
MCSA00038599 "FS103 Standard" KS 3.000 1
MCSA00038601 "FS103 Standard" KS 3.000 1
MCSA00038603 "FS103 Standard" KS 3.000 1
MCSA00038605 "FS103 Standard" KS 3.000 1
MCSA00038607 "FS103 Standard" KS 3.000 1
MCSA00038609 "FS103 Standard" KS 2.000 1
MCSA00038611 "FS103 Standard" KS 3.000 1
MCSA00038613 "FS103 Standard" KS 3.000 1
MCSA00038615 "FS103 Standard" KS 3.000 1
$ more standardeven
MCSA00038344 "FS27 Standard" KS 8.000 1
MCSA00038346 "FS27 Standard" KS 8.000 1
MCSA00038348 "FS27 Standard" KS 8.000 1
MCSA00038350 "FS27 Standard" KS 8.000 1
MCSA00038352 "FS27 Standard" KS 8.000 1
MCSA00038354 "FS27 Standard" KS 8.000 1
MCSA00038356 "FS27 Standard" KS 8.000 1
MCSA00038358 "FS103 Standard" KS 8.000 1
MCSA00038360 "FS103 Standard" KS 2.000 1
MCSA00038362 "FS103 Standard" KS 3.000 1
MCSA00038364 "FS103 Standard" KS 3.000 1
MCSA00038366 "FS103 Standard" KS 3.000 1
MCSA00038368 "FS103 Standard" KS 3.000 1
MCSA00038370 "FS103 Standard" KS 3.000 1
MCSA00038372 "FS103 Standard" KS 3.000 1
MCSA00038374 "FS103 Standard" KS 3.000 1
MCSA00038376 "FS103 Standard" KS 3.000 1
MCSA00038378 "FS103 Standard" KS 3.000 1
MCSA00038596 "FS103 Standard" KS 2.000 1
MCSA00038598 "FS103 Standard" KS 3.000 1
MCSA00038600 "FS103 Standard" KS 3.000 1
MCSA00038602 "FS103 Standard" KS 3.000 1
MCSA00038604 "FS103 Standard" KS 3.000 1
MCSA00038606 "FS103 Standard" KS 3.000 1
MCSA00038608 "FS103 Standard" KS 3.000 1
MCSA00038610 "FS103 Standard" KS 2.000 1
MCSA00038612 "FS103 Standard" KS 3.000 1
MCSA00038614 "FS103 Standard" KS 3.000 1
MCSA00038616 "FS103 Standard" KS 3.000 1
などで確認すれば分かるように標準星データは分類でき
FS27・・・38343〜38356
FS103・・・38361〜38378
FS103・・・38597〜38616
の３種類のデータある。ここでFS103についてexp timeが2秒や8秒のやつは試しどりなので省いた方が安全かも。
・３種類の標準光源リスト*2(odd or even)をつくる
$ ls flFSMCSA383[45][13579].fits > flFS27odd
$ vi FS27odd
（いらない行(38357,38359など)を消し(dd)、保存(ZZ)）
$ more flFS27odd
flFSMCSA38343.fits
flFSMCSA38345.fits
flFSMCSA38347.fits
flFSMCSA38349.fits
flFSMCSA38351.fits
flFSMCSA38353.fits
flFSMCSA38355.fits
$ ls flFSMCSA000383[45][02468].fits > flFS27even
$ vi flFS27even
（いらない行を消し、保存）
$ more flFS27even
flFSMCSA38344.fits
flFSMCSA38346.fits
flFSMCSA38348.fits
flFSMCSA38350.fits
flFSMCSA38352.fits
flFSMCSA38354.fits
flFSMCSA38356.fits
(同様の作業を繰り返し、flFS103_1odd,flFS103_1even,flFS103_2odd,flFS103_2evenというリストもつくる)
最終的な中身は
$ more flFS103_1odd
flFSMCSA38361.fits
flFSMCSA38363.fits
flFSMCSA38365.fits
flFSMCSA38367.fits
flFSMCSA38369.fits
flFSMCSA38371.fits
flFSMCSA38373.fits
flFSMCSA38375.fits
flFSMCSA38377.fits
$ more flFS103_1even
flFSMCSA38362.fits
flFSMCSA38364.fits
flFSMCSA38366.fits
flFSMCSA38368.fits
flFSMCSA38370.fits
flFSMCSA38372.fits
flFSMCSA38374.fits
flFSMCSA38376.fits
flFSMCSA38378.fits
$ more flFS103_2odd
flFSMCSA38597.fits
flFSMCSA38599.fits
flFSMCSA38601.fits
flFSMCSA38603.fits
flFSMCSA38605.fits
flFSMCSA38607.fits
flFSMCSA38609.fits
flFSMCSA38611.fits
flFSMCSA38613.fits
flFSMCSA38615.fits
$ more flFS103_2even
flFSMCSA38598.fits
flFSMCSA38600.fits
flFSMCSA38602.fits
flFSMCSA38604.fits
flFSMCSA38606.fits
flFSMCSA38608.fits
flFSMCSA38610.fits
flFSMCSA38612.fits
flFSMCSA38614.fits
flFSMCSA38616.fits
てな感じ
◎必要なフォーマットだけ切り出し
・標準星のどの部分が必要かdisplayで確認する
cl> !ds9 &
cl> reset stdimage=imt4096
cl> display flFSMCSA38343.fits 1
cl> display flFSMCSA38344.fits 2
cl> display flFSMCSA38355.fits 3
cl> display flFSMCSA38356.fits 4
→FS27(38343~38356)は513〜1536のフォーマットになってた
2048*2048の半分、1024*1024
(以下同様な感じで)
FS103(38361~38378)は769~1280のフォーマット（512*512）
FS103(38597~38616)は769~1280のフォーマット（512*512）
・flMCSA38???.fitsの内上記のフォーマットの部分だけimcopyしたいがそのための入力リストと出力リストを用意する
入力用リストの中身はflFSMCSA38???.fits[513:1536,513:1536](FS27)又はflFSMCSA38???.fits[769:1280,769:1280](FS103)
$ awk '{print $1"[513:1536,513:1536]"}' flFS27odd > flFS27oddin
$ awk '{print $1"[513:1536,513:1536]"}' flFS27even > flFS27evenin
$ awk '{print $1"[769:1280,769:1280]"}' flFS103_1odd >flFS103_1oddin
$ awk '{print $1"[769:1280,769:1280]"}' flFS103_1even >flFS103_1evenin
$ awk '{print $1"[769:1280,769:1280]"}' flFS103_2odd >flFS103_2oddin
$ awk '{print $1"[769:1280,769:1280]"}' flFS103_2even >flFS103_2evenin
出力用リストの中身はoutflFS38???.fitsとする
$ awk '{print "outflFS"substr($1,9,20)}' flFS27odd > flFS27oddout
$ awk '{print "outflFS"substr($1,9,20)}' flFS27even > flFS27evenout
$ awk '{print "outflFS"substr($1,9,20)}' flFS103_1odd > flFS103_1oddout
$ awk '{print "outflFS"substr($1,9,20)}' flFS103_1even > flFS103_1evenout
$ awk '{print "outflFS"substr($1,9,20)}' flFS103_2odd > flFS103_2oddout
$ awk '{print "outflFS"substr($1,9,20)}' flFS103_2even > flFS103_2evenout
・2048*2048の標準星ファイルを必要なフォーマットだけ切り出してimcopy
cl> imcopy @flFS27oddin @flFS27oddout
フォーマット小さくなってるか確認
cl> imhe @flFS27oddout
outflFS38343.fits[1024,1024][real]: FS27 Standard
outflFS38345.fits[1024,1024][real]: FS27 Standard
outflFS38347.fits[1024,1024][real]: FS27 Standard
outflFS38349.fits[1024,1024][real]: FS27 Standard
outflFS38351.fits[1024,1024][real]: FS27 Standard
outflFS38353.fits[1024,1024][real]: FS27 Standard
outflFS38355.fits[1024,1024][real]: FS27 Standard
→ok
あとは同様にしてflFS27evenout,flFS103_1oddout,flFS103_1evenout,flFS103_2oddout,flFS103_2evenoutの中身もつくる
cl> imcopy @flFS27evenin @flFS27evenout
cl> imcopy @flFS103_1oddin @flFS103_1oddout
cl> imcopy @flFS103_1evenin @flFS103_1evenout
cl> imcopy @flFS103_2oddin @flFS103_2oddout
cl> imcopy @flFS103_2evenin @flFS103_2evenout
◎切り出し行ったファイル(flFS〜outというリストに入ってる)に対してスカイ引き
・sky引きするスクリプト6つ、FS27odd.cl,FS27even.cl,FS103_1odd.cl,FS103_1even.cl,FS103_2odd.cl,FS103_2even.clをつくる
$ emacs sky1013odd.cl &
適当なclスクリプトを下のように書き換えてFS27odd.clという名前で保存
→FS27odd.cl
実行
cl> cl < FS27odd.cl
できたファイルをdisplayさせて確認
cl> display subsky38355.fits 1 →こんなの
一応よさそう？(画面左よりの２つの明るい星の左側が標準星)
(以上の作業を６つの切り出しリストflFS〜outにたいして行う。)
→FS27even.cl
cl> cl< FS27even.cl
→FS103_1odd.cl
cl> cl< FS103_1odd.cl
→FS103_1even.cl
cl> cl< FS103_1even.cl
→FS103_2odd.cl
cl> cl< FS103_2odd.cl
→FS103_2even.cl
cl> cl< FS103_2even.cl
これで一応標準星画像の一次処理できたはず。
詳しくはまた次週見る
☆其の二 SExtracter使って標準星の測光→zero点決めの練習
・まずはsextractor用のパラメータファイルつくり
$ emacs default.sex & #これを以下のように書き換えて1211.sexという名で保存#
→1211.sex
・さらに何を出力させるかの指定ファイルもdefault.paramを書き換えて1211.paramをつくる
$ emacs default.param & #これを書き換えて1211.paramとして別名保存#
(defaultでコメントアウト、#かかってる行を必要に応じて#外す)
→1211.param
・適当な一枚に対して1211.sexかけてやる
$ sex subsky38362.fits -c 1211.sex
・結果見てみる
$ more test.cat
# 1 NUMBER Running object number
# 2 FLUX_ISO Isophotal flux [count]
# 3 FLUXERR_ISO RMS error for isophotal flux [count]
# 4 FLUX_APER Flux vector within fixed circular aperture(s) [count]
# 5 FLUX_AUTO Flux within a Kron-like elliptical aperture [count]
# 6 FLUXERR_AUTO RMS error for AUTO flux [count]
# 7 X_IMAGE Object position along x [pixel]
# 8 Y_IMAGE Object position along y [pixel]
# 9 FLAGS Extraction flags
1 1599.016 115.7251 1713.862 1788.532 181.3064 57.842 33.720 0
2 430134.7 593.1713 174183.1 425662.9 559.3682 455.226 436.713 0
3 2788.279 181.3064 1810.942 3973.151 348.9243 290.534 367.545 0
4 6813.465 236.6519 3091.89 8407.732 397.8349 363.606 257.851 2
5 1801.845 135.1378 1372.988 2615.345 323.579 32.648 259.341 0
6 3997.791 115.7251 3856.722 3914.883 167.3382 350.971 96.165 0
7 10155.86 249.1818 5024.22 10164.2 406.9121 109.740 93.270 0
8 1297.903 110.3395 1350.296 1444.563 191.1137 493.098 92.082 0
二行目が一番明るいのでたぶん標準星
・http://www.jach.hawaii.edu/UKIRT/astronomy/calib/phot_cal/fs_izyjhklm.datからこの標準星が何等級なのかをしらべ、それとauto flux=425662.9をくらべてzero点を決定してやる
調べた結果、FS103は11.724mag(K band)。subsky38362.fitsに写ってるFS103のauto flux(count)=425662.9、積分時間=3sec。これをmagnitude=-2.5log10(count/sec)+zero点の式に入れてzero点を決定。awkで何のファイルにプリントさせるかはテキトー
$ awk '{print 11.724+2.5*log(425662.9/3.0)/log(10)}' test.cat
24.6039
24.6039
24.6039
24.6039
24.6039
24.6039
24.6039
24.6039
24.6039
24.6039
24.6039
24.6039
24.6039
24.6039
24.6039
24.6039
24.6039
とりあえずsubsky38362.fitsだけ使ってchip2のzero点を24.6039と決定した。次回以降、本格的にzero点を決めていく
