Dome Flat vs. Twilight Flatそしてアノマリパターン

2005-03-15

Ichi Tanaka

ドームフラットとトワイライトフラットを比較

1 日目データ。ドームフラットを 1 日目にまとめて取得している関係上、 1 日目の明け方に取った K バンドの Twilight データしか比較できるものがない。

ドームフラットはオン点、オフ点を単純にそれぞれメジアン合成後引き算して作成。

オフ点のカウント 6620 ADU (ch1:MCSA04045.fits), 8070 ADU (ch2:MCSA04046.fits)

オン点のカウント 7300 ADU (ch1:MCSA04063.fits), 8900 ADU (ch2:MCSA04064.fits)

Twilight フラットは、空が暗い時のデータを OFF 点としてメジアン合成。明るくなってきた段階データを On 点とするが、明るさに応じてアノマリパターンが明瞭に生じるため、オフ点、オン点ともアノマリ処理をしてから合成する。

•　オフ点データ処理。オフ点画像リストの一枚目の画像を基準として、 2 枚目以降の画像を1 枚目カウントにスケールし、 1 枚目データを引き算する。すると 2 枚目データに入っているアノマリパターンが見えてくる。それを imsurfit でモデル化し、規格化済み 2 枚目データから引く。同様に全部のデータを処理してからメジアン合成。→ OFF 点データとする。

•　オン点処理：オン点リストの各々の画像から OFF 点合成データを引く。オン点リストの一枚目のデータについて、 0) と同じ行程を行ってアノマリパターンを除き、合成。規格化してフラットとする。

オフ点の基準画像カウント　　 2191 ADU (ch1: MCSA03751.fits) 2097 ADU(ch2:MCSA03752.fits)

オン点の基準画像カウント　　 5912 ADU (ch1: MCSA03793.fits) 6887 ADU(ch2:MCSA03792.fits)

ドームフラットをトワイライトフラットで割った結果

どちらもきれいなアノマリパターンが支配的。これに伴う誤差は今の最大20% 程度。これを抑えないと、ドームフラット自身の特性は見えにくい。

よく見ると、若干迷光の兆候が、アノマリパターンの X 方向の傾きとして見える（矢印）。起源は不明。

ｃｈ１　　　　　　　　　　　　　　　　　　　　　　　　　ｃｈ２0.9 1.0 1.1 1.2 0.9 1.0 1.1 1.2

アノマリパターンを減らすにはアノマリはフラットを作った際のカウントの違いを反映して現れるグローバルパターン。実験として、 Twilight データから、ドームの On-Off のカウントに近いカウントのセットを選び、それで仮 twilight フラットを作り、ドームフラットと比較する事を考える。

Dome のオンオフ点と完全に同じカウントセットのデータは残念ながらないため、それに近い値の 2 画像を使って On ・ OFF 点とし、グローバルなパターンの振る舞いを見る。

Ch1: オフ点カウント 5912 ADU (MCSA03793) 、オン点 7070ADU (MCSA03795)

Ch2: オフ点カウント 6886 ADU (MCSA03794) 、オン点 8304ADU (MCSA03796)

前頁の Twilit フラット作成時に比べて、 On-OFF カウント差は

　 CH1: 3721 ADU （前頁 TW ）→ 1149 ADU （今回： dome は 680ADU ）

　 CH2: 4790 ADU （前頁 TW ）→ 1418 ADU （今回： dome は 830ADU ）

となった。まだ一致とは程遠いものの、カウント差は 3 分の 1 以下まで減っている。

アノマリを除く：割り算の結果3 ページ目と同じ表示レベル。アノマリパターンは弱まった。まだ残っているのは、ドームの OnOff 差と Twilight データの OnOff 差がまだ 2 倍近くズレていためかと思われる。

なおアノマリ以外の成分がよりはっきり見えてきたが、各象限で不連続に振る舞うため迷光あるいはフラット光の非一様性とは別起源。これによる誤差成分は最大で 5% 。

ｃｈ１　　　　　　　　　　　　　　　　　　　　　　　　　ｃｈ２0.9 1.0 1.1 1.2 0.9 1.0 1.1 1.2

比較

ｃｈ１　　　　　　　　　　　　　　　　　　　　　　　　　　　　　ｃｈ２

x=1050—1150pixels の範囲の列を平均しプロットしてみた。横軸は y 座標。1025pixel から立ち上がるパターンが第一象限のアノマリパターン断面。

通常の Twilight フラットで割った。

カウントを合わせた仮のTwilight フラットで割った。

　 On-Off カウント差をドームに近づけた方はアノマリが半分弱。通常の Twilightのカウント差より 1/3 以上カウント差を近づけたのだから、その程度には減って欲しいが、そうはならなかった。アノマリパターンはカウント差には比例していない。　もっと低照度の Twilight を使って、より Dome フラット作成時に近いカウント差を再現した画像も試したが、却ってアノマリは増えてしまった。当たり前だが、アノマリパターンの再現には同じ照度かつ同じカウント差が要求される。

カウントを合わせたドームフラット／ Twilight フラット

の結果

チャンネル 1 断面

アノマリパターン

謎のへたりパターン

アノマリパターン

カウントを合わせたドームフラット／ Twilight フラット

の結果

チャンネル 2 断面第一、第三象限外側のへたりパターンは、強度が違うがチャンネル 1 と共通で出ている。電気的なもの？ アノマリパターン

アノマリパターン

へたりパターン

アノマリをもう少し調べる。アノマリパターンを引く操作は、フラット作成の要となると思われる。しかし、アノマリパターンは 2 つのフレームを比較して初めて見えてくるため、あるカウントのフラットに入っているアノマリパターンの絶対強度を直接測定する方法はない。単位時間に来る入射光強度と、おそらく露出時間にも依存すると思われる。

その意味でアノマリパターンの相対強度と入射光との関係を見ておくのは、アノマリの絶対量を予測できるかどうか考える上で重要である。そこで K バンド Twilight フラットデータを使って、入射光とアノマリの関係を見た。

まず、 Twilight の「オン点」データ全てからオフ点データ（スカイカウントが上がり始める前のデータ平均）を引く。次に、オン点データリストの最初の一枚をアノマリ基準フレームとし、オン点データリストの 2 枚目以降の各画面中央のカウント値にアノマリ基準フレームのカウントをスケールした上で、それぞれから引く（この行程で初めてアノマリパターンが見えてくる）。結果画像の第 3 象限 [400:900,400:500] 付近のカウントを計り、それをアノマリのカウント値とする。さらに、その画像の生画像（オフ点を引く前）の中央平均カウントを出し、アノマリカウント値との関係をプロットする。

アノマリと画像カウント値との関係

横軸は生画像の中央付近平均値。アノマリの基準画像（ on データ１枚目）のアノマリカウントは解析行程上 0 になる。この基準画像からどの位アノマリが増加したか、が示されている。＊指数関数でフィットできそう。 f(x)=a*exp(b*x)+c で試しにフィットすると (gnuplot) 、ch1: a=70 b=1.70 c=-200ch2: a=30 b=1.70 c=-80カウントが 0 でアノマリが 0 になると仮定すると、ゲタが評価できるか？？

ch1

ch2

• ある照度レベルでのアノマリパターンをフィットし、別の低めの照度のアノマリパターンにスケールして引き算してみると、小さなレベルでパターンの違いが見えるが、基本的には引き算できた。ただし ch2 の額縁模様は違う振る舞いをするなど、それなりに残差は出た。

• アノマリパターンを消したといっても、 アノマリ引きの基準にした画像（オン点リストの最初の画像）自体に入っているアノマリパターンは評価できないため、アノマリ基準画像にどの位アノマリパターンが隠れているかが、この方法でのアノマリ引きの精度を決める。最終的には、同じ星を画面のあちこちに置いて測光し、実際に測光結果がどうばらつくかを確認する必要がある。

• 今回の結果を他の観測フレームにも適用できるかについては、まずアノマリ引きの基準にする画像のカウントレベルが違えば、当然結果も違ってくる（次ページ図）。さらに、露出時間などの違いに対してアノマリの振舞いがどう違ってくるかも見極める必要があり、まだまだ基礎データの蓄積が必要である。

アノマリ評価の基準画像が変わると・・・

ch1

K バンドで評価したもの

J バンドでの評価（１）

J バンドでの評価（ 2 ）

J バンドで同じ実験をした。ただし、 (1) はアノマリ基準画像を 2300ADU 程度の照度で評価し、 (2) では K バンドでの基準画像カウントに近い 5500ADU の画像で見ている。

評価 (2) と K バンドの値がずれているが、 Off 点として扱った画像が K の 2230ADU より低いカウントのため（同じカウントの画像を OFF 画像に選ぶと、ほぼ一致することは確認してある）。

「 Off 点」として扱った画像のカウント(1500ADU)

アノマリ基準とした画像のカウント (2280ADU)

•ドームフラットとトワイライトフラットの差は、アノマリパターンが支配的だった。両チップにほとんど性質の差は見られなかった。それゆえ、カウントを制御しやすいドームフラットが推奨フラットとなる。

•ただし、チップ１は残像が盛大で、オン点を取った後のオフ点は 5 枚以上影響が消えないケースもあった。取得は注意を要する。

•アノマリパターンをどうフラットフレームから除去するかが正確なフラットの作成の要となる、と考えているが、そのためのドームフラットの取得手順を次ページにまとめた。

•アノマリパターン抽出の基準にする画像に内在するアノマリパターンは除去できず、これが不定性となる。オフ点とアノマリ基準画像とのカウント差を小さく（ <数 100ADU ）保つ事で不定性を小さく抑える事が可能と考えている。

•これらの作業の上で、アノマリパターンは割り算でなく引き算するものである、というのは仮定だが、エンジニアリング的には合理的との事。

•最終的な確認は視野の多くの場所に星を置き測光をして確認する必要がある。

まとめ

アノマリを考えたドームの取り方• ライト OFF でたくさん取る (K) 。• ライトちょっとだけつけて、低照度で数枚取る。 OFF 点 + 数 100ADU 程度。安定するま

で時間が要る。• ライトをしっかりつけて、適当な枚数だけ取る。サチるかなり手前の、ピーク値 15000A

DU 以下に抑える方が安全。安定するまで時間が要る。• ライト On を取ったら 5 枚程度ダークを撮って残像消し。

解析• オフ点を合成する。• 低照度 On, およびオン点画像からオフ点画像を引く。• 合成した低照度 On 画像をスケールして、各オン点画像から引きアノマリだけを引き出す。• J バンドはオフ点は不要か。低照度 On 画像をスケールして引くだけでよいかもしれない

（枚数は要る）。迷光、バイアスパターンをキャンセルする上で、小さな影響かもしれないが、オフ点があると確実。ただカウントはほとんどない。あるいは、低照度を 2 種類のカウントレベルで撮るのが推奨される。

• アノマリパターンをフィットした面を作り、オン点画像から引く。タスク。• 全部加工したら On 合成。低照度 On 画像に含まれるレベルのアノマリは残ってしまうが、

これは仕方ない。これに起因する誤差は、低照度オン画像の照度レベル依存。• 様々な輝度レベルが取れたら（ Twiligh ｔ）それでゲタを推定できるかもしれない。

＊ここまでの解析や考察は、今後データがたまって来たら修正されることもあると思う。

Twilight vs. Self flats について

＊３月８日に投げたレポート「 Twilight flat と Self flat による標準星解析」において使用した K バンド Twilight フラット画像は、ここでのアノマリ除去フラット作りの行程に従って処理した画像です。図１にあるように Self Flat 画像と Twilight flat 画像とは感度比がだいぶ違います。ここで使用している SelfFlat の画像は標準星データで一枚あたりのカウントが 3000-4000ADU 程度である事から、アノマリもそれほど深刻ではないかも知れません（ 12 ページの図からの類推）。しかし、一般的な天体画像のカウントレベルである 10000ADU 超のデータの場合、アノマリも大きくなって 10% 近い値になっている事も想像されます。セルフフラットの誤差はこのように大変大きな (>0.1mag) ものになる可能性が大きいと、現在までの解析が間違えていなければいえると思います。