[論文レビュー] Approximation of length minimization problems among compact connected sets

Matthieu Bonnivard, Antoine Lemenant+1 more|arXiv (Cornell University)|Mar 9, 2014

Computational Geometry and Mesh Generation参考文献 38被引用数 44

ひとこと要約

本稿は、2次元におけるコンpact連結集合に制約を課した長さ最小化問題に対して、極限において接続性を強制する重み付き測地線距離項を用いた、アンブロシオ＝トルトレッリ型の新規近似を導入する。この手法は、Γ収束を用いてステイナー問題、平均距離問題、p適合エネルギー問題を成功裏に近似し、数値的検証により位相場のゼロレベル集合から連結最小化子が出現することを示している。

ABSTRACT

In this paper we provide an approximation \\`a la Ambrosio-Tortorelli of some classical minimization problems involving the length of an unknown one-dimensional set, with an additional connectedness constraint, in dimension two. We introduce a term of new type relying on a weighted geodesic distance that forces the minimizers to be connected at the limit. We apply this approach to approximate the so-called Steiner Problem, but also the average distance problem, and finally a problem relying on the p-compliance energy. The proof of convergence of the approximating functional, which is stated in terms of Gamma-convergence relies on technical tools from geometric measure theory, as for instance a uniform lower bound for a sort of average directional Minkowski content of a family of compact connected sets.

研究の動機と目的

• 2次元における連結性制約付き長さ最小化問題の数値近似手法の不足を解決すること。
• 新たな重み付き測地線距離項を用いて、極限において接続性を強制する位相場アプローチを開発すること。
• アンブロシオ＝トルトレッリフレームワークを、ステイナー問題、平均距離問題、p適合エネルギー最小化問題などの問題に拡張すること。
• 幾何測度論の道具を用いて、提案された汎関数が元の幾何学的問題にΓ収束することを証明すること。
• 反復最適化と部分勾配移動を用いて、連結最小化子を生成する数値的に実装可能なアルゴリズムを提供すること。

提案手法

• 標準のモディカ＝モルトラ項と重み付き測地線距離項を組み合わせた新しい汎関数を導入：$ \frac{1}{4\varepsilon}\int_{\Omega}(1-\varphi)^2dx + \varepsilon\int_{\Omega}|\nabla\varphi|^2dx + \frac{1}{c_\varepsilon}\sum_{i=1}^N d_\varphi(x_i,x_1) $、ここで$ d_\varphi $は$ \varphi $で重み付けされた測地線距離である。
• 重み付き測地線距離を$ d_\varphi(x,y) = \inf\left\{ \int_\gamma \varphi \, d\mathcal{H}^1 \,;\, \gamma \text{ が } x \text{ と } y \text{ を結ぶ} \right\} $として定義し、高速マーチング法による数値計算を可能にする。
• すべての$ x_i $ を含み、$ \{\varphi = 0\} $ が経路連結であることを保証するため、$ \sum_{i=1}^N d_\varphi(x_i,x_1) = 0 $ ならば、極限における接続性が強制されることを用いる。
• Γ収束理論を用いて、$ \varepsilon \to 0 $ のとき、近似汎関数が元の幾何学的問題に収束することを証明する。方向的ミンコフスキー内容の一様下界に依拠する。
• 部分勾配移動法を用いて$ d_\varphi $およびその上微分を計算し、最適化のための降下法を可能にする。
• 反復的最適化戦略を採用し、$ \varepsilon $ を徐々に小さくする。各ステップを直前の解から初期化することで、局所最小値を回避し、連結最小化子への収束を促進する。

実験結果

リサーチクエスチョン

• RQ1 連結性制約付き長さ最小化問題を近似する位相場法を構築可能か？
• RQ2 明示的な位相的制約なしに、拡散界面モデルに連結性制約をどのように組み込むことができるか？
• RQ3 重み付き測地線距離項の導入により、$ \varepsilon \to 0 $ の極限において、位相場のゼロレベル集合が連結のままである保証はあるか？
• RQ4 この枠組みは、平均距離問題やp適合エネルギー最小化といった他の幾何学的問題に拡張可能か？
• RQ5 提案された数値的手法は、異なる問題設定や初期化条件下でも連結最小化子を安定して生成できるか？

主な発見

• モディカ＝モルトラ項と重み付き測地線距離項を組み合わせた提案汎関数$ G^{h}_{\Lambda,\varepsilon} $は、$ \varepsilon \to 0 $ のとき、元の幾何学的問題にΓ収束し、最小化子の収束を保証する。
• 数値シミュレーション（$ \varepsilon = 0.05 $）では、$ \phi $ のゼロレベル集合が連結であり、$ y_0 $ を含むことが確認され、$ y_0 $ が中心からずれていても同様であった。
• 長さペナルティパラメータ$ \Lambda $が大きくなるにつれ、連結集合$ \{\phi = 0\} $ の長さが短くなることが観察され、期待される物理的挙動と整合的であった。
• 減少する$ \varepsilon $を用いた反復的最適化戦略により、局所最小値を効果的に回避し、複数回の実行において一貫した連結最小化子が得られた。
• 部分勾配移動法により、$ d_\varphi $ およびその上微分の効率的計算が可能となり、凹項の最適化が現実可能となった。
• 本手法は、連結性制約下で、ステイナー問題、平均距離問題、p適合エネルギー最小化問題という3つの異なる幾何学的問題を、成功裏に近似した。