delaunay三角剖分随机增量算法的演示性实现及球面扩展.pdfVIP

Delaunay 三角剖分随机增量算法的演示性实现及球面扩展 2010310573 袁文 Delaunay 三角剖分随机增量算法的演示性实现及球面扩展 2010310573 袁文 一、 问题背景 在计算几何课上，我们学习了Delaunay 三角剖分算法，并且详细学习了一个 ( ) 期望时间复杂度为O n log n 的随机增量算法。本次大实验，我的主要目标就是 实现该算法，并且分步实现算法中的每一步骤，使算法可以逐点、逐步执行，达 到较好的演示效果。并进一步将该二维平面算法扩展到三维球面上，达到类似的 演示效果。  程序的二维平面演示效果如下： 1 / 17 Delaunay 三角剖分随机增量算法的演示性实现及球面扩展 2010310573 袁文  程序的三维球面演示效果如下： 二、 算法及原理 Delaunay 三角剖分算法参照课程讲义上的随机增量算法来实现。 二维平面算法中具体实现了“无穷远点”概念，在“toLeft ”测试以及“inCircle” 测试上加入了关于无穷远点的判断。 在初始化三角形的选择上，即无穷远点的选择上，使用了如下三个无穷远点： ( ) ( ) ( ) InfiA 0, +∞ , InfiB −∞, −∞ , InfiC +∞, −∞ 以及无穷远三角形InfiA − InfiB − InfiC。 三维球面算法中，首先明确了球面顶点、球面边、球面三角形、“toLeft ”测 试以及“inCircle”测试的概念。 2 / 17 Delaunay 三角剖分随机增量算法的演示性实现及球面扩展 2010310573 袁文 球面顶点的概念较简单，就是在三维直角坐标系中任意一个到原点距离为1 的坐标点。 对于任意两个不重合的球面顶点A 、B，若它们加上坐标原点O 三点不共线， 则这三点确定了一个唯一的平面AOB 。这个平面与单位球R （原点O 为球心，1 为半径）相交成一个以原点O 为圆心，1 为半径的圆c 。圆c 上以A 、B 为端点 ̂ 的劣弧AB即为球面边AB 。当A 、B 与原点O 三点共线时，球面边AB 没有定义。 ⃑⃑⃑⃑⃑ ⃑⃑⃑⃑⃑ 同时，圆C 将单位球R 分为了两个部分Rc1和Rc2 。规定向量⃑n = OA × OB的方向 为球面边AB 的法向量方向。定义球面上球面边AB 的左侧为Rc1和Rc2 中法向量 所在的那一侧。即，在定点V 上，球面边AB 的toLeft 测试返回结果为是否满足 ⃑⃑⃑⃑⃑ OV ∙ ⃑n ≥ 0。 球面上的任意三个顶点A 、B、C，若球面边AB 、BC、CA 都有定义，且A 、B、 C 三点不在同一单位圆上，则由球面边AB 、BC、CA 首尾相连组成的封闭图形定 义为球面三角形。 与二维平面上的情形相比较，三维球面的情形较为复杂：任意给定两个顶点 A 、B，不一定存在合理的球面边AB （A 、B 可能与原点O 共线）。同样的，任意 给定三个顶点A 、B、C，也不一定存在球面三角形ABC 。但可以通过在初始化时， 在球体上插入一个“正四面体”作为初始三角形集合，来避免这些情况的出现。 初始化正四面体如下图所示：