[论文解读] Riemannian geometries on spaces of plane curves

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

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

一句话总结

本文在平面曲线空间上引入了一种黎曼度量，通过参数 $ 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, bR^2) = \operatorname{Emb}(S^1, bR^2)/\operatorname{Diff}(S^1) $ 的几何结构，以建模计算机视觉中的形状空间。
• 分析在所提度量下测地线与截面曲率的行为，特别是非负曲率与负曲率区域之间的过渡。
• 开发并验证用于计算曲线间测地线的数值方法，克服离散化与参数化带来的挑战。

提出的方法

• 提出一种 $ \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 进行数值延续，通过包含曲率与加速度估计的离散能量泛函来近似测地线。
• 采用带约束的离散能量最小化方法，实现均匀参数化与法向运动，使用一阶有限差分和调和平均数对线段长度进行加权。
• 通过在 $ \bbR^2 $ 中绘制路径上的中间曲线，以三角形布局可视化 $ B_e $ 中的测地线与测地三角形，提升可解释性。

实验结果

研究问题

• RQ1 如何在平面曲线空间上构造一种黎曼度量，使得任意两条不同曲线之间的测地距离严格为正？
• RQ2 在所提度量下，平面曲线空间的截面曲率行为如何？其与曲线光滑性及曲率的关系是怎样的？
• RQ3 当正则化参数 $ A $ 变化时，两个椭圆之间的测地线如何演化？这揭示了形状空间几何的哪些特征？
• RQ4 尽管存在离散化伪影与参数化歧义，是否仍可在曲线的无限维空间中实现稳定且精确的数值测地线计算？

主要发现

• 当 $ A > 0 $ 时，度量 $ G^A $ 确保任意两条不同曲线之间的测地距离严格为正，解决了 $ A = 0 $ 情况下距离为零的退化问题。
• 该空间表现出双重特性：对于大而光滑的曲线，截面曲率为非负；但对于高曲率或高频扰动的曲线，曲率为负。
• 数值模拟显示，当 $ A = 1 $ 时，两个方向相差 $ 60^\circ $ 的椭圆之间的测地线接近于椭圆的旋转运动；而当 $ A = 0.01 $ 时，测地线涉及不对称的凸起演化。
• 由三个分别位于 $ 0^\circ, 60^\circ, 120^\circ $ 方向的椭圆构成的测地三角形，当 $ A = 1 $ 时其内角和为 $ 102^\circ $，表明曲率为负；而当 $ A = 0.01 $ 时，其内角和为 $ 216^\circ $，表明曲率为正。
• 所提出的数值方法，结合一阶有限差分与包含曲率和加速度项的离散能量泛函，即使在离散化挑战下，仍能产生稳定且视觉上合理的测地线。
• 该模拟框架通过添加有利于恒定速度参数化与法向运动的项，成功处理了形状空间的重参数化不变性，但初始值的质量至关重要。