[論文レビュー] Riemannian geometries on spaces of plane curves

Peter W. Michor, David Mumford|ArXiv.org|Dec 19, 2003

Geometric Analysis and Curvature Flows参考文献 8被引用数 10

ひとこと要約

本稿では、パrameter $ A > 0 $ を用いて曲率を組み込むことで、測地線距離を正則化する平面曲線空間上のリーマン計量を導入する。これにより、$ A = 0 $ の退化した場合とは異なり、異なる曲線間の距離が正しく保証される。測地線方程式を導出し、断面曲率を計算し、滑らかな曲線では曲率が非負であるが、高曲率または高周波数の摂動では負になることを示し、数値シミュレーションによりこの二重性を確認する。

ABSTRACT

We study some Riemannian metrics on the space of regular smooth curves in the plane, viewed as the orbit space of maps from $S^1$ to the plane modulo the group of diffeomorphisms of $S^1$, acting as reparameterizations. In particular we investigate the metric for a constant $A> 0$: $$ G^A_c(h,k) := \int_{S^1}(1+A\ka_c(þ)^2)< h(þ),k(þ) > |c'(þ)| dþ$$ where $\ka_c$ is the curvature of the curve $c$ and $h,k$ are normal vector fields to $c$. The term $A\ka^2$ is a sort of geometric Tikhonov regularization because, for A=0, the geodesic distance between any 2 distinct curves is 0, while for $A>0$ the distance is always positive. We give some lower bounds for the distance function, derive the geodesic equation and the sectional curvature, solve the geodesic equation with simple endpoints numerically, and pose some open questions. The space has an interesting split personality: among large smooth curves, all its sectional curvatures are $\ge 0$, while for curves with high curvature or perturbations of high frequency, the curvatures are $\le 0$.

研究の動機と目的

• 滑らかで正則な平面曲線の空間に非退化なリーマン計量を定義し、$ A = 0 $ の場合に測地線距離が曲線を区別できないという問題を解決すること。
• コンピュータビジョンにおける形状空間をモデル化するための軌道空間 $ B_e(S^1, \bbR^2) = \operatorname{Emb}(S^1, \bbR^2)/\operatorname{Diff}(S^1) $ の幾何的構造を研究すること。
• 提案された計量下での測地線および断面曲率の挙動、特に非負曲率と負曲率の領域への遷移を分析すること。
• 離散化とパrametrizationの曖昧さに起因する課題を克服し、曲線の無限次元空間における測地線の計算に安定かつ正確な数値手法を開発・検証すること。

提案手法

• パrametrization $ \operatorname{Diff}(S^1) $ 不変なリーマン計量 $ G^A_c(h,k) = \int_{S^1} (1 + A\kappa_c(\theta)^2) \langle h(\theta), k(\theta) \rangle |c'(\theta)| \, d\theta $ を提案。ここで $ \kappa_c $ は曲率であり、$ A > 0 $ が正則化を提供する。
• 曲線空間上で変分法を用いて測地線方程式を導出し、法線ベクトル場を形状空間への接ベクトルとして扱う。
• 空間の断面曲率を計算し、大きな滑らかな曲線では $ \geq 0 $、高曲率または高周波数の摂動を伴う曲線では $ \leq 0 $ であることを示す。
• MATLABを用いた数値連続法により、曲率と加速度推定を含む離散的エネルギー関数を用いて測地線を近似する。
• 一階の有限差分とセグメント長の調和平均を用いた重み付けを用いて、均一なパラメトリゼーションと法線方向運動を制約する離散的エネルギー最小化を実装する。
• 形状空間 $ B_e $ における測地線および測地線三角形を、$ \bbR^2 $ 上の経路に沿った中間曲線をプロットすることで可視化し、解釈可能性を高めるために三角形配置を用いる。

実験結果

リサーチクエスチョン

• RQ1 平面曲線空間に、任意の異なる曲線間の測地線距離が正であるようなリーマン計量をどのように構築できるか？
• RQ2 提案された計量下での平面曲線空間における断面曲率の挙動は何か？また、曲線の滑らかさや曲率にどのように依存するか？
• RQ3 正則化パrameter $ A $ を変化させた場合、2つの楕円間の測地線はどのように変化するか？これにより形状空間の幾何的性質はどのような情報を示すか？
• RQ4 離散化の影響やパラメトリゼーションの曖昧さに起因する課題を克服して、曲線の無限次元空間において安定的かつ正確な測地線計算が可能か？

主な発見

• $ A > 0 $ の場合、計量 $ G^A $ により任意の異なる曲線間の測地線距離が正に保証され、$ A = 0 $ の場合に距離がゼロとなる退化を解消する。
• 空間は二面性を示す：大きな滑らかな曲線では断面曲率が非負であるが、高曲率または高周波数の摂動を伴う曲線では負になる。
• 数値シミュレーションにより、$ A = 1 $ の場合、$ 60^\circ $ の方向差を有する楕円間の測地線は回転する楕円に近く、一方 $ A = 0.01 $ の場合、非対称な隆起の進化を示す。
• $ 0^\circ, 60^\circ, 120^\circ $ の方向に配置された3つの楕円からなる測地線三角形の内角の和は、$ A = 1 $ の場合 $ 102^\circ $ であり負曲率を示し、$ A = 0.01 $ の場合 $ 216^\circ $ であり正曲率を示す。
• 一階の有限差分と曲率および加速度項を含む離散的エネルギー関数を用いた数値手法により、離散化の課題にかかわらず安定的かつ視覚的に妥当な測地線が得られた。
• 初期化の質が極めて重要であるが、再パラメトリゼーション不変性を確保するため、一定速度のパラメトリゼーションと法線方向運動を促進する項を追加することで、シミュレーションフレームワークは形状空間の性質をうまく扱えた。