引き続き独自のブラインド-デコンボリューションの話をしていきます。コスト関数の極値問題から得られる係数を代入して $\langle 1\rangle_k\mapsto\{\psi,\partial_r\psi\}$ が計算され、ヒルベルト変換により $\{\psi,\partial_r\psi\}\mapsto f$ が計算されます。それぞれの演算は FIR で表現できるので、これらを畳込んで合成したフィルタも FIR になります。この最終的な FIR のフィルタサイズは、ラプラシアンの冪級数での膨張分とヒルベルト変換分の和、すなわち目視で決まる径 $r$ 強となります。

このサイズが要因となり、実像から $r$ 強の位置にアーチファクトが発生します。このアーチファクトは遮断周波数で発生するアーチファクトを回避するためにトレードオフで発生していると考えることができます。

このアーチファクトを消すためにはどのようにすれば良いでしょうか。

$r$ に関する平均処理ができれば、アーチファクト部分はブラーが掛かったように目立たなくなります。$r$ は冪級数の次数で変化させることができるので、幾つかの冪数に渡りデコンボリューション結果若しくはフィルタを平均処理すれば $r$ に関する平均処理が達成できます。

以前、冪数の上限は決まっていると説明しました。したがって、どのような冪数をサンプリングするか迷ってしまうかもしれません。概ねの傾向としては、冪数とデコンボリューションの効果には相関があると考えれば良いです。

このような説明をすると、アーチファクトは平均化され目立たなくなっても、トレードオフで分解能も悪くなるのではないかと心配されるかもしれません。しかし、或るレベル以上の冪数であれば、殆ど分解能に悪影響はないのでご安心ください。

以上により、分解能へのデコンボリューションの効果を維持した状態でアーチファクトのみを低下させることができます。不思議に思われるかもしれませんが、アーチファクトの総和という視点では変わっていないので不可解な話ではありません。

アーチファクトの平均をノイズレベルまで下げてしまえば、実用上遮断周波数で発生するアーチファクトは回避されたことになります。(written by J.Nishigata)

研究100- 遮断周波数で発生するアーチファクトの回避