整約メモLast modified: Wed 6 Aug 2008
### 作成: 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