[論文レビュー] Singular solutions, momentum maps and computational anatomy

Colin J. Cotter, Darryl D. Holm|ArXiv.org|May 8, 2006

Nonlinear Waves and Solitons参考文献 24被引用数 15

ひとこと要約

この論文は、流体力学のオイラー＝ポワンカーテーリー理論を用いて計算解剖学におけるテンプレートマッチングを定式化し、微分同相群上の測地線流れとしてEPDiff方程式を導入する。Landmarks（特徴点）と画像の輪郭は、EPDiffの特異解であり、完全で重複のない表現と数値シミュレーションを可能にする運動量写像によって支配される。$H^1$正則化を伴うシンプレクティックな変分的粒子メッシュ法により、安定でソリトンに類似した変形が得られる。

ABSTRACT

This paper describes the variational formulation of template matching problems of computational anatomy (CA); introduces the EPDiff evolution equation in the context of an analogy between CA and fluid dynamics; discusses the singular solutions for the EPDiff equation and explains why these singular solutions exist (singular momentum map). Then it draws the consequences of EPDiff for outline matching problem in CA and gives numerical examples.

研究の動機と目的

• 流体力学の類似性を用いて計算解剖学の変分的枠組みを確立すること。
• 特徴点と画像の輪郭をEPDiff方程式の特異解として同定すること。
• 効率的で完全な形状表現のための運動量写像に基づくハミルトニアン定式化を開発すること。
• シンプレクティックな粒子メッシュ法を用いて数値的シミュレーションとデータ同化を可能にすること。
• 微分同相群の幾何的力学を通じて、多様な計算解剖学的手法を統合すること。

提案手法

• 速度場の$L^2$または$H^1$ノルムを用いた微分同相群上の測地線パスの最小化としてテンプレートマッチングを定式化する。
• 変分原理のオイラー＝ラグランジュ方程式としてEPDiff方程式を導出し、微分同相の測地線流れを支配する。
• ディラックのデルタ関数を用いた特異密度としての形状表現により、弱解と運動量写像の定式化を可能にする。
• 運動量写像を用いて、特徴点と輪郭をその正準的位置と運動量でパラメータ化する。
• 空間的離散化に三次スプラインを用い、$H^1$ノルムの効率的計算にFFTを用いた変分的粒子メッシュ（VPM）法を適用する。
• 離散的変分原理に基づいて導出されたシンプレクティックなEuler-A時間積分法を用い、幾何的構造を保存する。

実験結果

リサーチクエスチョン

• RQ1 EPDiff方程式は、計算解剖学における形状変形の測地線モデルとしてどのように機能するか？
• RQ2 運動量写像は、特徴点や輪郭のような特異解をパラメータ化する際に果たす役割は何か？
• RQ3 $H^1$正則化を施した変形は、画像の輪郭マッチングにおいてなぜソリトンに類似した挙動を示すのか？
• RQ4 シンプレクティックな粒子メッシュ法は、最小限の冗長性でテンプレートマッチングのための安定的かつ高精度な数値解を提供できるか？
• RQ5 随伴法を用いることで、最適な形状マッチングのための初期運動量をどのように計算できるか？

主な発見

• 計算解剖学における特徴点と画像の輪郭は、EPDiff方程式の特異解であり、1次元ではピークドソリトンに対応する。
• 運動量写像は、形状配置とその初期運動量の間の正準的同型写像を提供し、完全で重複のない表現を可能にする。
• 変分原理に$H^1$ノルムを用いることでカマッサ＝ホルム方程式が得られ、弾性的衝突特性を持つピークオン解を支持する。
• シンプレクティック時間積分を用いたVPM法は、ハミルトニアン系の幾何的構造を保存し、長期間にわたるシミュレーションの安定性を保証する。
• 数値結果では、$H^1$正則化作用と初期運動量に関する勾配降下法を用いた、形状間の最小コスト変形経路が成功裏に得られた。
• 点対点マッチングを避けるために、グリッド離散化された特異密度上でのカーネルベースの距離ノルムを用いることで、形状変動に対してよりロバストな性能が得られた。