《ArcGIS表面分析报告.docx》由会员分享，可在线阅读，更多相关《ArcGIS表面分析报告.docx（19页珍藏版）》请在三一文库上搜索。
1、实用文档外表分析来源：互联网外表分析主要通过生成新数据集,诸如等值线、坡度、坡向、山体阴影等派生数据,获得更多的反映原始数据集中所暗含的空间特征、空间格局等信息.在ArcGIS中,外表分析的主要功能有：查询外表值、从外表获取坡度和坡向信息、创立等值线、分析外表的可视性、从外表计算山体的阴影、确定坡面线的高度、寻找最陡路径、计算面积和体积、数据重分类、将外表转化为矢量数据等.一栅格插值一般情况下采集到的数据都是以离散点的形式存在的,只有在这些采样点上才有较为准确的数值,而 其它未采样点上都没有数值.然而,在实际应用中却很可能需要用到某些未采样点的值,这个时候就需要 通过已采样点的数值来推算未采样
2、点值.这样的一个过程也就是栅格插值过程.插值结果将生成一个连续 的外表,在这个连续外表上可以得到每一点的值.栅格插值包括简单栅格外表的生成和栅格数据重采样.通过栅格插值运算生成栅格外表是栅格空间分析中很重要的一局部,很多操作都将基于离散点插值生 成的外表来进行.用来插值的离散点可以是多种数据,例如,空气污染指数,土壤有机质含量,或离散高 程点.在ArcGIS 9栅格分析模块中,通过栅格插值运算生成外表主要有三种实现方式：反距离权重插值,Spat 161 fijialyEt *Siu fact Ajialyii iLcXiit cjrp olSplint . 於i粗ofBitonceDiihXi
3、 ly样条函数插值和克里格插值,如图1所示.Cell Statist lesOft we栅格插值方法下面以一组土壤元素 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、字段；(4) 在Power栏中填入进行插值计算的幂值；幂值就是距离的指数.如幂指数为2时那么进行反向距离平方插值.幂指数是一个正实数,其缺省值为2.(5) 在Search radius type栏中选择一种搜索半径设置类型；1) Variable :中选择此项时,搜索半径由下面两两个量来限制, ,Numberof points 和Maximumdistance首先在Number of points 中 输入搜索的最近点的个数(缺省值为12),然后 在Maximum distance 中输入一个限制距离.如果最近点的个数超岀限制距离,那么将会以限制距离为限制来选取较少的点；2) Fixed :中
5、选择此项时,搜索半径由下面两个量来限制,Distance 和Minimum number of points.首先在Distance 中输入搜索半径距离(缺省值是输出栅格大小的五倍),然后在Minimum number of points 中输入限制插值点个数的最小整数值.如果搜索半径距离内的点个数小于插值点个数的最小整数值,那么搜索半径自动增大.(6) Use barriers polyline为可选项,输入中断线文件.barriers是在插值中,如有某些地方岀现异常,(如某些断裂带),而要求插值时考虑到这样的因素,所设置的选项.它是一个打断外表的线 特征.这一线特征没有 Z值.悬崖,峭壁,
6、堤岸或某些障碍都是典型的barriers .barriers 限制了插值计算,它使得计算只在线的两侧各自进行.而落在线上的点那么会同时参与线两侧的计算.(7) Output cell size :指定输出结果的栅格大小；(8) Output raster:为输出结果指定目录及名称；(9) 点击0K按钮.生成结果如图 3所示.图3 IDW插值2 3 6 2 3Bu5 6 66677T口口2.样条函数插值(SPLINE)样条函数插值采用两种不同的计算方法,Regularized 和Tension.如选择Regularized ,它将生成一个平滑、渐变的外表,得岀的插值结果很可能会超岀样本点的取值范
7、围.如选择Tension,它会根据要生成的现象的特征生成一个比拟坚硬的外表,得岀结果的插值更接近限制在样本点的取值范围内.同时,计算过程中,除了需要选择不同的计算方法,还需要在每种方法中设定一个适宜的权重(weight )选择Regularized时,他决定了外表最小曲率三次导的权重.权重越高外表越光滑.可能用到的典型值有：0、0.001、0.01、1和5.选择Tension时,他决定了 Tension的权重.权重越高,外表越粗糙.可能用 到的典型值有：0、1、5和10.样条函数插值过程如下：(1) 在 Spatial Analyst下拉菜单中选择Interpolate to Raster,在
8、弹出的下一级菜单中点击Spline,出现Spline 对话框,如图4,以下所有设置如图中所示.图4 Spline 对话框(2) 在Input points的下拉菜单中选择被用来进行插值的离散点数据；在 Z value field 的下拉 菜单中选择要参加的字段；(3) 在Spline type中填入样条函数插值的类型；(4) 在Weight栏中填入一个影响插值的特征参数；(5) 在Number of points栏中输入每一个区域内用来估值点的个数.它的缺省值是12.(6) 指定输岀结果的栅格大小；(7) 为输岀结果指定目录及名称；8点击0K按钮.插值结果如图 5.图5 Spline 插值结果
9、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 在Sea
10、rch radius type栏中选择一种搜索半径设置类型；设置方法同IDW8指定输岀结果的栅格大小；9可选项,是否需要生成预测的标准误差；10为输出结果指定目录及名称；11点击0K按钮.结果如图 7.图6图7 Kriging 插值结果Eflfl TlS0 12 36 6 en- 6d- 5 6 76 6 A- 756 6 6 6 & B 6 L *4.数据重采样栅格插值除了包括简单栅格外表的生成还应包括栅格数据重采样.重采样是栅格数据空间分析中处理 栅格分辨率匹配问题的常用数据处理方法.进行空间分析时,用来分析的数据资料由于来源不同,经常会 出现栅格不同大小问题,这时为了便于分析,就需要将不
11、同的栅格大小转化为同样栅格大小,这样的一个 过程就是栅格数据的重采样过程.栅格数据的重采样主要基于三种方法：最邻近采样NEAREST,双线性ILINEAR 采样和三次卷积采样CUBIC.1最邻近采样：此重采样法用输入栅格数据中最邻近栅格值作为输岀值.因此,在重采样后的输岀栅格中的每个栅格值,都是输入栅格数据中真实存在而未加任何改变的值.这种方法简单易用,计算量小,重采样的速度最快.2双线性采样：此重采样法取待采样点x,y点周围四个邻点,在 y方向或x方向内插两 次,再在x方向或y方向内插一次,得到x, y点的栅格值.3三次卷积采样：这是进一步提升内插精度的一种方法.它的根本思想是增加邻点来获得
12、最正确插值函数.取待计算点周围相邻的16个点,与双线性采样类似,可先在某一方向上内插,如先在x方向上,每四个值依次内插四次,再根据四次的计算结果在y方向上内插,最终得到内插结果.在ArcGIS 中重采样功能是在 ArcToolbox 下来实现的.在 ArcToolbox 中,Data Management Tools 下的Raster子菜单中,选择 Resample命令图8.参数说明：1Input raster :选择输入栅格数据；2 Output raster:设置输出栅格数据路径及名称；3 Out put cell size:设置输出栅格单元大小；4Resampling technique
13、:选择重采样方法1NEAREST最邻近距离法；2BILINEAR :双线性内插法；3CUBIC:立方体法.图8重采样参数设置框二等值线绘制等值线是将外表上相邻的具有相同值的点连接起来的线,如地形图上的等高线、气温图上的等压线.等值线分布的疏密一定程度上说明了外表值的变化情况.值的变化越小的地方,等值线就越疏,反之越密.因此,通过研究等值线的疏密情况,可以获得对外表值变化的大致情况.ArcGIS中等值线的绘制方法如下：1.在Spatial Analyst 下拉菜单中选择 Surface analysis, 在弹出的下一级菜单中点击 Contour,出现Contour对话框,如图9,以下所有设置如
14、图中所示.图9等高线生成对话框2. 在Input Surface的下拉菜单中选择用来生成等高线的外表；3. 在Contour interval栏中设置等高线间距；4. 在Base contour栏中指定等高线基准高程；5. 在Z factor栏中设定高程变换系数：6在Output features栏中指定输出等高线的存放路径与文件名7.点击OK按钮.图10为某区域的等高线图,背景为该地区的地形光照晕渲图.图10某区域等高线图背景为光照晕渲图地形因子提取因子分析方法是GIS空间分析,尤其GIS数字地形分析常用的根本分析方法.不同的地形因子从不同侧面反映了地形特征性,实际应用人们提岀了各种各样的地
15、形因子.从其所描述的空间区域范围,常用的 地形因子可以划分为微观地形因子与宏观地形因子两种根本类型图11.根据提取地形因子差分计算的阶数,又可将地形因子分为一阶地形因子、二阶地形因子和高阶地形因子图12.坡向坡长坡形阖子览世起优度图11依据空间区域范围的坡面因子分类体系其中,坡度、坡向、平面曲率、剖面曲率在ArcGIS中可直接提取,其它因子的提取那么需要进行一系列的复合计算.后者的具体提取过程可以参阅相关资料,这里不作介绍.坡度-nJ图12基于提取算法的坡面因子分类体系1. 坡度的提取地外表任一点的坡度Slope 是指过该点的切平面与水平地面的夹角图13.坡度表示了地外表在该点的倾斜程度.实际
16、应用中,坡度有两种表示方式方法：(1) 坡度(degree of slope):既水平面与地形面之间夹角.(2) 坡度百分比(percent slope ):既高程增量与水平增量之比的百分数.(3) ArcGIS中坡度的提取过程为：Slope ,(4) 在Spatial Analyst 下拉菜单中选择 Surface analysis,在弹出的下一级菜单中点击 出现Slope对话框,如图14,以下所有设置如图中所示.图14坡度生成对话框(5) 在Input Surface的下拉菜单中选择用来生成坡度的外表；(6) 选择一种坡度表示方法,在此分别用两种方法做了坡度.(7) 在Z factor栏中
17、设定高程变换系数；(8) 在Output cell size栏中设定栅格大小；(9) 在Output raster栏中指定输出坡度的存放路径与文件名.(10) 点击OK按钮. 坡度结果如图15所示,坡度百分比方图16所示:图15坡度结果图 mt匚玄I图16坡度百分比图2. 坡向的提取坡向指地外表上一点的切平面的法线矢量在水平面的投影与过该点的正北方向的夹角.对于地面任何一点来说,坡向表征了该点高程值改变量的最大变化方向.在输岀的坡向数据中,坡向值有如下规定：正 北方向为0度,按顺时针方向计算,取值范围为0360 .ArcGIS中坡向的提取过程为：(1)在Spatial Analyst 下拉菜单
18、中选择 Surface analysis,在弹出的下一级菜单中点击Aspect ,出现Aspect对话框,如图17,以下所有设置如图中所示.图17坡向生成对话框(2)在 Input Surface的下拉菜单中选择用来生成坡向的外表；(3)在 Output cell size栏中设定栅格大小；(4)在 Output raster栏中指定输岀坡向的存放路径与文件名.坡向结果如图18所示图18坡向结果图up亡匚HMortK (O丝豹(2.5-67 S)East 7.5-ll2.5匸 JSouthet (1I3.5-J57 5) OfoatK (157 5-502,5) 5authft 阪 5-147
19、 3) Q47.5-292 砧 恤&92 5-3J7 5) Hkth 7 E-3E0)3. 平面曲率、剖面曲率的提取地面曲率是对地形外表一点扭曲变化程度的定量化度量因子,地面曲率在垂直和水平两个方向上分量 分别称为剖面曲率和平面曲率.剖面曲率是对地面坡度的沿最大坡降方向地面高程变化率的度量.平面曲 率指在地形外表上,具体到任何一点,指用过该点的水平面沿水平方向切地形外表所得的曲线在该点的曲 率值.平面曲率描述的是地表曲面沿水平方向的弯曲、变化情况,也就是该点所在的地面等高线的弯曲程 度.下面以一幅DEM图来说明平面曲率、剖面曲率的提取过程为：1 翻开 ArcGIS 的 Toolbox,在 Sp
20、atial 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,
21、平面曲率结果如图22：图19平面曲率、剖面曲率提取对话框图20总曲率结果图图21剖面曲率结果图|r*讷II图22平面曲率结果图四山体阴影山体阴影是根据假想的照明光源对高程栅格图的每个栅格单元计算照明值.山体阴影图不仅很好地表 达了地形的立体形态,而且可以方便的提取地形遮蔽信息.计算过程中包括三个重要参数：太阳方位角、 太阳高度角、外表灰度值.太阳方位角以正北方向为 0度,按顺时针方向度量,如 90度方向为正东方向图 23 由于人眼的视觉习惯,通常默认方位角为315度,即西北方向.太阳高度角为光线与水平面之间的夹角,同样以度为单位图为45度.ARCGIS中,默认情况下,光照产生的灰度外表的值的范
22、围为1.计算外表阴影在 Spatial Analyst 下拉菜单中选择 Surface analysis, 出现Hillshade 对话框,如图25,以下所有设置如图中所示.24.为符合人眼视觉习惯,通常默认0255.山体阴影的实现过程如下：在弹出的下一级菜单中点击HillshadeISO)图23太阳方位角度量示意图E (90)90度图24太阳高度角示意图*r*f41I刊D-： ulito 4xtrCdJL Sttf-iilacs.5Lpt.geilbxluGd St t ti 3 li cs.&临匚 Statistics.E枣包雲iiEy陆鼻1 和 Clculatri cnsJbfedL .
23、uH/F：IEGiwortO.lh 41. 图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 )方法生成二值图像,阴影区赋值为0,非阴影区赋值为1.图 27即为阴影二值图.图26外表阴影图图27阴影二值图