1、第七讲第七讲 DEMDEM与数字地形分析与数字地形分析第 2 页第 2 页第七讲第七讲 DEM和数字地形分析和数字地形分析数字地形模型（数字地形模型（DTMDTM）是）是2020世纪世纪5050年代由美国年代由美国MITMIT摄影测量试验室主任米勒（摄影测量试验室主任米勒（C.L.MILLERC.L.MILLER）首）首次提出，并用其成功地解决了道路工程中土方估算次提出，并用其成功地解决了道路工程中土方估算等问题。此后它被用于各种线路选线（铁路、公路等问题。此后它被用于各种线路选线（铁路、公路、输电线等）的设计以及各种工程的面积、体积、输电线等）的设计以及各种工程的面积、体积、坡度计算，任意两
2、点间的通视判断即任意断面图绘坡度计算，任意两点间的通视判断即任意断面图绘制。在测绘中被用于绘制等高线、坡度坡向图、立制。在测绘中被用于绘制等高线、坡度坡向图、立体透视图，制作正射影像图以及地图的修测。体透视图，制作正射影像图以及地图的修测。DTMDTM是地理信息系统的基础数据，可用于土地利是地理信息系统的基础数据，可用于土地利用现状的分析、合理规划和洪水险情预报等。在军用现状的分析、合理规划和洪水险情预报等。在军事上可用于导航、作战电子沙盘等。事上可用于导航、作战电子沙盘等。2022-6-1地理信息系统第 3 页第 3 页OUTLINE基本概念基本概念DEM的建立的建立数字地形分析数字地形分析
3、2022-6-1地理信息系统第 4 页第 4 页第一节第一节基本概念基本概念2022-6-1地理信息系统第 5 页第 5 页数字高程模型的概念数字高程模型的概念数字高程模型（数字高程模型（Digital Elevation Model，简称，简称DEM）是通过有限的地形高程数据实现对地形曲）是通过有限的地形高程数据实现对地形曲面的数字化模拟（即地形表面形态的数字化表示）面的数字化模拟（即地形表面形态的数字化表示），它是对二维地理空间上具有连续变化特征地理现，它是对二维地理空间上具有连续变化特征地理现象的模型化表达和过程模拟。象的模型化表达和过程模拟。2022-6-1地理信息系统第 6 页第 6
4、 页数字高程模型的数学概念数字高程模型的数学概念数学意义上的数字高程模型是定义在二维空间上的数学意义上的数字高程模型是定义在二维空间上的连续函数连续函数 。由于连续函数的无限性，。由于连续函数的无限性，DEM通常是通常是将有限的采样点用某种规则连接成一系列的曲面或将有限的采样点用某种规则连接成一系列的曲面或平面片来逼近原始曲面，因此平面片来逼近原始曲面，因此DEM的数学定义为的数学定义为区域区域D的采样点或内插点的采样点或内插点Pj按某种规则按某种规则 连接成的面连接成的面片片M的集合：的集合：2022-6-1地理信息系统, 1, 1,),()(minjDHyxPPMDEMjjjjjiDEM按
5、照其结构，可分为规则格网按照其结构，可分为规则格网DEM、TIN、基于点的基于点的DEM和基于等高线的和基于等高线的DEM等等。第 7 页第 7 页数字地形分析的概念数字地形分析的概念数字地形分析（数字地形分析（Digital Terrain Analysis, DTA），），是指在数字高程模型上进行地形属性计算和特征提是指在数字高程模型上进行地形属性计算和特征提取的数字信息处理技术。取的数字信息处理技术。DTA技术是各种与地形因技术是各种与地形因素相关空间模拟技术的基础。素相关空间模拟技术的基础。2022-6-1地理信息系统第 8 页第 8 页地形属性地形属性 具有明确的数学表达式和物理定义
6、，并可在具有明确的数学表达式和物理定义，并可在DEM上直接量算，如坡度、坡向、曲率上直接量算，如坡度、坡向、曲率等等地形曲面参数（地形曲面参数（parameters） 地表形态和特征的定性表达，可以在地表形态和特征的定性表达，可以在DEM上直接提取，其特点是定义明确，但边界上直接提取，其特点是定义明确，但边界条件有一定的模糊性，难以用数学表达式表达，如在实际的流域单元的划分中，往条件有一定的模糊性，难以用数学表达式表达，如在实际的流域单元的划分中，往往难于确定流域的边界往难于确定流域的边界地形形态特征（地形形态特征（features） 给定地表区域的统计学上的特征给定地表区域的统计学上的特征地
7、形统计特征（地形统计特征（statistics） 在地形曲面参数和地形形态特征的基础上，利用应用学科（如水文学、地貌学和土在地形曲面参数和地形形态特征的基础上，利用应用学科（如水文学、地貌学和土壤学）的应用模型而建立的环境变量，通常以指数形式表达壤学）的应用模型而建立的环境变量，通常以指数形式表达复合地形属性（复合地形属性（compound attributes）2022-6-1地理信息系统第 9 页第 9 页地形分析的主要内容地形分析的主要内容一是在复杂的现实世界地理过程中各影响因子和简单、高一是在复杂的现实世界地理过程中各影响因子和简单、高效、精确、易于理解的抽象与计算机实现中找到平衡。简
8、效、精确、易于理解的抽象与计算机实现中找到平衡。简单地说，就是提取描述地形属性和特征的因子，并利用各单地说，就是提取描述地形属性和特征的因子，并利用各种相关技术分析解释地貌形态、划分地貌形态等。种相关技术分析解释地貌形态、划分地貌形态等。二是二是DTM的可视化分析。数字地形分析中可视化分析的的可视化分析。数字地形分析中可视化分析的重点在于地形特征的可视化表达和信息增强，以帮助传达重点在于地形特征的可视化表达和信息增强，以帮助传达地形曲面参数、地表形态特征和复合地形属性的信息。地形曲面参数、地表形态特征和复合地形属性的信息。2022-6-1地理信息系统第 10 页第 10 页常用的数字地形分析的
9、方法常用的数字地形分析的方法提取坡面地形因子提取坡面地形因子 常用的坡面地形因子有坡度、坡向、平面曲率、坡面曲常用的坡面地形因子有坡度、坡向、平面曲率、坡面曲率、地形起伏度、粗糙度、切割深度等率、地形起伏度、粗糙度、切割深度等 提取特征地形要素提取特征地形要素 流域分析流域分析 可视域分析可视域分析地形统计特征分析地形统计特征分析 进行相关、回归、趋势面、聚类等统计分析进行相关、回归、趋势面、聚类等统计分析2022-6-1地理信息系统第 11 页第 11 页第二节第二节数字高程模型的建立数字高程模型的建立2022-6-1地理信息系统第 12 页第 12 页DEM建立的一般步骤建立的一般步骤从模
10、型论角度讲，就是将源域（地形）表现在另一从模型论角度讲，就是将源域（地形）表现在另一个域（目标域或个域（目标域或DEM）中的一种结构，建模的目）中的一种结构，建模的目的是对复杂的客体进行简化和抽象，并把对客体（的是对复杂的客体进行简化和抽象，并把对客体（源域，源域，DEM中为地形起伏）的研究转移到对模型中为地形起伏）的研究转移到对模型的研究上来。的研究上来。2022-6-1地理信息系统第 13 页第 13 页DEM建立的一般步骤建立的一般步骤首先要为模型构造一个合适的空间结构首先要为模型构造一个合适的空间结构 空间结构一般是规则的（如格网），或不规则的（如不规则三角网空间结构一般是规则的（如格
11、网），或不规则的（如不规则三角网TIN） 模型的可计算性要求有两点，一是空间域的数量、属性域和空间结构是有限的，二是域模型的可计算性要求有两点，一是空间域的数量、属性域和空间结构是有限的，二是域函数是可计算的。函数是可计算的。当空间结构为欧几里德平面，属性域是实数集合时，模型为一自然表当空间结构为欧几里德平面，属性域是实数集合时，模型为一自然表面。将欧几里德平面充当水平的面。将欧几里德平面充当水平的XY平面，属性域给出平面，属性域给出Z坐标（或高坐标（或高程），模型即为数字高程模程），模型即为数字高程模空间结构的构造过程即为空间结构的构造过程即为DEM的格网化过程（形成格网），属性值的格网化过
12、程（形成格网），属性值为高程，构造空间域函数即为内插函数的确定，利用空间域函数进行为高程，构造空间域函数即为内插函数的确定，利用空间域函数进行分析就是求取格网点的函数值。分析就是求取格网点的函数值。2022-6-1地理信息系统第 14 页第 14 页DEM的表示法的表示法区域地表高程模型可以采用多种方法表达，如下图所示。区域地表高程模型可以采用多种方法表达，如下图所示。2022-6-1地理信息系统DEM表示方法 数学方法 整体 局部 傅立叶级数 高次多项式 规则数学分块 不规则数学分块 图形法 点数据 线数据 规 则 不规则 典型特征 密度一致 密度不一致 三角网 邻近网 山峰、洼坑 隘口、边
13、界 水平线 垂直线 典型线 山脊线 谷底线 海岸线 坡度变换线 DEM的表示方法 第 15 页第 15 页DEM的图形表示方法（的图形表示方法（1）（1）线模式）线模式 等高线是表示地形最常见的形式，其地形特征线也是表达地面等高线是表示地形最常见的形式，其地形特征线也是表达地面高程的重要信息源，如山脊线、谷底线、海岸线及坡度变换线高程的重要信息源，如山脊线、谷底线、海岸线及坡度变换线等。等。（2）点模式）点模式 高程矩阵：高程矩阵：DEM最普通的形式是高程矩阵或规则矩形格网最普通的形式是高程矩阵或规则矩形格网（GRID），规则网格通常是正方形、矩形、三角形等规则网），规则网格通常是正方形、矩形
14、、三角形等规则网格。规则网格将区域空间切分为规则的格网单元，每一个格网格。规则网格将区域空间切分为规则的格网单元，每一个格网单元对应一个数值。数学上可以表示为一个矩阵，在计算机实单元对应一个数值。数学上可以表示为一个矩阵，在计算机实现中则是一个二维数组。每个格网单元或数组的一个元素对应现中则是一个二维数组。每个格网单元或数组的一个元素对应一个高程值。一个高程值。2022-6-1地理信息系统第 16 页第 16 页DEM的图形表示方法（的图形表示方法（2）2022-6-1地理信息系统（2）点模式）点模式 不规则三角网（不规则三角网（TIN） 不规则三角网是专为产生不规则三角网是专为产生DEM数据
15、而设计的一种采样表示数据而设计的一种采样表示系统。系统。TIN模型根据区域所有采样点取得的离散数据，按照模型根据区域所有采样点取得的离散数据，按照优化组合的原则，把这些离散点（各三角形的顶点）连接成优化组合的原则，把这些离散点（各三角形的顶点）连接成相互连续的三角面，在连接时尽可能地使每个三角形为锐角相互连续的三角面，在连接时尽可能地使每个三角形为锐角三角形或为三边的长度近似相等，如图所示。区域中任意点三角形或为三边的长度近似相等，如图所示。区域中任意点落在三角面的顶点、边上或三角形内。如果点不在顶点上，落在三角面的顶点、边上或三角形内。如果点不在顶点上，在边上用边的两个顶点的高程，在三角形内
16、的则用三个顶点在边上用边的两个顶点的高程，在三角形内的则用三个顶点的高程，该点的高程值通常通过线性插值的方法得到。所以的高程，该点的高程值通常通过线性插值的方法得到。所以TIN是一个三维空间的分段线性模型。是一个三维空间的分段线性模型。第 17 页第 17 页DEM的图形表示方法（的图形表示方法（3）2022-6-1地理信息系统（2）点模式）点模式不规则三角网示例 第 18 页第 18 页规则格网规则格网DEM的建立的建立2022-6-1地理信息系统x不规则分布点规则分布等高线分布 Y对每一格网点求取格网点高程第 19 页第 19 页DEM内插方法内插方法DEM建立过程中的关键环节是根据采样点
17、的值内插计算格网点上的高程值。内建立过程中的关键环节是根据采样点的值内插计算格网点上的高程值。内插是指根据分布在内插点周围的已知参考点的高程值求出未知点的高程值，它插是指根据分布在内插点周围的已知参考点的高程值求出未知点的高程值，它是是DEM的核心问题的核心问题2022-6-1地理信息系统DEM内插数据分布规则分布内插方法不规则分布内插方法等高线数据内插方法内插范围整体内插方法局部内插方法逐点内插方法内插曲面与参考点关系纯二维内插曲面拟合内插内插函数性质多项式内插线性插值双线性插值高次多项式插值样条内插有限元内插最小二乘配置内插地形特征理解克立金内插多层曲面叠加内插加权平均值内插分形内插傅立叶
18、级数内插第 20 页第 20 页DEM建立的方法（建立的方法（1）摄影测量法摄影测量法 摄影测量法以航空或航天遥感图像为数据源，利用遥感立体像对，摄影测量法以航空或航天遥感图像为数据源，利用遥感立体像对，采用摄影测量的方法建立空间地形立体模型，量取密集数字高程采用摄影测量的方法建立空间地形立体模型，量取密集数字高程数据，建立数据，建立DEM。采集数据的摄影测量仪器包括附有自动记录装。采集数据的摄影测量仪器包括附有自动记录装置的立体测图仪或立体坐标仪、解析测图仪及数字摄影测量系统。置的立体测图仪或立体坐标仪、解析测图仪及数字摄影测量系统。 涉及涉及DEM数据采集的摄影测量采样方法如下：数据采集的
19、摄影测量采样方法如下： （1）沿等高线采样 在地形复杂及陡峭地区，可采用沿等高线跟踪方式进行数据采集，而在平坦地区，则不宜沿等高线采样。 （2）规则格网采样 利用解析测图仪在立体模型中按规则矩形格网进行采样，直接构成规则格网DEM。 2022-6-1地理信息系统第 21 页第 21 页DEM建立的方法（建立的方法（2）摄影测量法摄影测量法 （3）选择采样）选择采样 为了准确反映地形，可根据地形特征进行选择采样，例如沿山脊线、山谷线、断裂线等为了准确反映地形，可根据地形特征进行选择采样，例如沿山脊线、山谷线、断裂线等进行采集以及离散碎部点（如山顶）的采集。这种方法获取的数据尤其适合于不规则三进行
20、采集以及离散碎部点（如山顶）的采集。这种方法获取的数据尤其适合于不规则三角网角网DEM的建立。的建立。 （4）渐进采样）渐进采样 这种采样方法的目的是使采样点分布合理，即平坦地区样点少，地形复杂区的样点较多。这种采样方法的目的是使采样点分布合理，即平坦地区样点少，地形复杂区的样点较多。 （5）混合采样）混合采样 为了同步考虑采样的合理性和效率，可将规则采样或渐进采样与选择采样结合进行混合为了同步考虑采样的合理性和效率，可将规则采样或渐进采样与选择采样结合进行混合采样，即在规则采样的基础上再进行沿特征线、点采样。利用混合采样可建立附加地形采样，即在规则采样的基础上再进行沿特征线、点采样。利用混合
21、采样可建立附加地形特征的规则格网特征的规则格网DEM，也可建立附加特征的不规则三角网，也可建立附加特征的不规则三角网DEM。 （6）自动化）自动化DEM数据采集数据采集 上述方法均是基于解析测图仪或机助制图系统利用半自动的方法进行上述方法均是基于解析测图仪或机助制图系统利用半自动的方法进行DEM数据采集，数据采集，现在已经可以利用自动化测图系统进行完全自动化的现在已经可以利用自动化测图系统进行完全自动化的DEM数据采集。此时可按相片上数据采集。此时可按相片上的规则格网利用数字影像匹配进行数据采集。的规则格网利用数字影像匹配进行数据采集。 最后摄影测量获取的最后摄影测量获取的DEM数据点都要按一
22、定插值方法转成规则格网数据点都要按一定插值方法转成规则格网DEM或规则三角网或规则三角网DEM格式数据。格式数据。2022-6-1地理信息系统第 22 页第 22 页DEM建立的方法（建立的方法（3）地形图数字化法地形图数字化法 这种方法以比例尺大于这种方法以比例尺大于1:1万的国家近期地形图为数据源，万的国家近期地形图为数据源，从中量取中等密度地面点集的高程数据，建立从中量取中等密度地面点集的高程数据，建立DEM。其方。其方法有下面几种：法有下面几种：（1）手工采集法）手工采集法 采用方格膜片、网点板或带刻划的平移角尺叠置在地形图上，采用方格膜片、网点板或带刻划的平移角尺叠置在地形图上，并使
23、地形图的格网与网点板或膜片的格网线逐格匹配定位，并使地形图的格网与网点板或膜片的格网线逐格匹配定位，自上而下，逐行从左到右量取高程。当格网交点落在相邻等自上而下，逐行从左到右量取高程。当格网交点落在相邻等高线之间时，用目视线性内插法估计高程值。它的优点是几高线之间时，用目视线性内插法估计高程值。它的优点是几乎不需要购置仪器设备，而且操作简便，但效率低。乎不需要购置仪器设备，而且操作简便，但效率低。2022-6-1地理信息系统第 23 页第 23 页DEM建立的方法（建立的方法（4）（2）手扶跟踪数字化仪采集法）手扶跟踪数字化仪采集法 可采用如下可采用如下3种方式进行采集：种方式进行采集： 逐条
24、等高线的线方式连续采集样点，并采集所有高程注记点作补充，逐条等高线的线方式连续采集样点，并采集所有高程注记点作补充，这种方式适合于等高线较稀疏的平坦地区。这种方式适合于等高线较稀疏的平坦地区。 沿主要等高线采集特征点，并选择采集高程注记点和线性加密点作补沿主要等高线采集特征点，并选择采集高程注记点和线性加密点作补充。充。 沿曲线和坡折线采集特征点，并补采峰沿曲线和坡折线采集特征点，并补采峰鞍线和水边线的支撑点，分鞍线和水边线的支撑点，分别以等高线，峰别以等高线，峰鞍链和边界格式存储。鞍链和边界格式存储。（3）扫描采集法）扫描采集法 这种方法利用扫描仪扫描地图，从地图扫描数据中自动地建立这种方法
25、利用扫描仪扫描地图，从地图扫描数据中自动地建立DEM。这种方式采集速度快。随着研究的不断深入，一些难点和瓶颈问题被这种方式采集速度快。随着研究的不断深入，一些难点和瓶颈问题被逐步解决，目前自动建立逐步解决，目前自动建立DEM技术已达到实用水平。技术已达到实用水平。 2022-6-1地理信息系统第 24 页第 24 页DEM建立的方法（建立的方法（5）地面测量法地面测量法 这种方法以地面实测记录为数据源，利用这种方法以地面实测记录为数据源，利用GPS、全站仪和电子手簿或测距经纬仪等设、全站仪和电子手簿或测距经纬仪等设备，在已知点位的测站上，观测到目标点的备，在已知点位的测站上，观测到目标点的方向
26、、距离和高差方向、距离和高差3个要素。计算出目标点的个要素。计算出目标点的（x,y,z）三维坐标，存储于电子手簿或袖珍）三维坐标，存储于电子手簿或袖珍计算机中，成为建立计算机中，成为建立DEM的原始数据。这种的原始数据。这种方法一般用于建立小范围大比例尺区域的方法一般用于建立小范围大比例尺区域的DEM，对高程的精度要求较高。，对高程的精度要求较高。 2022-6-1地理信息系统第 25 页第 25 页DEMDEM的空间插值方法（的空间插值方法（1 1）分块内插分块内插 是是把需要建立把需要建立DEM的地区，切割成一定大小的规则方块，的地区，切割成一定大小的规则方块，形状通常为正方形，它的尺寸应
27、根据地形复杂程度和数据源形状通常为正方形，它的尺寸应根据地形复杂程度和数据源的比例尺确定。在每一个分块上展铺一张数学面，相邻分块的比例尺确定。在每一个分块上展铺一张数学面，相邻分块之间有适当宽度的重叠带，以使重叠带内全部数据点成为相之间有适当宽度的重叠带，以使重叠带内全部数据点成为相邻块展铺数学面时的共用数据，保证一张数学面能够较平滑邻块展铺数学面时的共用数据，保证一张数学面能够较平滑地与相邻分块的数学面拼接。这种内插方法的优点是可以得地与相邻分块的数学面拼接。这种内插方法的优点是可以得到光滑连续的空间曲面。到光滑连续的空间曲面。 2022-6-1地理信息系统由于由于DEM采样的数据点呈离散分
28、布形式，或是数据点虽按格网排列，采样的数据点呈离散分布形式，或是数据点虽按格网排列，但格网的密度不能满足使用的要求，这就需要以数据点为基础进行插但格网的密度不能满足使用的要求，这就需要以数据点为基础进行插值运算。值运算。DEM内插按插点分布范围，可分为分块内插、剖分内插和内插按插点分布范围，可分为分块内插、剖分内插和单点移面内插三类。单点移面内插三类。第 26 页第 26 页DEMDEM的空间插值方法（的空间插值方法（2 2）剖分内插剖分内插 是把需要建立是把需要建立DTM的地区切割成大小和形状的地区切割成大小和形状不同的子区（剖分），子区间拥有公共边但不同的子区（剖分），子区间拥有公共边但不
29、重叠，在该区内展铺一个数学面，内插剖不重叠，在该区内展铺一个数学面，内插剖分区内任意点的高程。该法只在剖分间边界分区内任意点的高程。该法只在剖分间边界端点处重合，通常没有严格重合的边界，所端点处重合，通常没有严格重合的边界，所以既不连续，也不光滑。剖分多边形的顶点以既不连续，也不光滑。剖分多边形的顶点都是数据点，最常见的数据点个数为都是数据点，最常见的数据点个数为3，与，与TIN结构相同。结构相同。2022-6-1地理信息系统第 27 页第 27 页DEMDEM的生成流程的生成流程DEMDEM生成的全过程包括：生成的全过程包括：原始数据获取、原始数据获取、DEMDEM模模型构造、数据插值、在所
30、定数据结构支持下的数据型构造、数据插值、在所定数据结构支持下的数据存储和模型输出。存储和模型输出。如下图所示，以地形图数字化为如下图所示，以地形图数字化为例，说明例，说明DEMDEM的生成过程。的生成过程。2022-6-1地理信息系统第 28 页第 28 页DEMDEM的生成流程的生成流程2022-6-1地理信息系统全要素地形图地图扫描 人机交互等高线矢量化加测注记点 等高线赋值与检查 周边等高线地图的数据获取 构造三角网 内插DEM网格 DEM建库与刻盘 质量检测与元数据文件记录 第 29 页第 29 页第三节第三节 数字地面模型分析数字地面模型分析2022-6-1地理信息系统第 30 页第
31、 30 页 基本地形因子计算基本地形因子计算一、一、 1. 1. 极值高程和最大高差极值高程和最大高差 2. 2. 相对高程和平均高程相对高程和平均高程 3. 3. 坡度和坡向计算坡度和坡向计算 地面上某点的坡度是表示地面在该点倾斜程度的一个量，因此，它是一个既有大小又有方向的量，即矢量。坡度定义为水平面和地形表面之间夹角的正切值；坡向为坡面法线在水平面上的投影与正北方向的夹角如左下图所示。坡度和坡向的计算一般采用拟合曲面法，通常在33个DEM格网窗口中进行，如右下图所示。每个窗口中心为一个高程点。 ynoGxAze e1 e5 e8 e4 e2 e6 e3 e7 地表单元坡度和坡向示意图 3
32、3窗口计算点的坡度和坡向 第 31 页第 31 页 坡向分析坡向分析 坡向图是坡向的类别显示图，斜坡的倾斜方向可取方位角0360中的任意方向。坡向一般分为9类：东、南、西、北、东北、西北、东南、西南和平地。在实际应用中，可以综合为四种坡向，即平缓坡、阳坡、半阳坡和阴坡，分别用1、2、3和4表示。如图所示。324NWNESE3NESWSW坡向的综合表示第 32 页第 32 页 4. 4. 表面积的计算表面积的计算 根据数字高程模型很容易求得地表面积，其计算可看作是所包含各个网络的表面积之和。 5. 5. 投影面积的计算投影面积的计算 6. 体积的计算体积的计算 7. 7. 剖面积的计算剖面积的计
33、算 8. 8. 地表粗糙度计算地表粗糙度计算 地表粗糙度是反映地表的起伏变化与侵蚀程度的指标，一般定义为地表单元的曲面面积与投影面积之比： CZ = S表面积/S投影面积 显然，这种定义对于表面光滑的斜面也可求出不同的粗糙度，这是不合适的。这里用对顶点连线L1与L2中点的距离D来表示粗糙度，如下图所示，D值愈大，说明单元的四个顶点的起伏变化也愈大。YZZi, j-1i+1, j+1 ZDi, jZi+1, j2LZXO地表粗糙度计算坡向分析坡向分析第 33 页第 33 页在地表的基本形态中，山谷和山脊是常见的两种主要形态。它在区域地形研究和制图综合中具有重要的意义。利用数字高程模型可对谷脊特征
34、作概略分析。 1. 1. 谷点和脊点的判定谷点和脊点的判定 谷点是地势相对最低的点集，脊点为地势相对最高的点集如图711所示，要判定高 程为Z网格的形态特征，按照以下判别式可直接提取谷点和脊点。 如果（Z i,j-1 Z i,j ) （Z i,j+1 Z i,j )0 当Z i,j+1 Z i,j )时，则P(i,j)=1 当Z i,j+10 j 当Z i+1,j Z i,j时，则P(i,j)=1 当Z i+1,j Z i,j时，则P(i,j)=1 如果和或和同时成立,则P(i,j)= 2 如果以上条件均不成立,则P(i,j)= 0Zi,j Zi+1,j Zi-1,j Zi,j-1 Zi, +
35、1 P(i,j)=-1表示谷点 1表示脊点2表示鞍点0表示其它点谷脊特征分析谷脊特征分析 第 34 页第 34 页 2. 2. 沟谷密度分析沟谷密度分析 沟谷密度是表征地面破碎程度的一种指标,它是沟谷总长度(L)与地表单元总面积(A)之比.提取谷点和脊点,将地表单元内所有谷点在单元区域内的延伸长度累加,便获得单元的沟谷长度.沟谷密度为:nkmiikALVD11/nK,3,2, 1mi,3,2,1谷脊特征分析谷脊特征分析 第 35 页第 35 页 3. 3. 切割深度分析切割深度分析. 地表单元的谷点与最近脊点的平均高差为谷点的切割深度,区域平均切割深度为若干谷点切割深度的平均值。nKVKRKn
36、hhVH1/ )(nK,3,2,1式中 为距该谷点最近的脊点的平均高程值， 为谷点高程。 RhVh谷脊特征分析谷脊特征分析 第 36 页第 36 页地形剖面可以概括而直观地研究区域的地势、地质和水文特征，也可以与其他地理变量相叠加，提供农业生产布局的立体背影分析、土地利用规划以及工程选址选线等的参考依据。从DEM可以很方便地制作任一方向上的地形剖面图。只要知道所绘剖面线在数字高程模型中的起点位置和终点位置，就可以惟一地确定这条剖面线与DEM网格各个交点的平面位置及其高程，根据选定的垂直比例尺和水平比例尺，就能自动绘出所需要的剖面图。可用梯形法、辛普森法来计算剖面面积。当两个剖面面积相互平行时，
37、可用类似的方法计算该两个剖面之间所夹的体积。地形剖面地形剖面分析分析 第 37 页第 37 页DEM数据不仅直接用于各种地形因素的分析，而且还可以与有关信息进行复合，研究地形要素与其他要素之间的相互联系。其中洪水淹没损失估算，就是研究DEM与土地利用之间的关系。为了科学合理地估价，首先将数字高程的数据与土地利用数据进行匹配；其次是根据淹没高程，在DEM上确定淹没范围；最后统计淹没范围内的土地类型和面积，就能精确估算出淹没损失。 淹没损失估算淹没损失估算 第 38 页第 38 页2022-6-1地理信息系统在DEM数据文件的基础上，进行地貌类型的自动分类。首先根据区域的地形特点，拟定地形分类的高
38、程界值；然后，计算机根据确定的高程界值自动提取地形类型信息，便可获得区域的地形分类系统，如平原、丘陵、低山、中山和高山等；最后输出地貌类型图。 地貌类型的自动划分地貌类型的自动划分 第 39 页第 39 页高程矩阵没有存储山脊线、山谷线等地形特征线，或者在地形图数字化时，对这些地形特征线没有单独数字化，在这种情况下，用程序自动地将它们从高程矩阵中提取出来也许是必要的。例如从叠置到DEM的卫星图像上勾绘出集水范围线，使遥感图像与特殊地理景观联系在一起。高程矩阵用于其他数量分析如费用量、集水范围、旅行时间等时，应有一种方法来描述线、面特征。从从DEMDEM数据自动形成地形轮廓线数据自动形成地形轮廓
39、线第 40 页第 40 页 1. 1. 山脊线和谷底线的探测山脊线和谷底线的探测 为了自动探测山脊线和谷底线，设计了专门的运算算子。较为简单的算子是4个像元的局部算子。该算子在高程矩阵中移动并比较每一位置处4个像元的高程值，同时标出其中高程最大（探测谷底线）或最小（探测山脊线）的像元。标记过程完成后，剩下为标记的像元就是山脊线和山谷线所在的像元。下一步就是把它们连接成线模式，形成山脊线或山谷线。 2. 2. 集水范围的确定集水范围的确定 集水范围即流域范围的确定对流域分析十分必要。流域探测除确定边界线外，还要将整个范围从整个数据库中分离出来。探测方法是：首先需要交互式地确定河流流域的出口，并作
40、为搜索工作的起点。以33算子的中心像元置于起始点，比较中心像元相邻近的8个像元的坡向。如果坡向朝向中心像元，则认为它是中心像元的上游，算子的中心像元移至新的“上游”点，重复比较过程又能得到新的“上游”点。以由“上游”标志的不予比较，整个数据范围都运算完毕后，流域范围就全部标记出来了。用户可以对这些像元重新编码，形成某一流域的分布图。从从DEMDEM数据自动形成地形轮廓线数据自动形成地形轮廓线第 41 页第 41 页1、通视分析也称可视分析，它实质属于地形进行最优化处理的范畴，如设置雷达站、电视台的发射站、道路选择、航海导航等，在军事上如布设阵地（炮兵阵地、电子对抗阵地）、设置观察哨所、铺架通信
41、线路等。 数字高程模型的建立为这类分析提供了有利基础，能方便地算出一个观察点所能看到的各个部分。在DEM中辨认出观察点所在的位置，从这个位置引出一簇射线，比较射线通过的每个点（高程矩阵中即为像元）的高程，将不被物体隐藏的各点进行特殊编码，便可得到一幅简单的地图。HAvsTEODOCBhvs通视计算示意图 通视分析通视分析第 42 页第 42 页立体图是表现物体三维模型最直观形象的图形，它可以生动逼真地描述制图对象在平面和空间上分布的形态特征和构造关系。通过分析立体图可以了解地理模型表面的起伏状况，可以看出各个断面的形态，这对研究区域的轮廓形态、变化规律以及内部结构是非常有益的。计算机自动绘制透
42、视立体图的理论基础是透视原理，而DEM是其绘制的数据基础。制作透视立体图的基本步骤包括：建立透视变换基础、DEM高程阵列剖面布设、消除隐藏线，处理和粘贴表面影像与纹理。 根据根据DEMDEM制作透视立体图制作透视立体图第 43 页第 43 页2022-6-1地理信息系统晕渲图是通过模拟实际地面反映地形起伏特征的重要地图制图学方法，在各种小比例尺地形图、地理图，以及有关专题地图上得到非常广泛的应用。但是，传统的人工描绘晕渲图的方法费工、费时、成本较高。利用DEM数据作为信息源，以地面光照通量为依据，计算该栅格所输出的灰度值，由此产生的晕渲图具有相当逼真的立体效果。自动地貌晕渲图的计算步骤为：首先
43、是根据DEM数据计算坡度和坡向，然后将坡向数据与光源方向比较，向光源的斜坡得到浅色调灰度值，反方向的斜坡得到深色调灰度值，介于中间坡向的斜坡得到中间灰度值。灰度值的大小则按坡度进一步确定。根据根据DEMDEM制作晕渲图制作晕渲图第 44 页第 44 页2022年6月1日星期三44表面积yxffSAyxdd122采用Simpson法 近似计算：11121021222001020029491919136136191361361hhhhhhhhhaS三维几何参数计算三维几何参数计算第 45 页第 45 页2022年6月1日星期三45体积4)(43212hhhhaV第 46 页第 46 页2022年6
44、月1日星期三461、坡度和坡向分析 坡度：水平面与地形表面之间夹角的正切值； 坡向：坡面法线在水平面的投影与正北方向的夹角法线交线坡度坡向地形因子的自动提取地形因子的自动提取第 47 页第 47 页2022年6月1日星期三47坡度i-1,j-1i,j-1i+1,j-1i-1,ji-1,j+1i,ji,j+1i+1,ji+1,j+1yzzyzxzzxzjijijiji2)(2)(,1,11,1,22xzxztg坡向)(1xzyztg地形因子的自动提取地形因子的自动提取第 48 页第 48 页2022年6月1日星期三482、 地表粗糙度反映地表的起伏变化和侵蚀程度的指标，一般定义为地表单元的曲面面
45、积与其在水平面上的投影面积之比。xzyi,ji+1,ji,j+1i+1,j+1jijijijijijijijijizzzzzzzzDR, 11,1, 1, 11,1, 1,2122地形因子的自动提取地形因子的自动提取第 49 页第 49 页2022年6月1日星期三493、高程变异分析平均高程：格网4个顶点的高程平均值相对高程：格网平均高程与研究区最低点高程之差高程变异V：格网顶点高程的标准差与平均值之比21412min41)(41/41kkzkkzzSzSVzzDzz地形因子的自动提取地形因子的自动提取第 50 页第 50 页2022年6月1日星期三504、谷脊特征分析、谷脊特征分析谷：地势相
46、对最低点的集合脊：地势相对最高点的集合谷、脊线是汇水区提取的依据（1）当（zi，j-1zi，j）（zi，j+1zi，j）0时 若zi，j+1zi，j，则VR（i，j）=1 若zi，j+1zi，j，则VR（i，j）=1（2）当（zi-1，jzi，j）（zi+1，jzi，j）0时 若zi+1，jzi，j，则VR（i，j）=1 若zi+1，jzi，j，则VR（i，j）=1在其他情况下，VR（i，j）=0其中VR（i，j）=-1 表示谷点 1 表示脊点 0 表示其他点i-1,j-1i,j-1i+1,j-1i-1,ji-1,j+1i,ji,j+1i+1,ji+1,j+1地形因子的自动提取地形因子的自动提
47、取第 51 页第 51 页2022年6月1日星期三51地表地貌形态的自动分类地表地貌形态的自动分类DTM计算地有形态要素H、H、坡度、坡向等拟定地形分类决策表地形自动分类第 52 页第 52 页2022年6月1日星期三52地形剖面图的绘制地形剖面图的绘制1 2 3 4 5 6 7 8 9 10 nm87654321i2, j2i1, j1i2, j2dxdy第 53 页第 53 页2022年6月1日星期三53地形剖面图第 54 页第 54 页2022年6月1日星期三54通视性分析通视性分析通视问题可以分为五类：已知一个或一组观察点，找出某一地形的可见区域；欲观察到某一区域的全部地形表面，计算最
48、少观察点数量；在观察点数量一定的前提下，计算能获得的最大观察区域；以最小代价建造观察塔，要求全部区域可见；在给定建造代价的前提下，求最大可见区。通视分析一般有两种情况：v视点与待判定点之间的可视性（Intervisibility）；v根据已知视点计算可视域(Viewshed)，即视点能可视的地形表面区域集合。第 55 页第 55 页2022年6月1日星期三55可视域分析问题 点对区域的通视分析可在点对点通视分析算法的基础上扩展。判定点P沿DTM数据边缘顺时针移动，逐点检查视点至P点的直线上的点是否通视。 点对点通视问题 第 56 页第 56 页2022年6月1日星期三56考虑视点高度及地物高度的两点通视问题视点V视线判定点PhVhP第 57 页第 57 页2022-6-1地理信息系统
