实例与练习 栅格数据的空间分析 ArcGIS地理信息系统基础实验练习

2019-07-07 08:09 作者:产品案例 来源:环亚ag88

  第八章栅格数据的空间分析GIS 软件应用实习材料 实例与练习 练习1. 学校选址 背景随着城市的发展,城市人口不断膨胀,随之而来的即是教育资源的短缺。在此情形下, 考虑适当增加学校数量是较好的解决方案。学校的选址问题需要考虑地理位置、学生娱乐 场所配套、与现有学校的距离间隔等因素,从总体上把握这些因素能够确定出适宜性比较 好的学校选址区。 目的通过练习,熟悉ArcGIS 栅格数据距离制图、成本距离加权、数据重分类、多层面合并 等空间分析功能,熟练掌握利用 ArcGIS 上述空间分析功能分析和结果类似学校选址的实 际应用问题的基本流程和操作过程。 school(现有学校分布图)所有原始数据存放于随书光盘的..\Chp8\Ex1\目录下。 各数据层权重比为:距离娱乐设施占0.5,距离学校占 0.25,土地利用类型和地 势位置因素各占0.125。 实现过程运用ArcGIS的扩展模块(Extension)中的空间分析(Spatial Analyst) 部分功能,具体包括:坡度计算、直线距离制图功能、重分类及栅格计算器等功能完成。 实现流程图ArcGIS 中实现学校选址分析,首先利用现有学校数据集、现有娱乐场所数据集和高程 数据派生出坡度数据以及到现有学校、娱乐场所距离数据集。然后重分类数据集到相同的 等级范围,再按照上述数据集在学校选址中的影响率赋权重值,最后合并这些数据即可创 建显示新学校适宜位置分布的地图。其间用到的ArcGIS 扩展模块(Extensions)空间分析 第八章栅格数据的空间分析GIS 软件应用实习材料 功能包括距离制图中的直线距离制图、表面分析中的坡度计算、重分类及栅格计算器等。 学校选址的逻辑过程主要包括四个部分(图1,结果中深色部分为学校候选区)。 土地利用图 高程图 距离学校距离娱乐场所 计算坡度 0.1250.125 0.5 0.25 娱乐场所分布图 学校分布图 数据准备,确定需要哪些数据作为输入,包括高程数据(dem)、土地利用数据(landuse)、现有学校数据(school)娱乐场所数据(rec_sites)。 派生数据集,从现存数据派生出能提供学校选址的原始成本数据,包括坡度数据、到现有学校距离数据集和到娱乐场所数据集。 第八章栅格数据的空间分析GIS 软件应用实习材料 重分类各种数据集,消除各成本数据集的量纲影响,使各成本数据具有大致相同的可比分类体系。各成本数据均按等间距分类原则分为 10级,级数越高适宜性 越好。邦美居家电清洗公司官网创新诚意! 给各数据集赋权重。必要的话在适宜性模型中影响比较大的数据集赋比较高权重,然后合并各数据集合以寻找适宜位置 运行ArcMap,加载SpatialAnalyst 模块,如果Spatial Analyst 模块未能激活,点 击Tools 菜单下的Extensions,选择Spatial Analyst,点击Close 按钮。 单击File 菜单下的 Open 命令,打开加载地图文档对话框,选择 E:\Chp8\Ex1\school.mxd。 设置空间分析环境。点击SpatialAnalyst 模块的下拉箭头,打开Options 对话框, 设置相关参数: 打开Options 对话框中的 General 选项卡,设置默认工作路径为: “E:\Chp8\Ex1\result\”。 打开Options对话框中的Extent 选项卡,在Analysis Extent 下拉框中选择“Same Layerlanduse”。 打开Options 对话框中的 Cell Size 选项卡,在 Analyst Cell Size 下拉框中选择 “Same Layerlanduse”。 从DEM数据提取坡度数据集。选择DEM数据层,点击Spatial Analyst 模块的下拉 箭头,选择Surface Analysis 并点击Slope,生 成slope 数据集(图2)。 从娱乐场所数据“Rec_sites”提取娱乐场直线距离数据。选择 Rec_sites 数据层, 点击 Spatial Analyst 模块的下拉箭头,选择 Distance 并点击Straight line,生成dis_recsites 数据集(图3)。 从现有学校位置数据“School”提取学校直线距离数据库。选择School数据层, 第八章栅格数据的空间分析GIS 软件应用实习材料 重分类坡度数据点击 校直线距离数据集宜性好,采用分级分为 10 级,距离学校最远 重分类土地利用数据集易发现各种土 键,选择“water”、“wetland”、“grass”,点击“deleteentries”,删除“water”、 Sptial Analyst 模块的下拉箭头,选择Distance 并点击Straight line 命令创建数据集, 得到dis_School 数据集(图4)。 重分类坡度数据学校的位置在平坦地区比 方适宜性比较差。采用等间距分级分为 10 级,在平坦的地方适宜性好,赋以较大的 适宜性值;陡峭的地区赋比较小的值,得到坡 度适宜性数据recalssslope(图5)。 重分类娱乐场直线距离数据考虑到新学校距离娱乐场所比较近 好,采用等间距分级分为 10 级,距离娱 乐场所最近适宜性最高,赋值 10;距离最远 的地方赋值 。得到娱乐场所适宜性图reclassdisr(图6)。 重分类现有学考虑到新学校距离现有学校比较远时适 元赋值10,距离最近的单元赋值1。得到重分类学校距离图reclassdiss(图7)。 重分类娱乐场所距离图在考察土地利用数据时,容 地利用类型对学校适宜性也存在一定的影 如在有湿地、水体分布区建学校的适宜性极差,于是在重分类时删除这两个选项,实 现如下: 按Ctr etland” 、“grass”。然后 根据用地类型给各种类型赋值,得到reclassland(图8),深 色部分为比较适宜区,浅色部分表示适宜性比较差,白色表示该处不允许建学校。 第八章栅格数据的空间分析GIS 软件应用实习材料 适宜区分析分类后 然后合并数 ,各个数据集都统一到相同的等级体系之内,且每个数据集中那些被认为比较适宜性的属性都被赋以比较高的值,现在开始给四种因素赋以不同的权重, 以找出最适宜的位置。 点击Spatial Analyst 下拉列表框中的Raster Calculator 命令对各个重分类后数据集的合 并计算,最终适宜性数据集的 Suit(最终适宜性)= reclassdisr(娱乐场所)* 0.5 reclassdiss(现有学校)*0.25 reclassland(土地利用数据)*0.125 reclassslope(坡度数据)得到最终适宜性数据集(图9),适宜性较高区域(深色部分)为推荐学校选址区域。 第八章栅格数据的空间分析GIS 软件应用实习材料 练习2. 寻找最佳路径 背景随着社会经济发展需求,公路的重要性日益提高。在一些交通欠发达的地区,公路建 设迫在眉睫。如何根据实际地形情况设计出比较合理的公路规划,是一个值得研究的问题。 目的通过练习,熟悉ArcGIS 栅格数据距离制图、表面分析、成本权重距离、数据重分类、 最短路径等空间分析功能,熟练掌握利用 ArcGIS 上述空间分析功能分析和结果类似学校 选址的实际应用问题的基本流程和操作过程。 新建路径的成本数据计算时,考虑到河流成本(Reclass_river)是路径成本中较关键因素,先将坡度数据(reclass_slope)和起伏度数据(reclass_QFD)按照 0.6:0.4 重合并,然后与河流成本作等权重的加和合并,公式描述如下:cost 寻找最短路径的实现需要运用ArcGIS的空间分析(Spatial Analyst)中距离制图 中的成本路径及最短路径、表面分析中的坡度计算及起伏度计算、重分类及栅格计算器等 功能完成 startPot(路径源点数据) endPot(路径终点数据) river(小流域数据) 所有原始数据存放于随书光盘的..\Chp8\Ex2\目录下 实现流程图ArcGIS 中实现最佳路径分析,首先利用高程数据派生出坡度数据以及起伏度数据集。 然后重分类流域数据、坡度、起伏度数据集到相同的等级范围,再按照上述数据集在路径 选择中的影响率赋权重值,最后合并这些数据即可得到成本数据集。在得到成本数据集之 基础上,计算栅格数据中各单元到源点的成本距离与方向数据集。最后执行最短路径即得 到最佳路径。 具体逻辑过程如下 第八章栅格数据的空间分析GIS 软件应用实习材料 运行ArcMap,加载SpatialAnalyst 模块,如果Spatial Analyst 模块未能激活,点 10学校选址逻辑过程 计算起伏 计算坡度 数据准备派生数据 权重分配成本数据 计算距离方向 小流域 分布图 高程 第八章栅格数据的空间分析GIS 软件应用实习材料 击Tools 菜单下的Extensions,选择Spatial Analyst,点击Close 按钮。 单击File 菜单下的 Open 命令,打开加载地图文档对话框,选择 E:\Chp8\Ex2\road.mxd。 设置空间分析环境。点击SpatialAnalyst 模块的下拉箭头,打开Options 对话框, 设置相关参数: 打开Options 对话框中的 General 选项卡,设置默认工作路径为: “E:\Chp8\Ex2\result\”。 打开Options对话框中的Extent 选项卡,在Analysis Extent 下拉框中选择“Same Layerlanduse”。 打开Options 对话框中的 Cell Size 选项卡,在 Analyst Cell Size 下拉框中选择 “Same Layerlanduse”。 创建成本数据集要找到到学校的最佳路径,首先需要从适宜性地图创建源数据输入及成本数据集,把 它们作为成本加权函数输入。 考虑到山地坡度、起伏度对修建公路的成本影响比较大,其中尤其山地坡度更是人们 首先关注的对象,则在创建成本数据集时,可考虑分配其权重比为:0.6:0.4。但是在有 流域分布的情况下,河流对成本影响不可低估。在此情形下,成本数据集考虑为合并山地 坡度和起伏度之后的成本,加上河流对成本之影响即可。 坡度成本数据集选择 DEM 数据层,点击 Spatial Analyst 下拉列表框,选择 Surface Analysis 并点击 slope,生成坡度数据集,记为Slope。 11坡度成本数据 选择 Slope 数据层,点击 Spatial Analyst 下拉键头,选择Reclassify 命令实施重分类。 对坡度数据集实施重分类的基本原则是:采用 等间距分为 10 级,坡度最小一级赋值为 最大一级赋值为10,得到图11 所示坡度成本 数据(reclass_slope)。 起伏度成本数据集选择DEM数据层,点击Spatial Analyst 下拉列表框,选择Neighborhood Statistics,设 置如图12 所示参数设置,点击Ok 按钮,生成起伏度数据层,记为QFD。 选择QFD 数据层,点击Spatial Analyst 下拉键头,选择Reclassify 命令,按10 得到图13所示地形起伏成本数据(reclass_QFD)。 第八章栅格数据的空间分析GIS 软件应用实习材料 图12 生成起伏度 图13 起伏度成本数据图 河流成本数据集选择River 数据层,点击Spatial Analyst 下拉键头,选择 Reclassify 命令,按照河流 等级如下进行分类:4 级为10;如此依次为 8,5,2,1,生成图 14 所示河流成本 (reclass_river)。 14河流成本数据 加权合并单因素成本数据,生成最终成本数据集。 点击 Spatial Analyst 下拉箭头,选择 Raster Calculator 命令合并数据集,计算公式 如下: cost reclass_rough(重分类起伏度数据)*0.4) 图15 最终成本数据 16计算成本权重数据对线 所示最终成本数据集(cost),其中深色表示成本高的部分。 计算成本权重距离函数点击Spatial Analyst 模块下拉箭头,选择Distance 中的Cost Weighted,设置参数如图 第八章栅格数据的空间分析GIS 软件应用实习材料 19最短路径对线最佳路径图(黑色粗线, 尖点为源点。 Analyst下拉框,选择Distance 中的Shortest Path,设置参数如图19 所示, 点击 的最短路径图(图20,其中黑色粗线部分为确定的路径)。 熊猫分布密度制图大熊猫是我国国家级珍惜保护动物,熊猫的生存必须满足一定槽域(独占的猎食与活 。因此,科学准确的分析熊猫的分布情况,对合理制定保护措施和评价保护成效 点击OK 按钮。生成图17 所示成本距离图,其中浅色为源点;图18 所示成本方向图, 求取最短路径点击Spatia 17成本距离图 18成本方向图 OK 按钮,生成最终 条件具有重要的意义。 10 第八章栅格数据的空间分析 GIS 软件应用实习材料 目的通过练习,熟悉ArcGIS 密度制图函数的原理及差异性,掌握如何根据实际采样数据特 提供的密度制图功能和部分空间分析,灵活变通,制作符合实际需要的密度 猫活动具有一定的槽域范围,一个槽域范围只有一个或一对熊猫,在此练习中, 假设熊猫槽域半径为5km; 积。假定熊猫活动范围分布满足以采样点为中心的泰森多边 过,相同足迹只记载一次。数据存放于随书光盘的..\Chp8\Ex3\目录下。 提取熊猫的槽域范围,然后用理 果Spati能激活,点击Tools 菜单下的 Exte y.mxd。拉箭头, 结合 ArcG 虽然一个采样点代表一个熊猫,但由于熊猫的生存具有确定槽域特征,不同的采样点具有不同的空间控制面 如何将这一信息加入密度分布图是本练习的重点。 在野外实采的熊猫活动足迹数据的基础上,以每个熊猫槽域范围为权重,运用ArcGIS 中的区域分配功能和密度制图功能制作该地区熊猫 计算原理首先利用栅格数据空间分析模块提供的区域分配功能 域面积(际槽域面积,作为采样点的加权值(记为Power字段),生成熊猫分布密度图。 运行ArcMap,加载SpatialAnalyst 模块,如 al Analyst 模块未 nsio s,选择Spatial Analyst,点击Close 按钮。在 Options General页面中在设置默认工作路径,此 处假定为“E:\chp8\ex3\Result\”,图21。 单击File菜单下的Open 命令,打开加载地图 文档对线\ XMDensit 熊猫的生存必须满足一定的槽域范围。选择熊猫活动足迹数据图层,点击Spatial Analyst 模块的下 选择Distance 并点击Allocation,设置参数如 图22 所示;点击OK,生成熊猫槽域范围图(图23, 白色区域没有熊猫出现)。 图21 工作路径设置 11 第八章栅格数据的空间分析 GIS 软件应用实习材料 图22 区域分配对线 槽域分配图 选择FP数据层,点击鼠标右键并选择Open Attribute Table 命令,打开FB 属性表。 点击FB 属性表右下角Option 按钮的下拉箭头,选择Export 命令(图24),导出FB 属性 数据表,参数设置如图 25 所示,点击确定按钮;当出现提示是否需要加载该数据表,选 图24FB 属性表输出 图25 参数设置 12 第八章栅格数据的空间分析 GIS 软件应用实习材料 选择熊猫活动足迹数据图层(XMPoint),点击鼠标右键,选择Join Relates菜单下的Join,弹出数据连接对线 所示,点击OK 按钮,完成熊猫采 …,参数设置如图27所示,提取密 据与槽域范围数据的连接。图26 采样数据与槽域范围数据连接 选择熊猫活动足迹数据图层(XMPoint),点击鼠标右键并选择OpenAttribute Tabl 命令,打开XMPoint 属性表;点击属性表中的XMPoint.CaoYuArea 字段名,点击鼠标右键, 选择 Calculate Values…;在 Field Calculate 对话框计算公式输入视窗输入计算公式:nt]*500*500,500 为生成FB 数据层时设置的栅格大小。 图27密度计算 Poin .Power 字段名,点击鼠标右键,选择 Calculate Values…;在Field Calculate 对话框计 [XMPoint.CaoYuArea],3.1415926*5000*5000 为假定的最大槽域面积, 计算每个采样点的权重值,作为计算密度的样本 点击SpatialAnalyst 模块的下拉箭头, 选择Density 13 第八章栅格数据的空间分析 GIS 软件应用实习材料 上述密度以平方米为面积单位,数据值太小。点击SpatialAnalyst 模块的下拉箭 头,选择 练习 协调发展是保证国民经济健康持续稳定增长的关键。GDP是分析各地区经 目的供了三种空间插值方法,每种插值方法在原理上和应用上都大不相同,在 发展具有一定的连带效应和辐射作用。以该地区各区域年GDP 数据为依 数据(GDP.shp),数据范围:4601万元~132630 万元。数据存放于随书 光盘 Raster Caculate…,输入计算公式:XMDensity10 10000000,将面积单位换算为10 平方公里,结果如图28 所示。 图28 熊猫密度图(单位:10 平方公里) GDP区域分布图的生成与对比 背景 各地区经济 展状况的重要指标,科学准确分析各地区GDP 空间分布特征,对制定有效措施,指导 经济协调发展具有重要参考价值。 背景 ArcGIS中提 过具体实例练习如何利用IDW内插方法和Spline 内插方法进行GDP 空间分布特征的 分析,以此来引导读者对空间插值有一个更深刻的认识。 经济的,采用IDW和Spline 内插方法创建该地区GDP 空间分异栅格图; 分析两种内插方法生成的GDP 空间分析图的差异性,简单 数据:该地区的GD 的..\Chp8\Ex4\目录下 14 第八章栅格数据的空间分析 GIS 软件应用实习材料 析。选择ID ,加载SpatialAnalyst 模块,如果Spatial Analyst 模块未能激活,点 st下拉菜单中选择Options 选项,在Options 中的General 页面中 IDW内插方GDP 数据 内插方法 两种方法之间的对比,通过空间分析中 Raster Calculate 命令来进行分 W(Power=2),Spline(Regularized,Weight=0.01)。在Raster Calculate 中计算Abs(IDW 运行ArcMapools 菜单下的Extensions,选择Spatial Analyst,点击Close 按钮。 单击File 菜单下的 Open 命令,打开加载地图文档 hp8 x4\GDP.mx 在SpatialAnaly 置默认工作路径,此处假定为“E:\chp8\ex4\Result\”,并设置Analysis mask 为board.shp。 相减并取绝对值 IDW(Power=2) IDW(Power=5) 总结 图29 IDW内插试验流程图 相减并取绝对值 Weight=0 Weight=0.01 GDP 数据 Regularized Tension 相减并取绝对值Weight=0 Weight=5 图30Spline 内插试验流程图 15 第八章栅格数据的空间分析 GIS 软件应用实习材料 在SpatialAnalyst 下拉菜单 中选择 Interpolate Raster,在弹出 InverseDistance Weighted; 设置Zvalue field 为GDP; 设置Power 为2;设置Output cell size 500;其他参数不变,点击OK, 进行计算。 将Power值改为5,重复上 述步骤。 在SpatialAnalyst 下拉菜单 中选择 Raster Calculator Abs((Power=2)—(Power=5))。 IDW插值结果如图31 所示。 由于IDW 是一个加权距离 平均,平均值不能大于输入最高值或 是小于输入最低值,因此输出的结果 数据中,每一栅格值均处于采样数据 的最大值与最小值范围之内。的输出 值限制在用于插值的输入值的范围 幂指数不同,IDW的插值结果是不同的。幂指数越大,较远的点 对于输入的影响越小,即幂指数越 高,其局部影响的程度越高。 IDW属于精确性(生成表面通过采样点)插值方法。因此,即使 给出不同的幂指数值,在采样点上它 们差值的绝对值比较小,差值比较大 的地方出现在突然采样数据变化比 较剧烈和频繁的区域有。此外,如果 采用同样搜索半径和参加内插的最 小采样点数,生成表面在最近的几个 采样的几何中心区域具有比较稳定 的内插值,受幂指数的影响不是很明 Abs((Power=2)—(Power=5))图31 IDW插值结果 16 第八章栅格数据的空间分析 GIS 软件应用实习材料 Spline内插法 在SpatialAnalyst 下拉菜单 中选择Interpolate Raster,在弹出 的下一级菜单中点击Spline; 设置Zvalue field 为GDP; 设置 Spline type 0;设置Output cell size 为500;其他参数不变,点击OK, 进行计算; Weight值改为 0.01,重 复上述步骤; 设置Splinetype 为Tension, 并分别取Weight 在SpatialAnalyst 下拉菜单 中选择 Raster Calculator Abs((Weight=0)-(Weight=0.01))和 Tension Abs((Weight=0)—(Weight=0.01))。 图32 为采用Regularized 型时插值结果。修改Spline type Tension,Weight分别取0 和5,其它 参数不变,重复上述步骤,计算插值 结果。图 33 为采用 Tension Spline插值结果光滑,它的 插值结果范围将远远超出原值范围; ized选项,它 结果要Tension 对于Regular 次求导(坡度),二次求导(坡 度变化率)和三次求导(二次求导的 变化率)合并为其最低估计计算,它 决定了表面最小曲率表达的三次导 的权重。weight 值越高则生成的表面越光滑; Abs((Weight=0)—(Weight=5))图32 Spline Regularized 插值结果 17 第八章栅格数据的空间分析 GIS 软件应用实习材料 对于Tension选项,它 仅使用一次和二次导数。weight 值越高则生成的表面越粗糙,但 是表面的值越接近控制点; Tension选项中分别对两次权重 差值取绝对值的结果图中可以 看出,在采样点处差值的绝对值 都比较小,而在采样点之间比较 大,尤其是在采样点比较稀疏的 地方差值绝对值相对会大一些。 Abs((Weight=0)—(Weight 图33Spline Tension 插值结果 18 第八章栅格数据的空间分析 GIS 软件应用实习材料 IDW与Spline对比分析 SpatialAnalyst 下拉菜单中选择 Raster Calculator,输入公式:Abs(IDW Spline),点击Evaluate 来完成运算。结果如图34 所示。 图34 IDW与Spline 对比分析 从上一组图可以看出每两组值的差的绝对值在采样点处都比较小,在非采样点处相对较大; 在采样点比较密集的地方差的绝对值较小,在采样点分布比较稀疏的地方差的绝对值较大。 结合采样点处数据大小,可以看出在采样点变化比较大的地方差的绝对值较大,变化比较均衡的地方差的绝对值较小。 IDW与Spline插值结果之间的存在着比较大的区别。在采样点处变化都比较小, 在非采样点处相对较大;在采样点比较密集的地方内插值变化较小,在采样点分布比较稀 疏的地方变化较大;结合采样点处数据大小,可以看出在采样点分布比较规律或采样点值 比较均一的地方,内插结果的变化普遍较小,而在采样点分布零散并且采样值变化频繁、 剧烈的区域,内插结果的变化普遍较大。 IDW是一个加权距离平均,其每一栅格的输出值限制在采样点的输入值的范围 内,因此,对如山脊和沟谷这样的极端地形,如果没有采样点,IDW不会生成这些地形。 当取样点足够密时,IDW对局部变化具有非常好的效果。Spline 是基于生成具有连续的二 阶导数和最小的平方曲率的插值方法,所以它适合那些空间连续变且光滑的表面的生成。 19 第八章栅格数据的空间分析 GIS 软件应用实习材料 练习5. 提取山顶点、凹陷点 背景山顶点指哪些在特定邻域分析范围内,该点都比周围点高的区域。凹陷点指哪些在特 定邻域分析范围内,该点都比周围点低的区域区域。DEM 数据提供了丰富的地形信息, 基于DEM数据如何正确有效的提取这些地形特征信息,在数字地形分析中具有重要意义。 目的通过等高线、山顶点、凹陷点的提取和配置、引导读者熟练掌握利用ArcGIS 栅格数据 空间分析中等高线的提取、栅格数据邻域分析和窗口计算功能,完成栅格数据表面分析。 数据:黄土丘陵地区1:10000DEM数据。数据存放于随书光盘的..\Chp8\Ex5\目录中。 要求图35 设置工作路径 应用栅格数据空间分析模块中的等高线提取功能,分别提取等高距为15 米和75 米的等高线图, 并按标准地形图绘制等高线方法绘制等高线,作为山 顶点、凹陷点空间分布的背景图; 运行ArcMap,加载SpatialAnalyst 模块, 如果Spatial Analyst 模块未能激活,点击Tools 菜单下 的Extensions,选择Spatial Analyst,点击Close 按钮。 在Options 中的General 页面中在设置默认工 作路径,此处假定为“E:\chp8\ex5\Result 图36 等高线提取对线 所示,提取等高 米的等高线,输出文件名为Contour75; 35 单击File菜单下的Open 命令, hp8\HillTop.mxd。 点击SpatialAnalysis 下拉箭头, 选择 Sueface Analysis 子菜单并 点击 Contour,参数设置如图 15 等高线; 修改Contourinterval 为75 20第八章栅格数据的空间分析 GIS 软件应用实习材料 单击Contour15数据层图例,选择显示颜色为灰度60%(图37),点击OK 图37选择图例颜色 点击SpatialAnalysis 下拉箭头,选择Sueface Analysis 子菜单并点击Hillshade, 设置输出文件名为Hillshade,其它参数取默认值,提取该地区光照晕渲图,作为等高线三 维背景; 点击SpatialAnalysis 下拉箭头,选择Raster Calculator…,输入计算公式:Back 选择Back数据层,在弹出的属性对话框的Display 属性页设置透明度60%, 在Symbology 属性框中设置其显示颜色为Gray50%,点击OK 按钮; 图38 三维立体等高线、Back、Hillshade次序放置数据层,生成三维立体等 21 第八章栅格数据的空间分析 GIS 软件应用实习材料 高线) 点击Spatial Analysis 下拉箭头,选择 Neighborhood Statistics,设置参数如图39 所示, 点击OK 按钮,提取11*11 分析窗口最大值; 图39 计算最大值 (11) 修改图39 中Statistic type 为Mininum, Output Raster(输出文件名)为MinPoint,其它参数 不变,点击OK,提取11*11 分析窗口最小值; (12) 点击 Spatial Analysis 下拉箭头,选择 Raster Calculator…,输入计算公式:SD 0,提取山顶点区域;(13) 选择 SD 数据层,点击 Spatial Analysis 下拉箭头,选择Reclassify,设置参数如图40 所示, 重分类SD 数据。 (14) 选择 RE_SD 数据层,点击 Spatial Analysis 下拉箭头,选择Convert Features…,设置参数如图41 所示,输出矢量山顶 点数据,见图42。 图41 格式转换 图40 重分类SD 数据 图42 山顶点分布图 (15) 点击Spatial Analysis 拉箭头,选择RasterCalculator…, 输入计算公式:GD 0,提取凹陷点区域;(16) 选择 GD 数据层,点击 Spatial Analysis 下拉箭头,选择 Reclassify,设置参数同图40 所示, 22 第八章栅格数据的空间分析 GIS 软件应用实习材料 输出文件名改为Re_GD,重分类GD 数据。 (17) 选择RE_SD 数据层,点击Spatial Analysis 下拉箭头,选择Convert 子菜单并 点击Raster Features…,设置参数同图41所示,输出文件名改为SL_GD,输出矢量凹陷 点数据(图43)。 (18) Sl_SD、Sl_GD、Contour75、Contour15、Back、Hillshade从上到下依次放 置上述数据层,合并显示山顶点、凹陷点提取结果(图44)。 图43 凹陷点分布图 图44 山顶点、凹陷点分布图 23