不规则点建立TIN和等高线的方法.docVIP

[测绘]不规则点建立TIN和等高线的方法！TIN?????对于不规则分布的高程点，可以形式化地描述为平面的一个无序的点集P，点集中每个点p对应于它的高程值。将该点集转成TIN，最常用的方法是Delaunay三角剖分方法。生成TIN的关键是Delaunay三角网的产生算法，下面先对Delaunay三角网和它的偶图Voronoi图作简要的描述。Voronoi图，又叫泰森多边形或Dirichlet图，它由一组连续多边形组成，多边形的边界是由连接两邻点线段的垂直平分线组成。N个在平面上有区别的点，按照最近邻原则划分平面：每个点与它的最近邻区域相关联。Delaunay三角形是由与相邻Voronoi多边形共享一条边的相关点连接而成的三角形。Delaunay三角形的外接圆圆心是与三角形相关的Voronoi多边形的一个顶点。Delaunay三角形是Voronoi图的偶图，如图所示。对于给定的初始点集P，有多种三角网剖分方式，而Delaunay三角网有以下特性：1）其Delaunay三角网是唯一的；2）三角网的外边界构成了点集P的凸多边形“外壳”；3）没有任何点在三角形的外接圆内部，反之，如果一个三角网满足此条件，那么它就是Delaunay三角网。4）如果将三角网中的每个三角形的最小角进行升序排列，则Delaunay三角网的排列得到的数值最大，从这个意义上讲，Delaunay三角网是“最接近于规则化”的三角网。Delaunay三角形产生的基本准则：Delaunay三角形产生准则的最简明的形式是：任何一个Delaunay三角形的外接圆的内部不能包含其它任何点[Delaunay1934]。Lawson[1972]提出了最大化最小角原则：每两个相邻的三角形构成的凸四边形的对角线，在相互交换后，六个内角的最小角不再增大。Lawson[1977]又提出了一个局部优化过程LOP（LocalOptimization?Procedure）方法。如图所示。先求出包含新插入点p的外接圆的三角形，这种三角形称为影响三角形（InfluenceTriangulation）。删除影响三角形的公共边（图b中粗线），将p与全部影响三角形的顶点连接，完成p点在原Delaunay三角形中的插入。将该点集转成TIN，最常用的方法是Delaunay三角剖分方法，生成过程分两步完成：1）利用P中点集的平面坐标产生Delaunay三角网；2）给Delaunay三角形中的节点赋予高程值。 基于TIN绘制等高线直接利用原始观测数据，避免了DTM内插的精度损失，因而等高线精度较高；对高程注记点附近的较短封闭等高线也能绘制；绘制的等高线分布在采样区域内而并不要求采样区域有规则四边形边界。而同一高程的等高线只穿过一个三角形最多一次，因而程序设计也较简单。但是，由于TIN的存贮结构不同，等高线的具体跟踪算法跟踪也有所不同。对于记录了三角形表的TIN，按记录的三角形顺序搜索。其基本过程如下：1）对给定的等高线高程h，与所有网点高程zi（i=1,2,?，n），进行比较，若zi=h，则将zi加上（或减）一个微小正数ε0（如ε=10-4），以使程序设计简单而又不影响等高线的精度。2）设立三角形标志数组，其初始值为零，每一元素与一个三角形对应，凡处理过的三角形将标志置为1，以后不再处理，直至等高线高程改变。3）按顺序判断每一个三角形的三边中的两条边是否有等高线穿过。若三角形一边的两端点为P1（x1，y1，z1），P2（x2，y2，z2）则(z1-h)(z2-h)0表明该边有等高线点；(z1-h)(z2-h)0表明该边无等高线点。x，y）：4）搜索该等高线在该三角形的离去边，也就是相邻三角形的进人边，并内插其平面坐标。搜索与内插方法与上面的搜索起点相同，不同的只是仅对该三角形的另两边作处理。5）进入相邻三角形，重复第（4）步，直至离去边没有相邻三角形（此时等高线为开曲线）或相邻三角形即搜索起点所在的三角形（此时等高线为闭曲线）时为止。6）对于开曲线，将已搜索到的等高线点顺序倒过来，并回到搜索起点向另一方向搜索，直至到达边界（即离去边没有相邻三角形）。7）当一条等高线全部跟踪完后，将其光滑输出，方法与前面所述矩形格网等高线的绘制相同。然后继续三角形的搜索，直至全部三角形处理完，再改变等高线高程，重复以上过程，直到完成全部等高线的绘制为止。 DEM提取等高线 DEM生成等高线时，需要将其中的每个点视为一个几何点，而不是一个矩形区域，这样可以根据格网DEM中相邻四个点组成四边形进行等高线跟踪。其方法类似于后面描述的利用TIN提取等高线。实际上，也可以将每个矩形分割成为两个三角形，并应用TIN提取等高线算法，但是由于矩形有两种划分三角形的方