ArcGIS空间分析教程(汤国安)实习报告

1、精选优质文档-倾情为你奉上地理信息系统原理期末实习报告专 业： 资源环境与城乡规划专业 班 级： 1201班 姓 名： xxxxxx 学 号： xxxxxxxxxx 指导教师： xxxxxxxxxxxxxxxx 二零一五年一月专心-专注-专业目录1.重要概念····························

2、83;·················12.实验目的·······························&#

3、183;··············23.数据来源··································&

4、#183;···········24.要求·····································&#

5、183;············25.实验内容····································&

6、#183;·········26.附图·······································&#

7、183;··········3 实验一：数据处理白水县··························3实验二：寻找最佳路径·········&#

8、183;······················10实验三：土壤稳定性评估·························&

9、#183;····20 实验四：土壤侵蚀性分析建模··························32 实验五：水文分析··············

10、3;·····················44 实验六：找出某种珍贵药材的生长区域··················46实验七：地形鞍部的提取·····&

11、#183;························48实验八：沟谷网络的提取·······················

12、·······49实验九：TIN及DEM的生成与应用·······················50实验十：缓冲区分析的应用··············

13、83;·············577.实习心得···································&

14、#183;··········611.重要概念缓冲区：缓冲区是地理空间目标的一种影响范围或服务范围，具体指在点、线、面实体周围一定范围。空间叠置分析：指用来提取空间隐含信息的方法之一。它是将代表不同主题的各个数据层面进行叠置产生一个新的数据层面，叠置结果综合了原来两个或多个层面要素所具有的属性。通用土壤侵蚀方程：是结合了美国20世纪30年代起的8000多个土壤侵蚀试验观测点资料统计总结提出的一种用于计算土壤侵蚀强度的公式。水文分析：是DEM数字地形分析的一个重要方面，使用水文分析工具，基于DEM

15、对地形进行分析。地图：指依据一定的数学法则，使用制图语言，通过制图综合在一定载体上，表达地球或其他天体上各种事物的空间分布、联系及时间中的发展变化状态而绘制的图形。栅格数据空间分析方法：是指针对于栅格数据的空间分析方法，具有自动分析处理较简单，分析处理模式化很强的特征。TIN：是一种复合矢量模型，它采用一组互不叠置的三角形来近似表示地形。DEM：是一种数字模型，等间距高程数据以栅格格式排列。空间分析：是基于地理对象的空间布局的地理数据分析技术。土地信息系统：土地信息系统是综合应用地理信息系统和管理信息系统，对人类在土地利用过程中产生的土地数据进行采集、存储、检索、分析和管理的信息系统。2.实验

16、目的掌握ARCGIS10空间分析方法，熟练掌握Model Builder的创建和使用，掌握利用ARCGIS进行水文分析，缓冲分析，空间叠置分析，邻域分析、重分类等方法与技术。掌握利用ARCGIS进行DEM及TIN的创建，掌握地图的编制、整饰及输出。3.数据来源汤国安地理信息系统空间分析实验教程，西安科技大学测绘学院地理信息系统实习教案。4.要求（1）部分实验需要写出步骤，步骤要图文并茂的反映操作流程，截图并用文字说明。（2）所有附图均需添加姓名、学号、班级。（3）实验报告A4纸正反面打印，左侧装订。（4）重要概念查阅相关文献，电脑输入，不用手工填写。（5）附图黑白打印即可。（6）实习报告后附实

17、习心得体会，1000左右，手写。（7）此部分内容占总成绩25%。5.实验内容实验一：用Model Builder进行数据处理与分析，实例：汤国安P130。实验二：汤国安地理信息系统空间分析实验教程（P290寻找最佳路径）实验三：土壤稳定性评估6. 附图实验一：数据处理白水县1:25万的矢量数据vector一、流程图：2幅1:25万数据DEM1和DEM2 合并提取白水县行政区划范围1幅1:25万DEM数据白水县DEM数据裁接投影 投影数据图1 白水县数据更新变换流程图2、 步骤：（1）数据加载与环境设置：1）打开Arcmap，添加图层vector，如下图: 图2 2）打开Arctoolbox，右

18、击，点击“环境”，设置工作空间如图3： 图3（2）利用模型构建器构建模型： 1）模型环境设置：点击打开模型构建器，执行菜单命令模型属性，设置如图4。选择环境处理范围范围与图层vector相同，如图5所示： 图4 图5 2）白水县行政范围提取：将Arcmap目录下的vector拖至ModelBuilder中，在Arctoolbox中，选择分析工具提取筛选工具，将其拖至ModelBuilder窗口中，双击筛选框，在出现的对话框中选择输入栅格为vector，输出栅格命名为Vector_Select1.shp，如图6所示，确定，结果如图7： 图6 图7 3）白水县DEM数据的拼接：在ModelBuil

19、der中点击添加数据DEM1和DEM2，选择数据管理工具栅格栅格数据集镶嵌至新栅格工具，将其拖至ModelBuilder窗口，双击镶嵌至新栅格，在出现的对话框中设置如图8，确定，结果如图9所示： 图8 图9 4）利用白水县范围对DEM裁切：在Arctoolbox中选择spacial analyst提取分析按掩膜提取工具，将其拖至ModelBuilder中，如图10所示。双击按掩膜提取工具框，在出现的对话框中将DEM作为输入栅格，输入栅格数据或要素数据选为vector-select.shp，输出栅格为Extract-DEM1,如图11，裁切结果如图12： 图10 图11 图12 5）白水县DEM

20、的投影变换：选择数据管理工具投影与变化栅格投影栅格工具，将其拖至ModelBuilder窗口中，双击打开投影栅格对话框，将Extract-DEM1设为输入栅格，输出栅格命名为Ext-DEM1-Prj，点击输出栅格右边的按钮，进入空间参考属性对话框，点击选择按钮，浏览坐标系，选择Projected Coordinate SystemsGauss KrugerXian 1980Xian 1980 GK Zone 19.prj，重采样技术选择NEAREST,如图13所示，确定。 图136）运行模型:点击运行按钮运行模型，如图14所示： 图14三、实验结果 （1）最终模型： 图15白水县数据更新变换模

21、型图（2）专题地图： 图16 白水县DEM图实验二：寻找最佳路径一、流程图：小流域分布图river高程数据DEM起伏度数据 坡度数据起伏度重分类数据坡度重分类重分类图分配权重0.40.6成本数据回溯连接图成本距离图最佳路径图图1寻找最佳路径流程图二、实验步骤： （1）激活工具：打开Arctoolbox，激活spacial analyst空间分析和3D分析扩展模块（如图2）：自定义扩展模块 图2（2） 设置环境，加载数据：1）打开Arcmap，添加river.shp、startpot、endpot和DEM数据，如图3；打开Arctoolbox，右击选择环境，设置工作空间如图4： 图3 图42）点

22、击按钮打开模型构建器，执行菜单模型模型属性常规，设置模型名称及标签，如图5左所示，确定。菜单模型图属性符号系统样式2，如图5右，确定。将数据river.shp、startpot、endpot和DEM数据拖至模型构建器窗口中，如图6： 图5 图6（3） 利用modelbuilder构建模型：1）生成坡度成本数据集：选择Spacial Analyst工具表面分析坡度，将坡度工具拖拽至模型构建器窗口中，双击坡度工具框，在出现的对话框中，输入dem，输出栅格命名为slope，如图7左所示，确定。运行如图7右： 图6 选择Spacial Analyst工具重分类重分类，将其拖至模型构建器，双击重分类对话

23、框，输入slope，字段为Value，点击分类按钮，分类方法为相等间隔，类别为10，如图8所示，确定。 图8 2）生成起伏度成本数据集：选择Spacial Analyst工具邻域分析焦点统计，将其拖至模型构建器中，双击焦点统计工具框，输入dem输出栅格命名为QDF，具体设置如图9所示，确定。右击重分类工具框，选择运行。 图9 选择Spacial Analyst工具重分类重分类，将其拖至模型构建器中，双击重分类对话框，输入QDF数据层，点击分类按钮，选择相等间隔分类方法，类别选为10，如图10，确定。分类结果如图11，输出栅格为Reclass-QDF，确定。 图10 图11 3）生成河流成本数据

24、集：选择Spacial Analyst工具重分类重分类，将其拖至模型构建器中，双击重分类工具框，选择river数据层，其他设置如图12，确定。 图12 4）加权合并单因素成本数据，生成最终成本数据集：选择Spacial Analyst工具地图代数栅格计算器,并将其拖至模型构建器中，双击栅格计算器，计算公式为cost=Recriver（重分类流域数据）+（Reclassslop（重分类坡度数据）*0.6+ReclassQFD（重分类起伏度数据）*0.4），如图13所示，确定。 图13 5）计算成本权重距离函数：选择Spacial Analyst工具距离分析成本距离，将其拖至模型构建器中，双击成本

25、距离，输入栅格为startpot，输入成本栅格为cost，输出距离栅格为CostDis-star1，输出回溯连接栅格为blanklink，如图14所示： 图14 6）求取最短路径：选择Spacial Analyst工具距离分析成本距离，将其拖至模型构建器中，双击成本路径工具框，在出现的对话框中输入栅格为endpot，目标字段为Id，输入成本距离栅格为CostDis-star1，输入回溯连接栅格为blanklink，输出栅格为CostPat-endP1，如图15所示，确定。 图157）运行模型：点击运行按钮运行模型如图16： 图16（4）加载模型运行结果图：在Arcmap中加载模型运行后的结果，

26、如图18所示： 图173、 实验结果（1）模型： 图18 最佳路径模型图（2）专题图： 图19最佳路径专题图实验三：土壤稳定性评估一、流程图：DEM土地利用矢量图坡向图坡度图土地利用栅格图重分类重分类重分类阴、阳坡等级图阴阳坡等级图级土地利用等级图坡度等级图权重赋予 土壤稳定性数据（加权平均）土壤稳定性专题图图1土壤稳定性评估流程图2、 实验步骤：（1）激活工具：打开Arctoolbox，激活spacial analyst空间分析和3D分析扩展模块（如图2）：自定义扩展模块 图2（2）设置环境，加载数据： 1)打开Arcmap，添加DEM和landuse.shp数据，如图3。打开Arctool

27、box，右击选择环境，设置工作空间如图4所示： 图3 图42)点击按钮打开模型构建器，执行菜单模型模型属性常规，设置模型名称及标签，确定。菜单模型图属性符号系统样式2，确定，如图5所示；将Arcmap中的DEM及laduse.shp数据拖至模型构建器窗口中，如图6： 图5 图6（3） 构建模型： 1）提取坡度：打开spacial analyst工具表面分析坡度，将坡度工具拖拽至模型构建器窗口中，双击坡度工具框，在出现的对话框中将dem设为输入栅格，输出栅格为slope，如图7，确定。右击坡度工具框，选择运行，如图8所示： 图7 图82）坡度重分类：选择spacial analyst工具重分类重

28、分类工具并将其拖拽至模型构建器中，双击重分类工具框，在出现的对话框中选择slope作为输入栅格，值字段选择value，如图9，点击分类按钮，类别选为6,方法选为手动，修改中断值如图10，确定。 图9 图10 打开将重分类对话框中的新值一栏依次改为10、8、7、5、3、1，输出栅格命名为Re-slope，如图11所示，确定。模型如图12。 图11 图123）提取坡向：打开spaciaal anlyst工具表面分析坡向，将坡向工具拖至模型构建器窗口，双击坡向工具框，在工具对话框中输入l，dem，输出栅格命名为Aspect，如图13，确定。右击坡向工具框，选择运行，如图14所示： 图13 图144）

29、坡向重分类：选择spacial analyst工具重分类重分类工具并将其拖拽至模型构建器中，双击重分类工具框，在出现的对话框中选择Aspect作为输入栅格，值字段选为Value，点击分类，类别选为4，方法选为手动，中断值改为0、90、270、360，如图15，确定。将新值一栏改为5、10、1、10，输出栅格命名为Re-Aspect，如图16所示。模型如图17： 图15 图16 图17 5）土地利用图的转换：选择转换工具转为栅格面转栅格，将其拖至模型构建器，双击面转栅格工具，在其对话框中输入landuse，值字段LANDUSE,输出栅格命名为landusegrid，像元大小设为100，如图18，

30、确定。右击面转栅格，运行，结果如图19： 图18 图19 6）对landusegrid重分类：选择spacial analyst工具重分类重分类工具并将其拖拽至模型构建器中，双击重分类工具框，在出现的对话框中选择landusegrid作为输入栅格，值字段选为LANDUSE,新值分别改为2、8、6、10，输出栅格命名为Re_land，如图20，确定。模型如图21： 图20 图217)空间叠加分析：选择spacial analyst工具叠加分析加权总和，并将其拖至模型构建器窗口中，双击打开加权总和对话框，逐个添加Re-slope、Re-Aspect和Re_land，字段均选为VALUE,Re-sl

31、ope的权重为0.3,Re-Aspect的权重为0.3，Re_land的权重为0.4，输出栅格为stability，如图22所示，模型如图23，点击菜单中的按钮运行。 图22 图23(4) 添加模型运行结果数据：在Arcmap中加载叠加分析结果stability数据，如图24： 图24（5）土壤稳定性的分级显示：在Arcmap中右击stability打开属性对话框，选择符号系统，显示选为已分类，点击分类，设置分类方法为手动，类别为3，将中断值设为4、7、10，如图25所示；将每个范围对应的标注分别改为“不稳定”“较稳定”和“很稳定”，并选择合适的拉伸色带，如图26： 图25 图263、 实验结

32、果图27土壤稳定性专题图实验四：土壤侵蚀性分析建模一、流程：DEM矢量数据土壤类型栅格数据植被类型矢量数据DEM栅格数据植被类型栅格数据坡度图重分类坡度图50%25%25%加权叠加土壤侵蚀栅格图图1土壤侵蚀性分析建模流程图二、实验步骤： （1）加载数据：打开Arcmap并添加数据studyarea.shp 、vegetation.shp、soilsgrid.shp,如图2所示：图2（2） 激活工具：打开Arctoolbox，激活spacial analyst空间分析和3D分析扩展模块（如图3）：自定义扩展模块； 图3（3） 数据前期处理： 1）右击vegetation打开属性propertie

33、ssymbology Unique value，值字段选择vegtype,选择好色带，点击Add All Values，如图4所示，点击ok； 图4 2）调整图层顺序为studyarea，vegtation，soilsgrid如图5，保存地图文档； 图5（4） 利用ModelBuilder构建模型： 1）设置属性与环境：点击打开模型构建器ModelBuilder如图6所示，执行菜单命令modelmodelpropertis设置如图7，选择environment，设置如图8，点击value设置处理范围为Same As Layer StudyArea,如图9，确定； 图6 图7 图8 图9 在mo

34、delbuiler窗口中执行菜单ModelDiagram Propertis,选择style2，如图10； 图102）加载数据：将ArcMap中的vegtation，soilsgrid拖至ModelBuilder窗口中。3）将DEM数据转为栅格数据：从Arctoolbox中将Conversiontools to raster下的DEM to Raser拖至ModelBuilder窗口中，如图11所示。在modelbuilder中双击DEM to raser,在出现的工具设置对话框中指定输入USGS DEM文件为：elevation.dem,选择输出路径，确定，如图12： 图11 图12 4）生

35、成坡度数据：将Arctoolbox中3D Analyst ToolsRaster Surfaceslope拖至modelbuilder窗口中在 ModelBuilder 窗口中，点击添加连接按钮将派生数据图框DEMToRa_elev1与工具图框slope连接在一起。完成后效果如下图13所示： 图13右键点击图框output raster将其改名为: 坡度图5）坡度数据重分类：将Arctoolbox中Spatial Analyst Tools Reclass下的工具Reclassify 拖放到ModelBuilder窗口中，点击按钮将派生数据图框坡度图与工具图框Reclassify连接在一起，如

36、图14所示： 图14在 ModelBuilder窗口中，双击工具图框Reclassfy，对坡度图重分类，输入栅格为坡度图，字段为value，点击load装入数据表，将坡度分为10级，具体设置如图15所示： 图15 8）植被类型矢量数据转栅格：将Arctoolbox的Conversion ToolsTO Rasterfeature to raster拖放到ModelBuilder窗口中，在 ModelBuilder窗口中，点击按钮将数据图框vegetaion与工具图框feature to raster连接在一起，如图16所示： 图16 双击与Vegetaion相连的工具图框Feature to

37、Ratser，在出现的对话框中，设置字段为VEGTYPE,选择输入栅格Vegetaion和输出路径，确定，如图16；右击Feature to Raster，执行命令Run，如图17； 图16 图17 9）加权叠加分析： <1>添加工具：将Arctoolbox中Spatial Analyst Tools OverlayWeighted Overly拖至modelbuilder中，如图18； 图18 <2>添加数据：双击Weighted Overly，在出现的设置框中点击按钮，在出现的框中输入soilsgrid，字段选为S-value，如图19，确定；同理，添加植被栅格（字段值为vegtype，如图20）和Reclass-slop1（字段值为value，如图21），添加结果如图22所示： 图19 图20 图21 图22 <3>设置因子（如图23所示）： Soilsgrid：根据不同土壤类型对土壤侵蚀危险性的影响力，不同的土壤类型给定不同的的数值，数值1表示改天土壤侵蚀危险度较低，9表示较高； 植被栅格：根据不同植被类型对土壤侵蚀危险性的影响力，不同的植被类型给定不同的的数值，数值1表示改天土壤侵蚀危险度较低，9表示较高（wat