[论文解读] Computing distances and geodesics between manifold-valued curves in the SRV framework

Alice Le Brigant|arXiv (Cornell University)|Jan 11, 2016

Morphological variations and asymmetry参考文献 42被引用 19

一句话总结

本论文通过平方根速度函数（SRVF）框架，在流形值曲线空间上提出了一种重参数化不变的黎曼度量，实现了对流形几何结构和曲线重参数化均敏感的测地线与距离计算。该方法将切丛 TM 上的自然度量拉回到曲线空间，定义了一类一阶 Sobolev 度量，从而获得显式的测地线方程以及指数映射与 Jacobi 场的计算方法，并在双曲上半平面曲线上的雷达信号分析应用中得到验证。

ABSTRACT

This paper focuses on the study of open curves in a Riemannian manifold M, and proposes a reparametrization invariant metric on the space of such paths. We use the square root velocity function (SRVF) introduced by Srivastava et al. to define a Riemannian metric on the space of immersions M'=Imm([0,1],M) by pullback of a natural metric on the tangent bundle TM'. This induces a first-order Sobolev metric on M' and leads to a distance which takes into account the distance between the origins in M and the L2-distance between the SRV representations of the curves. The geodesic equations for this metric are given and exploited to define an exponential map on M'. The optimal deformation of one curve into another can then be constructed using geodesic shooting, which requires to characterize the Jacobi fields of M'. The particular case of curves lying in the hyperbolic half-plane is considered as an example, in the setting of radar signal processing.

研究动机与目标

• 在黎曼流形 M 上的开放、定向曲线空间上，发展一种重参数化不变的黎曼度量。
• 将 SRVF 框架从欧几里得空间推广至一般黎曼流形，同时保持曲线上几何保真性。
• 推导显式的测地线方程，并通过测地线射击法计算指数映射与 Jacobi 场，实现最优曲线形变。
• 将该框架应用于双曲上半平面，其动机源于其与雷达信号处理中高斯分布统计流形的等价性。
• 在该空间中计算曲线的 Fréchet 均值，以实现对局部平稳雷达信号的统计分析。

提出的方法

• 通过 SRVF 变换将切丛 TM 上的自然 L2 度量拉回到曲线空间 M = Imm([0,1], M)，定义其上的黎曼度量 G。
• 通过构造 G，使得对所有 φ ∈ Diff+([0,1])，均有 Gc◦φ(h◦φ, k◦φ) = Gc(h,k)，以确保重参数化不变性。
• 利用拉回结构推导 G 的测地线方程，从而可通过测地线射击法实现测地线的数值计算。
• 利用测地线方程在 M 上构造指数映射，以计算曲线之间的最优形变。
• 表征测地线上的 Jacobi 场，以支持测地线射击算法与形状空间分析。
• 将该框架应用于双曲上半平面 H，利用其已知的度量与测地线结构，在模拟雷达数据上验证算法。

实验结果

研究问题

• RQ1 如何在流形值曲线空间上定义一种重参数化不变的黎曼度量，以同时尊重基流形的内在几何结构？
• RQ2 该新度量下的测地线与指数映射具有何种结构？如何实现其高效计算？
• RQ3 与依赖于将所有曲线数据平行传输至单一切空间的先前方法相比，该 SRVF 方法有何优势？
• RQ4 该框架能否在非欧几里得空间（如双曲上半平面）中有效应用，特别是在信号处理场景中？
• RQ5 在该度量下，Fréchet 均值的行为如何？其在非平坦空间中对一组曲线的代表性如何？

主要发现

• 所提出的度量具有重参数化不变性，并在曲线空间上诱导出一种黎曼结构，该结构同时考虑了 M 中曲线起点之间的距离以及它们 SRV 表示之间的 L2 距离。
• 利用拉回结构显式推导出测地线方程，从而可通过指数映射与 Jacobi 场映射，实现测地线的数值计算。
• 该方法避免了将所有曲线数据平行传输至单一切空间的需要，从而在整条曲线上更完整地保留了流形曲率带来的几何信息。
• 在双曲上半平面中，该框架成功计算出 11 条雷达信号曲线的 Fréchet 均值，其中均值曲线（黑色）准确捕捉了在小转子转速变化下反射系数的平均演化过程。
• 该算法在模拟直升机雷达数据上得到验证，结果表明均值曲线可作为目标识别的稳健参考特征，其中红色与蓝色端点分别表示最高与最低转子转速。
• 理论证明确认了在双曲上半平面中对数映射、指数映射与平行传输映射的正确性，且通过 Möbius 变换与矩阵指数显式推导出测地线与向量传输公式。