[论文解读] Symmetries in LDDMM with higher order momentum distributions

Henry O. Jacobs|arXiv (Cornell University)|Jun 14, 2013

Quantum chaos and dynamical systems参考文献 19被引用 4

一句话总结

该论文通过引入微分同胚的k-喷射（k-jets）来扩展大变形微分同胚度量映射（LDDMM），从而实现对局部旋转、缩放和剪切等高阶动量分布的建模。它识别出这些高阶粒子中的内在对称性，利用诺特定理推导出守恒动量，为医学图像配准中的稳定、精确数值积分提供了几何框架。

ABSTRACT

In some implementations of the Large Deformation Diffeomorphic Metric Mapping formulation for image registration we consider the motion of particles which locally translate image data. We then lift the motion of the particles to obtain a motion on the entire image. However, it is certainly possible to consider particles which do more than translate, and this is what will be described in this paper. As the unreduced Lagrangian associated to EPDiff possesses $\Diff(M)$ symmetry, it must also exhibit $G \subset \Diff(M)$ symmetry, for any Lie subgroup. In this paper we will describe a tower of Lie groups $G^{(0)} \subseteq G^{(1)} \subseteq G^{(2)} \subseteq...$ which correspond to preserving $k$-th order jet-data. The reduced configuration spaces $Q^{(k)} := \Diff(M) / G^{(k)}$ will be finite-dimensional (in particular, $Q^{(0)}$ is the configuration manifold for $N$ particles in $M$). We will observe that $G^{(k)}$ is a normal subgroup of $G^{(0)}$ and so the quotient $G^{(0)} / G^{(k)}$ is itself a (finite dimensional) Lie group which acts on $Q^{(k)}$. This makes $Q^{(k)}$ a principle bundle over $Q^{(0)}$ and the reduced geodesic equations on $Q^{(k)}$ will possess $G^{(0)} / G^{(k)}$-symmetry. Noether's theorem implies the existence of conserved momenta for the reduced system on $T^{\ast}Q^{(k)}$.

研究动机与目标

• 通过赋予粒子旋转、缩放和剪切等高阶局部变换能力，推广地标LDDMM。
• 利用微分同胚的k阶喷射（k-jet）建模这些变换，捕捉超越简单平移的局部图像形变。
• 识别这些高阶粒子的内在对称性结构，类似于杨-米尔斯理论中的规范对称性。
• 通过由对称性引起的辛约化，利用动量映射和诺特定理推导守恒量，提升测地线计算的数值稳定性。
• 为一般流形（包括医学影像中非欧几里得空间如2-球面）提供一种无坐标、几何化的框架。

提出的方法

• 将粒子运动从基流形M提升至微分同胚的k阶喷射丛，形成主丛Q^(k) → Q。
• 利用微分同胚的k阶泰勒展开（即k-jet）建模高阶形变，编码局部线性和非线性变换。
• 将内部自由度的对称群G^(k)_z识别为GL(d)与对称张量空间S^1_k的半直积，表示其内部自由度。
• 应用由对称性引起的辛约化，结合动量映射理论与诺特定理，推导守恒量。
• 将约化配置空间Q^(k)构作为以地标为底的纤维丛，结构群为GL(d) ⋉ S^1_k。
• 利用Clebsch方法推导约化空间上的测地线方程，并验证其与先前工作（SNDP13）的一致性。

实验结果

研究问题

• RQ1 如何在LDDMM框架内系统地建模超越平移的局部图像形变（如旋转、缩放和剪切）？
• RQ2 编码微分同胚k阶喷射信息的粒子，其代数与几何结构为何？
• RQ3 这些高阶粒子中出现的内部对称性是什么？它们如何与测地流中的守恒动量相关联？
• RQ4 如何在一般流形M上，以无坐标、流形自适应的方式形式化k-喷射粒子的对称性结构？
• RQ5 由此产生的守恒动量能否提升LDDMM实现中的数值稳定性和精度？

主要发现

• 微分同胚的k-喷射为局部图像形变提供了有限维表示，使旋转和剪切等高阶变换得以建模。
• 具有k阶结构的粒子构成主丛Q^(k) → Q，其结构群为GL(d) ⋉ S^1_k，其中S^1_k为对称(1,k)-张量空间。
• 当k=1时，对称群为GL(d)；当k=2时，其为半直积GL(d) ⋉ S^1_2，可捕捉二阶局部形变。
• 这些粒子的内部对称性通过诺特定理导致守恒动量，可用于设计稳定的数值积分器。
• 当M = ℝ^d时，约化配置空间Q^(k)微分同胚于ℝ^d × GL(d) × S^1_2 × ⋯ × S^1_k，与先前公式一致。
• 推导出的速度场与测地线方程与SNDP13中提出的完全一致，验证了该几何框架作为地标LDDMM的自然延伸。