NIC画像 upper limitの測定(point source) のバックアップの現在との差分(No.1)

バックアップ一覧
差分を表示
ソースを表示
バックアップを表示
NIC画像 upper limitの測定(point source) へ行く。
- 1 (2014-07-07 (月) 11:51:31)

追加された行はこの色です。
削除された行はこの色です。

[[staff/soft tips]]

* NIC画像 upper limitの決定 [#a51b0153]

- WCSが正確に入った画像を用意する
- ds9 Search for Catalog をクリック
 1. "Name or Designation" に "2MASS"と入力し、以下のカタログを選択
 II/246 | 2MASS All-Sky Catalog of Point Sources (Cutri+ 2003)
 
 2. Retriveをクリックすると、画像中の点光源のリストが表示され、画像の2MASS sourceに印がつく
 3. File - Save - Starbase でカタログの値を保存(保存したファイルをds9.catとする)
 4. ra decのみを取り出して、"#"を消す
 cat ds9.cat | awk '($1 !~ "#" && $1 !~ "_" && $1 !~ "---"){print $1,$2}' > 2mass_radec.dat
 5. 各画像の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の測定値を使う)
- 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を計算
- 7. flu errorを計算
 error = sqrt (flux / epadu + area * stdev**2 + area**2 * stdev**2 / nsky)