[论文解读] Shape Splines and Stochastic Shape Evolutions: A Second Order Point of View

Alain Trouvé, François‐Xavier Vialard|arXiv (Cornell University)|Mar 19, 2010

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

一句话总结

本文提出了一种基于黎曼流形上哈密顿力学的二阶样条框架，用于在2D/3D特征点数据中建模平滑、时间连续的形状演化。通过将动量控制视为确定性变量，该方法将形状样条形式化为二阶能量的极小化问题，从而实现超越分段测地线的平滑插值。通过引入随机扰动，该模型可扩展至随机形状演化，为纵向形状数据的非参数统计分析提供了理论基础。

ABSTRACT

This article presents a new mathematical framework to perform statistical analysis on time-indexed sequences of 2D or 3D shapes. At the core of this statistical analysis is the task of time interpolation of such data. Current models in use can be compared to linear interpolation for one dimensional data. We develop a spline interpolation method which is directly related to cubic splines on a Riemannian manifold. Our strategy consists of introducing a control variable on the Hamiltonian equations of the geodesics. Motivated by statistical modeling of spatiotemporal data, we also design a stochastic model to deal with random shape evolutions. This model is closely related to the spline model since the control variable previously introduced is set as a random force perturbing the evolution. Although we focus on the finite dimensional case of landmarks, our models can be extended to infinite dimensional shape spaces, and they provide a first step for a non parametric growth model for shapes taking advantage of the widely developed framework of large deformations by diffeomorphisms.

研究动机与目标

• 解决一阶分段测地线插值在纵向形状分析中的局限性，即缺乏平滑性与统计稳健性。
• 基于黎曼流形理论，开发一种用于时间索引的2D/3D形状平滑插值的确定性二阶样条模型。
• 通过在哈密顿测地线方程中引入随机力，构建随机形状演化模型，实现形状变异性的概率建模。
• 在确定性与随机设定下，建立解的存在性与唯一性的理论基础。
• 提供一个统一的时空形状数据统计分析框架，可推广至无穷维形状空间。

提出的方法

• 基于变分原理，将特征点匹配问题形式化为微分同胚流形上的测地线问题，利用哈密顿方程求解。
• 在哈密顿系统中引入控制变量，以建模二阶形状样条，其中控制变量调控特征点轨迹的加速度。
• 通过流形上的变分法推导样条估计问题的欧拉-拉格朗日方程，确保插值路径的C2连续性。
• 利用动量映射与李代数结构，将动力学表达为伴随作用对偶形式，导出欧拉-庞加莱方程。
• 通过将确定性控制替换为随机力，将模型扩展至随机演化，得到具有良好定义动力学的随机哈密顿系统。
• 利用流形上随机微分方程的理论，证明随机模型的全局存在性与适定性，且噪声作用于动量变量。

实验结果

研究问题

• RQ1 如何利用黎曼几何与哈密顿力学，形式化时间索引的2D/3D形状的二阶平滑插值？
• RQ2 在特征点空间中，确定性形状样条解的存在性与唯一性的理论基础是什么？
• RQ3 如何通过哈密顿方程的随机扰动来建模真实、含噪声的生物形状演化？
• RQ4 所提出的样条模型与黎曼流形上经典三次样条之间有何关系？
• RQ5 该框架能否推广至无穷维形状空间（如曲线与曲面）？

主要发现

• 所提出的二阶样条模型在时间索引的特征点构型之间实现了平滑、C2连续的插值，克服了分段测地线模型的非光滑性缺陷。
• 在哈密顿框架下推导出形状样条的欧拉-拉格朗日方程，表明控制变量调控形状轨迹的加速度。
• 模型的随机扩展在动量上引入随机力，得到具有全局时间存在性的适定随机微分方程。
• 当控制变量设为零时，该模型在数学上等价于黎曼流形上的三次样条，验证了与既有理论的一致性。
• 通过相同的哈密顿与动量映射形式化，该框架可推广至无穷维形状空间（如曲线与曲面）。
• 数值实验表明，该模型能够生成平滑且具有生物学合理性的形状轨迹，在视觉与定量评估中均优于一阶插值方法。