Delaunay三角剖分算法.docxVIP

[1] [1] Delaunay 三角剖分算法 点集的三角剖分（Triangulation），对数值分析（比如有限元分析）以及图形学来说，都是 极为重要的一项预处理技术。尤其是 Delaunay 三角剖分，由于其独特性，关于点集的很 多种几何图都和 Delaunay 三角剖分相关，如 Voronoi 图，EMST 树，Gabriel 图等。 Delaunay 三角剖分有最大化最小角，“最接近于规则化的“的三角网和唯一性（任意四点不 能共圆）两个特点。 定义 编 辑 【定义】三角剖分 ：假设 V 是二维实数域上的有限点集，边 e 是由点集中的点作为 端点构成的封闭线段, E 为 e 的集合。那么该点集 V 的一个三角剖分 T=(V,E)是一个平面图 G，该平面图满足条件： 除了 端 点，平面图中的边不包含点集中的任何点。 没有相交边。 平面图中所有的面都是三角面，且所有三角面的合集是散点集 V 的 凸 包。 在实际中运用的最多的三角剖分是 Delaunay 三角剖分，它是一种特殊的 三 角 剖 分。 先从 Delaunay 边说起： 【定义】Delaunay 边：假设 E 中的一条边 e（两个端点为 a,b），e 若满足下列条件，则称之为 Delaunay 边：存在一个圆经过 a,b 两点，圆内(注意是圆 内，圆上最多三点共圆)不含点集 V 中任何其他的点，这一特性又称空圆特性。 【定义】Delaunay 三角剖分：如果点集 V 的一个三角剖分 T 只包含 Delaunay 边， 那么该三角剖分称为 Delaunay 三角剖分。 优化处理: 理论上为了构造 Delaunay 三 角 网 [2] ，Lawson 提出的局部优化过程 LOP(Local Optimization Procedure)，一般三角网经过 LOP 处理，即可确保成为 Delaunay 三角网， 其基本做法如下所示： 将两个具有共同边的三角形合成一个多边形。 以最大空圆准则作检查，看其第四个顶点是否在三角形的 外 接 圆之内。 如果在，修正对角线即将对角线对调，即完成局部优化过程的处理。 LOP 处理过程如下图所示： Delaunay 剖分的算法 Delaunay 剖分是一种三角剖分的标准，实现它有多种算法。 准则特性 编 辑 准则: 要满足 Delaunay 三角剖分的定义，必须符合两个重要的准则： 空圆特性：Delaunay 三 角 网是唯一的（任意四点不能共圆），在 Delaunay 三角 形网中任一三角形的外接圆范围内不会有其它点存在。如下图所示： 最大化最小角特性：在 散 点集可能形成的三角剖分中，Delaunay 三角剖分所形成 的三角形的最小角最大。从这个意义上讲，Delaunay 三角网是“最接近于规则化的“的三角 网。具体的说是指在两个相邻的三角形构成 凸 四 边 形的对角线，在 相 互 交 换后，六个内角 的最小角不再增大。如下图所示： 特性: 以下是 Delaunay 剖分所具备的优异特性： 最接近：以最近的三点形成三角形，且各 线 段(三角形的边)皆不相交。 唯一性：不论从区域何处开始构建，最终都将得到一致的结果。 最优性：任意两个相邻三角形形成的 凸 四 边 形的对角线如果可以互换的话，那么两 个三角形六个内角中最小的角度不会变大。 最规则：如果将 三 角 网中的每个三角形的最小角进行升序排列，则 Delaunay 三角 网的排列得到的数值最大。 区域性：新增、删除、移动某一个顶点时只会影响临近的三角形。 　　　Lawson 算 　　　 Lawson 算法产生的三角剖分 [3] 6.具有 凸 多 边 形的外壳：三角网最外层的边界形成一个凸多边形的外壳。 计算方法 编 辑 a)Lawson 算法 逐点插入的 Lawson 算法 是 Lawson 在 1977 年提出的，该算法思路简单，易于编 程实现。基本原理为：首先建立一个大的三角形或多边形，把所有数据点包围起来，向其 中插入一点，该点与包含它的三角形三个顶点相连，形成三个新的三角形，然后逐个对它 们进行空外接圆检测，同时用 Lawson 设计的局部优化过程 LOP 进行优化，即通过交换对 角线的方法来保证所形成的 三 角 网为 Delaunay 三角网。 上述基于 散 点的构网算法理论严密、唯一性好，网格满足空圆特性，较为理想。由其 逐点插入的构网过程可知，遇到非 Delaunay 边时，通过删除调整，可以构造形成新的 Delaunay 边。在完成构网后，增加新点时，无需对所有的点进行重新构网，只需对新点的 影响三角形范围进行局部联网，且局部联网的方法简单易行。同样，点的