第6章-数字地形模型与地形分析课件

数字地形模型（DTM）与地形分析数字地形模型（DTM）与地形分析第6章-数字地形模型与地形分析课件1．概述数字地形模型（DTM,DigitalTerrainModel）最初诞生于高速公路的自动设计。应用领域：用于各种线路选线（铁路、公路、输电线）的设计以及各种工程的面积、体积、坡度计算，任意两点间的通视判断及任意断面图绘制。测绘：用于绘制等高线、坡度坡向图、立体透视图，制作正射影像图以及地图的修测。遥感：作为分类的辅助数据。GIS：基础数据，可用于土地利用现状的分析、合理规划及洪水险情预报等。军事：用于导航及导弹制导、作战电子沙盘等。1．概述数字地形模型（DTM,DigitalTerrai写景法写景法根据等高线素描的地貌写景图根据等高线素描的地貌写景图根据等高线作密集而平行的地形剖面，然后按一定的方法叠加根据等高线作密集而平行的地形剖面，晕滃法是沿斜坡方向布置晕线表示地貌的一种方法。

晕滃法是沿斜坡方向布置晕线表示地貌的一种方法。

1.1DTM和DEMDTM：数字地面模型（DigitalTerrainModel，简称DTM）是定义于二维区域上的一个有限项的向量序列，它以离散分布的平面点来模拟连续分布的地形。地形模型不仅包含高程属性，还包含其它的地表形态属性，如坡度、坡向等。DEM：数字地形模型中地形属性为高程时称为数字高程模型（DigitalElevationModel,简称DEM）广义的DEM还包括等高线、三角网等所有表达地面高程的数字表示。在地理信息系统中，DEM是建立DTM的基础数据。1.1DTM和DEMDTM：数字地面模型（Digital1.2DEM的表示法1）数学方法整体拟合方法：根据区域所有的高程点数据，用傅立叶级数和高次多项式拟合统一的地面高程曲面。局部拟合方法：将地表复杂表面分成正方形规则区域或面积大致相等的不规则区域进行分块搜索，根据有限个点进行拟合形成高程曲面。1.2DEM的表示法1）数学方法1.2DEM的表示法2）图形方法（1）线模式等高线：山脊线、谷底线、海岸线及坡度变换线等。（2）点模式：离散采样数据点规则格网采样：密度一致的或不一致的；不规则采样：不规则三角网、邻近网模型等；选择性地采样：山峰、洼坑、隘口、边界等重要特征点。1.2DEM的表示法2）图形方法DEM的表示方法DEM的表示方法DEM的主要表示模型2.1规则格网模型通常是正方形，每个格网单元对应一个高程值DEM的主要表示模型2.1规则格网模型DEM的主要表示模型2.1规则格网模型通常是正方形，每个格网单元对应一个高程值DEM的主要表示模型2.1规则格网模型规则格网模型的特点很容易用计算机进行处理，如计算等高线、坡度坡向、山坡阴影和自动提取流域地形，是DEM最广泛使用的格式。缺点：格网DEM的缺点是不能准确表示地形的结构和细部。解决办法：附加地形特征数据，如地形特征点、山脊线、谷底线、断裂线，以描述地形结构。规则格网模型的特点很容易用计算机进行处理，如计算等高线、坡度2.2等高线模型每一条等高线对应一个已知的高程值，等高线集合和它们的高程值一起就构成地面高程模型。2.2等高线模型每一条等高线对应一个已知的高程值，等高线集第6章-数字地形模型与地形分析课件第6章-数字地形模型与地形分析课件第6章-数字地形模型与地形分析课件陡石崖崩崖冲沟陡崖陡石崖崩崖冲沟陡崖等高线模型等高线模型只表达了区域的部分高程值，要用插值方法来计算落在等高线外的其它点的高程，因为这些点是落在两条等高线包围的区域内，通常只使用外包的两条等高线的高程进行插值。等高线模型等高线模型只表达了区域的部分高程值，要用插值方法来等高线模型存储方法二维链表存储坐标用图来表示等高线的拓扑关系，将等高线之间的区域表示成图的节点，用边表示等高线本身。这种拓扑关系图可以改造成一种无圈自由树等高线模型存储方法等高线和相应的自由树等高线和相应的自由树2.3不规则三角网模型（TIN）TIN的产生原因——规则格网DEM存在的缺陷：1）在地形平坦的地方，存在大量的数据冗余；2）在不改变格网大小的情况下，难以表达复杂地形的突变现象；3）在某些计算，如通视问题，过分强调网格的轴方向。2.3不规则三角网模型（TIN）TIN的产生原因——规则TINTIN由Peuker和他的同事于1978年设计的一个系统，表示连续面状实体；

——把一系列不规则分布的数据点以直线相连形成一个不规则三角网，网中所有三角形相互邻接，互不交叉、互不重叠。数据点连接为三角网的方法：

Delaunay三角构网法——形成的每一个三角形，它的外接圆不含有除三个顶点以外的其它数据点，外接圆的圆心正是与该三角形三个顶点相对应的Voronoi(Thiessen)多边形的公共顶点。一个Voronoi多边形仅包含一个数据点。不规则三角网（TIN）数据结构由Peuker和他的同事于1978年设计的一个系统，表示连续Delauay三角形和Voronoi多边形Delauay三角形和Voronoi多边形Delaunay三角形产生的基本准则任何一个Delaunay三角形的外接圆的内部不能包含其它任何点最大化最小角原则：每两个相邻的三角形构成的凸四边形的对角线，在相互交换后，六个内角的最小角不再增大局部优化过程LOP（LocalOptimizationProcedure）方法：求出包含新插入点p的外接圆的三角形——影响三角形（InfluenceTriangulation）删除影响三角形的公共边（粗线）将p与全部影响三角形的顶点连接，完成p点在原Delaunay三角形中的插入Delaunay三角形产生的基本准则任何一个Delaunay狄洛尼三角网的构建：（三角网生长法）

设有参考点集Pi（i＝1,2,…,n），首先找出离散点集中相距最短的两点，连线成为D-TIN的初始基线，然后在这附近找第三点。在找第三点的过程中要逐点比较，一般取第三点到前两点的距离平方和最小的参考点作为候选点，以这三点作一外接圆，然后判断周围是否有落入该外接圆的点，如果有，则该三角形不是Delaunay三角形，用周围的其他点作为候选点，重新作外接圆，重新判断周围是否有点落入该外接圆。一直找到没有其他参考点落入外接圆内，该三角形就是Delaunay三角形。以该三角形的任意一条边作为基边，用同样的方法形成其他三角形。一直到所有的参考点都参与构造了Delaunay三角网为止。

狄洛尼三角网的构建：（三角网生长法）设有参考点集Pi（用狄洛尼三角网构建泰森多边形给定一个D-TIN，对于它的所有内边，连接共有每条内边的两个三角形的外接圆的圆心，即构成该TIN的平面点集Voronoi图。

（1）首先构建离散平面点集的D-TIN；

（2）然后求取各三角形的外接圆心；

（3）对每一个离散点，按顺时针或逆时针方向连接与其关联的三角形的外接圆心，即得到该离散点的泰森多边形；

（4）将各离散点的泰森多边形形成集合，即得到本平面点集的泰森多边形。

用狄洛尼三角网构建泰森多边形给定一个D-TIN，对于它的所不规则三角网模型(TIN)

用不规则三角形拟合地表，主要描述地表高程不规则三角网模型(TIN) 用不规则三角形拟合地表，主要描述

将地理空间的表面表示成一系列相连接的三角形。它是在一组结点之中，按照“一定的规则”连接相邻结点形成的边组成的。 将地理空间的表面表示成一系列相连接的三角形。它是在一组结

在TIN模型中，点的位置控制着三角形的顶点，这些三角形尽可能接近等边，地表地形就可由一组三角形很好地表示出来。

三角网的一个优点是：其三角形大小随点密度变化而自动变化，当数据点密集时生成的三角形小，表示地形陡峭，数据点较稀时生成的三角形较大，表示地形平缓。TIN也能表示不连续对象，如悬崖、断层、海岸线和山谷谷底。把TIN转化为栅格，可用线性内插方法，也可用非线内插方法。

在TIN模型中，点的位置控制着三角形的顶点，这些三角形第6章-数字地形模型与地形分析课件第6章-数字地形模型与地形分析课件TIN的数据存储方式存储内容每个顶点的高程、平面坐标，节点连接的拓扑关系，三角形及邻接三角形等关系。存储方式对于每一个三角形、边和节点都对应一个记录，三角形的记录包括三个指向它三个边的记录的指针；边的记录有四个指针字段，包括两个指向相邻三角形记录的指针和它的两个顶点的记录的指针；也可以直接对每个三角形记录其顶点和相邻三角形。每个节点包括三个坐标值的字段，分别存储X，Y，Z坐标。TIN的数据存储方式存储内容三角网存储方式三角网存储方式TIN与规则格网的比较TIN随地形起伏变化的复杂性而改变采样点的密度和决定采样点的位置，能够避免地形平坦时的数据冗余，又能按地形特征点如山脊、山谷线、地形变化线等表示数字高程特征。TIN与规则格网的比较TIN随地形起伏变化的复杂性而改变采样2.4层次模型层次地形模型（LayerofDetails，LOD）用于表达多种不同精度水平的数字高程模型。层次地形模型允许根据不同的任务要求选择不同精度的地形模型。2.4层次模型层次地形模型（LayerofDetai3.DEM模型之间的相互转换为满足不同需要常需要在DEM和TIN之间进行转换。3.DEM模型之间的相互转换为满足不同需要常需要在DEM和3.1等高线转成格网DEM局部插值算法，估计未知格网点的高程要在一个半径范围内搜索落在其中的已知点数据，再计算它的加权平均值。面临问题：如果搜索到的点都具有相同的高程，那待插值点的高程也同为此高程值。结果导致在每条等高线周围的狭长区域内具有与等高线相同的高程，出现“阶梯”地形。当低海拔平原地区等高线距离更远时，搜索到一条等高线上的数据的可能性就越大，问题更严重。以带“阶梯”地形的DEM为基础，计算坡度往往会出现不自然的条斑状分布模式。3.1等高线转成格网DEM局部插值算法，估计未知格网点的3.2利用格网DEM提取等高线方法：根据格网DEM中相邻四个点组成四边形进行等高线跟踪。在格网DEM提取等高线中，除了划分为三角形之外，也可以直接使用四边形跟踪等高线。但是在图9-11所示的情形中，仍会出现等高线跟踪的二义性，即对于每个四边形，有两条等高线的离去边。进行取舍判断的方法一般是计算距离，距离近的连线方式优于距离远的连线方式。在图9-11种，就要采用（b）图所示的跟踪方式。格网DEM提取等高线另一个值得注意的问题是，如果一些网格点的数值恰好等于要提取的等高线的数值，会使判断过程变得复杂，并且会生成不闭合的等高线，一般的解决办法是将这些网格点的数值增加一个小的偏移量。3.2利用格网DEM提取等高线方法：根据格网DEM中相邻3.2格网DEM转成TIN可以看作是一种规则分布的采样点生成TIN。目的：尽量减少TIN的顶点数目，并尽可能多地保留地形信息，如山峰、山脊、谷底和坡度突变处。算法实现的两个重要的特征：1）筛选要保留或丢弃的格网点；2）判断停止筛选的条件。保留重要点法启发丢弃法。3.2格网DEM转成TIN可以看作是一种规则分布的采样点生3.2.1保留重要点法通过比较计算格网点的重要性，保留重要的格网点。重要点（VIP，VeryImportantPoint）通过3×3的模板来确定，根据八邻点的高程值决定模板中心是否为重要点。比较格网点的高程值与8邻点的高程值，将差分超过某个阈值的格网点保留下来。作为三角网顶点生成Delaunay三角网3.2.1保留重要点法通过比较计算格网点的重要性，保留重要3.2.2启发丢弃法（DH，DropHeuristic）算法：给定一个格网DEM和转换后TIN中节点的数量限制，寻求一个TIN与规则格网DEM的最佳拟合。首先输入整个格网DEM，迭代进行计算，逐渐将那些不太重要的点删除，处理过程直到满足数量限制条件或满足一定精度为止。1）输入TIN，每次去掉一个节点进行迭代，得到节点越来越少的TIN。2）取TIN的一个节点O及与其相邻的其它节点，将O的邻点进行Delaunay三角形重构。3）判断节点O位于哪个新生成的Delaunay三角形中（图为BCE）。计算O点的高程和过O点与三角形BCE交点O’的高程差d。若高程差d大于阈值de，则O点为重要点，保留，否则，删除（de为阈值）。4）对TIN中所有的节点，重复进行上述判断过程。5）直到TIN中所有的节点满足条件d>de，结束。3.2.2启发丢弃法（DH，DropHeuristic）保留重要点法(VIP)和启发丢弃法(DH)比较VIP方法在保留关键网格点方面（顶点、凹点）最好；DH方法在每次丢弃数据点时确保关键信息丢失最少，但要求计算量大。两种方法各有利弊，实际应用中根据不同的需要，可以选择使用不同的方法。检测极值点高效存储最小误差保留重要点法(VIP)和启发丢弃法(DH)比较VIP方法在保3.3TIN转成格网DEM等同于不规则点生成格网DEM的过程。按要求的分辨率大小和方向生成规则格网对每一个格网搜索最近的TIN数据点，按线性或非线性插值函数计算格网点高程3.3TIN转成格网DEM等同于不规则点生成格网DEM的4.DEM模型在GIS中的应用1、由TIN获取任意点P的高程已知Q1、Q2、Q3三个顶点，高程线性内插求P的Zx y z 1x1 y1 z1 1x2 y2 z2 1x3 y3 z3 1 =0平面方程Z=Z1-（x-x1)(y21z31-y31z21)+(y-y1)(z21x31-z31x21)/(x21y31-x31y21)4.DEM模型在GIS中的应用1、由TIN获取任意点P的高2、由TIN进行曲面拟合由于在TIN中可以获取任意点高程,进行密集插值可获得拟合的地形曲面.2、由TIN进行曲面拟合