[论文解读] The Square Root Velocity Framework for Curves in a Homogeneous Space

Zhe Su, Eric Klassen+1 more|arXiv (Cornell University)|Jun 9, 2017

Morphological variations and asymmetry被引用 4

一句话总结

该论文将平方根速度（SRV）框架扩展至齐性空间 $M = G/K$ 中的曲线，实现了形状空间中曲线的重参数化不变黎曼度量。通过将曲线提升至李群 $G$ 中的水平路径，该方法在商空间中计算测地线、Karcher均值并执行主成分分析，实验基于飓风路径数据，实现了高效的数值计算，并清晰地按区域聚类。

ABSTRACT

In this paper we study the shape space of curves with values in a homogeneous space $M = G/K$, where $G$ is a Lie group and $K$ is a compact Lie subgroup. We generalize the square root velocity framework to obtain a reparametrization invariant metric on the space of curves in $M$. By identifying curves in $M$ with their horizontal lifts in $G$, geodesics then can be computed. We can also mod out by reparametrizations and by rigid motions of $M$. In each of these quotient spaces, we can compute Karcher means, geodesics, and perform principal component analysis. We present numerical examples including the analysis of a set of hurricane paths.

研究动机与目标

• 开发齐性空间 $M = G/K$ 中无参数化曲线空间上的重参数化不变黎曼度量。
• 将平方根速度（SRV）框架推广至非欧几里得、弯曲的齐性空间。
• 在模去重参数化和刚性运动的商空间中，实现测地线、Karcher均值和主成分分析的计算。
• 将该框架应用于真实世界数据，特别是大西洋和太平洋飓风路径，以分析形状差异。
• 克服先前SRV方法的局限性，这些方法受限于欧几里得空间或需要曲率加权度量。

提出的方法

• 使用基于联络的水平提升方法，将 $M = G/K$ 中的曲线提升为李群 $G$ 中的水平曲线，以保持度量结构。
• 通过将 $AC([0,1],G)$ 上的 $L^2$ 度量拉回到水平提升空间，定义曲线空间上的重参数化不变度量。
• 采用离散测地线积分方案计算水平提升：给定 $\beta(t) \in S^n$，通过旋转 $R_{\beta(t),\beta(t+\Delta t)}$ 计算 $\alpha(t)$，以保持水平条件。
• 通过适配齐性空间几何结构的改进SRV变换实现度量，从而高效计算距离和测地线。
• 对 $\operatorname{Diff}^+([0,1])$ 和 $\operatorname{SO}(n+1)$ 的作用进行商化，得到用于形状分析的商形状空间。
• 使用多维缩放和Karcher均值计算，可视化并分析数据中的形状可变性。

实验结果

研究问题

• RQ1 平方根速度框架能否推广至非欧几里得齐性空间 $M = G/K$ 中的曲线？
• RQ2 所提出的基于提升的方法是否能在 $M$ 中曲线的形状空间上生成定义良好、重参数化不变的黎曼度量？
• RQ3 该框架中能否实现测地线、Karcher均值和主分量的高效数值计算？
• RQ4 在模去刚性运动的商空间中分析时，飓风路径的形状差异是否反映出区域模式？
• RQ5 绝对连续曲线空间在 $M$ 中的度量完备性是否表现良好且可计算？

主要发现

• 该框架成功将SRV方法推广至齐性空间中的曲线，实现了重参数化不变的形状分析。
• 水平提升方法确保 $G$ 中的测地线投影到形状空间 $\mathcal{S}^{AC}([0,1],M)$ 中的测地线，保持了度量结构。
• 在标准笔记本电脑上，150条飓风路径的成对距离计算耗时不足两分钟，证明了其高效性。
• 多维缩放揭示了在未模去刚性运动的形状空间中，大西洋和太平洋飓风路径存在明显的聚类，表明路径形状存在区域差异。
• 在模去刚性运动（$SO(3)$）后，聚类消失，表明原始数据中分离的主要原因是位置而非内在形状。
• 第一主成分方向捕捉了形状可变性，而第二主成分主要反映了路径长度的变化。