WFCAM Reduction - Making Mask to Reject Crosstalk

Crosstalk 除去

クロストーク除去は以下の手順で行っている。

特定のセット/検出器のフレーム(たとえば sdf1-?_1.fits とか)すべて
```
imcombine offsets=wcs
```
で足し合わせる
ピクセル値がある一定 X 以上になった場合にクロストークが生じると考えて、
```
imreplace [image] 0 upper=X
```
のようにして、ピーク以外をゼロにしてしまう
で、daofindでピーク検出。パラメータは以下の通り。
```
daofind fwhmpsf=5 thresho=20 verify- sigma=2
```
じつは、この手法では
- 低いピーク（Xをわずかしか越えていない）を取りこぼす
- 検出器のコドラントによって X が微妙に変わるらしく、低いレベルのクロストークの除去率がコドラントによって違うだろう
- そもそも、露出時間の違う(10sec/5sec)フレームを足し合わせているので 5secだとクロストークが生じていない場合も除去してしまう
という欠点がある。特に、2,3番目についてはちゃんと切り分けて処理するべきかもしれないが、今のところはこのまま。
この座標をregion fileに変換して、ds9の画面上で確認。
明らかに取りこぼしているもの、あるいはクロストークがないのにチェックされているものを編集して直す

完全に取りきれない/取りすぎていることはあるだろうが、明らかなものは確実に取れているはず。
まあ、天体の数が多いので統計的には問題ないだろうと言うことでこのまま進める。
（逆に、特殊な天体を探す際には十分な注意が必要となる。）

ポイントは以下の通り

クロストークは、採取画像では親ピクセルからx or y方向に +/-256ピクセル（生データだと128ピクセル）離れた位置に生じる。
また、6pix離れたpositive/negativeのピークとして発生するので、マスクの位置は 256+/-3pixの2個所となる。
ただし、強いピークはそれだけではなく、コドラントの mod 256 が同じになる行/列にも出る。
クロストークのサイズは、おおむね直径10pix以下であるが、これもやはり強いピークについては 20pix以上になることも多い。

ピーク強度の評価はなかなか難しい（特にピーク値がサチッてしまうとピーク値はあまり信用できない）が、とりあえず daofind が返してくる簡易 magnitude測定結果を使う。
この結果を主に、あとは目で見ながら

をmagnitudeで切り分ける。
で、マスクのサイズごとにマスクの場所の座標テーブルをつくって、それぞれの露出セットごとに、露出時間マップを

imedit [exposure map file] "" cursor=mask.coo display- radius=5 default="e" search=0

のようにして、クロストークのある場所の露出時間を0にしてしまう。
(なので、この時点では画像データにマスクをかけるわけではない)

最後に足し合わせる際に、exposure map の露出時間を参照にマスクをつくるので、その際に反映される。

以下に、クロストークを除去した実例を示す。

クロストークを除去前。緑丸はクロストークの原因となるピークの両隣、1次のクロストークが出現すると思われる場所

クロストークを除去後。緑丸の位置以外にもクロストークが出現していて、それもちゃんとマスクされているのがわかる。

Last Updated 2007/7/9
kmotohara@ioa.s.u-tokyo.ac.jp