NIC解析マニュアルのバックアップ(No.6)

(まだまだ編集中)

概要 †

上図は各段階での画像の例。(a) 生画像、 (b) 4までを行った画像、 (c) 5で生成したスカイフレーム (d) 6のスカイ差し引き処理後の画像、 (e) 縦縞パターン差し引き (7) 後の画像、 (f) 最終画像

ダークフレームの取得手順に従って取得したダークフレームを用いる。
オブジェクトフレームと積分時間を揃えたダークフレームを10-15枚程度使い、重ね合わせる。
camera マシン (制御室) にインストールされている自動スクリプトでは、3sigma clipping をした後で平均 (average) をしている。
ダークフレーム取得の際に、コールドシャッターの位置が適切でないと、光が入る。Pol 窓や Open 窓が入ってしまったフレームは、前もって削除しておく。
重ね合わせたフレームを、オブジェクトフレーム (含・フラットフレーム) から引き算する。
- 自動スクリプトで処理すると、*.da.fits という名前のファイルができる。

検出器上の位置による感度のむらを補正する。
通常、トワイライトフラットを用いる。フラットフレームの取得手順を参照。
マスターフラットを用いる場合には、camera の /home/nhao/nic/nicred.latest/common 以下にあるもののうち、日付の近いものを選んで用いる。
フラットの作成に関しては特別なツールを作ってはいない。上記以外は通常の手順をふめば良い。 - 手順としては以下の通り。
1. カウントの高いフレームからダークを引く。カウントがほぼ同じフレームが複数あれば重ね合わせる。
2. カウントの低いフレームにも同様の処理をする。
3. (高カウントのフレーム) - (低カウントのフレーム) の引き算をする。
4. 全体の平均が1になるよう、規格化する。
几帳面にやりたい場合は、(3)-(6) の処理 (ホットピクセル・ダークピクセル補正)までをやってもいい。ただしどのみち後の段階でこれらはやることになるので、ここでやる必要は必ずしもない。
規格化したフレームで、ダーク処理後のフレームを割り算する。
- 自動スクリプトで処理すると、*.fl.fits という名前のファイルができる。

ここでいう「スカイ」は、いわゆるスカイというよりも、ダーク、フラットなどで取り切れない「うねうね」のパターンを差し引くといった方がが正しい。NIC のデータに対する対処療法的な操作。
5. スカイ画像作成は、ディザー観測で得られた画像を重ね合わせることで行う。位置合わせはせずに median で重ね合わせる。(IRAF/imcombine 等)
- 星の少ない天域であれば、オブジェクトフレームをそのまま用いることができる。
- 星の混んだ天域を観測している場合には、なるべく近傍の星の少ない天域をスカイ用に別途撮像する。積分時間・モードは同じにする。
- 視野内には数個程度の星があっても OK。
6. スカイ差し引きは、基本的には単純にスカイ画像を差し引けばよい。(IRAF/imarith 等)
- ただし、差し引く場合にはスカイ画像を各天体画像のスカイレベルを合わせるために、スケーリングしておく。

NIC のデータには検出器由来の縦縞パターンが存在する。これも現状では対処療法的に引く必要がある。
- ダーク、フラットなどで取り切れない。読み出し時の電圧の変化が原因と考えられている。
- 前述の「スカイ」と違い、列ごとにほぼ一定の値を示す。
- 下記の2つの方法がある。

星の少ない天域の場合 (Sampling 法)
- 星の写っていない行 (典型的に 20-100行) を取り出して各列を縦方向に平均し、縦縞パターンのプロファイルを得る。
- 同様の操作を上半分、下半分に対してそれぞれ行う。
星の多い(混んだ)天域の場合 (Median 法)
- 画像の上下半分で、各行のカウントのメディアン (中央値) を取り、縦縞パターンのプロファイルを得る。

基準となる星の中心位置を測る。(IRAF/imexamine、SExtractor 等)
その基準星がディザリングセット内ですべて同一のピクセルに来るように、画像を平行移動する。
- NIC はさほど視野が広くないため、ゆがみ補正などはひとまず考えなくてよい。
- 当然、ローテーターの位置角を一定にしておけば、回転なども考えなくてよい。
平行移動によって位置合わせが出来たら、重ね合わせる。(IRAF/imcombine 等)