一种快速生成平面Delaunay三角网的横向扩张法.PDFVIP

第 10卷 第 1期 地 球 信 息 科 学 Vol10, No1 2008年 2 月 GEO IN FORMA TION SC IENCE Feb. , 2008 一种快速生成平面 D elaunay三角网的横向扩张法 1 2 1 刘永和 , 王燕平 , 齐永安 ( 1河南理工大学 , 焦作 　454000; 2 中国科学院大气物理所东亚区域气候 - 环境重点实验室 , 全球变化东亚区域研究中心 , 北京 　100029) 摘要 : 目前已有多种基于平面上离散点集构造 D elaunay三角网的算法 , 其中三角网扩张法 、逐点插入法的平均 时间复杂度为 O ( n2 ) , 分治算法和其他分块合并算法能使平均时间复杂度接近线性 , 但增加了算法的复杂性 , 从而使浮点计算误差错误发生的机率增大 。本文作者提出了一种新算法 : 将用于构网的离散点集先按横坐标从 小到大排序 , 在空间上表现为从左到右排列 ; 然后先以点序列中的前三个点作为初始三角网 , 每次将剩余点集 中最左边的点联入三角网 , 最终得到一个三角剖分 , 再用 LOP法优化三角剖分 。该算法的优势是具有快速的三 角剖分过程 , 使整体的平均时间复杂度为 O ( n) , 并且构网效率高 , 算法简单 。 关键词 : D elaunay三角网; 快速算法 ; LOP优化 ; 不规则三角网; 时间复杂度 1　引言 似 O ( n) 的算法相比, 该算法思路更简单 , 效率更 高 , 避免了现有算法在上面介绍的三个方面的不 一个好的构网算法应该具有 以下特点 : ①可 足 , 因此 , 适合于海量数据的构网。因为该算法 靠 , 即保证算法能够生成正确的三角网及三角形 中每次从左向右扩张一个点 , 所 以, 笔者将其命 之间的拓扑关系 ; ②高效 , 即算法在生成三角网 名为 “横向扩张法 ”。 时要尽量节省时间 ; ③简单 , 即算法要尽量简单 易行 , 因为复杂的算法需要 占用大量编程和调试 2　横向扩张法构网思路 的时间 , 而且会产生更多不可预料的错误 。学者 们已经提 出了很多种生成 D elaunay 三角 网的算 本算法的基本思路是三角剖分 - LO P 优化构 法 , 其中影响较大的有三角网扩张法 [ 1～3 ] 、逐点 成 , 即先在点集构造一个所有三角形互不相交的 插入法 [ 4 ] 和分治算法 [ 5 ～8 ] , 它们的平均时间复杂 三角剖分 , 然后采用 L aw son 提出的 LO P 方法 [ 14 ] ( 2 ) ( 2 ) ( ) 优化所有三角形 , 使优化后的三角形满足空外接 度分别是 O n 、O n 和 O n 。经笔者测试 , 逐点插入法的效率要比三角网扩张法快三分之一 圆、最小角最大化的属性 。由两步法思想可以派 左右 。分治算法是采用递归分割的方式生成三角 生出很多种构网法 , 如逐点插入法 、凸包环切边 网子块 , 最后合并各子块 , 虽然每个子块的生成 界法等 , 这些方法的差别在于第一步即三角剖分 效率仍为 O ( n2 ) , 但因为子块的分割使每块的平 方法的不同 , 各种不同三角剖分方法的剖分效率 均点数减少 , 避免 了消耗时间随 n 的指数级增 有很大差别 。对于构