地形初步分析.docx
- 文档编号:11716369
- 上传时间:2023-06-02
- 格式:DOCX
- 页数:25
- 大小:2.89MB
地形初步分析.docx
《地形初步分析.docx》由会员分享,可在线阅读,更多相关《地形初步分析.docx(25页珍藏版)》请在冰点文库上搜索。
地形初步分析
练习一、地形指标提取
一、题目:
地形指标提取
二、目的:
通过本实验,加深对各基本地形指标的概念及其应用意义的理解。
熟练掌握使用
ArcGIS软件提取这些地形指标的方法和步骤。
三、意义:
地形指标是最基本的自然地理要素,也是对人类的生产和生活影响最大的自然要素。
地形特征广泛应用于诸多研究和应用领域。
地形指标的提取对水土流失、土地利用、土地资源评价、城市规划等方面的研究起着重要的作用。
根据研究区域尺度的不同,地形指标有许多因子。
基于ArcGIS的地形指标的提取,大多均是基于DEM
数据完成。
四、实验数据:
本实验采用某区域栅格DEM。
图1是某区域分辨率为5米的DEM数据。
下文中关于地形指标的提取都是以这个数据为基础。
五、实验软件:
ArcGIS9.2
六、要求:
利用所提供DEM数据,提取该区域坡度变率、坡向变率、地形起伏度、地面粗糙度等四个基本地形指标的栅格图层。
七、实验步骤:
1.坡度变率
地面坡度变率,是地面坡度在微分空间的变化率,是依据坡度的求算原理,在所提取的坡度值的基础上对地面每一点再求算一次坡度。
即坡度之坡度(Slopeof
Slope,简称SOS)。
坡度是地面高程的变化率的求解,因此,坡度变率表征了地表面高程相对于水平面变化的二阶导数。
坡度变率在一定程度上可以很好的反映剖面曲率信息,其提取方法如下:
1)选中DEM图层数据,选择表面分析中的坡度(Slope)工具,提取坡度,得
到坡度数据层,命名为Slope。
计算过程如下:
A.选择表面分析的坡度工具(Slope),见图2,弹出Slope对话框(图3)
图2打开坡度计算工具
图3坡度计算对话框
B.在InputSurface中选择用来生成坡度图的DEM。
C.在Outputmeasurement中选择坡度单位(度/Degree或百分数/Percent)。
D.在Zfactor中设定高程转换系数(当输入数据所定义的空间参考具有高程
单位时,自动进行转换计算,这里采取默认值1)。
E.在Outputcellsize中指定输出图的栅格单元大小(默认与输入DEM相同)。
F.在Outputraster中指定输出路径与文件名。
得到坡度数据层,如图所示:
图4坡度提取
2)选中坡度数据层Slope,对其再用上述的方法提取坡度,得到坡度变率数据,
命名为SOS(图5)
图5坡度变率
2.坡向变率
地面坡向变率,是指在提取坡向基础上,提取坡向的变化率,亦即坡向之坡度
(SlopeofAspect,SOA)。
它可以很好地反映等高线弯曲程度。
地面坡向变率在所提取的地表坡向矩阵的基础上沿袭坡度的求算原理,提取地表局部微小范围内坡向的最大变化情况。
需要注意:
SOA在提取过程中在北面坡将会有误差产生。
北面坡坡向值范围为0~90度和270~360度,在正北方向附近,如15度和345度两个坡向之间坡向差值只是30度,而计算结果却是330度(如图6所示)。
所以要将北坡地区的坡向变率误差进行纠正,具体的操作方法为:
图6
1)求取原始DEM数据层的最大高程值,记为H;通过SpatialAnalyst下的
栅格计算器RasterCalculator。
公式为(H-DEM),得到与原来地形相反
的DEM数据层,即反地形DEM数据。
A.加载SpatialAnalyst模块,打开SpatialAnalyst下的栅格计算器(图
7)。
图7栅格计算器
B.提取DEM数据层的最大高程值,约为1153.79,将公式redem=1153.79-
[dem]输入公式编辑器中(图8)。
注意,在公式编辑器中如果引用Layers
选择框的数据层,数据层名必须用[]括起来,如,[dem]。
图8用栅格计算器计算反地形DEM数据
C.点击Evaluate进行计算,得到反地形DEM数据,记为redem(图9)
图9反地形DEM数据
2)基于反地形DEM数据求算坡向值。
A.在3DAnalyst下选择表面分析的坡向工具(Aspect),见图10,弹出
“Aspect”对话框(图11)。
图10坡向工具
图11计算坡向对话框
B.选择输入表面数据(Inputsurface),即第一步做的反地形DEM数据。
C.在Outputcellsize中指定输出栅格单元大小。
D.在Outputraster中指定输出路径和文件名。
E.点击OK完成反地形DEM数据求算坡向值,记为redemaspect(图12)
图12反地形DEM数据求算坡向值
3)利用SOA方法求算反地形的坡向变率,记为SOA2,由原始DEM数据求算出的
坡向变率值为SOA1。
A.对上一步求出的反地形DEM坡向值再求算坡度值,即得到反地形的坡向变率
SOA2(图13)。
图13反地形DEM数据坡向变率
B.对原始DEM数据,先求算坡向值记为aspect(图14)
图14原始DEM数据坡向值
C.然后再对aspect求算坡度值,即得到原始DEM数据的坡向变率SOA1(图15)
图15原始DEM数据坡向变率
4)在SpatialAnalyst下使用栅格计算器RasterCalculator,公式为SOA=
(([SOA1]+[SOA2])-Abs([SOA1]-[SOA2]))/2,即可求出没有误差的DEM
的坡向变率,如图16、17所示。
图16栅格计算器计算
图17计算出没有误差的坡向变率
3.地形起伏度
地形起伏度是指特定的区域内,最高点海拔高度与最低点海拔高度的差值。
它是描述一个区域地形特征的宏观性指标。
地形起伏度的计算,可先求出一定范围内海拔高度的最大值和最小值,然后对其求差值即可。
可以使用SpatialAnalyst中的栅格邻域计算工具Neighborhood
Statistics求得最大值和最小值,邻域的设置可以为圆,也可以为矩形,邻域的大小可以根据自己的要求来确定。
地形起伏度的提取方法如下:
1)选中DEM数据,在SpatialAnalyst下使用栅格邻域计算工具Neighborhood
Statistics(图18)。
设置Statistictype为最大值,邻域的类型为矩形(也
可以为圆),邻域的大小为11×11(这个值也可以根据自己的需要进行改变),
则可得到一个邻域为11×11的矩形的最大值层面,记为A(图19)。
图18栅格邻域计算工具
图19邻域为11×11的矩形的最大值层面A
2)重复第1步,只是把Statistictype值设置为最小值,即可得到DEM数据的
最小值层面,记为B(图20)
图20邻域为11×11的矩形的最小值层面B
3)在SpatialAnalyst下使用栅格计算器RasterCalculator,公式为[A]-[B],
如图21所示,即可得到一个新层面,其每个栅格的值是以这个栅格为中心的
确定邻域的地形起伏度。
提取的结果如图22。
图21用栅格计算器进行计算
图22地形起伏度
4.地面粗糙度
地面粗糙度是特定的区域内地球表面积与其投影面积之比。
它也是反映地表形态的一个宏观指标。
根据地面粗糙度的定义,求每个栅格单元的表面积与其投影面积之比,可以用如下方法来完成。
如图23,假如ABC是一个栅格单元的纵剖面,α为此栅格单元的坡度,则AB面的面积为此栅格的表面积,AC面为此栅格的投影面积(也即是此栅格的面积),根据公式:
图23
则可得出此栅格单元的地面粗糙度M为:
M=AB面的面积/AC栅格单元的面积=(AC*AB)/(AC*AC)=1/Cosα
地面粗糙度的提取步骤如下:
1)点击DEM数据层,选择表面分析中的坡度(Slope)工具,提取得到坡度数据
层,命名为Slope(注意:
这次要选择计算度数,而不是百分比)。
2)点击Slope数据层,在SpatialAnalyst下使用栅格计算器RasterCalculator
公式为:
1/Cos([Slope]*3.14159/180),即可得到地面粗糙度数据层,如图24。
图24表面粗糙度
需要注意的是,在ArcGIS中,Cos使用弧度值作为角度单位,而利用表面分析
工具提取得到的坡度是角度值,所以在计算时必须把角度转为弧度。
此外,地形指标还包括一些常用的水文因子,如坡长、沟壑密度等,该类因子的
提取一般通过水文方法实现。
练习二、地形特征信息提取
一、题目:
地形特征信息提取
二、目的:
通过本实例,掌握山脊线和山谷线这两个基本地形特征信息的理论及其基于DEM
的提取方法与原理。
同时,熟练掌握利用ArcGIS软件对这两个地形特征信息的提取
方法。
三、意义:
特征地形要素,主要是指对地形在地表的空间分布特征具有控制作用的点、线或面状要素。
特征地形要素构成地表地形与起伏变化的基本框架。
与地形指标的提取主要采用小范围的邻域分析不同的是,特征地形要素的提取更多地应用较为复杂的技术方法,如山谷线、山脊线、沟沿线等的提取采用了全局分析法(globalprocess),成为栅格数据地学分析中很具特色的数据处理内容。
特征地形要素从表示的内容上可分为地形特征点和特征线两大类。
地形特征点主要包括山顶点(peak)、凹陷点(pit)、脊点(ridge)、谷点(channel)、鞍点(pass)、平地点(plane)等。
基于DEM提取地形特征点,可利用一个3×3或更大的栅格窗口,通过中心格网点与8个邻域格网点的高程关系来进行判断获取。
山脊线和山谷线构成了地形起伏变化的分界线(骨架线),因此它对于地形地貌研究具有重要的意义。
另一方面,对于水文物理过程研究而言,由于山脊、山谷分别表示分水性与汇水性,山脊线和山谷线的提取实质上也是分水线与汇水线的提取。
下文通过山脊线和山谷线的提取,介绍了如何基于ArcGIS完成地形特征信息的
提取。
自动提取山脊线和山谷线的主要方法都是基于规则格网DEM数据的,从算法设计原理上来分,大致可以分为以下五种:
1)基于图像处理技术的方法。
2)基于地形表面几何形态分析的方法。
3)基于地形表面流水物理模拟分析方法。
4)基于地形表面几何形态分析和流水物理模拟分析相结合的方法。
5)平面曲率与坡形组合法。
其中,平面曲率与坡形组合法提取的山脊、山谷的宽度可由选取平面曲率的大小来调节,方法简便,效果好。
该方法基本处理过程为:
首先利用DEM数据提取地面的平面曲率及地面的正负地形,取正地形上平面曲率的提取比较繁琐,而坡向变率(SOA)在一定程度上可以很好地表征平面曲率。
因此,下面的提取过程以SOA代替平面曲率。
四、实验数据:
某区域栅格DEM
五、实验软件:
ArcGIS9.2
六、要求:
利用所给区域DEM数据,提取该区域山脊线、山谷线栅格数据层
七、实验步骤:
1.山脊线计算
具体提取过程为:
1)点击DEM数据,使用表面分析中的坡向(Aspect)工具,提取DEM的坡向数
据层,命名为A。
2)点击数据层A,使用表面分析中的坡度(Slope)工具,提取数据层A的坡度
数据,命名为SOA1。
3)求取原始DEM数据层的最大高程值,记为H;使用空间分析工具集中的栅格
计算器(RasterCalculator),公式为(H-DEM),得到与原来地形相反的数据层,
即反地形DEM数据。
4)基于反地形DEM数据求算坡向值。
5)利用SOA方法求算反地形的坡向变率,记为SOA2。
6)使用空间分析工具集中的栅格计算器(RasterCalculator),公式为
SOA=(([SOA1]+[SOA2])-Abs([SOA1]-[SOA2]))/2
即可求出没有误差的DEM的坡向变率SOA。
7)再次点击初始DEM数据,使用空间分析工具集中的栅格邻域计算工具
(NeighborhoodStatistics);设置统计类型(Statistictype)为平均值(Mean),
邻域的类型为矩形(也可以为圆),邻域的大小为11×11(这个值也可以根据自己的需要进行改变),则可得到一个邻域为11×11的矩形的平均值数据层,记为B。
8)使用空间分析工具集中的栅格计算器(RasterCalculator),公式为
C=[DEM]-[B],即可求出正负地形分布区域(图25)
图25正负地形分布示意图
9)使用空间分析工具集中的栅格计算器(RasterCalculator),公式为
shanji=[C]>0&SOA>70,即可求出山脊线(图26)
图26山脊线
2.山谷线计算
10)同理,在栅格计算器(RasterCalculator)中,键入公式为
shangu=[C]<0&SOA>70,即可求出山谷线(图27)。
图27山谷线
盛年不重来,一日难再晨。
及时宜自勉,岁月不待人。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 地形 初步 分析