在不可定向的流形网格曲面上计算测地距离的一般方法.pdfVIP

陈双敏等：在不可定向的流形网格曲面上计算测地距离的一般方法 尽管所有的二流形曲面都具有 局“部与圆盘同胚”的性质，但其离散的网格表示却依赖于曲面的 全局可定向性：可定向的二流形 曲面常用半边结构 【，。】(halfedgedatastructure)来表达，而不可定向的 二流形曲面大多表达成若干多边形的集合 (polygonsoup1．在 以可定向网格 曲面作为主要研究对象的 数字几何处理领域，绝大多数算法[6一la]依赖于半边结构，因而无法直接应用于不可定向的情形． 从经典数学的观点来看 由不可定向的流形 曲面得到可定向曲面有两种方式 [，12]，一种是考虑 由 坐标和定向构成的二元组，计算不可定向流形的二重覆盖；另一种方式是考虑基本多边形的边界走 向， 找出逆转定向的环路，沿着这些环路把原始 曲面剪开，并将两份拷贝在断开处重新接在一起，保持定 向的连贯性．这两种方式得到同样的结果． 本文从离散化的角度分别考虑了上述两种方法 提出了把不可定向的二流形网格曲面上的测地距 离问题转化到可定向曲面上进行处理的一般算法框架．由此我们在不可定向的二流形网格 曲面与传统 数字几何处理方法之间搭起一座桥梁．为了展示该算法框架的普适应性，我们将其应用于不可定向曲 面上的三个重要场合，包括测地距离的求解、离散指数映射和最远点采样． 以下是本文的结构安排：第 2节叙述相关的研究工作，并同时介绍基本多边形和同伦基底的定义 以及半边结构的表示；第 3节阐述把不可定向的二流形网格 曲面转化到可定向 面的一般步骤；第 4 节展示不可定向曲面上的三个应用场合：最后给出本文的结论并提出未来的研究方向． 2 相关研究工作介绍及准备工作 2．1 相关研究工作介绍 测地线可以看作直线在弯 曲空间中的推广，是 Riemann几何中的重要概念．测地距离，即测地 线的长度，是曲面的内蕴度量．1984年，Sharir和 SchorrIs]率先提出一个时间复杂为 O(n3logn)的 算法，用于求解凸多面体上一点到其他任何点的测地距离．1987年，Mitchell等人f9】把该算法改善至 0(n0logn)，改进后的算法适用于一般的可定向网格 曲面．1990年，Chen和 Han[1o]进一步把算法的 时间复杂度和空间复杂度分别降至 o(n)和 o(n1．2009年，Xin和 Wang[11】提出了 过“滤定理”对 Chen和Han1【0】的算法加以改进，得到了时空效率更高的测地线算法．最近，Liu[13】把每条边正面和 反面的窗元组织在一个数据结构中，从而使Mitchell等人9【1的算法在效率上有显著的提高．有意思的 是，这种创新性改进脱离了半边结构，完全可以用于求解不可定向曲面上的测地距离．但该算法 占用 的实际空间仍然较多，而且与传统的精确算法有较大差异，实现起来 比较困难．此外，Liu等人_14】还基 于精确的测地度量构建了曲面上的等距线、平分线和 Voronoi图． 2．2 准备工作 本节主要介绍基本多边形、同伦基底和半边结构，为叙述我们的算法作准备．对于每个闭曲面 S， 存在一个由偶数条边组成的有向多边形P，通过成对地粘合边，由P构造出s．我们称多边形 P为 基本多边形 (fundamentalpolygon)．图 1为环状 曲面和 Klein瓶的基本多边形．事实上，这个构造可 以表示成一个长度为2 的字符串，它包含g个不同的符号，每个符号出现两次且带有指数 +1或 一1． 指数 一1对应于该边的定向与基本多边形 P 的定向相反，在这里，g实际上是 曲面的亏格，它与Euler 示性数 (Eulercharacteristic)的关系是 2—29，可定向， )=(： 一 l2一g， 不可定向． 780 陈双敏等：在不可定向的流形网格曲面上计算测地距离的一般方法 第 1步 将顶点唯一化，计算 由S诱导出来的图结构 G； 第 2步 在图G上求同伦基底和基本多边形； 第 3步 根据基本多边形中各边界的行进方向找出造成定向逆转的那些边界，并沿着这些环路将 原始曲面 s剪开，并根据定向的连贯性将两份拷贝拼合得到可定向的曲面 s． 3．1．1 将不可定向的网格曲面转化为图结构 如果输入的不可定向曲面 尚没有进行顶点唯一化，我们有必要先把坐标相同的顶点赋予同一个ID， 这是将不可定向的网格 曲面转化为图结构的