ArcGIS表面分析.doc

《ArcGIS表面分析.doc》由会员分享，可在线阅读，更多相关《ArcGIS表面分析.doc（13页珍藏版）》请在金锄头文库上搜索。

1、表面分析 来源：互联网表面分析主要通过生成新数据集，诸如等值线、坡度、坡向、山体阴影等派生数据，获得更多的反映原始数据集中所暗含的空间特征、空间格局等信息。在ArcGIS 中，表面分析的主要功能有：查询表面值、从表面获取坡度和坡向信息、创建等值线、分析表面的 可视性、从表面计算山体的阴影、确定坡面线的高度、寻找最陡路径、计算面积和体积、 数据重分类、将表面转化为矢量数据等。一 栅格插值 一般情况下采集到的数据都是以离散点的形式存在的，只有在这些采样点上才有较为准确的数值，而其它未采样点上都没有数值。然而，在实际应用中却很可能需要用到某些未采样点的值，这个时候就需要通过已采样点的数值来推算未采样

2、点值。这样的一个过程也就是栅格插值过程。插值结果将生成一个连续的表面，在这个连续表面上可以得到每一点的值。栅格插值包括简单栅格表面的生成和栅格数据重采样。 通过栅格插值运算生成栅格表面是栅格空间分析中很重要的一部分，很多操作都将基于离散点插值生成的表面来进行。用来插值的离散点可以是多种数据，例如，空气污染指数，土壤有机质含量，或离散高程点。 在ArcGIS 9 栅格分析模块中，通过栅格 插值运算生成表面主要有三种实现方式：反距离权重插值，样条函数插值和克里格插值， 如图1 所示。图1 栅格插值方法下面以一组土壤元素PH 值的插值来逐一说明在Spatial Analyst 中三种表面生成插值的实

3、现过程。 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 rad

4、ius type 栏中选择一种搜索半径设置类型； 1) Variable：当选择此项时,搜索半径由下面两两个量来控制，Number of points 和Maximum distance。首先在Number of points 中 输入搜索的最近点的个数（缺省值为12），然后 在Maximum distance 中输入一个控制距离。如果最近点的个数超出控制距离，则将会以控制距离为限制来选取较少的点； 2) Fixed：当选择此项时，搜索半径由下面两个量来控制，Distance 和Minimum number of points。首先在Distance 中输入搜索半径距离（缺省值是输出栅格大小的

5、五倍），然后在 Minimum number of points 中输入控制插值点个数的最小整数值。如果搜索半径距离内的点个数小于插值点个数的最小整数值，则搜索半径自动增大。 （6） Use barriers polyline 为可选项，输入 中断线文件。barriers 是在插值中，如有某些地方出现异常，（如某些断裂带），而要求插值时考虑到这样的因素，所设置的选项。它是一个打断表面的线特征。这一线特征没有 Z 值。悬崖，峭壁，堤岸或某些障碍都是典型的barriers。barriers 限制了插值计算，它使得计算只在线的两侧各自进行。而落在线上的点则会同时参与线两侧的计算。 （7） Outpu

6、t cell size：指定输出结果的栅格大小； （8） Output raster：为输出结果指定目录及名称； （9） 点击OK 按钮。生成结果如图3 所示。图3 IDW 插值2. 样条函数插值（SPLINE） 样条函数插值采用两种不同的计算方法，Regularized 和Tension。如选择Regularized， 它将生成一个平滑、渐变的表面，得出的插值结果很可能会超出样本点的取值范围。如选择Tension，它会根据要生成的现象的特征生成一个比较坚硬的表面，得出结果的插值更接近限制在样本点的取值范围内。 同时，计算过程中，除了需要选择不同的计算方法，还需要在每种方法中设定一个合适的权重

7、（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 的下拉菜单中选择被用来进行插值

8、的离散点数据；在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, 在弹出的下一级菜单中

9、点击Kriging，出现Kriging 对话框，如图6，以下所有设置如图中所示。 （2） 在Input points 的下拉菜单中选择被用来进行插值的离散点数据； （3） 在Z value field 的下拉菜单中选择要加入的字段； （4） 选择你所需要的克里格方法； （5） 在Semivariogram model 的下拉菜单中选择你所需要使用的模型； （6） 点击Advanced Parameters 按钮，如果知道可以指定这些参数，另外空间分析模块也将为你估算这些参数； （7） 在Search radius type栏中选择一种搜索半径设置类型；（设置方法同IDW）（8） 指定输出结果的

10、栅格大小； （9） 可选项，是否需要生成预测的标准误差； （10） 为输出结果指定目录及名称； （11） 点击OK 按钮。结果如图7。图6图7 Kriging 插值结果4. 数据重采样 栅格插值除了包括简单栅格表面的生成还应包括栅格数据重采样。重采样是栅格数据空间分析中处理栅格分辨率匹配问题的常用数据处理方法。进行空间分析时，用来分析的数据资料由于来源不同，经常会出现栅格不同大小问题，这时为了便于分析，就需要将不同的栅格大小转化为同样栅格大小，这样的一个过程就是栅格数据的重采样过程。 栅格数据的重采样主要基于三种方法：最邻近采样（NEAREST），双线性 ILINEAR）采样和三次卷积采样（C

11、UBIC）。 （1) 最邻近采样：此重采样法用输入栅格数据中最邻近栅格值作为输出值。因此, 在重采样后的输出栅格中的每个栅格值, 都是输入栅格数据中真实存在而未加任何改变的 值。这种方法简单易用，计算量小，重采样的速度最快。 （2) 双线性采样：此重采样法取待采样点（x，y）点周围四个邻点，在y方向（或x 方向 ）内插两次，再在x方向（或y方向）内插一次，得到（x，y）点的栅格值。 （3) 三次卷积采样：这是进一步提高内插精度的一种方法。它的基本思想是增加邻点来获得最佳插值函数。取待计算点周围相邻的16个点，与双线性采样类似，可先在某一方向上内插，如先在x方向上，每四个值依次内插四次，再根据四

12、次的计算结果在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）

13、 CUBIC：立方体法。 图8 重采样参数设置框 二 等值线绘制 等值线是将表面上相邻的具有相同值的点连接起来的线，如地形图上的等高线、气温图上的等压线。等值线分布的疏密一定程度上表明了表面值的变化情况。值的变化越小的地方，等值线就越疏，反之越密。因此，通过研究等值线的疏密情况，可以获得对表面值变化的大致情况。 ArcGIS 中等值线的绘制方法如下： 1. 在Spatial Analyst 下拉菜单中选择Surface analysis, 在弹出的下一级菜单中点击Contour，出现Contour 对话框，如图9，以下所有设置如图中所示。图9 等高线生成对话框2. 在Input Surface

14、 的下拉菜单中选择用来生成等高线的表面； 3. 在Contour interval 栏中设置等高线间距； 4. 在Base contour 栏中指定等高线基准高程； 5. 在Z factor栏中设定高程变换系数： 6 在Output features 栏中指定输出等高线的存放路径与文件名。 7. 点击OK 按钮。 图10 为某区域的等高线图,背景为该地区的地形光照晕渲图。图10 某区域等高线图（背景为光照晕渲图）三 地形因子提取 因子分析方法是GIS 空间分析，尤其GIS 数字地形分析常用的基本分析方法。不同的地形因子从不同侧面反映了地形特征性，实际应用人们提出了各种各样的地形因子。从其所描述

15、的空间区域范围，常用的地形因子可以划分为微观地形因子与宏观地形因子两种基本类型（图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） 选择一种坡度表示方