[論文レビュー] Total-variation minimization with bound constraints

Musa Maharramov, Stewart A. Levin|arXiv (Cornell University)|May 21, 2015

Sparse and Compressive Sensing Techniques参考文献 5被引用数 5

ひとこと要約

本稿では、境界制約付き全 variation 正則化のための効率的な ADMM アルゴリズムを提示する。投影を split-Bregman フレームワークに統合することで、追加の反復を回避する。ノイズの多い表面隆上データから、高レベルのコンpartimentalized 圧力変化を物理的境界を満たす形で正確に逆推定可能であり、制約なしおよび単純な投影手法に比べ、鋭い対比を保持し、制約を満たす点で優れている。

ABSTRACT

We present a powerful and easy-to-implement algorithm for solving constrained optimization problems that involve $L_1$/total-variation regularization terms, and both equality and inequality constraints. We discuss the relationship of our method to earlier works of Goldstein and Osher (2009) and Chartrand and Wohlberg (2010), and demonstrate that our approach is a combination of the augmented Lagrangian method with splitting and model projection. We test the method on a geomechanical problem and invert highly compartmentalized pressure change from noisy surface uplift observations. We conclude the paper with a discussion of possible extension to a wide class of regularized optimization problems with bound and equality constraints.

研究の動機と目的

• 大規模な全変動量正則化最適化問題に、境界制約および等式制約を組み合わせた、ロバストでスケーラブルな手法を開発すること。
• 地質力学的逆問題において、制約なし TV 正則化が物理的に不適切な負圧力値を生じるという限界を解消すること。
• 最適性条件を破らずに、追加の計算コストを増加させることなく、境界制約を split-Bregman/ADMM フレームワークに自然に統合すること。
• ノイズの多い表面隆上観測と高レベルのコンパートメント化された圧力変化を含む地質力学的逆問題において、本手法の有効性を示すこと。
• 大規模非線形ソルバと互換性があり、追加の等式制約および不等式制約への拡張が可能な汎用フレームワークを提供すること。

提案手法

• ADMM を用いて、モデル、勾配、および投影モデル変数の三重分割により、TV 正則化境界制約問題を再定式化する。
• 勾配（x）および投影モデル（y）のための補助変数を導入し、それらを元のモデル（m）と等式制約で結ぶ。
• 制約を双対化するための増分ラグランジュ法を用い、部分問題を反復的に解く。TV 最小化と境界投影を統合する。
• 境界制約を直接 ADMM ループ内に組み込むために、投影演算子 Π(m) = max{min{m, m₂}, m₁} を用いる。これにより、別個の近位項を導入する必要がなくなる。
• TV 部分問題にはソフトしきい値処理を、データ適合部分問題には非線形共役勾配法などの反復的勾配ベース手法を用いる。
• 境界制約の強制を split-Bregman 迴帰に統合し、追加の外回り反復やペナルティ項の必要性を排除する。

実験結果

リサーチクエスチョン

• RQ1 制約なし TV 正則化逆問題において、収束性を低下させず、計算コストを増加させずに、境界制約を効率的かつ正確に強制できるか。
• RQ2 提案手法は、制約なし TV 正則化および単純な投影ベース手法と比較して、ブロック型でコンパートメント化された圧力モデルをどれほど正確に回復できるか。
• RQ3 split-Bregman フレームワークに境界制約を統合することで、解の品質および最適性がどの程度向上するか。
• RQ4 本手法は、ノイズの多い観測と複雑な制約を伴う大規模な地球物理的逆問題に適用可能か。
• RQ5 アルゴリズムは、収束が速く保たれるとともに、鋭い対比を維持し、有効な境界を満たすことができるか。

主な発見

• 境界制約付きアルゴリズムにより、制約なし TV 逆問題に見られた物理的に不適切な負圧力値が完全に除去された。これらは真のモデルに存在しない。
• 制約なし解で生じた負圧力の補償として発生した、x ≈ 1 km 付近の不自然な圧力スパイクが、本手法により除去された。
• 収束は速く、初期の数回の反復以外は、最終解の近傍に非常に短時間で収束した。
• 単純なモデルの境界への反復的投影手法は、KKT 条件を破り、定性的および定量的に劣った結果を生じた。
• 15% のノイズレベルでさえも、ノイズの多い表面隆上データから、ブロック型でコンパートメント化された圧力構造を正確に回復できた。
• ADMM ループ内に境界制約を統合することで、追加の計算コストが発生せず、split-Bregman フレームワークの効率性が維持された。