三维分析.pdfVIP

三维分析 本章介绍几个实例与练习，加深对三维分析理论与方法的理解，熟练各种分析工具的 操作和使用。 一 地形指标提取 1. 背景 地形指标是最基本的自然地理要素，也是对人类的生产和生活影响最大的自然要素。 地形特征广泛应用于诸多研究和应用领域。地形指标的提取对水土流失、土地利用、土地资 源评价、城市规划等方面的研究起着重要的作用。根据研究区域尺度的不同，地形指标有许 多因子。基于ArcGIS 的地形指标的提取，大多均是基于DEM 数据完成。 2. 目的 通过本实验，使读者加深对各基本地形指标的概念及其应用意义的理解。熟练掌握使 用ArcGIS 软件提取这些地形指标的方法和步骤。 3. 实验数据 本实验采用某区域栅格数据DEM ，存放在数据/exe1 内。 4. 要求 利用所提供 DEM 数据，提取该区域坡度变率、坡向变率、地形起伏度、地面粗糙度 等四个基本地形指标的栅格图层。 5. 实验步骤 （1）坡度变率 地面坡度变率，是地面坡度在微分空间的变化率，是依据坡度的求算原理，在所提取 的坡度值的基础上对地面每一点再求算一次坡度（slope of slope，简称SOS）。坡度是地面 高程的变化率的求解，因此，坡度变率表征了地表面高程相对于水平面变化的二阶导数。 坡度变率在一定程度上可以很好的反映剖面曲率信息，其提取方法如下： 1）选中DEM 图层数据，选择表面分析中的坡度（slope ）工具（图2.1 ），提取坡度，得到 坡度数据层，命名为slope 。 图2.1 用surface analysis 提取DEM 的坡度 2 ）选中坡度数据层slope，对其再用上述的方法提取坡度，得到坡度变率数据，命名为SOS （图2.2 ）。 图2.2 坡度变率 （2 ）坡向变率 地面坡向变率，是指在提取坡向基础上，提取坡向的变化率，亦即坡向之坡度（slope of aspect，SOA）。它可以很好地反映等高线弯曲程度。 地面坡向变率在所提取的地表坡向矩阵的基础上沿袭坡度的求算原理，提取地表局部 微小范围内坡向的最大变化情况。需要注意：SOA 在提取过程中在北面坡将会有误差产生。 北面坡坡向值范围为0~90 度和270~360 度，在正北方向附近，如 15 度和345 度两个坡向 之间坡向差值只是30 度，而计算结果却是330 度（如图2.3 ）。所以要将北坡地区的坡向变 率误差进行纠正，具体操作方法为： 1）求取原始 DEM 数据层的最大高程值，记 为H ；通过spatial analyst 下的栅格计算器 raster calculator ，公式为（H-DEM ），得到 与原理地形相反的 DEM 数据层，即反地 形DEM 数据。 2 ）基于反地形DEM 数据求算坡向值。 3 ）利用 SOA 方法求算反地形的坡向变率， 记为SOA2，由原始DEM 数据求算出的坡 向变率值为SOA1。 图2.3 北坡坡向差值示意图 4 ）在spatial analyst 下使用栅格计算器raster calculator，公式为SOA= （（[SOA1]+[SOA2]）-Abs （[SOA1]-[SOA2]））/2 ，即可求出没 有误差的DEM 的坡向变率，如图2.4 。 图2.4 坡向变率 （3 ）地形起伏度 地形起伏度是指特定的区域内，最高点海拔高度与最低点海拔高度的差值。它是描述 一个区域地形特征的宏观性指标。 地形起伏度的计算，可先求出一定范围 内海拔高度的最大值和最小值，然后对其求差 值即可。可以使用spatial analyst 中的栅格邻域 计算工具neighborhood statistics 求得最大值和 最小值。 地形起伏度的提取方法如下： 1）选中DEM 数据，在spatial analyst 下使用 栅格邻域计算工具 neighborhood statistics （图2.5 ）。设置statistic type 为最大值， 邻域的类型为矩形（也可为圆