ArcGIS表面分析 - 图文

表面分析

来源：互联网

表面分析主要通过生成新数据集，诸如等值线、坡度、坡向、山体阴影等派生数据，获得更多的反映原始数据集中所暗含的空间特征、空间格局等信息。在ArcGIS 中，表面分析的主要功能有：查询表面值、从表面获取坡度和坡向信息、创建等值线、分析表面的 可视性、从表面计算山体的阴影、确定坡面线的高度、寻找最陡路径、计算面积和体积、 数据重分类、将表面转化为矢量数据等。

一 栅格插值

一般情况下采集到的数据都是以离散点的形式存在的，只有在这些采样点上才有较为准确的数值，而其它未采样点上都没有数值。然而，在实际应用中却很可能需要用到某些未采样点的值，这个时候就需要通过已采样点的数值来推算未采样点值。这样的一个过程也就是栅格插值过程。插值结果将生成一个连续的表面，在这个连续表面上可以得到每一点的值。栅格插值包括简单栅格表面的生成和栅格数据重采样。

通过栅格插值运算生成栅格表面是栅格空间分析中很重要的一部分，很多操作都将基于离散点插值生成的表面来进行。用来插值的离散点可以是多种数据，例如，空气污染指数，土壤有机质含量，或离散高程点。

在ArcGIS 9 栅格分析模块中，通过栅格 插值运算生成表面主要有三种实现方式：反距离权重插值，样条函数插值和克里格插值， 如图1 所示。

图1 栅格插值方法

下面以一组土壤元素PH 值的插值来逐一说明在Spatial Analyst 中三种表面生成插值的实现过程。 1. 反距离权重插值（IDW）

（1） 在Spatial Analyst 下拉菜单中选择Interpolate to Raster，在弹出的下一级菜单中点击Inverse Distance Weighted 命令, 弹出IDW 对话框，如图2。

图2 IDW 对话框

（2） 在Input points 的下拉菜单中选择被用来进行插值的离散点数据； （3） 在Z value field 的下拉菜单中选择要加入的字段；

（4） 在Power 栏中填入进行插值计算的幂值；幂值就是距离的指数。如幂指数为2时则进行反向距离平方插值。幂指数是一个正实数，其缺省值为2。

（5） 在Search radius type 栏中选择一种搜索半径设置类型；

1) Variable：当选择此项时,搜索半径由下面两两个量来控制，，Number of points 和Maximum distance。首先在Number of points 中 输入搜索的最近点的个数（缺省值为12），然后 在Maximum distance 中输入一个控制距离。如果最近点的个数超出控制距离，则将会以控制距离为限制来选取较少的点；

2) Fixed：当选择此项时，搜索半径由下面两个量来控制，Distance 和Minimum number of points。首先在Distance 中输入搜索半径距离（缺省值是输出栅格大小的五倍），然后在 Minimum number of points 中输入控制插值点个数的最小整数值。如果搜索半径距离内的点个数小于插值点个数的最小整数值，则搜索半径自动增大。

（6） Use barriers polyline 为可选项，输入 中断线文件。barriers 是在插值中，如有某些地方出现异常，（如某些断裂带），而要求插值时考虑到这样的因素，所设置的选项。它是一个打断表面的线特征。这一线特征没有 Z 值。悬崖，峭壁，堤岸或某些障碍都是典型的barriers。barriers 限制了插值计算，它使得计算只在线的两侧各自进行。而落在线上的点则会同时参与线两侧的计算。

（7） Output cell size：指定输出结果的栅格大小； （8） Output raster：为输出结果指定目录及名称； （9） 点击OK 按钮。生成结果如图3 所示。

图3 IDW 插值

2. 样条函数插值（SPLINE）

样条函数插值采用两种不同的计算方法，Regularized 和Tension。如选择Regularized， 它将生成一个平滑、渐变的表面，得出的插值结果很可能会超出样本点的取值范围。如选择Tension，它会根据要生成的现象的特征生成一个比较坚硬的表面，得出结果的插值更接近限制在样本点的取值范围内。

同时，计算过程中，除了需要选择不同的计算方法，还需要在每种方法中设定一个合适的权重（weight）。选择Regularized 时，他决定了表面最小曲率三次导的权重。权重越高表面越光滑。可能用到的典型值有：0、0.001、0.01、1 和5。选择Tension 时，他决定了Tension 的权重。权重越高，表面越粗糙。可能用到的典型值有：0、1、5 和10。 样条函数插值过程如下：

（1） 在Spatial Analyst 下拉菜单中选择 Interpolate to Raster, 在弹出的下一级菜单中点击 Spline，出现Spline 对话框，如图4，以下所有设置如图中所示。

图4 Spline 对话框

（2） 在Input points 的下拉菜单中选择被用来进行插值的离散点数据；在Z value field 的下拉菜单中选择要加入的字段；

（3） 在Spline type 中填入样条函数插值的类型 ； （4） 在Weight 栏中填入一个影响插值的特 征参数；

（5） 在Number of points 栏中输入每一个区域内用来估值点的个数。它的缺省值是 12。 （6） 指定输出结果的栅格大小； （7） 为输出结果指定目录及名称；

（8） 点击OK 按钮。插值结果如图5。

图5 Spline 插值结果

3. 克里格插值（Kriging）

（1） 在Spatial Analyst 下拉菜单中选择Interpolate to Raster, 在弹出的下一级菜单中点击Kriging，出现Kriging 对话框，如图6，以下所有设置如图中所示。

（2） 在Input points 的下拉菜单中选择被用来进行插值的离散点数据； （3） 在Z value field 的下拉菜单中选择要加入的字段； （4） 选择你所需要的克里格方法；

（5） 在Semivariogram model 的下拉菜单中选择你所需要使用的模型；

（6） 点击Advanced Parameters 按钮，如果知道可以指定这些参数，另外空间分析模块也将为你估算这些参数；

（7） 在Search radius type栏中选择一种搜索半径设置类型；（设置方法同IDW） （8） 指定输出结果的栅格大小；

（9） 可选项，是否需要生成预测的标准误差； （10） 为输出结果指定目录及名称； （11） 点击OK 按钮。结果如图7。

图6

图7 Kriging 插值结果

4. 数据重采样

栅格插值除了包括简单栅格表面的生成还应包括栅格数据重采样。重采样是栅格数据空间分析中处理栅格分辨率匹配问题的常用数据处理方法。进行空间分析时，用来分析的数据资料由于来源不同，经常会出现栅格不同大小问题，这时为了便于分析，就需要将不同的栅格大小转化为同样栅格大小，这样的一个过程就是栅格数据的重采样过程。

栅格数据的重采样主要基于三种方法：最邻近采样（NEAREST），双线性 ILINEAR）采样和三次卷积采样（CUBIC）。

（1) 最邻近采样：此重采样法用输入栅格数据中最邻近栅格值作为输出值。因此, 在重采样后的输出栅格中的每个栅格值, 都是输入栅格数据中真实存在而未加任何改变的 值。这种方法简单易用，计算量小，重采样的速度最快。

（2) 双线性采样：此重采样法取待采样点（x，y）点周围四个邻点，在y方向（或x 方向 ）内插

两次，再在x方向（或y方向）内插一次，得到（x，y）点的栅格值。

（3) 三次卷积采样：这是进一步提高内插精度的一种方法。它的基本思想是增加邻点来获得最佳插值函数。取待计算点周围相邻的16个点，与双线性采样类似，可先在某一方向上内插，如先在x方向上，每四个值依次内插四次，再根据四次的计算结果在y方向上内插，最终得到内插结果。

在ArcGIS 中重采样 功能是在ArcToolbox 下来实现的。在ArcToolbox 中，Data Management Tools 下的Raster 子菜单中，选择Resample 命令（图8）。

参数说明：

（1) Input raster：选择输入栅格数据；

（2) Output raster：设置输出栅格数据路径 及名称； （3) Out put cell size：设置输出栅格单元大小； （4)Resampling technique：选择重采样方法 1）NEAREST：最邻近距离法； 2） BILINEAR：双线性内插法； 3） CUBIC：立方体法。

图8 重采样参数设置框

二 等值线绘制

等值线是将表面上相邻的具有相同值的点连接起来的线，如地形图上的等高线、气温图上的等压线。等值线分布的疏密一定程度上表明了表面值的变化情况。值的变化越小的地方，等值线就越疏，反之越密。因此，通过研究等值线的疏密情况，可以获得对表面值变化的大致情况。

ArcGIS 中等值线的绘制方法如下：

1. 在Spatial Analyst 下拉菜单中选择Surface analysis, 在弹出的下一级菜单中点击Contour，出现Contour 对话框，如图9，以下所有设置如图中所示。

图9 等高线生成对话框

2. 在Input Surface 的下拉菜单中选择用来生成等高线的表面； 3. 在Contour interval 栏中设置等高线间距； 4. 在Base contour 栏中指定等高线基准高程； 5. 在Z factor栏中设定高程变换系数：

6 在Output features 栏中指定输出等高线的存放路径与文件名。 7. 点击OK 按钮。

图10 为某区域的等高线图,背景为该地区的地形光照晕渲图。

图10 某区域等高线图（背景为光照晕渲图）

三 地形因子提取

因子分析方法是GIS 空间分析，尤其GIS 数字地形分析常用的基本分析方法。不同的地形因子从不同侧面反映了地形特征性，实际应用人们提出了各种各样的地形因子。从其所描述的空间区域范围，常用的地形因子可以划分为微观地形因子与宏观地形因子两种基本类型（图11）。按照提取地形因子差分计算的阶数，又可将地形因子分为一阶地形因子、二阶地形因子和高阶地形因子（图12）。

图11依据空间区域范围的坡面因子分类体系

其中，坡度、坡向、平面曲率、剖面曲率在ArcGIS 中可直接提取，其它因子的提取则需要进行一系列的复合计算。后者的具体提取过程可以参阅相关资料，这里不作介绍。

图12 基于提取算法的坡面因子分类体系

1. 坡度的提取

地表面任一点的坡度（Slope）是指过该点的切平面与水平地面的夹角（图13）。坡度表示了地表面在该点的倾斜程度。

图13地表单元坡度示意图

实际应用中，坡度有两种表示方式方法：

（1） 坡度(degree of slope)：既水平面与地形面之间夹角。

（2） 坡度百分比（percent slope）：既高程增量与水平增量之比的百分数。 （3） ArcGIS 中坡度的提取过程为：

（4） 在Spatial Analyst 下拉菜单中选择Surface analysis, 在弹出的下一级菜单中点击 Slope，出现Slope 对话框，如图14，以下所有设置如图中所示。

图14 坡度生成对话框

（5） 在Input Surface 的下拉菜单中选择用来生成坡度的表面； （6） 选择一种坡度表示方法，在此分别用两种方法做了坡度。 （7） 在Z factor 栏中设定高程变换系数； （8） 在Output cell size 栏中设定栅格大小；

（9） 在Output raster 栏中指定输出坡度的存放路径与文件名。

（10） 点击OK 按钮。 坡度结果如图15所示，坡度百分比如图16 所示：

图15坡度结果图 图16坡度百分比图

2. 坡向的提取

坡向指地表面上一点的切平面的法线矢量在水平面的投影与过该点的正北方向的夹角。对于地面任何一点来说，坡向表征了该点高程值改变量的最大变化方向。在输出的坡向数据中，坡向值有如下规定：正北方向为0 度，按顺时针方向计算，取值范围为0°～ 360°。

ArcGIS 中坡向的提取过程为：

（1） 在Spatial Analyst 下拉菜单中选择Surface analysis, 在弹出的下一级菜单中点击 Aspect，出现Aspect 对话框，如图17，以下所有设置如图中所示。

图17 坡向生成对话框

（2） 在Input Surface 的下拉菜单中选择用来生成坡向的表面； （3） 在Output cell size 栏中设定栅格大小；

（4） 在Output raster 栏中指定输出坡向的存放路径与文件名。 坡向结果如图18 所示。

图18 坡向结果图

3. 平面曲率、剖面曲率的提取

地面曲率是对地形表面一点扭曲变化程度的定量化度量因子，地面曲率在垂直和水平两个方向上分量分别称为剖面曲率和平面曲率。剖面曲率是对地面坡度的沿最大坡降方向地面高程变化率的度量。平面曲率指在地形表面上，具体到任何一点，指用过该点的水平面沿水平方向切地形表面所得的曲线在该点的曲率值。平面曲率描述的是地表曲面沿水平方向的弯曲、变化情况，也就是该点所在的地面等高线的弯曲程度。

下面以一幅DEM 图来说明平面曲率、剖面曲率的提取过程为：

（1） 打开ArcGIS 的Toolbox，在Spatial Analyst Tools 底下选择Surface，在Surface 的下一级菜单中选择Curvature。打开Curvature 对话框，如图19。

（2） 在Input raster 的下拉菜单中选择用来计算曲率的栅格数据；

（3） 在Output curvature raster 栏中设定输出总曲率的存放路径与文件名； （4） 在Z factor 栏中设定高程变换系数；

（5） 在Output profile curve raster 栏中指定输出剖面曲率的存放路径与文件名。

（6） 在Output plan curve raster 栏中指定输出平面曲率的存放路径与文件名。 总曲率结果如图20，剖面曲率结果如图21，平面曲率结果如图22：

图19 平面曲率、剖面曲率提取对话框 图20 总曲率结果图

图21剖面曲率结果图 图22 平面曲率结果图

四 山体阴影

山体阴影是根据假想的照明光源对高程栅格图的每个栅格单元计算照明值。山体阴影图不仅很好地表达了地形的立体形态，而且可以方便的提取地形遮蔽信息。计算过程中包括三个重要参数：太阳方位角、太阳高度角、表面灰度值。

太阳方位角以正北方向为0 度，按顺时针方向度量，如90 度方向为正东方向（图23）。 由于人眼的视觉习惯，通常默认方位角为315 度，即西北方向。

太阳高度角为光线与水平面之间的夹角，同样以度为单位（图24）。为符合人眼视觉习惯，通常默认为45 度。

ARCGIS 中，默认情况下，光照产生的灰度表面的值的范围为0~255。山体阴影的实现过程如下： 1. 计算表面阴影

(1) 在Spatial Analyst 下拉菜单中选择Surface analysis, 在弹出的下一级菜单中点击 Hillshade，出现Hillshade 对话框，如图25，以下所有设置如图中所示。

图23 太阳方位角度量示意图 图24太阳高度角示意图

图25 计算表面阴影

(2) 在Input surface 下拉菜单中选择用来计算阴影的表面； (3) 在Azimuth 栏中设置太阳方位角； (4) 在Altitude 栏中设置太阳高度角； (5) 在Z factor 栏中设定高程变换系数；

(6) 在Output cell size 栏中指定输出栅格单元大小；

(7) 在Output raster 栏中指定输出表面阴影的存放路径与文件名。

(8) 点击OK 按钮。 图26为某区域设定入射方位角为315 度，高度角为45 度所生成的表面阴影栅格图像。

2. 阴影化

通过阴影建模工具（在图24 对话框中选中该复选框），当给定太阳高度角和方位角 时，可以快速提取被阴影遮蔽的栅格单元，它们会被赋值为0。再通过空间分析的重分类 （Reclassify）方法生成二值图像，阴影区赋值为０，非阴影区赋值为１。图27 即为阴影二值图。

图26 表面阴影图 图27阴影二值图

本文来源：https://www.bwwdw.com/article/3mk3.html