第7章 数字地面模型

7.1数字地面模型概述7.2DEM常用数据模型7.3DEM数据的获取7.4数字地形可视化7.5数字地形分析7.6DEM的应用,第七章数字高程模型及其应用,地理空间实质是三维的，只是人们通常在二维地理空间上描述并进行分析。如在土地利用，土地分级等问题上，都用平面专题图来描述。数字地面模型的提出，从时间上实际上早于GIS，但GIS的发展大大促进人们对数字地面模型的研究。,目前，数字地面模型已成为GIS的重要内容，GIS的很多功能以数字地面模型为基础。,7.1数字地面模型概述,DTM(DigitalTerrainModel)数字地形模型50年代由MIT摄影测量实验室提出，是用数字形式描述地形表面的模型。实质上这是对地面形态和属性信息的数字表达。DEM(DigitalElevationModel)数字高程模型当DTM模型中数字属性为高程时称DEM模型，即数字高程模型。DEM模型是DTM模型的一种特例。,DEM模型是新一代的地形图，它通过存储在介质上的大量地面点空间数据和地形属性数据，以数字形式来描述地形地貌。,为了表示地形起伏必需存储三维数据，这首先必需研究三维数字地面模型。,从测绘的角度看,什么是4D产品（DEM,DLG,DRG,DOM）,数字高程模型（DigitalElevationModel，缩写DEM）是在某一投影平面（如高斯投影平面）上规则格网点的平面坐标（X，Y）及高程（Z）的数据集。DEM的格网间隔应与其高程精度相适配，并形成有规则的格网系列。根据不同的高程精度，可分为不同类型。为完整反映地表形态，还可增加离散高程点数据。数字线划地图（DigitalLineGraphic，缩写DLG）是现有地形图要素的矢量数据集，保存各要素间的空间关系和相关的属性信息，全面地描述地表目标。数字栅格地图（DigitalRasterGraphic，缩写DRG）是现有纸质地形图经计算机处理后得到的栅格数据文件。每一幅地形图在扫描数字化后，经几何纠正，并进行内容更新和数据压缩处理，彩色地形图还应经色彩校正，使每幅图像的色彩基本一致。数字栅格地图在内容上、几何精度和色彩上与国家基本比例尺地形图保持一致。数字正射影像图（DigitalOrthophotoMap，缩写DOM）是利用数字高程模型（DEM）对经扫描处理的数字化航空像片，经逐像元进行投影差改正、镶嵌，按国家基本比例尺地形图图幅范围剪裁生成的数字正射影像数据集。它是同时具有地图几何精度和影像特征的图像，具有精度高、信息丰富、直观真实等优点。,DEM的表示方法,DEM表示方法,数学方法,图形法,整体,局部,傅立叶级数,高次多项式,规则数学分块,不规则数学分块,点数据,线数据,规则,不规则,水平线,典型线,典型特征,密度一致,密度不一致,三角网,邻近网,山峰、洼坑,隘口、边界,垂直线,山脊线,谷底线,海岸线,坡度变换线,7.2DEM的常用数据模型,规则格网模型将空间区域分成规则的等距离单元，每个单元对应一个数值，通常在数学上表示为一个矩阵，在计算机中表现为一个二维数组，每个格网单元或数组元素对应一个高程值。用规则采样点数据（或把不规则采样点数据内插成规则点数据），而后，以矩阵形式来表地面形状。它已成为栅格数据结构中DEM的通用形式。,一、规则格网(grid)模型,按平面上等间距规则采样，或内插所建立的数字地面模型，称为基于栅格的数字地面模型，可以写成以下形式:DTMZi,j,i=l,2,m;j=1,2,n式中，Z为栅格结点（i,j）上的地面属性数据，包括土地权属、土壤类型、土地利用等。当该属性为海拔高程时，则该模型即为数字高程模型。,1.模型的表示,规则格网（grid)模型,1）格网面元组成格网的四个相邻格点在水平面上所包含的面积单元。2）格网面元的趋势面格网面元的四个角点高程支撑的数学面，通常该数学面用三种形式表示。按最小二乘法将格网面元四个角点高程拟合为一个平面，称格网面元的平面趋势面。将格网面元四个角点高程拟合为双线性趋势面。将格网面元四个角点高程拟合为双三次趋势面。3)格网点针对格网点的值,2.格网的含义,格网数据结构是典型的栅格数据结构，可采用栅格矩阵及其压缩编码的方法表示。其数据包括三部分：1）元数据描述DEM数据的数据，如数据表示的时间、边界、测量单位、投影参数。2）数据头DEM数据的起点坐标、坐标类型、格网大小、行列数等。3）数据体行列数分布的数据阵列。,3.格网数据结构,规则格网数据模型的优点：1）数据结构简单，算法实现容易，便于空间操作和存储。尤其适合在栅格数据结构的GIS系统中。2）容易计算等高线、坡度、坡向、自动提取地域地形等。规则格网是DEM最广泛使用的格式。目前，很多国家都以规则格网的数据矩阵作为DEM提供方式。,4、规则格网的优缺点,规则格网数据模型的缺点：1）数据量大，通常采用压缩存储无损压缩存储，如游程编码、链码、四叉树编码；有损压缩存储，如离散余弦（DCTDiscreteCosineTransformation),小波变换（WaveletTransformation)2）不规则的地面特性与规则的数据表示之间本身就不协调。它对不同地形采用一律平等的规则格网，不利于表示复杂地形。,由于受观测手段所限，或专业要求，在实际中获取的数据常不是规则格网数据，大多为不规则的离散数据。如地震观测中观测的地层结构数据，水利中观测的地下水资源数据等。,二、不规则三角网(TIN)(TriangulatedIrregularNetwork),不规则三角网（TIN）模型通过不规则分布的数据点，生成连续的三角形面来逼近地面的地形表面。三角形面的形状和大小取决于不规则分布的观测点数据的位置和密度。TIN将区域内有限的点集，连成相互连续的三角面，使离散点为三角面的顶点。,1、不规则三角网表示,不规则三角网模型是一种三维空间的分段线性模型，其数据格式在概念上类似二维数据结构中带拓扑结构的矢量数据结构，只是TIN模型中不定义岛和洞的拓扑关系建立。,不规则三角网（TIN）模型把不规则分布的数据点，按优化组合的方法，生成连续的三角形面来逼近地面的地形表面。使每个离散点为三角面的顶点。即TIN将区域内有限的点集划分为相连的三角面网。,2、模型的表示,不规则三角网（TIN）模型,生成不规则三角网的算法要求：1）连接时尽可能使三角形的三个边长接近，成锐角三角形；2）为保证从最邻近点生成不规则三角形，使三角形边长之和为最小。为此，已知三角形一边之后，选择到该一边两端点距离之和最小的一点作为三角形的另一个点；3）生成三角形的任意一边不能和地性线相交，不能和已生成的三角形的任意一边相交。为此，采样点要合理，它应能较好的反映地形真实情况。如对特殊的地形线，如山谷线、山瘠线、断裂线处，出现不完全反映真实情况时，要剔除这类三角形进行调整处理。在众多不规则三角网生成算法中，Delaunay算法用得最多。,3、不规则三角网的生成,1）逐步生成法以任一点为起始点；找出离起始点最近的数据点连线，作为起始基线；按联三角形法则找出第三点连成三角形；再以连成三角形的新边作基线重复上述步骤，直到终止。,Delaunay算法的三种类型,2）分割归并法（分而治之法）将点集分到足够小，使其成为易生成三角网的子集，然后将子集合并。3）逐点插入法先在包含所有数据点的多边形中建立初始三角网，然后将余下的点逐一插入。建立初始三角网的方法可能不同。,点文件三角形文件点号坐标点三角形号顶点邻接三角形N1X1Y1Z1T1N1N5N6T2T5/N2X2Y2Z1T2N1N4N3T1T3T6N3X3Y3Z1T3N1N2N4/T4T2N4X4Y4Z4T4N2N3N4T3/T8N5X5Y5Z5T5N8N5N6T1/T6N6X6Y6Z6T6N4N5N8T2T5T7N7X7Y7Z7T7N4N7N8T6T8/N8X8Y8Z8T8N3N4N7T4T7/,T1T2T3T4,T5T6T7T8,N1,N6,N2,N8,N7,N3,N4,N5,4、模型的数据结构,不规则三角网模型的优点：1）克服栅格数据中的数据冗余问题；2）表示地面形态效率高，数据精度高。它能较好地表示地性线，充分表示复杂的地形特征，适应起伏不同的地形。不规则三角网模型的缺点：1）算法实现复杂，由于形成三角网方法不同有不同算法；2）对特殊的地形线要调整。,5、不规则三角网模型的优缺点,大比例尺数据高程模型通常采用能表示地性线的不规则三角网，以便较精确地显示小区域地形特性。,小比例尺数据高程模型通常可采用规则格网模型，以显示大区域宏观地形特性。,等值线是DEM模型的平面表示形式，是地形表示中广泛使用的一种表示方法。等值线图以符号化的模型来表示空间立体地形形态。它由一系列数值相等的点，连成的曲线反映连续递变的面状分布地形特征。它反映的是连续递变的面状分布特征，如等高线、等温线等。,三、等值线模型,等值线是指X-Y平面上f(x,y)=c的轨迹分布线。这里的c为某一常数值，该值所表示的物理意义可以是地形高程数据，温度场中的温度数据，气象上气压的数据等等。等值线符合下述要求：给定值的等值线在相应域内不能互相交错；一根等值线通常是一条连续曲线；给定值后，相应域上等值线不限于一条；等值线可以是闭合曲线，也可以和域外连续。,1、等值线图的概念,等高线是通过对地球表面进行一系列不同高程的水平切割后，切割面和地球表面产生的一系列交线。它由一系列等高线和其高程值构成，其上的每根等高线有一已知的高程值。每条等高线是一组有序的坐标点序列，也可认为是带有高程属性的多边形或弧段。由于生成等高线时丢失了大量地表信息，所以用等高线重建地貌形态时只能近似表示。,2、等高线,等高线模型由一系列等高线和其高程值构成；每条等高线是一组有序的坐标点序列，也可认为是带有高程属性的多边形或弧段。等高线模型只能表达区域内部分高程值，等高线外的高程值用其外包等高线的高程值插值得到。,等高线离散化，等高线内插可生成格网数据。,等高线可生成TIN数据,3、等高线的数据组织结构,等高线的数据组织结构例,DEM的常用数据模型的比较,1）航空航天遥感影像为数据源大比例尺数字高程模型用航空影像；小比例尺数字高程模型用航天影像。通过立体象对，用摄影测量的方法建立空间立体模型量取密集的高程数据。2）以地形图为数据源主要用比例尺不大于己于1:1万的国家近期地形图为数据源从中量取等密度地面点集的高程数据，建立DEM。数据采集通常用手工、数字化仪几及扫描仪。,7.3DEM数据的获取,1、DEM数据源的种类,3）以地面实测记录为数据源对小范围的大比例尺（如1:5000）的DEM数据可用GPS，电子测速仪（全站仪），测距经纬仪等获取数据。4）其它数据源数字摄影测量航空测空仪可获取精度要求不很高的DEM数据，近景摄影测量，在地面摄取立体象对，构造解析模型，也可获的小区域的DEM数据。,由于实际地形无一定数学规律可循，因此，影响DEM精度的主要因素是原始数据的获取，其中主要包括:数据采集的密度;数据采集位置（选点）。注意:任何一种内插方法均不能弥补由于取样不当所造成的信息损失。,2、DEM数据源的采集,1)数据采集的精度,2)数据采集的方法,(1)人工网格法将地形图蒙上格网，逐格读取中心或角点的高程值、构成数字高程模型。(2)立体像对分析通过遥感立体像对，根据视差模型、自动选配左右影像的同名点，可建立数字高程模型。,(3)三角网方法（TIN）,(4)曲面拟合,根据有限个离散点的高程、采用多项式或样条函数求得拟合公式，再逐一计算各点的高程，可得到拟合的DEM。,对有限个离散点，每三个邻近点联接成三角形，每个三角形代表一个局部平面，再根据每个平面方程，可计算各格网点高程，生成DEM。,三角网上的任意数据获取,点数据的生成与叠加,(5)等值线插值,根据各局部等值线上的高程点，通过插值公式计算各点的高程，得到DEM。等值线插值法是比较常用的方法，输入等值线后，可在矢量格式的等值线数据基础上进行，插值效果较好。,（1）采集前根据建立DEM模型的精度要求确立合理的取样密度，如对单调地形可均匀采集，密度不必过大；（2）地形变化处采用密集采样，确保地形转折处的数据如山谷线，山脊线，谷缘线，崖线，山坡转折线等；（3）不应出现大的空白区，即使大片平坦地区也应有最底的采集密度。,3）数据采集的原则,采集的DEM数据必需进行后处理方可应用，后处理包括：(1)格式转换由于采集数据所存储格式包括数据内容、数据类型等各异，因此必须进行数据格式转换，使数据成为自己所要的格式。(2)坐标转换有时采集数据要进行坐标系转换，如将像片坐标转换成大地坐标。,4)数据采集的后处理,(3)数据编辑处理任何数据获取后总要编辑，编辑过程通常是一种交互过程，主要包括剔除错误的、过密的、重复的点；加密要加密区域的点。(4)数据分块由于采集数据的不同，数据常具有不同排列顺序，如利用地形图采集的等高线数据通常按条采集和存储查寻到，而对等高线进行区域插值时要用待插点周周的数据为依据，为迅速查询到待插点周周的数据，需要将数据重新分块存储，分块方法通常是将整个区域分成等间隔的格网，每个格网之间有一定的重叠度，相互之间可用链指针连接。,1)网格尺寸的确定网格尺寸的确定是DEM数据网格化的重要问题，它关系到数据模型的精度。网格尺寸的确定首先应分析地形的形态特征。一般地说，采样点的密度基本上确定了网格点的密度。在采集优选点的情况下，网格点数应大于或接近于样点数，取n0ANDZi,j+1Zi,jANDZi+1,jZi,j2)脊点满足下面4个条件(Zi,j-1-Zi,j)*(Zi,j+1-Zi,j)0AND(Zi,j-1-Zi,j)*(Zi+1,j-Zi,j)0ANDZi,j+1Zi,jANDZi+1,jZi,j,i,j-1i,ji,j+1,i-1,j-1i-1,ji-1,j+1,i+1,j-1i+1,ji+1,j+1,1、山谷点、山脊点的计算,2、沟谷密度,沟谷密度是表征地面破碎程度的一种指标，它是沟谷总长度(L)与地表单元总面积(A)之比。提取谷点和脊点，将地表单元内所有谷点在单元区域内的延伸长度累加，便获得单元的沟谷长度。沟谷密度为:,反映某一面积单元内地势起伏变化与侵蚀程度的复杂程度，是地表面积与投影面积之比：,3、地表粗糙度,计算辐照度需考虑日照条件（太阳赤纬、高度角、时角及大气状况）与坡面几何条件的相互关系由下式决定：其中:大气透过率，与太阳高度和大气状况有关；Sc为太阳常数；Sa为太阳高度角可由球面三角公式求出；t是时角；a、b为坡面方程系数；为坡度。,4、地表辐照度,1、地形特征的提取地形特征包括地形特征点、地形特征线、地形特征面。地形特征点：山峰点、谷低点、鞍部点地形特征线：山谷线、山瘠线、地形特征面：凹凸面2、水系特征的提取水系特征如分水岭、汇水流域,大部分地形特征的提取算法是基于格网的DEM,四、特征的提取,四、DEM模型的可视化分析,可视性分析也称通视分析，被广泛地应用在军事上、森林火灾监测点的设置，无线发射塔如微波站、广播电台、电视塔的设置等方面。,A,B,AB,1、剖面分析,地形的剖面分析可用来表示区域的地貌形态、地势变化、地质构造等等。地形的剖面分析中常画剖面图，画剖面图的关键是求剖面线。首先通过直线AB向地平面作一垂直于XOY的面S，而后再求S面同地形模型面交线，即剖面线，最后用剖面线来表示剖面,剖面线分析法在观察点A和目标点B之间画线，并作于XY平面的垂直面S，求出地形模型面与S面的交线判别AB直线与剖面线是否相交，如不相交，AB两点之间可视；如相交，AB两点之间不可视。,A点和B点之间的通视性（intervisibility)，指在A点是否能看到B点。,2、通视性分析,可视域指某观察点可视的区域范围。可视域分析的应用极为广泛。1）基于格网数据结构的可视域分析可视域分析中用得最多的是用格网数据结构，即用“可视