整約メモ

###  作成: 2008 08/05

 高速カメラのCCDチップが1.99°回転しているため、スリット像およびスペクトルも
 1.99°傾いている
  (補足: 天体のスペクトル画像は大気分散効果で1.99°とは異なる傾きになる)
 
 データ整約では、2種類の方法がある
    ・flat画像やスペクトル画像全てに対して傾き補正をして整約。
    ・傾き補正をせずに整約。

 以下、 傾き補正をして整約を行う手順について説明:
  ・1次処理のための準備:
      ・ダークフレームの用意
      ・ピクセル感度ムラ補正フレームの用意
      ・スリット長方向グローバルパターン補正フレームの用意
  ・1次処理
      ・ダーク引き
      ・ピクセル感度ムラ補正
      ・傾き補正(回転: rotate)
           rotate @inlist @outlist -1.99 interp=sinc
           (補間はsinc関数を推奨)
      ・スリット長方向グローバルパターン補正
  ・スペクトル一次元化
  ・波長較正


具体的な手順:
  ・ピクセル感度ムラ補正フレーム作成:
      ・boxcarでsmoothingをかけたフレームを作成し、オリジナルフレームを割る
          ことで作成。
          smoothingのbox sizeは 使用slitによって変える方が良い:
 
             slit    box-size
              なし     3x3
            0.11mm     9x3
             0.2mm     9x5

          <= 画像回転によるカウントの変化を少しでも抑えるために、
             slit使用時には、box-sizeを大きくしている。
             詳細な調査結果はこちら。

         ・例: slitなし観測での感度ムラ補正フレーム作成:
             boxcar flat.d.fits flat.bc3x3.fits 3 3 
             imarith flat.d.fits / flat.bc3x3.fits flat.px.fits

  ・スリット長方向グローバルパターン補正フレーム作成:
       ・boxcarでsmoothingをかけたフレームを元に作成。
       ・傾き補正(rotate、補間関数はsincを使用):
            rotate flat.bc3x3.fits flat.bc3x3.cw199s.fits -1.99 interp=sinc
       ・波長分散方向をヘッダーに明記
            hedit flat.bc3x3.cw199s.fits DISPAXIS 2 add+
       ・スリット長方向で平均した値でカウントを規格化(response)。
            twodspec
            longslit
            response flat.bc3x3.cw199s.fits flat.bc3x3.cw199s.fits[40:470,*] resp.bc3x3.cw199s.fits functio=spline3 order=14 low_rej=3 high_re=3 niterate=2


  ・スペクトル1次元化
       ・apallを使用(twodspec.apextract)
       ・オプションは以下の通り:
          ・プリズム分光
             apall HR5501_noslit_2prism_250ms_0001.cw199s.n.fits output=hr5501_1d_01.cw199s format=multispec extras+ line=INDEF nsum=5 lower=-12 upper=12 b_order=2 b_sampl="-140:-30,30:140" b_naver=-150 b_niter=1 b_low_r=3 b_high_=3 width=20 nfind=1 llimit=-15 ulimit=15 ylevel=INDEF t_nsum=5 t_step=5 t_order=4 t_sampl="171:271" t_niter=1 backgro=fit

          ・グリズム分光
             apall HR5501_noslit_grism_8s_0001.cw199s.n.fits output=hr5501_1d_01.cw199s format=multispec extras+ line=INDEF nsum=5 lower=-12 upper=12 b_order=2 b_sampl="-90:-30,30:90" b_naver=-150 b_niter=1 b_low_r=3 b_high_=3 width=20 nfind=1 llimit=-15 ulimit=15 ylevel=INDEF t_nsum=5 t_step=5 t_order=4 t_niter=1 backgro=fit

---