约束条件下不规则Delaunay三角网构建方法.pdfVIP

崔铁军 昊正升 李爱光 5 2 解放军信．怠工程大学测绘学院 郑州4005 摘要：车文详细论述了约束条件下不规则D。laonay三角网构建的矢量方法和橱格方法．在等高 线转换成三角网数字地面模型时，将地形信息作为约束条件下建立不规则DeIaunay三角同，可以解 决等高线结点三角网化时出现的两个问题(水平三角形和不希望的三角形边)，大大改善了三角网 BTM的质量． 关键词：数字地面模型Delaunay三角构网TIN 许多应用中，除了对三角形的形状有要求外，还对离散数据本身作出某种约束。例如， 某两点的连线对形状有决定性大的影响，必须考虑把它们加入到离散数据中，作为一种特征 约束来考虑，才能得到高质量的模型。在二维情况下，这种特征约束可归纳为两种形态：一 种是离散数据带有若干组“有向折线”；一种是带有若干组“封闭多边形环”。通常．我们 把上述问题称作为在约束条件限制下离散数据的不规则三角形构造问题。该问题的解决基于 两种方法：矢量方法和栅格方法。 1 基于矢量方式约束条件下TlN的构建 不相交(端点除外)的线段(约束条件)成为三角网的边。严格地讲，带有约束条件的三角 质。它对于描述地形特征及地表现象才其有重要意义。基于矢量方式约束条件下不规则三角 形的生成算法在国内外早有研究，朱庆和陈楚江(1998)提到了一种顾及地形特征的针对 带状TIN的算法，它在每构建一个三角网对均要判段三角形的边是否与地形特征线相交， 程序计算量较大。另一种有效算法是：首先将地形特征点、地形特征线、等高线看作离散点 进行标准三角剖分，然后强行嵌入不在剖分中的约束边，最后删除域外三角形。卢朝阳和吴 成柯(1997)介绍了一种“(m+2)边形对角线交换技术”，从而解决了影响域任意复杂的 情形，但该方法要比较点到直线的距离，因而程序实现量较大。与之相比，周晓云和刘慎权 (1996)提出的“轨迹生成”法则易于实现．但有局限性，会出现交换不下去的情形。而 李立新和谭建荣(1998)的“递减算法”和“循环算法”正可以解决此问题，但“循环算 法”简单，程序易于实现。为了解决上述问题，我们在进行标准三角剖分时采用前沿边推进 算法，在嵌入约束边时采用“循环算法”，保证在任何情况下构建三角网的正确性。 前沿边推进法是以包括所有离散点凸包的每一条边作为新生成三角形的边(所有三角 形的边按照逆时针存储)，向凸包内“推进”，找到一点满足Delaunay构网法则，对于每 条新生成的边，最多参与两次构网，直到构网结束。 (1)建立凸包 凸包(convex hull)是在二维平面上给定点集内的点构成的，是包含给定点集所有点 的凸多边形。这里我们使用Larkin改进和完善的凸壳生成算法，采用数据分块的方法提高 查询效率。其过程如下(图5)： 凸包，并按逆时针方向以链表的形式存储。 ②对于链表中相邻的两点v；和v…，查询位于其右侧距离最大的点(计算三角形面 积即可)。如果点存在或撮大距离为0，且点在线段VjV{+1之间则将该点插入Vi和V…之 间；否则不插入。重复上述过程，直到对链表中每一条边判断一次。 ③统计每次没有插入点的操作的次数，如果等于链表中顶点的总数，则终止判断，凸 包建构完成。整个凸包计算过程中，关键是查询到有向边距离最大的点。采用点集分块的方 法，没有必要在所有的点内查询，只需要查询有向边一侧我们“感兴趣”的区域即可。这样 可大大减二E五余计算。 图l凸包建立过程 (2)凸包“推进” 凸包“推进”过程如下(图2)： 角形(求出外接圆，判断外接圆内是否有点)。从边的中点所在的格网依次向边的左侧的格 网外扩查询，直到找到这样一点。 ②判断边与点构成的两条新边是否存在。如果存在且已使用两次(凸包上的边只使用 一次)，则将该边标识使用过两次，不生成新的三角形；如果存在但边只使用了一次．则将 改变标识使用过两次．生成新的三角形；如果不存在，将新生成的边加入elist中．则生成 新的三角形。 ③对新生成的三角形及边建立编号，几何拓扑关系及边的分块索引。 ④循环链表elist直到生成的边满足点、边、三角形的数量关系，结束构网过程。 圆黟零