[论文解读] Geometric dynamics of optimization

François Gay–Balmaz, Darryl D. Holm+1 more|arXiv (Cornell University)|Dec 16, 2009

Mathematical Dynamics and Fractals参考文献 62被引用 6

一句话总结

本文提出优化动力学作为基于变质（metamorphosis）的图像配准的几何框架，其中形状变换通过惩罚性非微分同胚变形进行建模。通过将拉格朗日-泊松和哈密顿-泊松约化应用于最优控制问题，作者推导出广义EPDiff的哈密顿与拉格朗日运动方程，并得到新的可积系统，如两分量Camassa-Holm方程。

ABSTRACT

This paper investigates a family of dynamical systems arising from an evolutionary re-interpretation of certain optimal control and optimization problems. We focus particularly on the application in image registration of the theory of \emph{metamorphosis}. Metamorphosis is a means of tracking the optimal changes of shape that are necessary for registration of images with various types of data structures, without requiring that the transformations of shape be diffeomorphisms. This is a rich field whose possibilities are just beginning to be developed. In particular, metamorphosis and its related variants in the geometric approach to control and optimization can be expected to produce many exciting opportunities for new applications and analysis in geometric dynamics.

研究动机与目标

• 开发用于图像配准中优化问题的几何力学框架，特别是利用变质来放松微分同胚约束。
• 将图像配准重新解释为由变分原理与动量映射控制的演化动力系统。
• 通过最优控制与约化理论，将EPDiff方程推广至非微分同胚、带惩罚项的形状变换。
• 通过拉格朗日-泊松与哈密顿-泊松约化，为计算解剖学中的优化动力学建立统一的几何表述。
• 在多种系统（包括流体动力学、刚体与可积PDE）上展示该框架，揭示新的动力学结构。

提出的方法

• 将拉格朗日-泊松约化适配于带有惩罚项的最优控制问题，以建模形状变形中的不精确或软约束。
• 利用动量映射与李-泊松括号，从半直积李代数上的变分原理推导哈密顿表述。
• 构造勒朗德变换，将拉格朗日系统转换为以(μ, n, β)为变量的哈密顿系统，分别代表动量、密度与对偶变量。
• 利用斜对称矩阵𝒪推导泊松括号结构，编码约化系统的动力学。
• 将该框架应用于具体实例：重摆、刚体、流体方程与Camassa-Holm方程，展示最优约化如何导出新的可积系统。
• 使用回标签映射与微分同胚的左作用，以几何且坐标无关的方式建模流体与形状演化。

实验结果

研究问题

• RQ1 如何通过几何力学将图像配准中的最优控制问题重新解释为动力系统？
• RQ2 动量映射与李-泊松结构在连接数据结构（如特征点、轮廓）与演化PDE解之间起到何种作用？
• RQ3 引入非微分同胚变换的二次惩罚项后，如何导致超越EPDiff的新动力系统？
• RQ4 拉格朗日-泊松与哈密顿-泊松约化能否系统性地应用于推导优化动力学的运动方程？
• RQ5 从经典力学与流体模型的最优约化中，会涌现出哪些新的可积系统，例如Camassa-Holm方程？

主要发现

• 本文通过最优约化推导出两分量Camassa-Holm方程，作为一维Camassa-Holm方程的可积哈密顿推广。
• SE(2)上的变质方程被作为一般框架的特例推导得出，展示了形状与强度变化如何通过惩罚项耦合。
• 勒朗德变换后的变量(μ, n, β)满足由斜对称矩阵𝒪定义的李-泊松括号结构，该结构编码了系统的完整动力学。
• 泊松括号(8.61)在μ上线性，在n与β上二次，证实了该系统在半直积李代数对偶上的李-泊松性质。
• 该框架成功将EPDiff方程推广至非微分同胚、带惩罚项的变形，使计算解剖学中的应用范围更广。
• 动量映射在图像数据（如特征点、轮廓）与演化PDE的典型变量之间提供了典范同构，从而实现最优路径计算。