[論文レビュー] A class of second-order geometric quasilinear hyperbolic PDEs and their application in imaging science

Guozhi Dong, Michael Hintermüller+1 more|arXiv (Cornell University)|Jan 1, 2019

Numerical methods in inverse problems参考文献 31被引用数 1

ひとこと要約

本稿は、画像処理のための第二階微分幾何的準線形双曲型偏微分方程式のクラスを導入し、減衰付き第二階全 Variation (TV) 流および平均曲率流に焦点を当てる。解の存在と一意性を確立し、簡単な初期データに対して解析的解を提示し、特にノイズ除去およびデジッタリング応用において、第一階流と比較して優れた収束性とエッジ保持性能を示す数値比較により、その有効性を検証する。

ABSTRACT

In this paper, we study damped second-order dynamics, which are quasilinear hyperbolic partial differential equations (PDEs). This is inspired by the recent development of second-order damping systems for accelerating energy decay of gradient flows. We concentrate on two equations: one is a damped second-order total variation flow, which is primarily motivated by the application of image denoising; the other is a damped second-order mean curvature flow for level sets of scalar functions, which is related to a non-convex variational model capable of correcting displacement errors in image data (e.g. dejittering). For the former equation, we prove the existence and uniqueness of the solution. For the latter, we draw a connection between the equation and some second-order geometric PDEs evolving the hypersurfaces which are described by level sets of scalar functions, and show the existence and uniqueness of the solution for a regularized version of the equation. The latter is used in our algorithmic development. A general algorithm for numerical discretization of the two nonlinear PDEs is proposed and analyzed. Its efficiency is demonstrated by various numerical examples, where simulations on the behavior of solutions of the new equations and comparisons with first-order flows are also documented.

研究の動機と目的

• 第二階微分幾何的準線形双曲型偏微分方程式を、最適化における第二階動的法の優れた数値的性能に着想を得て、画像修復のための改善を図る。
• 減衰付き第二階流を、全 Variation 流および平均曲率流に拡張し、画像のノイズ除去および変位誤差補正（例：デジッタリング）を実現する。
• これらの新しい偏微分方程式の解析的および数値的基盤を確立する。具体的には、解の存在、一意性、収束性を含む。
• 数値実験を通じて、第二階流がエッジをより効果的に保持し、振動を低減する点で、第一階勾配流を上回ることを示す。

提案手法

• 減衰付き第二階全 Variation 流（TVF）を、双曲型偏微分方程式として提案する：$\ddot{w} + \eta(t)\dot{w} = -\partial\Phi(w)$、ここで$\Phi(w)$は全 Variation 関数である。
• レベル集合のための減衰付き第二階平均曲率流を導入し、超曲面の進化をモデル化することで、周囲積分を最小化し、変位誤差を是正する。
• エネルギー安定性と数値収束性を保証するため、シンプレクティック・オイラー法を時間離散化に用いる。
• 時間ステップ演算子を表す行列$\mathbf{A}_k$を用いて、行列形式の半離散系$\mathbf{z}_{k+1} = \mathbf{A}_k \mathbf{z}_k$を導出する。
• Gershgorin 円定理および固有値解析を用いて、CFL 型条件のもとで離散スキームが収縮的であることを証明する。
• 正則化された平均曲率流偏微分方程式を実装し、適切な定義と数値的安定性を保証する。

実験結果

リサーチクエスチョン

• RQ1 ノイズ除去やデジッタリングなどの画像修復タスクに適した第二階微分幾何的偏微分方程式を定式化・分析することは可能か？
• RQ2 これらの新しい第二階双曲型偏微分方程式の解に対して、存在性および一意性の性質はどのようなものか？
• RQ3 第二階流の解は、第一階勾配流と比較して、収束速度およびエッジ保持性においてどのように異なるか？
• RQ4 第二階 TVF フレームワークにおいて、簡単な初期データに対して解析的解を導出することは可能か？
• RQ5 半離散数値スキームの安定性および収束性を保証するための時間ステップの条件は何か？

主な発見

• 減衰付き第二階全 Variation 流は一意解を有し、簡単な初期データに対して明示的な解析的解が導出される。
• 正則化された第二階平均曲率流に対して、解の存在および一意性が厳密に証明される。
• シンプレクティック・オイラー法に基づく数値スキームは、減衰係数$\eta$およびシステム行列の最大固有値に依存する時間ステップ条件のもとで収束性が保証される。
• 離散反復行列$\mathbf{A}_k$の固有値$\mu_{i,\pm}^k$が$|\mu_{i,\pm}^k| \leq 1$を満たすことが示され、収縮性および収束性が保証される。
• 数値比較により、第二階流が第一階流よりも収束が速く、特にエッジが豊富な領域でエッジをより効果的に保持することが明らかになった。
• レベル集合を平均曲率流で進化させることで、画像の変位誤差（例：デジッタリング）を効果的に是正し、元の画像構造を回復することができる。