gis第7章教案.docx
- 文档编号:4279552
- 上传时间:2022-11-28
- 格式:DOCX
- 页数:6
- 大小:25.77KB
gis第7章教案.docx
《gis第7章教案.docx》由会员分享,可在线阅读,更多相关《gis第7章教案.docx(6页珍藏版)》请在冰豆网上搜索。
gis第7章教案
第七章数字地形模型及其应用学习目标〃理解DTM和DEM概念〃掌握DEM的主要表示方法,并能分析其优缺点〃了解DEM的数据采集方法〃了解DEM的插值方法〃掌握DEM建立的一般作业过程〃了解DEM在地学分析中的主要用途〃了解地形可视化的方法数字地形模型(DTM)是20世纪50年代由美国MIT摄影测量试验室主任米勒(C.L.MILLER)首次提出,并用其成功地解决了道路工程中土方估算等问题。
此后它被用于各种线路选线(铁路、公路、输电线等)的设计以及各种工程的面积、体积、坡度计算,任意两点间的通视判断即任意断面图绘制。
在测绘中被用于绘制等高线、坡度坡向图、立体透视图,制作正射影像图以及地图的修测。
DTM是地理信息系统的基础数据,可用于土地利用现状的分析、合理规划和洪水险情预报等。
在军事上可用于导航、作战电子沙盘等。
第一节数字高程模型的概念和表示一、DTM和DEM的概念DTM是描述地形表面形态空间位置和地形属性分布的有序数值阵列。
DTM以离散分布的平面点来模拟连续分布的地形。
数字地形模型中地形属性为高程时称为数字高程模型(DEM)。
不过,被描述的地形属性也可以是地理空间上的地价、污染负荷量、绿化率、降雨量、气温、人口密度等。
显然,DEM是DTM的一个子集,是DTM的一个特例。
二、DEM的表示法区域地表高程模型可以采用多种方法表达,如图7-1所示。
1.数学方法用数学方法来表达,可采用整体拟合表面方法,即根据区域所有的高程点数据,用傅立叶级数和高次多项式拟合统一的地面高程曲面。
也可用局部拟合方法,将地表复杂的表面分成正方形规则区域或面积大致相等的不规则区域进行分块搜索,根据有限个点进行拟合形成高程曲面。
126
傅立叶级数整体高次多项式数学方法规则数学分块局部不规则数学分块密度一致规则DEM表示方法密度不一致三角网点数据不规则邻近网山峰、洼坑图形法典型特征隘口、边界水平线垂直线山脊线线数据谷底线典型线海岸线127
坡度变换线图7-1DEM的表示方法(据邬伦等,2001年)2.图形方法
(1)线模式等高线是表示地形最常见的形式,其地形特征线也是表达地面高程的重要信息源,如山脊线、谷底线、海岸线及坡度变换线等。
(2)点模式①高程矩阵:
DEM最普通的形式是高程矩阵或规则矩形格网(GRID),规则网格通常是正方形、矩形、三角形等规则网格。
规则网格将区域空间切分为规则的格网单元,每一个格网单元对应一个数值。
数学上可以表示为一个矩阵,在计算机实现中则是一个二维数组。
每个格网单元或数组的一个元素对应一个高程值。
规则格网的高程矩阵,其优点是用计算机处理很容易;缺点是:
一是在地形简单区存在大量冗余数据;二是如不改变格网大小,无法适用于起伏复杂程度不同的地区;三是对于某些特种计算如通视计算,过分依赖格网轴线;四是不能准确表示地形的结构和细部。
不规则的离散采样点可以按照两种方法产生高程矩阵:
一是将规则格网覆盖在这些数据的分布图上,然后用内插技术产生高程矩阵,当然,内插技术也可用来从一个粗略的高程矩阵产生更精细的高程矩阵;二是把离散采样点作为点模式中不规则三角网系统的基础。
②不规则三角网(TIN)不规则三角网是专为产生DEM数据而设计的一种采样表示系统。
TIN模型根据区域所有采样点取得的离散数据,按照优化组合的原则,把这些离散点(各三角形的顶点)连接成相互连续的三角面,在连接时尽可能地使每个三角形为锐角三角形或为三边的长度近似相等,如图7-3所示。
区域中任意点落在三角面的顶点、边上或三角形内。
如果点不在顶点上,在边上用边的两个顶点128
的高程,在三角形内的则用三个顶点的高程,该点的高程值通常通过线性插值的方法得到。
所以TIN是一个三维空间的分段线性模型。
图7-3不规则三角网示例TIN的数据存储方式比GRID复杂,它不仅要存储每个点的高程,还要存储其平面坐标、节点连接的拓扑关系,三角形及邻接三角形等关系。
TIN模型在概念上类似于多边形网络的矢量拓扑结构,只是TIN模型不需要定义“岛”和“洞”的拓扑关系。
不规则三角网数字高程由连续的三角面组成,三角面的形状和大小取决于不规则分布的测点,或节点的位臵和密度。
不规则三角网与高程矩阵方法不同之处是随地形起伏变化的复杂性而改变采样点的密度和决定采样点的位臵,因而它能够避免地形平坦时的数据冗余,又能按地形特征点如山脊、山谷线等表示数字高程特征。
第二节数字高程模型的建立为了建立DEM,必须量测一些点的三维坐标,这就是DEM数据采集。
一、DEM的数据源与采集方法1.摄影测量法摄影测量法以航空或航天遥感图像为数据源,利用遥感立体像对,采用摄影测量的方法建立空间地形立体模型,量取密集数字高程数据,建立DEM。
采集数据的摄影测量仪器包括附有自动记录装臵的立体测图仪或立体坐标仪、解析测图仪及数字摄影测量系统。
用摄影测量方法生产DEM,数据点的采样方法因产品的要求不同而异。
沿等高线、断面线、地性线进行采样往往是有目的的采样,而许多产品要求高程矩阵形式,所以基于规则网格或不规则网格点的面采样是必须的,这种方式与其它空间属性的采样方式一样,只是采样密度高一些。
涉及DEM数据采集的摄影测量采样方法如下:
(1)沿等高线采样在地形复杂及陡峭地区,可采用沿等高线跟踪方式进行数据采集,而在平坦地区,则不宜沿等高线采样。
129
(2)规则格网采样利用解析测图仪在立体模型中按规则矩形格网进行采样,直接构成规则格网DEM。
(3)选择采样为了准确反映地形,可根据地形特征进行选择采样,例如沿山脊线、山谷线、断裂线等进行采集以及离散碎部点(如山顶)的采集。
这种方法获取的数据尤其适合于不规则三角网DEM的建立。
(4)渐进采样这种采样方法的目的是使采样点分布合理,即平坦地区样点少,地形复杂区的样点较多。
(5)混合采样为了同步考虑采样的合理性和效率,可将规则采样或渐进采样与选择采样结合进行混合采样,即在规则采样的基础上再进行沿特征线、点采样。
为了区别一般的数据点和特征点,应当给不同的点以不同的特征码,以便于处理时可按不同的方式进行。
利用混合采样可建立附加地形特征的规则格网DEM,也可建立附加特征的不规则三角网DEM。
(6)自动化DEM数据采集上述方法均是基于解析测图仪或机助制图系统利用半自动的方法进行DEM数据采集,现在已经可以利用自动化测图系统进行完全自动化的DEM数据采集。
此时可按相片上的规则格网利用数字影像匹配进行数据采集。
最后摄影测量获取的DEM数据点都要按一定插值方法转成规则格网DEM或规则三角网DEM格式数据。
摄影测量方法是空间数据采集最有效的手段,它具有劳动强度低、效率高的优点。
数据采样可以全部由人工操作,通常费时且容易出错;半自动采样可以辅助操作人员进行采样,以加快速度和改善精度,通常是由人工控制高程Z,由机器自动控制平面坐标X,Y的驱动;全自动方法利用计算机视觉代替人眼的立体观测,速度虽然快,但精度较差。
2.地形图数字化法这种方法以比例尺大于1:
1万的国家近期地形图为数据源,从中量取中等密度地面点集的高程数据,建立DEM。
其方法有下面几种:
(1)手工采集法采用方格膜片、网点板或带刻划的平移角尺叠臵在地形图上,并使地形图的格网与网点板或膜片的格网线逐格匹配定位,自上而下,逐行从左到右量取高程。
当格网交点落在相邻等高线之间时,用目视线性内插法估计高程值。
它的优点是几乎不需要购臵仪器设备,而且操作简便。
(2)手扶跟踪数字化仪采集法可采用如下3种方式进行采集:
•逐条等高线的线方式连续采集样点,并采集所有高程注记点作补充,这种方式适合于等高线较稀疏的平坦地区。
•沿主要等高线采集特征点,并选择采集高程注记点和线性加密点作补充。
130
•沿曲线和坡折线采集特征点,并补采峰—鞍线和水边线的支撑点,分别以等高线,峰—鞍链和边界格式存储。
(3)扫描采集法这种方法利用扫描仪扫描地图,从地图扫描数据中自动地建立DEM。
这种方式采集速度快。
随着研究的不断深入,一些难点和瓶颈问题被逐步解决,目前自动建立DEM技术已达到实用水平。
3.地面测量法这种方法以地面实测记录为数据源,利用GPS、全站仪和电子手簿或测距经纬仪等设备,在已知点位的测站上,观测到目标点的方向、距离和高差3个要素。
计算出目标点的(x,y,z)三维坐标,存储于电子手簿或袖珍计算机中,成为建立DEM的原始数据。
这种方法一般用于建立小范围大比例尺区域的DEM,对高程的精度要求较高。
二、DEM的空间插值方法DEM数据的内插就是在一个由平面坐标(x,y)构成的二维空间中,由已知若干离散点P的高程,估算待内插点的高程值。
由于DEM采样的数据点呈i离散分布形式,或是数据点虽按格网排列,但格网的密度不能满足使用的要求,这就需要以数据点为基础进行插值运算。
DEM内插按插点分布范围,可分为分块内插、剖分内插和单点移面内插三类,如图7-4所示。
趋势面拟趋势面拟合内插傅立叶级数内插多次多项式内插二元样条函数内插分块内插孔斯曲面内插特殊曲面内插B样条曲面内多层叠加插值面函数内插塞尔曲面内插最小二乘配臵法内插131
有限元法(加权残差法)内插DEM空间内插多项式内插剖分内插多项式拟合内插移动拟合法单点移面内插加权平均法Kriging内插法图7-4DEM空间内插方法分类分块内插,是把需要建立DEM的地区,切割成一定大小的规则方块,形状通常为正方形,它的尺寸应根据地形复杂程度和数据源的比例尺确定。
在每一个分块上展铺一张数学面,相邻分块之间有适当宽度的重叠带,以使重叠带内全部数据点成为相邻块展铺数学面时的共用数据,保证一张数学面能够较平滑地与相邻分块的数学面拼接。
这种内插方法的优点是可以得到光滑连续的空间曲面。
剖分内插是把需要建立DTM的地区切割成大小和形状不同的子区(剖分),子区间拥有公共边但不重叠,在该区内展铺一个数学面,内插剖分区内任意点的高程。
该法只在剖分间边界端点处重合,通常没有严格重合的边界,所以既不连续,也不光滑。
剖分多边形的顶点都是数据点,最常见的数据点个数为3,与TIN结构相同。
单点移面内插是以待插点为中心,以适当半径或边长的圆或正方形作为移动面去捕捉适当数目的数据点,并以此展播一张数学面,内插该中心的高程。
如图7-5(a)所示分块内插区域;如图7-5(b)所示剖分内插区域;如图7-5(c)所示单点移面插值区域。
132
(a)(b)(c)图7-5DTM空间内插三、DEM的生成流程DEM生成的全过程包括:
原始数据获取、DEM模型构造、数据插值、在所定数据结构支持下的数据存储和模型输出。
如图7-6所示,以地形图数字化为例,说明DEM的生成过程。
全要素地形图地图扫描人机交互等高线矢量化加测注记点等高线赋值与检查133
周边等高线地图的数据获取构造三角网内插DEM网格DEM建库与刻盘质量检测与元数据文件记录从等高线数据可以直接生成TIN,也可直接生成格网DEM,另一方面,格网DEM也可由等高线先生成TIN再内插获得。
经过实践证明,由等高线先生成TIN再内插格网DEM的精度和效率都是很好的。
第三节DEM分析与应用数字高程模型是各种地学分析、工程设计和辅助决策的重要基础性数据,有着广泛的应用领域。
在地学分析中,用于自动提取各种地形因子,制作地形剖面图和划分地表形态类型;在工程设计中,可用于各种线路的自动选线、水库大坝的选址,以及土方、库容和淹没损失的自动估算等;通过与各专题数据的匹配分析,还可进行遥感影像地形畸变的自动校正,以及进行农业、林业、土地等的规划研究。
一、基本地形因子计算1.极值高程和最大高差2.相对高程和平均高程3.坡度和坡向计算地面上某点的坡度是表示地面在该点倾斜程度的一个量,因此,它是一个既有大小又有方向的量,即矢量。
坡度定义为水平面和地形表面之间夹角的正134z
切值;坡向为坡面法线在水平面上的投影与正北方向的夹角如图7-7所示。
坡度和坡向的计算一般采用拟合曲面法,通常在3×3个DEM格网窗口中进行,如图7-8所示。
每个窗口中心为一个高程点。
eee526eee13eee847图7-7地表单元坡度和坡向示意图图7-83×3窗口计算点的坡度和坡向设栅格的边长为ΔS,点e的x、y方向的坡度分别为:
Gx=(e-e)/(2×ΔS)13Gy=(e-e)/(2×ΔS)42点e的坡度和坡向求解公式如下:
坡度的计算公式:
G=arctan(Gx+Gy)221/2坡向的计算公式:
A=arctan(Gy/Gx)坡向图是坡向的类别显示图,斜坡的倾斜方向可取方位角0º~360º中的任意方向。
坡向一般分为9类:
东、南、西、北、东北、西北、东南、西南和平地。
在实际应用中,可以综合为四种坡向,即平缓坡、阳坡、半阳坡和阴坡,分别用1、2、3和4表示。
如图7-9所示。
NNWNE43EW32SWSES图7-9坡向的综合表示4.表面积的计算根据数字高程模型很容易求得地表面积,其计算可看作是所包含各个网络的表面积之和。
5.投影面积的计算135
6.体积的计算7.剖面积的计算8.地表粗糙度计算地表粗糙度是反映地表的起伏变化与侵蚀程度的指标,一般定义为地表单元的曲面面积与投影面积之比:
CZ=S/S表面积投影面积显然,这种定义对于表面光滑的斜面也可求出不同的粗糙度,这是不合适的。
这里用对顶点连线L与L中点的距离D来表示粗糙度,如图7-10所示,12D值愈大,说明单元的四个顶点的起伏变化也愈大。
计算公式为:
1RDZZZZi,ji1,j1i,ji,j1i1,j2ZZZi,ji,j-1L2DZi+1,jZi+1,j+1OYX图7-10地表粗糙度计算(据黄杏元等,2001年)二、谷脊特征分析在地表的基本形态中,山谷和山脊是常见的两种主要形态。
它在区域地形研究和制图综合中具有重要的意义。
利用数字高程模型可对谷脊特征作概略分析。
1.谷点和脊点的判定谷点是地势相对最低的点集,脊点为地势相对最高的点集如图7-11所示,要判定高程为Z网格的形态特征,按照以下判别式可直接提取谷点和脊i,j点。
)()0如果(ZZZZ>i,j1i,ji,j1i,jZ当Z时,则>i,j1i,j1=………………………………………………………①P(i,j)Z当Z时,则>i,j1i,j1=………………………………………………………②P(i,j)如果(ZZ)(ZZ)>0i,ji1,ji,ji1,jZ当Z时,则>i1,ji,j1P(i,j)=………………………………………………………③Z当Z时,则>i1,ji,j1P(i,j)=…………………………………………………………④136
如果①和④或②和③同时成立,则=P(i,j)2如果以上条件均不成立,则=P(i,j)0其中:
1表示谷点1表示脊点=P(i,j)2表示鞍点0表示其他点Zi+1,jZZZi,j-1i,ji,j+1Zi-1,j图7-11谷点和脊点的判定若对谷脊特征作精确分析时,可以建立地表单元的拟合曲面方程,然后确定曲面上各插值点的极值,以及当插值点在两个相互垂直方向上,分别为极大值和极小值时,确定出谷点、脊点和鞍点。
2.沟谷密度分析沟谷密度是表征地面破碎程度的一种指标,它是沟谷总长度(∑L)与地表单元总面积(∑A)之比.提取谷点和脊点,将地表单元内所有谷点在单元区域内的延伸长度累加,便获得单元的沟谷长度.沟谷密度为:
nmK1,2,3,,ni1,2,3,,mVDL/Akik1i1沟谷密度的精度与插值间隔有关.提取谷脊信息时给定的插值间距小,则一般沟谷密度的计算精度高.3.切割深度分析.地表单元的谷点与最近脊点的平均高差为谷点的切割深度,区域平均切割深度为若干谷点切割深度的平均值。
nK1,2,3,,nVH(hh)/nRKVKK1式中hh为距该谷点最近的脊点的平均高程值,为谷点高程。
RV三、地形剖面分析地形剖面可以概括而直观地研究区域的地势、地质和水文特征,也可以与其他地理变量相叠加,提供农业生产布局的立体背影分析、土地利用规划以及工程选址选线等的参考依据。
从DEM可以很方便地制作任一方向上的地形剖面图。
只要知道所绘剖面线在数字高程模型中的起点位臵和终点位臵,就可以137
惟一地确定这条剖面线与DEM网格各个交点的平面位臵及其高程,根据选定的垂直比例尺和水平比例尺,就能自动绘出所需要的剖面图。
可用梯形法、辛普森法来计算剖面面积。
当两个剖面面积相互平行时,可用类似的方法计算该两个剖面之间所夹的体积。
四、淹没损失估算DEM数据不仅直接用于各种地形因素的分析,而且还可以与有关信息进行复合,研究地形要素与其他要素之间的相互联系。
其中洪水淹没损失估算,就是研究DEM与土地利用之间的关系。
为了科学合理地估价,首先将数字高程的数据与土地利用数据进行匹配;其次是根据淹没高程,在DEM上确定淹没范围;最后统计淹没范围内的土地类型和面积,就能精确估算出淹没损失。
五、地貌类型的自动划分在DEM数据文件的基础上,进行地貌类型的自动分类。
首先根据区域的地形特点,拟定地形分类的高程界值;然后,计算机根据确定的高程界值自动提取地形类型信息,便可获得区域的地形分类系统,如平原、丘陵、低山、中山和高山等;最后输出地貌类型图。
六、从DEM数据自动形成地形轮廓线高程矩阵没有存储山脊线、山谷线等地形特征线,或者在地形图数字化时,对这些地形特征线没有单独数字化,在这种情况下,用程序自动地将它们从高程矩阵中提取出来也许是必要的。
例如从叠臵到DEM的卫星图像上勾绘出集水范围线,使遥感图像与特殊地理景观联系在一起。
高程矩阵用于其他数量分析如费用量、集水范围、旅行时间等时,应有一种方法来描述线、面特征。
1.山脊线和谷底线的探测为了自动探测山脊线和谷底线,设计了专门的运算算子。
较为简单的算子是4个像元的局部算子。
该算子在高程矩阵中移动并比较每一位臵处4个像元的高程值,同时标出其中高程最大(探测谷底线)或最小(探测山脊线)的像元。
标记过程完成后,剩下为标记的像元就是山脊线和山谷线所在的像元。
下一步就是把它们连接成线模式,形成山脊线或山谷线。
2.集水范围的确定集水范围即流域范围的确定对流域分析十分必要。
流域探测除确定边界线外,还要将整个范围从整个数据库中分离出来。
探测方法是:
首先需要交互式地确定河流流域的出口,并作为搜索工作的起点。
以3×3算子的中心像元臵于起始点,比较中心像元相邻近的8个像元的坡向。
如果坡向朝向中心像元,则认为它是中心像元的上游,算子的中心像元移至新的“上游”点,重复比较过程又能得到新的“上游”点。
以由“上游”标志的不予比较,整个数据范围都运算完毕后,流域范围就全部标记出来了。
用户可以对这些像元重新编码,形成某一流域的分布图。
七、通视分析通视分析也称可视分析,它实质属于地形进行最优化处理的范畴,如设臵138
雷达站、电视台的发射站、道路选择、航海导航等,在军事上如布设阵地(炮兵阵地、电子对抗阵地)、设臵观察哨所、铺架通信线路等。
数字高程模型的建立为这类分析提供了有利基础,能方便地算出一个观察点所能看到的各个部分。
在DEM中辨认出观察点所在的位臵,从这个位臵引出一簇射线,比较射线通过的每个点(高程矩阵中即为像元)的高程,将不被物体隐藏的各点进行特殊编码,便可得到一幅简单的地图。
ABHDhCOETOvvss′′图7-12通视计算示意图(距李志林等,2001年)八、地形可视化可视化是指人脑中形成对某种事物(人物)的图像,是一个心智处理过程,促进对事物的观察及概念的建立。
地形可视化是根据地表高程变化或加上地表覆盖物种类而建立的反映地物空间起伏分布的视觉模型。
地形的可视化方法主要有:
1.根据DEM制作等高线图利用DEM绘制等高线图是以格网点高程数据或者将离散的高程数据由栅格追踪法原理转换为矢量等值线所产生的如图7-13所示。
该方法可以适用于所有的利用格网数据方法绘制等值线图。
322328262435244338343435403241524348473130442735594726362230373143482423341627431737112519112328312.根据DEM制作透视立体图立体图是表现物体三维模型最直观形象的图形,它可以生动逼真地描述制图对象在平面和空间上分布的形态特征和构造关系。
通过分析立体图可以了解139
地理模型表面的起伏状况(图7-14所示),可以看出各个断面的形态,这对研究区域的轮廓形态、变化规律以及内部结构是非常有益的。
计算机自动绘制透视立体图的理论基础是透视原理,而DEM是其绘制的数据基础。
制作透视立体图的基本步骤包括:
建立透视变换基础、DEM高程阵列剖面布设、消除隐藏线处理和粘贴表面影像与纹理如图7-14所示。
3.根据DEM制作晕渲图晕渲图是通过模拟实际地面反映地形起伏特征的重要地图制图学方法,在各种小比例尺地形图、地理图,以及有关专题地图上得到非常广泛的应用。
但是,传统的人工描绘晕渲图的方法费工、费时、成本较高。
利用DEM数据作为信息源,以地面光照通量为依据,计算该栅格所输出的灰度值,由此产生的晕渲图具有相当逼真的立体效果如图7-15所示。
自动地貌晕渲图的计算步骤为:
首先是根据DEM数据计算坡度和坡向,然后将坡向数据与光源方向比较,向光源的斜坡得到浅色调灰度值,反方向的斜坡得到深色调灰度值,介于中间坡向的斜坡得到中间灰度值。
灰度值的大小则按坡度进一步确定。
140
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- gis 教案
![提示](https://static.bdocx.com/images/bang_tan.gif)