NIC画像 upper limitの測定(point source)

Last-modified: 2014-07-07 (月) 22:50:00 (3580d)

Top / NIC画像 upper limitの測定(point source)

NIC画像 upper limitの決定 †

WCSが正確に入った画像を用意する

1.ds9 Search for Catalog をクリック

"Name or Designation" に "2MASS"と入力し、以下のカタログを選択
 II/246 | 2MASS All-Sky Catalog of Point Sources (Cutri+ 2003)
 Retriveをクリックすると、画像中の点光源のリストが表示され、画像の2MASS sourceに印がつく
 File - Save - Starbase でカタログの値を保存(保存したファイルをds9.catとする)

4. ra decのみを取り出して、"#"を消す

 cat ds9.cat | awk '($1 !~ "#" && $1 !~ "_" && $1 !~ "---"){print $1,$2}' > 2mass_radec.dat

5. 各画像のwcsと座標を対応付けして、座標ファイルを作る

sky2xy jGRB140705A_01.cmw.fits @2mass_radec.dat | awk '{print $5,$6}'> 2mass_xy.j01.dat
すべての画像にたいして行う

6. IRAF/photで測光する(ap, skyvalueはimstatの測定値を使う)

phot jGRB140705A_01.cmw.fits coord=2mass_xy.h01.dat ap=23.4 annulus=5 dannulus=5 cbox=5 readnoise=50 epadu=9.2 skyvalue=145 interac- verify- 
pdump jGRB140705A_01.cmw.fits.mag.1  id,flux,area,epadu,stdev,nsky,mag,merr yes > j01.mag

7. flu errorを計算

error = sqrt (flux / epadu + area * stdev**2 + area**2 * stdev**2 / nsky)