2019-09-20カラーCMOSのフラット補正とホワイトバランス画像処理覚書2019/11/11追記:以下の記述は不正確な部分があります。フラット補正後の周辺のG浮きは、著者のダーク画像の取得ミスが原因だったと思われます。詳細はASI294MC-pro:ポカミス顛末報告とか、正しい処理フローのまとめなど。先日、蔵王で撮影したIC1396は、SharpCapのホワイトバランス設定の問題でライトフレームのホワイトバランスが緑に偏っていた。これをステライメージのマニュアルに記載されている手続きに従ってフラット補正・ダーク補正したところ、色ズレが起こって補正に失敗した。この結果から、カラーCMOSのFITS画像の場合、背景がニュートラルになるようにライトフレームのホワイトバランスを調整しておかないとフラット・ダーク補正で失敗するのだろうと推測した。これは、ホワイトバランスを後から調整できるデジカメのRAW画像と比較して対照的(だなあ)。ホワイトバランスが崩れたFITS画像をフラット補正・ダーク補正するための別の手続きも考えました、という話。FITS画像の通常のフラット・ダーク補正の手続きはじめにカメラはZWOのASI294MC-Proを使っています。SharpCapのフラット画像・ダーク画像の取得機能を使って得たマスターダークとマスターフラットを使って補正を行った場合のはなしです。マスターフラットはモノクロのFITS画像です*1*2。いって今回の話の一般性は、あまり損なわれないと考えています。フラット補正とダーク補正の手続きとその結果で、そのようにして得られたマスターフラットとマスターダークを使い、ステライメージのマニュアルに従って使って次のように前処理を行いました:FITS画像に対して、ダーク・フラット補正を行うRGGBのパターン(ASI294MC)でベイヤー変換を行い、カラー化オートストレッチ機能で各ライトフレームのWB調整コンポジット上記の方法でコンポジット後に得られた画像がこんな感じです:一見すると、フラット補正はうまくいっているように見えますが、実はこの画像は失敗しております。いいますのも、このまま彩度だけを強調しますとこんな風に周辺が緑になっています。RとBのフラットは正しく補正できていて、Gだけが補正不足過剰補正になってしまっているようなのです。ひとつ気になるのは、ガーネットスターの右のあたりの緑が強く出ているところ。ここってアンプグローが生じる位置ですよね。念のためマスターダーク画像を見てみますとまさにアンプグローの位置で、緑が強く出ています。ダーク補正も上手くできていないようです。ちなみに、DSSをつかった場合でも結果はうえと変わらず同じでした。原因は、ライトフレームのカラーバランスと思われるASI294MCの我々のファーストライトで、Crescent星雲を撮影したときは、上記のような問題は発生しませんでした。実は、あのときは偶然、撮影したライトフレームのWBがある程度ニュートラルに整っていたのです。それに対して、今回のIC1396の撮って出し画像はこのように思い切りGに偏っていました。撮影時はデジカメの感覚で「ホワイトバランスは後から調整すればよい」と思い込んでいて、このG偏りはまったく気にしていなかったのです。元の画像がGに偏っているのにも関わらず、ベイヤー変換前にモノクロのフラットと(ほぼ)モノクロのダークで補正を行うものだから、RとBは補正が合致しても、Gだけが補正不足過剰補正になる。いうのが、我々の推測です。(9/20追記:ここの記述,すこし話がアベコベですね。)(9/20追記:上の赤で強調した件,ちょっと数式とグラフで考えます。中心の輝度を10とし,周辺減光が二次関数で表されると仮定すればはという風に書けるはずです(係数2は式を簡単にするためで,意味なし)。5だけ大きいとします:グラフで表すと,下のような状態です。フラット補正というのは,この分布と同じ割合の輝度をもつ画像でそれぞれを割り算するわけです。割ってあげると当然ながらと定数になる。フラットが合うわけです。はとなって定数にはならない。グラフにするとこのように周辺が浮き上がってしまうわけです。(どうでも良いことですが,はてなのtex機能,昔に比べると数式が綺麗になりましたね。天文のブログで数式使いたい人多いと思うのですが,yahooブログが潰れたとき,はてなに引っ越してくれば良かったのにと思います。)一方、デジカメのフラット補正では、液晶ディスプレイを使用して、ライトフレームと同じカラーバランスになるようにフラットフレームを作成していたせいか、あるいはRStackerやDeepSkyStackerに,WBを調整してからフラット補正を行う演算機能があらかじめ備わっているのか不明ですが,ともかく上のような問題は生じていませんでした。以上の推測が正しければ、まずはベイヤー変換でカラー化してから、ホワイトバランス調整で背景をニュートラルにしてから、そのあとフラット・ダーク補正を行えばよいはずです。具体的には各ライトフレームをベイヤー変換でカラー化。各ライトフレームにオートストレッチ機能でWB調整いっぽう、マスターダークとマスターフラットもベイヤー変換してカラー化(マスターフラットはもともとSharpCapでモノクロ画像として取得しているので、ベイヤー変換してもモノクロのまま)。各ライトフレームに対して、ベイヤー変換後のマスターダークとマスターフラットを使って、ダーク・フラット補正コンポジットという流れです。その結果得られた画像を、再び彩度強調してみますと、こんなふうに、周辺のG偏りは解消しました。今度はBに偏っているように見えますが、これは画面全体のカラーバランス調整で解消します。厳密に言いますと、若干アンプグローの影響は残っており、周辺の輝度はすこし過剰補正になっています。この理由はわからないのですが、あらたにライセンス購入したフラットエイドプロに助けてもらいました。終わりにステライメージのマニュアルによると、ダーク減算はベイヤー変換によるカラー化の前に行うべきとあり、なぜならRGGBのベイヤーのうち、ひとつがノイズであった場合、カラー変換後にその影響が周辺のピクセルに波及してしまうから、とのことです。すると、上の手続きは、変換後にダーク減算してるので、うまくないです。できればダーク減算だけは、ベイヤー変換前に行うようにしたいですね。DSSだったのですが、今回の画像処理でステライメージの自由度の高さが気に入ってしまいました。離れることができなそうです。「カラーのフラットで補正を成功させるのは難しいよ」とソフトに怒られるので、このアドバイスに従っています。*2:ちなみに、SharpCapのこの機能はフラットダークもバイアスも自動で引いてくれるので便利です。01:30カラーCMOSのフラット補正とホワイトバランスTweet広告を非表示にする関連記事2019-11-12ASI294MC-pro:ポカミス顛末報告とか、正しい処理フローのまとめなど。今年の夏から天文部に導入したASI294MC-pro。RASA11inchのファーストライト(失敗)東北地方の夏休みは8月25日で終了。sky仙台市天文台スタッフブログはてなブログをはじめよう!snct-astroさんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。