第七章数字地面模型（DTM）及其应用7.1数字地形模型（DTM）与地形分析7.1.1 概述数字地形模型（DTM，DigitalTerrainModel）最初是为了高速公路的自动设计提出来的（Mi11er，1956）。此后，它被用于各种线路选线（铁路、公路、输电线）的设计以及各种工程的面积、体积、坡度计算，任意两点间的通视判断及任意断面图绘制。在测绘中被用于绘制等高线、坡度坡向图、立体透视图，制作正射影像图以及地图的修测。在遥感应用中可作为分类的辅助数据。它还是地理信息系统的基础数据，可用于土地利用现状的分析、合理规划及洪水险情预报等。在军事上可用于导航及导弹制导、作战电子沙盘等。对DTM的研究包括DTM的精度问题、地形分类、数据采集、DTM的粗差探测、质量控制、数据压缩、DTM应用以及不规则三角网DTM的建立与应用等。7.1.2DTM和DEM从数学的角度，高程模型是高程Z关于平面坐标X，Y两个自变量的连续函数，数字高程模型（DEM）只是它的一个有限的离散表示。高程模型最常见的表达是相对于海平面的海拔高度，或某个参考平面的相对高度，所以高程模型又叫地形模型。实际上地形模型不仅包含高程属性，还包含其它的地表形态属性，如坡度、坡向等。数字地形模型是地形表面形态属性信息的数字表达，是带有空间位置特征和地形属性特征的数字描述。数字地形模型中地形属性为高程时称为数字高程模型（DigitalElevationModel，简称DEM）。高程是地理空间中的第三维坐标。由于传统的地理信息系统的数据结构都是二维的，数字高程模型的建立是一个必要的补充。DEM通常用地表规则网格单元构成的高程矩阵表示，广义的DEM还包括等高线、三角网等所有表达地面高程的数字表示。在地理信息系统中，DEM是建立DTM的基础数据，其它的地形要素可由DEM直接或间接导出，称为“派生数据”，如坡度、坡向。7.1.3DEM的表示法一个地区的地表高程的变化可以采用多种方法表达，用数学定义的表面或点、线、影像都可用来表示DEM，如图1所示
第七章 数字地面模型（DTM）及其应用 7.1 数字地形模型（DTM）与地形分析 7.1.1 概述 数字地形模型（DTM, Digital Terrain Model）最初是为了高速公路的自动 设计提出来的（Miller，1956）。此后，它被用于各种线路选线（铁路、公路、 输电线）的设计以及各种工程的面积、体积、坡度计算，任意两点间的通视判断 及任意断面图绘制。在测绘中被用于绘制等高线、坡度坡向图、立体透视图，制 作正射影像图以及地图的修测。在遥感应用中可作为分类的辅助数据。它还是地 理信息系统的基础数据，可用于土地利用现状的分析、合理规划及洪水险情预报 等。在军事上可用于导航及导弹制导、作战电子沙盘等。对 DTM 的研究包括 DTM 的精度问题、地形分类、数据采集、DTM 的粗差探测、质量控制、数据压缩、DTM 应用以及不规则三角网 DTM 的建立与应用等。 7.1.2 DTM 和 DEM 从数学的角度，高程模型是高程 Z 关于平面坐标 X，Y 两个自变量的连续函 数，数字高程模型（DEM）只是它的一个有限的离散表示。高程模型最常见的表 达是相对于海平面的海拔高度，或某个参考平面的相对高度，所以高程模型又叫 地形模型。实际上地形模型不仅包含高程属性，还包含其它的地表形态属性，如 坡度、坡向等。 数字地形模型是地形表面形态属性信息的数字表达，是带有空间位置特征和 地形属性特征的数字描述。数字地形模型中地形属性为高程时称为数字高程模型 （Digital Elevation Model，简称 DEM）。高程是地理空间中的第三维坐标。 由于传统的地理信息系统的数据结构都是二维的，数字高程模型的建立是一个必 要的补充。DEM 通常用地表规则网格单元构成的高程矩阵表示，广义的 DEM 还包 括等高线、三角网等所有表达地面高程的数字表示。在地理信息系统中，DEM 是 建立 DTM 的基础数据，其它的地形要素可由 DEM 直接或间接导出，称为“派生数 据”，如坡度、坡向。 7.1.3 DEM 的表示法 一个地区的地表高程的变化可以采用多种方法表达，用数学定义的表面或 点、线、影像都可用来表示 DEM，如图 1 所示
1）数学方法用数学方法来表达，可以采用整体拟合方法，即根据区域所有的高程点数据，用傅立叶级数和高次多项式拟合统一的地面高程曲面。也可用局部拟合方法，将地表复杂表面分成正方形规则区域或面积大致相等的不规则区域进行分块搜索，根据有限个点进行拟合形成高程曲面。2）图形方法1.线模式等高线是表示地形最常见的形式。其它的地形特征线也是表达地面高程的重要信息源，如山脊线、谷底线、海岸线及坡度变换线等。2.．点模式用离散采样数据点建立DEM是DEM建立常用的方法之一。数据采样可以按规则格网采样，可以是密度一致的或不一致的：可以是不规则采样，如不规则三角网、邻近网模型等；也可以有选择性地采样，采集山峰、洼坑、隘口、边界等重要特征点。厂傅立叶级数整体高次多项式数学方法规则数学分块局部不规则数学分块厂密度一致规则L密度不一致DEM表示方法一厂三角网一点数据+不规则L邻近网山峰、洼坑L典型特征口、边界L图形法一水平线一山脊线一线数据十垂直线谷底线L典型线海岸线L坡度变换线图1：DEM的表示方法在地理信息系统中，DEM最主要的三种表示模型是：规则格网模型，等高线模型和不规则三角网模型。7.2DEM的主要表示模型
1）数学方法 用数学方法来表达，可以采用整体拟合方法，即根据区域所有的高程点数据， 用傅立叶级数和高次多项式拟合统一的地面高程曲面。也可用局部拟合方法，将 地表复杂表面分成正方形规则区域或面积大致相等的不规则区域进行分块搜索， 根据有限个点进行拟合形成高程曲面。 2）图形方法 1. 线模式 等高线是表示地形最常见的形式。其它的地形特征线也是表达地面高程的重 要信息源，如山脊线、谷底线、海岸线及坡度变换线等。 2. 点模式 用离散采样数据点建立 DEM 是 DEM 建立常用的方法之一。数据采样可以按规 则格网采样，可以是密度一致的或不一致的；可以是不规则采样，如不规则三角 网、邻近网模型等；也可以有选择性地采样，采集山峰、洼坑、隘口、边界等重 要特征点。 图 1：DEM 的表示方法 在地理信息系统中，DEM 最主要的三种表示模型是：规则格网模型，等高线 模型和不规则三角网模型。 7.2 DEM 的主要表示模型
7.2.1规则格网模型规则网格，通常是正方形，也可以是矩形、三角形等规则网格。规则网格将区域空间切分为规则的格网单元，每个格网单元对应一个数值。数学上可以表示为一个矩阵，在计算机实现中则是一个二维数组。每个格网单元或数组的一个元素，对应一个高程值，如图1所示。2591786350536344554035945816402505151OU10084.66425564665465511034784665665968249662D915619o6St8678570829828056155374b16962668388731458634570566274577174图1：格网DEM对于每个格网的数值有两种不同的解释。第一种是格网栅格观点，认为该格网单元的数值是其中所有点的高程值，即格网单元对应的地面面积内高程是均一的高度，这种数字高程模型是一个不连续的函数。第二种是点栅格观点，认为该网格单元的数值是网格中心点的高程或该网格单元的平均高程值，这样就需要用一种插值方法来计算每个点的高程。计算任何不是网格中心的数据点的高程值，使用周围4个中心点的高程值，采用距离加权平均方法进行计算，当然也可使用样条函数和克里金插值方法。规则格网的高程矩阵，可以很容易地用计算机进行处理，特别是栅格数据结构的地理信息系统。它还可以很容易地计算等高线、坡度坡向、山坡阴影和自动提取流域地形，使得它成为DEM最广泛使用的格式，目前许多国家提供的DEM数据都是以规则格网的数据矩阵形式提供的。格网DEM的缺点是不能准确表示地形的结构和细部，为避免这些问题，可采用附加地形特征数据，如地形特征点、山脊线、谷底线、断裂线，以描述地形结构。格网DEM的另一个缺点是数据量过大，给数据管理带来了不方便，通常要进
7.2.1 规则格网模型 规则网格，通常是正方形，也可以是矩形、三角形等规则网格。规则网格将 区域空间切分为规则的格网单元，每个格网单元对应一个数值。数学上可以表示 为一个矩阵，在计算机实现中则是一个二维数组。每个格网单元或数组的一个元 素，对应一个高程值，如图 1 所示。 图 1：格网 DEM 对于每个格网的数值有两种不同的解释。第一种是格网栅格观点，认为该格 网单元的数值是其中所有点的高程值，即格网单元对应的地面面积内高程是均一 的高度，这种数字高程模型是一个不连续的函数。第二种是点栅格观点，认为该 网格单元的数值是网格中心点的高程或该网格单元的平均高程值，这样就需要用 一种插值方法来计算每个点的高程。计算任何不是网格中心的数据点的高程值， 使用周围 4 个中心点的高程值，采用距离加权平均方法进行计算，当然也可使用 样条函数和克里金插值方法。 规则格网的高程矩阵，可以很容易地用计算机进行处理，特别是栅格数据结 构的地理信息系统。它还可以很容易地计算等高线、坡度坡向、山坡阴影和自动 提取流域地形，使得它成为 DEM 最广泛使用的格式，目前许多国家提供的 DEM 数据都是以规则格网的数据矩阵形式提供的。格网 DEM 的缺点是不能准确表示地 形的结构和细部，为避免这些问题，可采用附加地形特征数据，如地形特征点、 山脊线、谷底线、断裂线，以描述地形结构。 格网 DEM 的另一个缺点是数据量过大，给数据管理带来了不方便，通常要进
行压缩存储。DEM数据的无损压缩可以采用普通的栅格数据压缩方式，如游程编码、块码等，但是由于DEM数据反映了地形的连续起伏变化，通常比较“破碎”，普通压缩方式难以达到很好的效果；因此对于网格DEM数据，可以采用哈夫曼编码进行无损压缩；有时，在牺牲细节信息的前提下，可以对网格DEM进行有损压缩，通常的有损压缩大都是基于离散余弦变换（DiscreteCosineTransformation，DcT）或小波变换（WaveletTransformation）的，由于小波变换具有较好的保持细节的特性，近年来将小波变换应用于DEM数据处理的研究较多。7.2.2等高线模型等高线模型表示高程，高程值的集合是已知的，每一条等高线对应一个已知的高程值，这样一系列等高线集合和它们的高程值一起就构成了一种地面高程模型。如图2所示。图2：等高线等高线通常被存成一个有序的坐标点对序列，可以认为是一条带有高程值属性的简单多边形或多边形弧段。由于等高线模型只表达了区域的部分高程值，往往需要一种插值方法来计算落在等高线外的其它点的高程，又因为这些点是落在两条等高线包围的区域内，所以，通常只使用外包的两条等高线的高程进行插值。等高线通常可以用二维的链表来存储。另外的一种方法是用图来表示等高线的拓扑关系，将等高线之间的区域表示成图的节点，用边表示等高线本身。此方法满足等高线闭合或与边界闭合、等高线互不相交两条拓扑约束。这类图可以改
行压缩存储。DEM 数据的无损压缩可以采用普通的栅格数据压缩方式，如游程编 码、块码等，但是由于 DEM 数据反映了地形的连续起伏变化，通常比较“破碎”， 普通压缩方式难以达到很好的效果；因此对于网格 DEM 数据，可以采用哈夫曼编 码进行无损压缩；有时，在牺牲细节信息的前提下，可以对网格 DEM 进行有损压 缩 ， 通 常 的 有 损 压 缩 大 都 是 基 于 离 散 余 弦 变 换 （ Discrete Cosine Transformation，DCT）或小波变换（Wavelet Transformation）的，由于小波 变换具有较好的保持细节的特性，近年来将小波变换应用于 DEM 数据处理的研究 较多。 7.2.2 等高线模型 等高线模型表示高程，高程值的集合是已知的，每一条等高线对应一个已知 的高程值，这样一系列等高线集合和它们的高程值一起就构成了一种地面高程模 型。如图 2 所示。 图 2：等高线 等高线通常被存成一个有序的坐标点对序列，可以认为是一条带有高程值属 性的简单多边形或多边形弧段。由于等高线模型只表达了区域的部分高程值，往 往需要一种插值方法来计算落在等高线外的其它点的高程，又因为这些点是落在 两条等高线包围的区域内，所以，通常只使用外包的两条等高线的高程进行插值。 等高线通常可以用二维的链表来存储。另外的一种方法是用图来表示等高线 的拓扑关系，将等高线之间的区域表示成图的节点，用边表示等高线本身。此方 法满足等高线闭合或与边界闭合、等高线互不相交两条拓扑约束。这类图可以改
造成一种无圈的自由树。下图为一个等高线图和它相应的自由树（图9-4）。其它还有多种基于图论的表示方法。图3：等高线和相应的自由树7.2.3不规则三角网（TIN）模型尽管规则格网DEM在计算和应用方面有许多优点，但也存在许多难以克服的缺陷：1）在地形平坦的地方，存在大量的数据穴余：2）在不改变格网大小的情况下，难以表达复杂地形的突变现象；3）在某些计算，如通视问题，过分强调网格的轴方向。不规则三角网（TriangulatedIrregularNetwork，TIN）是另外一种表示数字高程模型的方法[Peuker等，1978]，它既减少规则格网方法带来的数据穴余，同时在计算（如坡度）效率方面又优于纯粹基于等高线的方法。TIN模型根据区域有限个点集将区域划分为相连的三角面网络，区域中任意点落在三角面的顶点、边上或三角形内。如果点不在顶点上，该点的高程值通常通过线性插值的方法得到（在边上用边的两个顶点的高程，在三角形内则用三个顶点的高程）。所以TIN是一个三维空间的分段线性模型，在整个区域内连续但不可微。TIN的数据存储方式比格网DEM复杂，它不仅要存储每个点的高程，还要存储其平面坐标、节点连接的拓扑关系，三角形及邻接三角形等关系。TIN模型在概念上类似于多边形网络的矢量拓扑结构，只是TIN模型不需要定义“岛”和“洞”的拓扑关系
造成一种无圈的自由树。下图为一个等高线图和它相应的自由树（图 9-4）。其 它还有多种基于图论的表示方法。 B A F C G E H D 图 3：等高线和相应的自由树 7.2.3 不规则三角网（TIN）模型 尽管规则格网 DEM 在计算和应用方面有许多优点，但也存在许多难以克服的 缺陷： 1）在地形平坦的地方，存在大量的数据冗余； 2）在不改变格网大小的情况下，难以表达复杂地形的突变现象； 3）在某些计算，如通视问题，过分强调网格的轴方向。 不规则三角网（Triangulated Irregular Network, TIN）是另外一种表示 数字高程模型的方法[Peuker 等，1978]，它既减少规则格网方法带来的数据冗 余，同时在计算（如坡度）效率方面又优于纯粹基于等高线的方法。 TIN 模型根据区域有限个点集将区域划分为相连的三角面网络，区域中任意 点落在三角面的顶点、边上或三角形内。如果点不在顶点上，该点的高程值通常 通过线性插值的方法得到（在边上用边的两个顶点的高程，在三角形内则用三个 顶点的高程）。所以 TIN 是一个三维空间的分段线性模型，在整个区域内连续但 不可微。 TIN 的数据存储方式比格网 DEM 复杂，它不仅要存储每个点的高程，还要存 储其平面坐标、节点连接的拓扑关系，三角形及邻接三角形等关系。TIN 模型在 概念上类似于多边形网络的矢量拓扑结构，只是 TIN 模型不需要定义“岛”和“洞” 的拓扑关系