[論文レビュー] Inverse Problems with Invariant Multiscale Statistics

Ivan Dokmanić, Joan Bruna+2 more|arXiv (Cornell University)|Sep 18, 2016

Photoacoustic and Ultrasonic Imaging被引用数 4

ひとこと要約

本稿では、測定空間におけるデータ整合性と非線形マルチスケール散乱変換領域における統計的相関構造の2つの制約に交互に射影することで、不適切に定義された線形逆問題を解くための新規なアルゴリズムを提案する。位相情報を破棄することで強力なスペクトル相関が露わになるため、欠損スペクトルの安定的かつ線形な推定が可能となり、特に重度のデータ損失が生じる超解像およびトモグラフィーのタスクにおいて、凸正則化より優れた性能を発揮する。

ABSTRACT

We propose a new approach to linear ill-posed inverse problems. Our algorithm alternates between enforcing two constraints: the measurements and the statistical correlation structure in some transformed space. We use a non-linear multiscale scattering transform which discards the phase and thus exposes strong spectral correlations otherwise hidden beneath the phase fluctuations. As a result, both constraints may be put into effect by linear projections in their respective spaces. We apply the algorithm to super-resolution and tomography and show that it outperforms ad hoc convex regularizers and stably recovers the missing spectrum.

研究の動機と目的

• 凸正則化の限界を克服すること。特に、凸性が平均への回帰を引き起こし、高周波数成分の詳細が失われるような、重度に不適切に定義された逆問題において。
• フーリエ変換やその他の変換領域における位相ゆらぎに隠された信号の統計的相関を特定すること。
• 測定値の整合性とマルチスケール統計的構造を同時に強制するフレームワークを構築し、再構成の忠実度を向上させること。
• 極度のデータ損失下でも、従来の手法が失敗する超解像およびトモグラフィーのタスクにおいて、本手法の有効性を示すこと。

提案手法

• 位相情報を破棄することで強力なスペクトル相関を露わにする非線形マルチスケール散乱変換を適用し、反復的なウェーブレット係数の絶対値を計算する。
• 散乱係数を、スケール間の統計的相関が保持され、線形に推定可能な表現として用いる。
• 測定データと一致する信号の空間への線形射影により、データ整合性を強制する（例：超解像における低周波成分の置き換え）。
• 散乱領域における非線形射影により、正しい統計的構造を強制し、散乱係数が元の信号のそれと一致するようにする。
• データ整合性と統計的構造の2つの射影を交互に繰り返し、両方の制約を満たす信号に収束させる。
• 散乱係数が局所的に並進不変であり、非常に相関が強いことを利用し、欠損スペクトル成分のロバストな線形推定を実現する。

実験結果

リサーチクエスチョン

• RQ1 位相を破棄したマルチスケール変換が、元来隠れており、逆問題の正則化に有用な統計的相関を露わにできるか？
• RQ2 測定値の整合性と散乱に基づく統計的構造の両方を強制することで、凸正則化よりも重度に不適切に定義された問題における再構成性能が向上するか？
• RQ3 提案された反復的射影法が、極度のデータ損失下でも超解像およびトモグラフィーにおいて、安定的に欠損スペクトル成分を回復できるか？
• RQ4 再構成結果の高次モーメント（例：尖度）は元の信号と比較してどの程度か？また、それらは評価のための信頼できる指標として機能するか？

主な発見

• 本手法は、特に超解像における256倍のデータ損失下でも、ℓ1最小化や全変動正則化を上回り、高周波数の構造的詳細の再構成に優れた性能を示す。
• 20回の反復後、Ising超解像の再構成では過剰尖度が1290に達し、元の信号の1760に近く、ℓ1最小化ではわずか3164にとどまるため、高次モーメントの回復が著しく劣っていることが示された。
• 半角度ラドン測定を用いたIsingトモグラフィーでは、最初の反復で視覚的品質が向上し、最終的な反復ではTV法やℓ1ベースの手法よりも元の信号の共尖度テンソル構造をよりよく再現した。
• 前方作用素が非常にコherentlyな場合でも、本手法は欠損スペクトルの回復を安定化させ、ℓ1最小化が一意な解を導けない場合でも有効である。
• Cox過程超解像において、本手法の過剰尖度は514であり、ℓ1法の3164より低い結果となり、MSEが高めでも信号構造の保存が優れていることが示された。
• 極度のデータ損失（例：超解像における256倍）下でも、構造的詳細および高次統計量の回復に成功し、不適切に定義された問題に対するロバスト性を示した。