|
上次课程，我们讲解了拓扑圆盘间的调和映照。这次，我们讨论拓扑球面间的调和映照。 图1. 共形脑图。 拓扑球面间的调和映照在医学图像领域被频繁使用，特别是构造大脑皮层曲面到单位球面间的映射，我们称之为“共形脑图”，如图1所示。
非线性热流方法 外蕴方法 我们将目标曲面单位球面嵌在三维欧式空间中，这样曲面间的映射可以被表示成从源曲面到三维欧式空间的映射，并且像集被限制在单位球面上，
传统的热流方法将一个函数经过“热力学扩散”成一个调和函数，使得其调和能量随时间单调下降：
这样我们可以保证对任意时间和任意点，映射像一直在单位球面上，
非线性热流方法可以确保映射的调和能量单调递降，但是由于调和映射不唯一，计算并不稳定，我们需要添加更多的条件来确保解的唯一性。事实上，如果存在两个映射都是调和映射，则它们彼此相差一个球面到自身的共形变换。首先，我们用球极投影将球面映到扩展复平面上，
图4. 球极投影。 如图4所示，我们在单位球的北极放置一盏灯，从北极发射的光线将球面上所有点映到过南极的切平面上，同时将北极点映成无穷远点。直接计算表明，球极投影是共形映射。扩展复平面到自身的所有共形变换构成所谓的莫比乌斯变换群，每一个自映射都具有如下形式：
为了去掉莫比乌斯变换群带来的歧义性，我们添加一些归一化条件，例如我们要求映射满足
这个限制去掉了3个自由度，余下的莫比乌斯变换是单位球面到自身的旋转。非线性热流方法是稳定的，不会诱发旋转，因此这一归一化条件足够确保调和映射解的唯一性。 非线性热流方法将一个初始映射“扩散”成一个调和映射，那么我们如何选取初始映射？理论上，任何一个映射度为1的映射都可以，例如最为常见的高斯映射。但在实践中，光滑曲面被离散的多面体网格所逼近，优化过程有可能落入局部最优的陷阱。一种有效的避免局部最优陷阱的方法如下：我们将源曲面一分为二，将每一半用拓扑圆盘的调和映射映到单位圆盘，并且边界映射相互一致。然后我们用球极投影将两个单位圆盘映到上半球面和下半球面，如此得到曲面到单位球面的初始映射。然后，我们再用非线性热流方法将初始映射进一步扩散。 内蕴法 内蕴方法只需要曲面的黎曼度量，不需要曲面在欧式空间中的等距嵌入。外蕴法引导我们得到计算方法，内蕴法使我们能够洞察调和映射更为深刻的性质。假设曲面间的映射给定，我们取曲面的等温坐标：
我们记
映射的调和能量密度为
映射的调和能量为
进一步，我们得到调和能量的欧拉-拉格朗日为
从而，内蕴的非线性热流方程为
图6. 共形映射 调和映照和共形映照的关系 通常情况下，共形映射一定是调和的，调和映射不一定是共形的。在上图情形，当我们固定边界映射，优化调和能量，我们得到调和映射。如果，我们放开边界，令边界的像在单位圆上可以自由滑动，从而进一步减小调和能量，则我们得到共形映射。换言之，共形映射是所有调和映射中调和能量最小者。 但是，对于球面而言，拓扑球面间的调和映射一定是共形的。直观来说，拓扑球面间的映射可以在球面上自由滑动，从而使调和能量无障碍地达到最优，从而得到共形映射。下面，我们将这一直觉述诸严格的证明。我们定义曲面间映射所诱导的二次微分，霍普夫微分：
我们可以证明，如果霍普夫微分全纯，则映射必为调和；如果霍普夫微分为0，则映射必为共形。 曲面上所有的全纯二次微分构成一个群，根据黎曼-罗赫定理，亏格为g>1的曲面，这个群的维数为6g-6。黎曼-罗赫定理是指标定理的特殊形式，本质上是说流形上椭圆型偏微分方程解的空间维数被流形拓扑决定。 拓扑球面上所有的全纯二次微分必然为0。假如球面上存在一个全纯二次微分，那么它诱导一个带有奇异点的平直度量，奇异点对应着全纯二次微分的零点，因此奇异点处的曲率测度为-180度。根据高斯博内定理，总曲率等于球面的欧拉示性数乘以360度，从而奇异点的个数为0，全纯二次微分没有零点。我们能够取其平方根的一个分支，得到一个全纯1-形式。球面上的全纯1-形式的实部为调和1-形式。因为球面的一阶上同调群为0，所以调和1-形式必为0，进而全纯二次微分必为0。 因此，拓扑球面间的所有的调和映射，其霍普夫微分为全纯二次微分，必然为0。我们得到拓扑球面间的所有调和映射必为共形映射。 黎曼映照的算法 共形映射和调和能量 我们直接计算调和能量：
这里雅克比矩阵
因此我们得到
这里
是映射诱导的Beltrami系数，
图2. Beltrami系数和伸缩商的几何意义。 我们考察从人脸曲面到平面单位圆盘的映射，如图2所示，映射将曲面上的无穷小椭圆映到平面的无穷小圆。无穷小椭圆的偏心率，就是伸缩商K，偏心率和椭圆长轴方向构成了Beltrami系数。从上面的推导我们可以看出，调和能量的下界等于目标曲面的面积；当所有的伸缩商K都是1时，调和能量的下界能够被达到，这时，映射为共形映射。 换言之，在所有可能的调和映射中，如果存在共形映射，则共形映射的调和能量达到全局最优。 调和映照唯一性
当映射为调和映射时，其像为闭测地线。如果目标曲面的曲率有正有负，那么在同一同伦类中，测地线并不唯一。如果目标曲面的曲率处处为负，那么同伦的测地线一定重合。直观上讲，假如曲面上有两条彼此同伦的闭测地线，它们之间的同伦在曲面上扫过一个拓扑圆柱面。根据高斯-博内定理，我们有
等式左边，因为圆柱面的边界为测地线，所有第二项为0，第一项非正。等式右边为0，所以第一项为0，圆柱面面积为0。两条测地线彼此重合。 封闭曲面间的调和映射，如果目标曲面上的高斯曲率处处为负，且M在N上的像不是一条闭测地线，则同伦的调和映射必然重合。一种想法是基于调和能量的凸性，设定二阶光滑同伦连接着两个调和映射，
我们计算调和能量的二阶导数
拉回联络算子，源曲面上的标准正交基底和高斯曲率分别为
我们选取测地同伦，亦即
为测地线，那么最后一项恒为0，因为 因此调和能量的二阶导数恒正，调和能量为严格凸。 同时在时间为0和1点处，映射为调和映射，调和能量关于时间的一阶导数为0，由此我们有 调和映照存在性 我们首先定义两个辅助函数
那么，映射的雅克比行列式为
对于调和映射，我们得到Bochner公式
调和映射诱导的辅助函数的零点是孤立零点，我们可以定义零点的阶数。 如果辅助函数不恒为0，那么辅助函数零点的总阶数和曲面的拓扑及映射的映射度之间有着整体的关系
由此，我们能够在特定条件下，讨论调和映射的存在性。例如，不存在从拓扑环面到拓扑球面，度为一的调和映射。 如果调和映射存在，并且非共形，那么零点的总阶数应该非负，
矛盾。因此调和映射必然为共形映射，进一步映射必为分支覆盖映射。因为映射的度为1，所以映射为同胚。但是我们知道拓扑环面和拓扑球面之间不存在同胚，矛盾。由此，拓扑环面和拓扑球面之间不存在度为1的调和映射。 调和映照正则性 我们用反证法给出简略证明。假设在某一内点，雅克比行列式为负，因此区域
非空。研究函数
则它在D的边界上为0， 在D的内部处处为负，因此函数
为上调和函数(super harmonic), 在D的内部处处为正，那么雅克比行列式
在D的内部处处为正，这和D的定义相矛盾。 References 1. Xianfeng Gu, Yalin Wang, Tony F. Chan, Paul M. Thompson and Shing-Tung Yau. Genus Zero Surface Conformal Mapping and Its Application to Brain Surface Mapping. IEEE Transaction on Medical Imaging (TMI), 23(8):949-958, August 2004. 2. Craig Gotsman, Xianfeng Gu and Alla Sheffer. Fundamentals of Spherical Parameterization for 3D Meshes. ACM Transaction on Graphics (TOG), 22(3):358-363,2003.
|
|
|
来自： 九亩荷塘一书屋 > 《机器学习及人工智能》
进一步，这一映射由三个坐标函数来表示，
映射的拉普拉斯由坐标函数的拉普拉斯给出，
但是，在我们目前的情况下，映射的像被限制在单位球面上，因此每一点的像只能沿着单位球面的切方向移动，而无法沿着球面的法方向移动。因此，我们需要将传统的热流方法修改成非线性热流：
这里等式右侧是映射拉普拉斯的切向分量，
由此可见，所谓的“非线性”是指在每一点，每一步我们都需要向切空间投影，因此能量泛函的欧拉-拉格朗日方程不是线性偏微分方程，其离散化的近似不是线性方程组：
