数字高程模型及其应用.ppt

《数字高程模型及其应用.ppt》由会员分享，可在线阅读，更多相关《数字高程模型及其应用.ppt（45页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第六章第六章数字高程模型数字高程模型及其应用及其应用DTM与与DEM的概念的概念 数字地面模型（Digital Terrain Model,DTM)是定义于二维区域上的一个有限项的向量序列，它以离散分布的平面点来模拟连续分布的地形。数学的角度 Kpfk(up,vp)(kl，2，3，m;p1，2，3，n)或：DTM=Z i,j,i=1,2,3,m；j=1,2,3,n DTM是地形表面形态属性信息的数字表达，是带有空间位置特征和地形属性特征的数字描述。如地面温度、降雨、地球磁力、重力、土地利用、土壤类型等其他地面诸特征。数字地形模型中地形属性为高程时称为数字高程模型（DEM，Digital Ele

2、vation Model）。高程是地理空间中的第三维坐标。数学表达为：数学表达为：z=fz=f（x x，y y）DEM是DTM的一个子集，是DTM的基础数据，最核心部分，可以从中提取出各种地形信息，如高度、坡度、坡向、粗糙度，并进行通视分析，流域结构生成等应用分析。数字高程模型的用途数字高程模型的用途 谷脊特征分析 淹没边界的计算 剖面图的自动绘制 坡度、坡向分析 地表粗糙度计算 地表形态的自动分类 地形曲面拟合 立体透视图 通视分析 地貌晕渲图及其与专题地图叠置 交通线路选择 工程土方量估算 项目选址 土地利用规划数字地面模型数据采集DTM的数据采集的数据采集数据源决定采集方法：（1）航空或

3、航天遥感图像为数据源（2）以地形图为数据源（3）以地面实测记录为数据源（4）其它数据源 DEM表示方法表示方法1、数学方法 用数学方法来表达，可以采用整体拟合方法，即根据区域所有的高程点数据，用傅立叶级数和高次多项式拟合统一的地面高程曲面。也可用局部拟合方法，将地表复杂表面分成正方形规则区域或面积大致相等的不规则区域进行分块搜索，根据有限个点进行拟合形成高程曲面。2.图形方法1）线模式 等高线是表示地形最常见的形式。其它的地形特征线也是表达地面高程的重要信息源，如山脊线、谷底线、海岸线及坡度变换线等。等高等高线线通常被存通常被存储储成一个成一个有序有序的的坐坐标标点序列点序列，可以，可以认为认

4、为是一条是一条带带有高程有高程值值属性的属性的简单简单多多边边形形或多或多边边形弧段。形弧段。由于等高由于等高线线模型只是表达了区域的模型只是表达了区域的部部分高程分高程值值，往往需要一种，往往需要一种插插值值方法方法来来计计算落在算落在等高等高线线以外以外的的其他点其他点的高程，的高程，又因又因为这为这些点是落在两条等高些点是落在两条等高线线包包围围的区域内，所以，通常只要使用的区域内，所以，通常只要使用外包外包的的两条等高两条等高线线的高程的高程进进行行插插值值。等高线通常被存成一个有序的坐标点对序列，可以认为是简单多边形或多边形弧段。由于等高线模型只表达了区域的部分高程值，往往需要一种插

5、值方法来计算落在等高线外的其它点的高程，又因为这些点是落在两条等高线包围的区域内，所以，通常只使用外包的两条等高线的高程进行插值。等高线通常可以用二维的链表来存储。另外的一种方法是用图来表示等高线的拓扑关系，将等高线之间的区域表示成图的节点，用边表示等高线本身。此方法满足等高线闭合或与边界闭合、等高线互不相交两条拓扑约束。这类图可以改造成一种无圈的自由树。下图为一个等高线图和它相应的自由树。等高线和相应的自由树 2.图形方法2）点模式 用离散采样数据点建立DEM是DEM建立常用的方法之一。数据采样可以按规则格网采样，可以是密度一致的或不一致的；可以是不规则采样，如不规则三角网、邻近网模型等；也

6、可以有选择性地采样，采集山峰、洼坑、隘口、边界等重要特征点。3）规则格网模型规则格网模型规则网格，通常是正方形，也可以是矩形、三角形等规则网格。规则网格将区域空间切分为规则的格网单元，每个格网单元对应一个数值。数学上可以表示为一个矩阵，在计算机实现中则是一个二维数组。每个格网单元或数组的一个元素，对应一个高程值，如图所示。规则规则格网法是把格网法是把DEMDEM表示成表示成高程矩高程矩阵阵，此，此时时，DEMDEM来源于直接来源于直接规则规则矩形格网采矩形格网采样样点或由不点或由不规则规则离散数据点内插离散数据点内插产产生。生。结构简单结构简单，计计算机算机对对矩矩阵阵的的处处理理比比较较方便

7、方便，高程矩，高程矩阵阵已成已成为为DEMDEM最通用最通用的形式。高程矩的形式。高程矩阵阵特特别别有利于有利于各种各种应应用用。但但GridGrid系系统统也有下列也有下列缺点缺点：a)a)地形地形简单简单的地区存在的地区存在大量冗余数大量冗余数据据；b)b)如不改如不改变变格网大小格网大小,则则无法适用无法适用于起伏程度不同的地区；于起伏程度不同的地区；c)c)对对于某些特殊于某些特殊计计算如算如视线视线计计算算时时，格网的，格网的轴线轴线方向被夸大方向被夸大；d)d)由于由于栅栅格格过过于粗略，于粗略，不能精确不能精确表示地形的关表示地形的关键键特征特征,如山峰、如山峰、洼坑、山脊等；洼

8、坑、山脊等；对于每个格网的数值有两种不同的解释。第一种是格网栅格观点，认为该格网单元的数值是其中所有点的高程值，即格网单元对应的地面面积内高程是均一的高度，这种数字高程模型是一个不连续的函数。第二种是点栅格观点，认为该网格单元的数值是网格中心点的高程或该网格单元的平均高程值，这样就需要用一种插值方法来计算每个点的高程。4）不规则三角网）不规则三角网(TIN)模型模型 不规则三角网（Triangulated Irregular Network,TIN）是另外一种表示数字高程模型的方法Peuker等，1978，它既减少规则格网方法带来的数据冗余，同时在计算（如坡度）效率方面又优于纯粹基于等高线的方

9、法。TIN模型根据区域有限个点集将区域划分为相连的三角面网络，区域中任意点落在三角面的顶点、边上或三角形内。如果点不在顶点上，该点的高程值通常通过线性插值的方法得到（在边上用边的两个顶点的高程，在三角形内则用三个顶点的高程）。所以TIN是一个三维空间的分段线性模型，在整个区域内连续但不可微。4）不规则三角网）不规则三角网(TIN)模型模型 不规则三角网(TIN)表示法克服了高程矩阵中冗余数据的问题，而且能更加有效地用于各类以DTM为基础的计算。但其结构复杂。TIN的数据存储方式比格网DEM复杂，它不仅要存储每个点的高程，还要存储其平面坐标、节点连接的拓扑关系，三角形及邻接三角形等关系。TIN模

10、型在概念上类似于多边形网络的矢量拓扑结构，只是TIN模型不需要定义“岛”和“洞”的拓扑关系。TIN存储方式：存储方式：129758103641 17 76 63 32 24 45 59 98 811111010TIN存储方式：存储方式：6388379:817369554433221NoNoP1P2P3112321343451:：11678坐标与高程值表三角形表NoXYZP190.010.043.51250.710.067.35367.223.962.68:1010.090.081.036 不规则三角网数字高程由连续的三角面组成，三角面的形状和大小取决于不规则分布的测点，或节点的位置和密度。不规

11、则三角网与高程矩阵方法不同之处是随地形起伏变化的复杂性而改变采样点的密度和决定采样点的位置，因而它能够避免地形平坦时的数据冗余，又能按地形特征点如山脊、山谷线、地形变化线等表示数字高程特征。DEM DEM 特点特点1 1）容容易易以以多多种种形形式式显显示示地地形形信信息息。地地形形数数据据经经过过计计算算机机软软件件处处理理过过后后，产产生生多多种种比比例例尺尺的的地地形形图图、纵纵横横断断面面图图和和立立体体图图。而而常常规规地地形形图图一一经经制制作作完完成成后后，比比例例尺尺不容易不容易改变或需要人工处理。改变或需要人工处理。2 2）精精度度不不会会损损失失。常常规规地地图图随随着着时

12、时间间的的推推移移，图图纸纸将将会会变变形形，失失掉掉原原有有的的精精度度。而而DEMDEM采采用用数数字字媒媒介介，因因而而能能保保持持精精度度不不变变。另另外外，由由常常规规的的地地图图用用人人工工的的方方法法制作其他种类的地图，精度会制作其他种类的地图，精度会受到损失受到损失，而由，而由DEMDEM直接输出，精度可直接输出，精度可得到控制得到控制。3 3）容容易易实实现现自自动动化化、实实时时化化。常常规规地地图图要要增增加加和和修修改改都都必必须须重重复复相相同同的的工工序序，劳劳动动强强度度大大而而且且周周期期长长，而而DEMDEM由由于于是是数数字字形形式式的的，所所以以增增加加和

13、和修修改改地地形形信信息息只只需需将修改信息直接将修改信息直接输入计算机输入计算机，经软件经软件处理后即可得各种地形图。处理后即可得各种地形图。与与传统传统地形地形图图比比较较，DEMDEM作作为为地形表面的一种数字表达形式有如下特点：地形表面的一种数字表达形式有如下特点：DEMDEM应用应用 1 1）作为国家地理信息的基础数据；）作为国家地理信息的基础数据；2 2）土木工程、景观建筑与矿山工程规划与设计；）土木工程、景观建筑与矿山工程规划与设计；3 3）为军事目的而进行的三维显示；）为军事目的而进行的三维显示；4 4）景观设计与城市规划；）景观设计与城市规划；5 5）流水线分析、可视性分析；

14、）流水线分析、可视性分析；6 6）交通路线的规划与大坝选址；）交通路线的规划与大坝选址；7 7）不同地表的统计分析与比较；）不同地表的统计分析与比较；8 8）生成坡度图、坡向图、剖面图、辅助地貌分析、估计侵蚀和径流等；）生成坡度图、坡向图、剖面图、辅助地貌分析、估计侵蚀和径流等；9 9）作作为为背背景景叠叠加加各各种种专专题题信信息息如如土土壤壤、土土地地利利用用及及植植被被覆覆盖盖数数据据等等，以以进进行行显显示示与与分析；分析；10 10）与）与GISGIS联合进行空间分析；联合进行空间分析；11 11）虚拟现实）虚拟现实(Virtual Reality)(Virtual Reality)

15、；此外，从此外，从DEMDEM还还能派生以下主要能派生以下主要产产品：平面等高品：平面等高线图线图、立体等高、立体等高线图线图、等坡度、等坡度图图、晕晕渲渲图图、通、通视图视图、纵纵横断面横断面图图、三、三维维立体透立体透视图视图、三、三维维立体彩色立体彩色图图等。等。返回返回DEMDEM建立建立1 1、数据获取与处理、数据获取与处理沿断面采集沿断面采集选点采集选点采集1 1）数据采集数据采集 为了建立DEM，必需量测一些点的三维坐标，这就是DEM数据采集。DEM数据采集方法数据采集方法1）地面测量）地面测量利用自动记录的测距经纬仪（常用电子速测经纬仪或全站经纬仪）在野外实测。这种速测经纬仪一

16、般都有微处理器，可以自动记录和显示有关数据，还能进行多种测站上的计算工作。其记录的数据可以通过串行通讯，输入计算机中进行处理。2）现有地图数字化）现有地图数字化利用数字化仪对已有地图上的信息（如等高线）进行数字化的方法，目前常用的数字化仪有手扶跟踪数字化仪和扫描数字化仪。3）空间传感器）空间传感器利用全球定位系统GPS，结合雷达和激光测高仪等进行数据采集。4）数字摄影测量方法）数字摄影测量方法这是DEM数据采集最常用的方法之一。利用附有的自动记录装置（接口）的立体测图仪或立体坐标仪、解析测图仪及数字摄影测量系统，进行人工、半自动或全自动的量测来获取数据。DEM DEM 生成生成1 1）人工网格

17、法）人工网格法 在地形图上蒙上格网，在地形图上蒙上格网，逐格逐格读取读取中心点中心点或或交点交点的高程值。的高程值。对有限个离散点，对有限个离散点，每三每三个个邻近点邻近点联结成联结成三角形，三角形，每个三角形代表一个每个三角形代表一个局部局部平面平面，再根据每个平面方，再根据每个平面方程，可计算程，可计算各格网点各格网点高程，高程，生成生成DEMDEM。2 2）三角网法）三角网法3 3）立体像对法）立体像对法 资料来源于张超主编的地理信息系统教程所配光盘资料来源于张超主编的地理信息系统教程所配光盘4 4）曲面拟合法）曲面拟合法 根据根据有限个有限个离散点的高程，采用离散点的高程，采用多项式多

18、项式或或样条函数样条函数求得求得拟合公式拟合公式，再逐个计算各点的高程，得到拟合的再逐个计算各点的高程，得到拟合的DEMDEM。可反映。可反映总的地势总的地势，但，但局部局部误误差差较大较大。可分为：可分为：整体拟合整体拟合：根据研究区域内：根据研究区域内所有采样点所有采样点的观测值建立的观测值建立趋势面模型趋势面模型。特。特点是不能反映内插区域内的局部特征。点是不能反映内插区域内的局部特征。局部拟合局部拟合：利用：利用邻近的邻近的数据点数据点估计估计未知点的值，能未知点的值，能反映局部反映局部特征。特征。5 5）等值线插值法）等值线插值法 三、三、DEM DEM 的应用的应用（一）（一）基于

19、基于DEMDEM的信息提取的信息提取1 1、坡度、坡度定义为地表单元的法向与定义为地表单元的法向与Z Z轴的夹角，即切平面与水平面的夹角。轴的夹角，即切平面与水平面的夹角。在在计计算算出出各各地地表表单单元元的的坡坡度度后后，可可对对不不同同的的坡坡度度设设定定不不同同的的灰灰度度级级，可得到坡度图。可得到坡度图。2 2、坡向、坡向坡向是地表单元的法向量在水平面上的投影与坡向是地表单元的法向量在水平面上的投影与X X轴之间的夹角，轴之间的夹角，在在计计算算出出每每个个地地表表单单元元的的坡坡向向后后，可可制制作作坡坡向向图图，通通常常把把坡坡向向分分为为东东、南南、西西、北北、东东北北、西西北

20、北、东东南南、西西南南8 8类类，再再加加上上平平地地，共共9 9类类，用用不同的色彩显示，即可得到坡向图。不同的色彩显示，即可得到坡向图。（Grid DEMGrid DEM上制作坡度、坡向图上制作坡度、坡向图）3 3、地表粗造度（破碎度）、地表粗造度（破碎度）是是反反映映地地表表的的起起伏伏变变化化和和侵侵蚀蚀程程度度的的指指标标，一一般般定定义义为为地地表表单单元元的的曲曲面面积与其水平面上的投影面积之比。面面积与其水平面上的投影面积之比。坡度坡向的计算坡度坡向的计算又：又：所以所以(0,0)(1,1)(1,0)Z10Z11Z01POQRSGrid DEMGrid DEM上制上制作坡度、坡

21、向作坡度、坡向图图 通常用通常用3*33*3的格的格网窗口在网窗口在DEMDEM数据数据矩阵中连续移动矩阵中连续移动计算完成。计算完成。数字地面模型应用数字地面模型应用坡度坡向的计算坡度坡向的计算4 4、高程变异分析、高程变异分析平地平地岗丘岗丘丘陵丘陵低山低山高山高山绝对绝对h h400600400600相对相对h h100100100-200100-200坡向坡向33DEMDEM计算地表形态要素计算地表形态要素H,H,H,H,坡度、坡向等坡度、坡向等地形分类标准表地形分类标准表地形自动分类地形自动分类5 5、地貌形态的自动分类地貌形态的自动分类包括平均高程、相对高程、高程标准差，高程变异。

22、包括平均高程、相对高程、高程标准差，高程变异。高程变异：为格网顶点的高程变异：为格网顶点的高程标准差高程标准差与与平均高程平均高程的的比值比值。（二）等高线的绘制（二）等高线的绘制 在格网在格网DEMDEM上上自动绘制自动绘制等高线主要包括两个步骤：等高线主要包括两个步骤：1 1、等高线追踪、等高线追踪，利用，利用DEMDEM矩形格网点的高程内插出矩形格网点的高程内插出格网边上格网边上的等高的等高线点，并将这些等高线点线点，并将这些等高线点排序排序；2 2、等高线光滑、等高线光滑，进一步，进一步加密加密等高线点并绘制光滑曲线。等高线点并绘制光滑曲线。数字地面模型应用数字地面模型应用等高线绘制等

23、高线绘制（三）（三）基于基于DEMDEM的可视化分析的可视化分析1 1、剖面分析、剖面分析1 1）意义）意义：可以：可以以线代面以线代面，研究区域的地貌形态、轮廓形状、地势变化、地质，研究区域的地貌形态、轮廓形状、地势变化、地质构造、斜坡特征、地表切割强度等。构造、斜坡特征、地表切割强度等。如果在地形剖面上如果在地形剖面上叠加叠加其它地理变量，例如坡度、土壤、植被、土地利用其它地理变量，例如坡度、土壤、植被、土地利用现状等，可以提供土地利用规划、工程选线和选址等的决策依据。现状等，可以提供土地利用规划、工程选线和选址等的决策依据。2 2）绘制）绘制 可在可在格网格网DEMDEM或或三角网三角网

24、DEMDEM上进行。上进行。已知两点的坐标已知两点的坐标A(xA(x1 1，y y1 1)，B(xB(x2 2，y y2 2)，则可求出，则可求出两点两点连线连线与与格网格网或或三角网三角网的的交点交点，并并内插内插交点上的高程，以及交点上的高程，以及各交点之间的距离。然后按各交点之间的距离。然后按选定选定的垂直比例尺和水平比的垂直比例尺和水平比例尺，按例尺，按距离距离和和高程高程绘出剖绘出剖面图。面图。剖面图剖面图不一定不一定必须沿直线必须沿直线绘制，也可沿一条绘制，也可沿一条曲线绘制曲线绘制。2 2、通视分析通视分析通视分析是指以通视分析是指以某一点某一点为观察点，研究为观察点，研究某一区

25、域某一区域通视情况通视情况的地形的地形分析。分析。1）方法）方法：a a、以以O O为为观观察察点点，对对格格网网DEMDEM或或三三角角网网DEMDEM上上的的每每个个点点判判断断通通视视与与否否，通通视视赋赋值值为为1 1，不不通通视视赋赋值值为为0 0。由由此此可可形形成成属属性性值值为为0 0和和1 1的的格格网网或或三三角网。对此以角网。对此以0.50.5为值追踪等值线，即得到以为值追踪等值线，即得到以O O为观察点的通视图。为观察点的通视图。b b、以以观察点观察点O O为轴，以一定的为轴，以一定的方位角间隔方位角间隔算出算出00360360的的所有方所有方位线上位线上的通视情况。

26、对于每条方位线，通视的地方绘线，不通视的地的通视情况。对于每条方位线，通视的地方绘线，不通视的地方断开，或相反。这样可得出射线状的通视图。方断开，或相反。这样可得出射线状的通视图。2）关键算法）关键算法均是判断格网或三角网上的某一点是否通视均是判断格网或三角网上的某一点是否通视（即（即两点是否可见两点是否可见）。）。(两点是否可见的算法两点是否可见的算法)两点是否可见的算法：两点是否可见的算法：a a）倾角法）倾角法格网格网DEMDEM为例，为例，O(xO(xo o，y yo o，z zo o)为为观察点观察点，P(xP(xp p,y,yp p,z,zp p)为某一为某一格网格网点点，OPOP

27、与格网的与格网的交点交点为为A A、B B、C C。OP OP的的倾角倾角为为 观察点与观察点与各交点各交点的的倾角倾角为为 i i(i(iA,B,C)A,B,C)若若tgtgmax(tgmax(tgi i,i,iA A、B B、C)C)，则则OPOP通视通视,否则，不通视。否则，不通视。b)b)剖面图剖面图 两点连线是否与剖面两点连线是否与剖面相交。相交。A AB BA A返回返回3）通视分析示例）通视分析示例观察点观察点不通视不通视通视通视雷雷达达盲盲区区飞飞行行 可可视视域域的的三三维维显显示示资料来源于张超主编的地理信息系统教程所配光盘资料来源于张超主编的地理信息系统教程所配光盘3 3

28、、地形三维图绘制、地形三维图绘制 DEMDEM高程点高程点建立几建立几何模型何模型透视透视变换变换隐藏线、隐藏线、面的消除面的消除光照光照模型模型计算计算贴贴纹纹理理图形图形输出输出开放式图形标准库（开放式图形标准库（OpenGLOpenGL）。）。部分图片摘自张超主部分图片摘自张超主编的地理信息系统编的地理信息系统教程所配光盘教程所配光盘4 4、地貌晕渲图绘制、地貌晕渲图绘制资料来源于张超主编的地理信息系统教程所配光盘资料来源于张超主编的地理信息系统教程所配光盘数字地面模型应用数字地面模型应用晕渲图晕渲图的制作的制作5 5、模拟飞行、模拟飞行资料来源于资料来源于张超主编的张超主编的地理信息地理信息系统教程系统教程所配光盘所配光盘数字地面模型应用数字地面模型应用透视图的绘制透视图的绘制（四）流域水文特征及土木工程（四）流域水文特征及土木工程 水淹示例水淹示例 三维规划设计三维规划设计 用于工程项目中的开挖填方、线路勘用于工程项目中的开挖填方、线路勘测设计、水利建设工程等。测设计、水利建设工程等。资料来源于资料来源于张超主编的张超主编的地理信息地理信息系统教程系统教程所配光盘所配光盘作业作业说明DEM的概念及建立方法。说明坡度，坡向，坡长的概念与基于DEM的提取方法。