[論文レビュー] Left-invariant Stochastic Evolution Equations on SE(2) and its Applications to Contour Enhancement and Contour Completion via Invertible Orientation Scores

Remco Duits, Erik Franken|ArXiv.org|Nov 6, 2007

Image and Signal Denoising Methods被引用数 4

ひとこと要約

本稿では、ユークリッド運動群 SE(2) 上の左不変確率的拡散方程式に対して、熱核として得られる正確なグリーン関数を用いて明示的解を提示する。これらの解により、輪郭の強調および補完が可能な可逆的方位スコアが可能となり、一様な曲率適合拡散および測地線に基づくモルフォロジーを介して、医療画像解析への応用が可能となる。

ABSTRACT

We provide the explicit solutions of linear, left-invariant, (convection)-diffusion equations and the corresponding resolvent equations on the 2D-Euclidean motion group SE(2). These diffusion equations are forward Kolmogorov equations for stochastic processes for contour enhancement and completion. The solutions are group-convolutions with the corresponding Green's function, which we derive in explicit form. We mainly focus on the Kolmogorov equations for contour enhancement processes which, in contrast to the Kolmogorov equations for contour completion, do not include convection. The Green's functions of these left-invariant partial differential equations coincide with the heat-kernels on SE(2), which we explicitly derive. Then we compute completion distributions on SE(2) which are the product of a forward and a backward resolvent evolved from resp. source and sink distribution on SE(2). On the one hand, the modes of Mumford's direction process for contour completion coincide with elastica curves minimizing $\int κ^{2} + εds$, related to zero-crossings of 2 left-invariant derivatives of the completion distribution. On the other hand, the completion measure for the contour enhancement concentrates on geodesics minimizing $\int \sqrt{κ^{2} + ε} ds$. This motivates a comparison between geodesics and elastica, which are quite similar. However, we derive more practical analytic solutions for the geodesics. The theory is motivated by medical image analysis applications where enhancement of elongated structures in noisy images is required. We use left-invariant (non)-linear evolution processes for automated contour enhancement on invertible orientation scores, obtained from an image by means of a special type of unitary wavelet transform.

研究の動機と目的

• 2次元ユークリッド運動群 SE(2) 上の左不変拡散過程の明示的解を導出すること。これにより、輪郭の強調および補完が可能となる。
• 補完分布のモードと、曲率に基づくエネルギー関数を最小化するエラスティカ／測地線の間の関係を確立すること。
• ユニタリウェーブレット変換を用いた可逆的方位スコアのフレームワークを構築し、補完プロセスにおけるソースおよびシンクの両方の分布を可能とすること。
• SE(2) 上の微分幾何学を用いて、方位スコアにおける曲率適合拡散を組み込むことで、一様性を高める拡散の改善。
• ハミルトン＝ジャコビ方程式および測地線に基づく、グラフィカルなスリムニングおよびエロージョンのためのモルフォロジー理論の導入。

提案手法

• 対流項なしの前向きコルモゴロフ方程式の解として、SE(2) 上の左不変拡散に対する正確なグリーン関数（熱核）を導出する。
• グリーン関数との群畳み込みを用いて、前向きおよび後向きのリゾルベント進化により、補完分布を計算する。
• ホルマンダー条件の下で、SE(2) 上の熱核の取り扱い可能な近似を得るために、ヘイセンベルク型群の近似を用いる。
• カルタン接続および主バンドル理論を用いて水平曲線を定義し、方位スコアから空間的曲率を抽出する。
• ブライアントとグリフィスの還元法をオイラー＝ラグランジュ方程式に適用し、特殊関数を用いずに解析的で非特殊関数の解を得る。
• ゲージ座標を用いて曲率を組み込むことで、方位スコア上での非線形的適合拡散を導入し、標準的一致性を高める拡散を改善する。

実験結果

リサーチクエスチョン

• RQ1 左不変拡散過程を、熱核およびグリーン関数を用いて、SE(2) 上でどのように明示的に解くことができるか？
• RQ2 補完分布のモードと、曲率エネルギーを最小化するエラスティカまたは測地線の間にはどのような関係があるか？
• RQ3 SE(2) 上の測地線は、形状的特徴および解析的取り扱いやすさの観点で、エラスティカ曲線とどのように比較できるか？
• RQ4 可逆的方位スコアを、輪郭の強調および補完におけるソースおよびシンクの両方として使用できるか？
• RQ5 交差する輪郭が存在する状況下で、方位スコア上での曲率適合拡散は、一貫性を高める拡散をどのように改善できるか？

主な発見

• SE(2) 上の正確な熱核が閉形式で導出され、左不変拡散過程の明示的グリーン関数が得られた。
• 補完分布のモードは、∫(κ² + ε)ds を最小化するエラスティカ曲線と一致し、2つの左不変微分のゼロクロスを介して検出可能である。
• 予想される寿命が小さい場合（1/α → 0）、輪郭強調のための補完測度は、∫√(κ² + ε)ds を最小化する測地線に集中する。
• ノエーターの定理を用いた Pfaffian 系の還元により、特殊関数を用いずに SE(2) 上の測地線が解析的に解ける。
• ヘイセンベルク近似により、熱核に対してガウス上界が得られ、鋭い定数 c = √[4]{2} ≈ 1.19 を有する。
• 方位スコア上での非線形的適合拡散は、交差する輪郭を効果的に処理でき、標準的一致性を高める拡散に対して定性的な改善を示した。