地形指标提取.docx

地形指标提取1. 背景:地形指标是最基本的自然地理要素，也是对人类的生产和生活影响最大的自然要素，地 形特征广泛应用于诸多研究领域和应用领域地形指标的提取对水土流失、土地利用、 土地资源评价、城市规划等方面的研究起着重要的作用根据研究区域尺度的不同，地 形指标有许多因子基于 ArcGIS 的地形指标的提取，大多均是基于 DEM 数据完成的2. 目的： 通过本实验，使读者加深对各基本地形指标的概念及其应用意义的理解，熟练掌握使用 ArcGIS 软件提取这些地形指标的方法和步骤3. 要求：利用所提供的DEM数据，提取该区域坡度变率、坡向变率、地形起伏度、地面粗糙度 等四个基本地形指标的栅格图层4. 实验步骤:（1）坡度变率 地面坡度变率，是地面坡度在微分空间的变化率，是依据坡度的求算原则，在所提取的 坡度值的基础上对地面每一点再求算一次坡度即坡度之坡度坡度是地面高程的变化 率的求解，因此，坡度变率表征了地面高程相对于水平面变化的二阶导数坡度变率在一定程度上可以很好的反映剖面曲率信息，其提取方法如下：1）选中DEM图层数据，选择表面分析中的坡度工具，提取坡度，得到坡度数据层， 命名为 Slope2）选中坡度数据层Slope，对其再用上述的方法提取坡度，得到坡度变率数据，命名 为 SOS（2）坡向变率地面坡向变率，是指在提取坡向基础上，提取坡向的变化率，亦即坡向之坡度（S0A）, 它可以很好的反应等高线弯曲程度。

地面坡向变率在所提取的地表坡向矩阵的基础上沿袭坡度的求算原理，提取地表局部微 小范围内的最大变化情况需要注意：SOA在提取过程中在背面坡将会有误差产生北面坡 坡向值范围为 0~90度和 270~360度，在正北方向附近，如15 度和345 度两个坡向之间坡向 差只是 30度，而计算结果却是 330度所以要将北坡地区的坡向变率误差进行纠正，具体 方法如下：1）求取原始DEM数据层的最大高程值，记为H,通过SpatialAnalyst下的栅格计算器 Raster Calculator公式为H=0-DEM，得到与原来地形相反的DEM数据层，即反地形 DEM 数据2）基于反地形 DEM 数据求算坡向值3）利用SOA方法求算反地形的坡向变率，即为S0A2，由原始DEM数据求算出的坡向变率 值为 SOA14）在SpatialAnalyst下使用栅格计算器RasterCalculator，公式为（3）地形起伏度 地形起伏度是指特定地区域内，最高点海拔高度与最低点海拔高度的差值，它是描述一 个区域地形特征的宏观性指标地形起伏度的计算，可先求出一定范围内海拔高度的最大值和最小值，然后对其求差值 即可可以使用 SpatialAnalyst 中的栅格邻域计算工具 NeihborhoodStatistics 求得最大 值和最小值，邻域的设置可以为圆，也可以为矩形，邻域的大小可根据自己的要求来确定。

地形起伏的提取方法如下：1）选中 DEM 数据，在 SpatialAanlyst 下使用栅格领域工具 NeighborhoodStatistics 设置 Statistictype 为最大值，邻域的类型为矩形（也可以为圆），邻域的大小为 11*11,，则可得到一个邻域为 11*11 的矩形的最大值层面，记为 A2）重复第一步，只是把StatisticType值设置为最小值，即可得到DEM数据的最小值层面,记为 B3）在SpatialAnalyst下使用栅格计算器Raster Calculator,公式为A-B即可得到一 个新的层面，其每个栅格的值是以这个栅格为中心确定邻域的地形起伏度提取结果 如图4）地面粗糙度 地面粗糙度是特定的区域内地球表面积与其投影面积之比，它也是反映地表形态的一个 宏观指标地面粗糙度的提取步骤如下：1）点击DEM数据层，选择表面分析中的坡度（slope）工具，提取得到坡度数据层，命名 为 Slope2) 点击 Slope 数据层，在 SpatialAnalyst 下使用栅格计算器 RasterCalculator 公式: l/cos(slope*3.14159/180),即可得到地面粗糙度数据层。

About Building Expressions Evaluate Cancel需要注意的是，在ArcGIS中，Cos使用弧度值作为角度单位，而利用表面分析工具提取 得到的坡度是角度值，所以在计算时必须把角度转换为弧度此外，地形指标还包括一些常用的水文因子，如坡长、沟壑密度等，该因子的提取一般通 过水文方法实现。