第六章 数字地面模型

地理信息系统原理及应用GIS第六章数字高程模型（DEM）学习目标掌握数字地面模型及数字高程模型的概念，数字地形分析与应用。重点：数字地形分析和应用。难点：空间数据的内插，数字高程模型应用。地理信息系统原理及应用GIS空间数据库返回地形信息地理信息系统原理及应用GIS第一节概述（一）地形表达的方法（二）数字高程模型的概念（三）DEM与传统地形图比较（四）DEM的数学表达（五）DEM数据的分布特征(六)DEM数据获取和采样方法地理信息系统原理及应用GIS（一）地形表达的方法1.绘图2.地图3.摄影4.摄像5.三维图6.DEM与三维表达地理信息系统原理及应用GIS什么是DEM？——DigitalElevationModel地理信息系统原理及应用GIS（二）数字高程模型的概念数字高程模型是新一代的地形图，地貌和地物不再用直观的等高线和图例符号在纸上表达，而是通过储存在磁性介质上的大量密集的地面点的空间坐标和地形属性编码，以数字的形式描述。DEM以数字的形式按一定的组织结构组织在一起，表示实体地形特征空间分布的模型，是地形形状大小和起伏的数字描述。是国家基础空间数据的重要组成部分。地理信息系统原理及应用GIS从研究对象与应用范畴角度出发，DEM可以归纳为狭义和广义两种定义。•从狭义角度定义，DEM是区域表面海拔高程的数字化表达。•从广义角度定义，DEM是地理空间中地理对象表面海拔高度的数字化表达。这是随着DEM的应用不断向海底、地下岩层以及某些不可见的地理现象（如空中的等气压面等）延伸，而提出的更广义的概念。地理信息系统原理及应用GIS（三）DEM与DTM:数字地面模型就是对某一种或多种地面特性空间分布的数字描述，是叠加在二维地理空间上的一维或多维地面特性向量空间，是是对二维地理空间上具有连续变化特征地理现象的模型化表达和过程模拟。数字地面模型的本质共性是二维地理空间定位和数字描述。地理信息系统原理及应用GISDEM是DTM的一个分支。实际上，DEM是DTM中最基本的部分，它通过有限的地形高程数据实现对地形曲面的数字化模拟，即地形表面形态的数字化表示。DTM的另外两个分支是各种非地貌特性的以矩阵形式表示的数字模型，包括自然地理要素以及与地面有关的社会经济及人文要素，如土壤类型、土地利用类型、岩层深度、地价、商业优势区等等。地理信息系统原理及应用GIS（四）DEM与传统地形图比较地球表面的的高低起伏变化是一种连续变化的曲面，是无法用平面地图来确切表示的。1）DEM容易以多种形式显示地形信息，地形数据经计算机处理后能产生比例尺、纵横断面图与立体图，而常规地图一旦制作形成，比例尺不容易改变，在绘制其他的地形图需要人工处理；2）DEM精度不会损失，没有载体变形的问题；3）DEM形象逼真。地理信息系统原理及应用GIS（五）DEM的数学表达DEM的核心是地形表面特征点的三维坐标数据和一套对地表提供连续描述的算法最基本的DEM由一系列地面x,y位置及其相联系的高程z所组成，它表示地表区域上地形的三维向量的有限序列，即地表单元上高程的集合。数学表达:z=f(x,y)，x,y属于DEM所在区域。地理信息系统原理及应用GIS（六）DEM数据的分布特征按其空间分布特征分为：格网状数据和离散数据。把DEM覆盖区划分为规则格网，每个网格大小和形状都相同，用相应矩阵元素的行列号来实现网格点的二维地理空间定位，第三维为特征值，可以是高程和属性。网格大小代表数据精度。地理信息系统原理及应用GIS不可能用规格网获取数据时，则获取离散数据，离散数据DEM的平面二维地理空间定位由不规则分布的离散样点平面坐标实现，第三维仍为高程或属性特征值。如气象、水文与其他地理抽样条调查等呈不规则分布，需取离散数据。地理信息系统原理及应用GIS(七)DEM数据采样方法数据源决定采集方法:⑴摄影测量数据采集方法(以航空或航天遥感影像为数据源)⑵从现有地图获取数据(以地形图为数据源)⑶地面实测记录(地面实测记录为数据源)(4)基于不规则三角网的方法(5)其他地理信息系统原理及应用GIS第二节DEM的表示方法及其转换一表示方法1.数学方法2.图形法地理信息系统原理及应用GIS1.数学方法用数学方法拟合表面时，在整体表示时，需依靠连续三维函数，连续的三维函数能以高平滑度表示复杂表面，一般用到傅里叶级数或高次多项式。表示局部使用规则块或不规则块，是将地表分为正方形像元，或面积大致相等的不规则形状的小块。据部分块模拟广泛应用于复杂表面模拟的机助设计系统，现在地下水、土壤特征或其他环境数据的表面内插。地理信息系统原理及应用GIS2.图形法分点模式与线模式，点模式分规则与不规则，规则两类，处理典型特征如山峰，洼地，边界等。线数据可处理水平线、垂直线与典型线，如山脊线、谷底线、海岸线坡度变化线等。⑴点模式：DEM的最普通形式是高程矩阵或规则矩形网格，高程数据直接由解析立体测量仪从立体航空相片上定量测量。⑵线模式：线模式是一系列描述高程测量曲线的等高线。地理信息系统原理及应用GIS1）等高线法等高线通常被存储成一个有序的坐标点序列，可以认为是一条带有高程值属性的简单多边形或多边形弧段。地理信息系统原理及应用GIS2）TIN法TIN（TriangulatedIrregularNetwork)表示法利用所有采样点取得的离散数据，按照优化组合的原则，把这些离散点（各三角形的顶点）连接成相互连续的三角面（在连接时，尽可能地确保每个三角形都是锐角三角形或是三边的长度近似相等--Delaunay）。因为TIN可根据地形的复杂程度来确定采样点的密度和位置，从而减少了地形较平坦地区的数据冗余。但存储方式较DEM复杂。地理信息系统原理及应用GIS邻接三角形1XYZ2XYZ3XYZ4XYZ5XYZ6XYZ7XYZ8XYZ12345678123456781234顶点5687542365746544888721X31264534XX587X62867XX点文件三角形文件11125443TIN存储方式地理信息系统原理及应用GIS3）规则格网法(Grid)规则格网法是把DEM表示成高程矩阵，此时，DEM来源于直接规则矩形格网采样点或由不规则离散数据点内插产生。结构简单，计算机对矩阵的处理比较方便，高程矩阵已成为DEM最通用的形式。高程矩阵特别有利于各种应用。但Grid系统也有下列缺点：a)地形简单的地区存在大量冗余数据；b)如不改变格网大小,则无法适用于起伏程度不同的地区；c)对于某些特殊计算如视线计算时，格网的轴线方向被夸大；d)由于栅格过于粗略，不能精确表示地形的关键特征,如山峰、洼坑、山脊等；地理信息系统原理及应用GIS二DEM建立的一般步骤构建DEM的整体思路是首先在二维平面上对研究区域进行格网划分（格网大小取决于DEM的应用目的），形成覆盖整个区域的格网空间结构，然后利用分布在格网点周围的地形采样点内插计算格网点的高程值，最后按一定的格式输出，形成该地区的格网DEM。地理信息系统原理及应用GISx不规则分布点规则分布等高线分布Y对每一格网点求取格网点高程格网DEM建立流程xXyY地理信息系统原理及应用GIS三DEM模型之间的相互转换在实际应用中，DEM模型之间可以相互转换。1不规则点集生成TIN对于不规则分布的高程点，可以形式化地描述为平面的一个无序的点集P，点集中每个点p对应于它的高程值。将该点集转成TIN，最常用的方法是Delaunay三角剖分方法。生成TIN的关键是Delaunay三角网的产生算法。地理信息系统原理及应用GIS对于给定的初始点集P，有多种三角网剖分方式，而Delaunay三角网有以下特性：（1）其Delaunay三角网是惟一的；（2）三角网的外边界构成点集P的凸多边形“外壳”（或群点的凸壳）；（3）外接圆规则：没有任何点在三角形的外接圆内部（即三角形的外接圆中没有其他点），反之，如果一个三角网满足此条件，那么它就是Delaunay三角网。地理信息系统原理及应用GIS（4）如果将三角网中的每个三角形的最小角进行升序排列，则Delaunay三角网的排列得到的数值最大，从这个意义上讲，Delaunay三角网是“最接近于规则化”的三角网。即三角形最大限度地保持均衡，避免狭长三角形出现————最大最小角规则）。（5）遵守平面图形的欧拉定理：（面＋点）的个数减去边的个数等于2地理信息系统原理及应用GIS2TIN转成等高线可以直接利用原始观测数据，避免了DEM内插的精度损失，因而等高线精度较高；对高程注记点附近的较短封闭等高线也能绘制；绘制的等高线分布在采样区域内而并不要求采样区域有规则四边形边界。而同一高程的等高线只穿过一个三角形最多一次，因而程序设计也较简单。地理信息系统原理及应用GISTIN直接由不规则数据点连成三角形网。重要点赛选－最优化组合－构网。3等高线转成TIN地理信息系统原理及应用GIS4格网DEM转成TIN是一种规则分布的采样点生成TIN的特例，其目的是尽量减少TIN的顶点数目，同时尽可能多地保留地形信息，如山峰、山脊、谷底和坡度突变处。绝大多数算法都有两个重要的特征：（1）筛选要保留或丢弃的格网点；（2）判断停止筛选的条件。其中两个代表性的方法算法是保留重要点法和启发丢弃法。地理信息系统原理及应用GIS5TIN转成格网DEMTIN转成格网DEM可以看作普通的不规则点生成格网DEM的过程。方法是按要求的分辨率大小和方向生成规则格网，对每一个格网搜索最近的TIN数据点，按线性或非线性插值函数计算格网点高程。地理信息系统原理及应用GIS表示地形的最常见的线模式是一系列描述高程曲线的等高线。由于现有地图大多数都绘有等高线，这些地图便是数字高程模型的现成数据源，可以将纸面等高线图扫描后，自动获取DEM数据。由于数字化的等高线不适合于计算坡度或制作地貌渲染图等地形分析，因此，必须要把数字化等高线转为格网高程矩阵。6等高线转成格网DEM地理信息系统原理及应用GIS使用局部插值算法，如距离倒数加权平均或克里金插值算法，可以将数字化等高线数据转为规则格网的DEM数据，但插值的结果往往会出现一些许多不令人满意的结果，而且数字化等高线时越小心，采样点越多，问题越严重。使用针对等高线插值的专用方法，把等高线数据点减少到最少，增加标识山峰、山脊、谷底和坡度突变的数据点，同时使用一个较大的搜索窗口。地理信息系统原理及应用GIS需要将其中的每个点视为一个几何点，而不是一个矩形区域，根据格网DEM中相邻四个点组成四边形进行等高线跟踪。其方法类似于利用TIN提取等高线。实际上，也可以将每个矩形分割成为两个三角形，并应用TIN提取等高线算法。7利用格网DEM生成等高线地理信息系统原理及应用GIS第三节空间数据的内插DEM空间内插的概念十分简单，即在一个由x、y坐标平面构成的二维空间中，由已知若干离散点Pi的高程，估算待内插点的高程值。DEM内插按插点分布范围，可分为整体内插、分块内插和逐点内插三类。地理信息系统原理及应用GISDEM内插方法分类地理信息系统原理及应用GIS1整体内插某种地理属性在空间的连续变化，可以用一个平滑的数学平面加以描述。其数学模型多用二元高次多项式来拟合。地理信息系统原理及应用GIS2分块内插把需要建立DTM的地区，切割成一定大小的规则方块，它的尺寸应根据地形复杂程度和数据源的比例尺确定。在每一个分块上展铺一张数学面，一般要求相邻分块之间有适当宽度的重叠带，以使重叠带内全部数据点成为相邻块展铺数学面时的共用数据，保证一张数学面能够较平滑地与相邻分块的数学面拼接，这在表达地形变化特征的数字高程模型(DEM)内插中应用尤为广泛。地理信息系统原理及应用GIS使用局部插值方法需要注意的几个方面是：所使用的插值函数；邻域的大小、形状和方向；数据点的个数；数据点的分布方式是规则的还是不规则的等。地理信息系统原理及应用GIS3逐点内插以待插点为中心，以适当半径或边长的圆或正方形作为移动面去捕捉适当数目的数据点，并以此展铺一张数学面，内插该中心的高程。地理信息系统原理及应