基于标准差地形三维表面模型建立方法.docVIP

基于标准差地形三维表面模型建立方法 　　摘要：针对二维数据域的Delaunay三角剖分推广到空间时存在的关键问题，提出一种基于标准差的地形三维表面模型建立算法。与空外接圆准则相比，该算法在计算时考虑了附加高程信息，引入标准差作为构网判断准则。剖析了标准差的含义，为什么使用标准差以及基于标准差的剖分准则，最后给出基于标准差的地形三维表面模型建立算法的具体步骤，并通过一个具体实例对算法进行了验证。 　　关键词：标准差；算法；表面模型；三角剖分；地形 　　中图分类号：TP311文献标志码：A 　　文章编号：1001-3695(2007)11-0295-03 　　 　　TIN（triangulated irregular network，不规则三角网）就是用一系列互不交叉、互不重叠的连接在一起的三角形来表示地形表面。TIN作为地表的数字化表现手段和分析工具，在诸多领域，如公（铁）路勘测设计一体化、地理信息系统等得到了广泛的应用[1]。其核心技术就是散点的三角剖分算法。三角剖分过程也就是TIN的建立过程，TIN的三角剖分就是按照三角剖分准则，将地形采样点用互不相交的直线段连接起来，并按一定的结构进行存储[2]。目前的三角剖分均是在二维平面上进行的，然后在三角形的顶点上叠加所对应的高程值，从而形成空间三角形平面。三角剖分中使用最广泛的算法是Delaunay直接三角剖分算法。Tsai在1993年根据其实现过程将DT三角剖分分成分割合并算法、三角网增长算法和逐点插入算法三类。其中，分割合并算法和三角网增长算法是静态过程，已形成的三角网不会因为新点的插入而破坏。而逐点插入算法是动态构网过程，新点的插入会导致已有三角网的改变，而且过程简单、容易编程及实现[3]。?? 　　然而这些传统的算法是在进行散点域的三角剖分时，依据散点的平面分布先确定一个三角形网络，然后将第三维数据如高程附加于平面三角网络上，从而形成包含连续的三角形面元的地表模型，因此它是一种面向曲面（三角形面元）的表面重构技术。面向曲面的建模方法仅考虑了曲面元素（三角形面元）的平面几何形状，而没有考虑其在空间的分布情况，故对一些特殊的约束条件如陡坎、虎口、悬崖等特殊地形需进行特殊处理[4]。这些基于二维平面上的三角剖分由于忽略了点的高程值，只考虑了x，y坐标，就会导致出现如下投影变形，形成的表面与实际地形表面相差甚远，如图1、2所示。具体来说，就是：a）有可能将正三角形投影成狭长三角形；b）有可能将狭长三角形投影成正三角形。?? 　　要解决上述问题，可以通过直接在空间上对采样点进行剖分，就是通常所说的体剖分。体剖分的关键是确定剖分原则、数据结构和算法设计。本文在逐点插入算法的基础上，研究了一种基于标准差的地形三维表面模型建立方法，理论和实验均证实了该算法简单有效，并且能解决平面三角剖分由于投影导致的剖分结果不准确的问题。该方法利用标准差的特性，在确定剖分准则时考虑了附加高程信息，即计算标准差时用到了空间点的x、y、z坐标。将标准差作为一种剖分准则，生成的表面模型更符合自然地形表面。?? 　　 　　1三角剖分准则分析?? 　　 　　三角剖分准则是指三角形网络的形成法则，它决定着三角形的几何形状和生成的TIN质量。应用不同的准则将会得到不同的三角形网络。目前常见的三角剖分准则有六种，分别是空外接圆准则、最大最小角准则、最短距离和准则、张角最大准则、面积比准则、对角线准则。其中，空外接圆准则、张角最大准则和最大最小角准则能做到三角网的惟一性，而这也正是Delaunay三角网所要求的[2]。?? 　　实际上，Lawson在1977年提出的LOP法则能优化任何三角剖分准则下得到的TIN，从而得到惟一的DT三角网络[5]。其基本思想是利用Delaunay三角形的空外接圆特性[1]，即一个三角形为 Delaunay三角形，该三角形外接圆中不含有其余任何数据点 (图3)。方法为：对具有公共边的两个三角形组成的四边形进行判断，如果其中一个三角形的外接圆中包含有另一三角形除公顶点外的第三顶点，则交换公共边，图4表示了该过程。?? 　　容易看出，现有的LOP法则考虑的是二维平面上的三角剖分，在LOP过程中，每个三角形都要经过空外接圆检测，而空外接圆检测公式中只涉及空间点的x，y坐标，忽略了空间点的高程信息，也就是在运算过程中没有考虑空间点的z坐标。这导致的后果是剖分是平面上的剖分，其剖分结果不能准确地反映空间真实地形。比如考虑空间四个点，其三维坐标分别为A(-2,0,0)，B(0,-23,42)，C(2,0,0)，D(0, 23,42)。经分析可以得出，在空间上，三角形ABD和BCD都是等边三角形，其边长为43（可通过空间两点的距离公式求各边边长）；而三角形ABC和ACD都是等腰