基于VS平台D―TIN生长算法优化.docVIP

基于VS平台D―TIN生长算法优化 　　【摘 要】本文对传统不规则三角网生长算法实现的研究，并对其进行优化，利用了Visual Studio 2012平台强大的可视化用户界面及其编程语言的灵活性及简单易懂特点，基于各行业对于DEM的需要，从而开发出一种构网速度快，运行效率高的一种三角网构建算法。 　　【关键词】不规则三角网；Delaunay三角网；VS环境；算法优化 　　0 引言 　　地球表面高低起伏，呈现为一种连续变化的曲面，这种曲面无法用平面地图来确切表示，于是我们就利用数字高程模型来表达，这种方法已经被广泛使用。数字高程模型即DEM（Digital Elevation Model），是以数字形式按一定结构组织在一起，表示实际地形特征空间分布的模型，也是地形形状大小和起伏的数字描述。 　　由于地理信息系统的普及，DEM作为数字地形模拟的重要成果已经成为国家空间数据基础设施（NSDI）的基本内容之一，并被纳入数字化空间框架（DGDF）进行规模化生产，已经成为独立的标准基础产品[1]。DEM有三种主要的表示模型：规则格网模型，等高线模型和不规则三角网。格网（即GRID）DEM在地形平坦的地方，存在大量的数据冗余，在不改变格网大小情况下，难以表达复杂地形的突变现象，在某些计算，如通视问题，过分强调网格的轴方向。不规则三角网（简称TIN，即Triangulated Irregular Network）是另外一种表示数字高程模型的的方法（Peuker等，1978），它既减少了规则格网带来的数据冗余，同时在计算（如坡度）效率方面又优于纯粹基于等高线的方法。不规则三角网能随地形起伏变化的复杂性而改变采样点的密度和决定采样点的位置，因而它能够避免地形起伏平坦时的数据冗余，又能按地形特征点如山脊，山谷线，地形变化线等表示数字高程特征。 　　基于三角形的表面建模可适合所有的数据结构，且三角形在形状和大小方面有很大灵活性，能很容易地融合断裂线，生成线或其他任何数据，因此基于三角形的方法在地形表面建模中得到了越来越多的注意，已经成为表面建模的主要方法之一。C#语言简洁易学并且开发效率高，易于移植，对于学习GIS的学生来说无疑是接受很容易而且较快的一门计算机编程和开发语言，也是大多数学生最熟悉和了解的语言。正是基于对生成不规则三角网算法的研究和满足学GIS的学生对C#语言熟悉的情况下，本文就主要介绍用三角网生长算法生成不规则三角网及其在Visual Studio 2012环境下的实现。 　　1 TIN的算法种类及各算法特点 　　在介绍构成TIN各种算法之前我们要来了解认识一下一个重要法则――Delaunay三角网法则。通常构建三角网并不考虑地性线（山脊线，山谷线）的骨架作用，但是，由于用等高线数据构建三角网时，由于地形的复杂多样，有的地区存在因地形突变而形成的断裂线等特殊地貌。另外一些地区存在大面积水域等内部不需要构网的区域，因此，在精度要求较高的TIN中，必须考虑以上问题。因此此时应顾及地性线，断裂线，水域线等特殊情况，也就是应构建约束―Delaunay三角网。约束法是基于约束图计算约束D―三角剖分[2，4]（简称CDT，即Constrained Delaunay Triangulation）构造算法，这种Delaunay三角网满足这样的法则：Delaunay三角网为相互邻接且互不重叠的三角形的集合，每一个三角形的外接圆内不包含其他点。Delaunay三角网由对应Voronoi多边形的点连接而成。Delaunay三角形有三个相邻点连接而成，这三个相邻顶点对应的Voronoi多边形有一个公共的顶点，此顶点是Delaunay三角形外接圆的圆心。 　　1.1 三角网生成算法的分类 　　根据构建三角网的步骤，可将三角网生成算法分为三类：（1）分而治之算法（由Shmaos和Hoey提出），其基本思路是使问题简化，把点集划分到足够小，使其易于生成三角网，然后把子集中的三角网合并生成最终的三角网，用局部优化（LOP，即Local Optimization Procedure）算法保证其成为Delaunay三角网[3]，它的优点是时间效率高，但需要大量递归运算，因此占用内存空间较多，如果计算机没有足够的内存，这一方法就无法使用；（2）数据点渐次插入算法（由Lawson提出），其思路很简单，先在包含所有数据点的一个多边形中建立初始三角网，然后将余下的点逐一插入，用LOP算法保证其成为Delaunay三角网[3]。此算法虽然容易实现，但效率极低；（3）三角网生长算法，在这三种算法中，三角网生长算法在20世纪80年代以后的文献中已很少见，较多的是前两种算法[3]，三角网生长算法目前较少人研究，笔者在这里讨论的就是这一算法，该算法是由Michael