[论文解读] Computing Quasiconformal Maps on Riemann surfaces using Discrete Curvature Flow

Wei Zeng, Lok Ming Lui+5 more|arXiv (Cornell University)|May 25, 2010

Advanced Numerical Analysis Techniques参考文献 29被引用 3

一句话总结

本文提出了一种新颖的数值方法，通过离散曲率流方法求解贝尔特拉米方程，以在一般黎曼曲面上计算拟共形映射。该方法引入一个辅助度量，使拟共形映射在新度量下变为共形映射，从而可通过离散亚梅贝流求解；随着网格分辨率的提高，该方法收敛至连续解，且在具有不同拓扑结构的真实扫描曲面上表现出高度的准确性和通用性。

ABSTRACT

Surface mapping plays an important role in geometric processing. They induce both area and angular distortions. If the angular distortion is bounded, the mapping is called a {\it quasi-conformal} map. Many surface maps in our physical world are quasi-conformal. The angular distortion of a quasi-conformal map can be represented by Beltrami differentials. According to quasi-conformal Teichmüller theory, there is an 1-1 correspondence between the set of Beltrami differentials and the set of quasi-conformal surface maps. Therefore, every quasi-conformal surface map can be fully determined by the Beltrami differential and can be reconstructed by solving the so-called Beltrami equation. In this work, we propose an effective method to solve the Beltrami equation on general Riemann surfaces. The solution is a quasi-conformal map associated with the prescribed Beltrami differential. We firstly formulate a discrete analog of quasi-conformal maps on triangular meshes. Then, we propose an algorithm to compute discrete quasi-conformal maps. The main strategy is to define a discrete auxiliary metric of the source surface, such that the original quasi-conformal map becomes conformal under the newly defined discrete metric. The associated map can then be obtained by using the discrete Yamabe flow method. Numerically, the discrete quasi-conformal map converges to the continuous real solution as the mesh size approaches to 0. We tested our algorithm on surfaces scanned from real life with different topologies. Experimental results demonstrate the generality and accuracy of our auxiliary metric method.

研究动机与目标

• 开发一种有效的数值算法，用于在一般黎曼曲面上计算拟共形映射，包括任意亏格的曲面。
• 解决在复平面上简单区域之外的曲面上数值求解贝尔特拉米方程这一长期存在的挑战。
• 在三角形网格上建立拟共形映射的离散类比，以支持计算几何处理。
• 通过基于度量的变换方法，实现从给定的贝尔特拉米微分重构拟共形映射。

提出的方法

• 通过定义一个使映射在新度量下变为共形的辅助度量，在三角形网格上构建离散拟共形映射。
• 基于给定的贝尔特拉米微分构造辅助度量，确保原始拟共形映射在新度量下变为共形映射。
• 应用离散亚梅贝流，将度量演化至常曲率，从而在辅助度量下计算共形映射。
• 利用离散亚梅贝流的解，恢复原始度量下的拟共形映射。
• 借助Teichmüller理论中贝尔特拉米微分与拟共形映射之间的1对1对应关系，确保解的唯一性与正确性。
• 通过分析能量泛函的Hessian矩阵，证明离散亚梅贝流在双曲设定下的局部凸性与收敛性，且其Hessian为正定。

实验结果

研究问题

• RQ1 能否通过离散曲率流方法在一般黎曼曲面上有效求解贝尔特拉米方程？
• RQ2 如何通过离散辅助度量构造，将拟共形映射转化为共形映射？
• RQ3 随着网格分辨率的提高，离散拟共形映射向连续解的收敛行为如何？
• RQ4 该方法在不同拓扑结构与真实几何形态的曲面上是否具有鲁棒性与高精度？
• RQ5 能否证明离散亚梅贝流收敛至满足辅助度量下拟共形映射条件的解？

主要发现

• 所提方法成功在任意亏格的黎曼曲面上计算了拟共形映射，包括来自真实扫描数据的复杂拓扑结构。
• 当网格尺寸趋近于零时，离散拟共形映射收敛至连续解，表现出良好的数值一致性。
• 辅助度量变换确保了拟共形映射在新度量下变为共形映射，从而可应用离散亚梅贝流。
• 证明了离散亚梅贝能最的Hessian在相关定义域上为正定，保证了流的局部凸性与收敛性。
• 离散亚梅贝流的解表现出曲率偏差的指数衰减，满足 |K_i(t) - K̄_i| ≤ c₁e^(-c₂t)，表明其收敛稳定且迅速。
• 在真实扫描曲面上的实验结果证实了该方法在多样拓扑与几何复杂性下的通用性、准确度与鲁棒性。