ArcGIS空间分析教程实习报告.docx

ArcGIS空间分析教程实习报告.docx

《ArcGIS空间分析教程实习报告.docx》由会员分享，可在线阅读，更多相关《ArcGIS空间分析教程实习报告.docx（21页珍藏版）》请在冰豆网上搜索。

ArcGIS空间分析教程实习报告

"地理信息系统原理"期末实习报告

专业：

资源环境与城乡规划专业

班级：

1201班

姓名：

******

学号：

**********

指导教师：

****************

二零一五年一月

1.重要概念··············································1

2.实验目的··············································2

3.数据来源··············································2

4.要求··················································2

5.实验容··············································2

6.附图··················································3

实验一：

数据处理——白水县··························3

实验二：

寻找最正确路径································10

实验三：

土壤稳定性评估······························20

实验四：

土壤侵蚀性分析建模··························32

实验五：

水文分析····································44

实验六：

找出某种珍贵药材的生长区域··················46

实验七：

地形鞍部的提取······························48

实验八：

沟谷网络的提取······························49

实验九：

TIN及DEM的生成与应用·······················50

实验十：

缓冲区分析的应用····························57

7.实习心得··············································61

1.重要概念

缓冲区：

缓冲区是地理空间目标的一种影响围或效劳围，具体指在点、线、面实体周围一定围。

空间叠置分析：

指用来提取空间隐含信息的方法之一。

它是将代表不同主题的各个数据层面进展叠置产生一个新的数据层面，叠置结果综合了原来两个或多个层面要素所具有的属性。

通用土壤侵蚀方程：

是结合了美国20世纪30年代起的8000多个土壤侵蚀试验观测点资料统计总结提出的一种用于计算土壤侵蚀强度的公式。

水文分析：

是DEM数字地形分析的一个重要方面，使用水文分析工具，基于DEM对地形进展分析。

地图：

指依据一定的数学法那么，使用制图语言，通过制图综合在一定载体上，表达地球或其他天体上各种事物的空间分布、联系及时间中的开展变化状态而绘制的图形。

栅格数据空间分析方法：

是指针对于栅格数据的空间分析方法，具有自动分析处理较简单，分析处理模式化很强的特征。

TIN：

是一种复合矢量模型，它采用一组互不叠置的三角形来近似表示地形。

DEM：

是一种数字模型，等间距高程数据以栅格格式排列。

空间分析：

是基于地理对象的空间布局的地理数据分析技术。

土地信息系统：

土地信息系统是综合应用地理信息系统和管理信息系统，对人类在土地利用过程中产生的土地数据进展采集、存储、检索、分析和管理的信息系统。

2.实验目的

掌握ARCGIS10空间分析方法，熟练掌握ModelBuilder的创立和使用，掌握利用ARCGIS进展水文分析，缓冲分析，空间叠置分析，邻域分析、重分类等方法与技术。

掌握利用ARCGIS进展DEM及TIN的创立，掌握地图的编制、整饰及输出。

3.数据来源

汤地理信息系统空间分析实验教程，科技大学测绘学院地理信息系统实习教案。

4.要求

〔1〕局部实验需要写出步骤，步骤要图文并茂的反映操作流程，截图并用文字说明。

〔2〕所有附图均需添加、学号、班级。

〔3〕实验报告A4纸正反面打印，左侧装订。

〔4〕重要概念查阅相关文献，电脑输入，不用手工填写。

〔5〕附图黑白打印即可。

〔6〕实习报告后附实习心得体会，1000左右，手写。

〔7〕此局部容占总成绩25%。

5.实验容

实验一：

用ModelBuilder进展数据处理与分析，实例：

汤P130。

实验二：

汤地理信息系统空间分析实验教程〔P290寻找最正确路径〕

实验三：

土壤稳定性评估

6.附图

实验一：

数据处理——白水县

一、流程图：

图1白水县数据更新变换流程图

二、步骤：

〔1〕数据加载与环境设置：

1〕翻开Arcmap，添加图层vector，如下列图:

图2

2〕翻开Arctoolbo*，右击，点击“环境〞，设置工作空间如图3：

图3

〔2〕利用模型构建器构建模型：

1〕模型环境设置：

点击

翻开模型构建器，执行菜单命令模型属性，设置如图4。

选择环境→处理围→围→与图层vector一样，如图5所示：

图4

图5

2〕白水县行政围提取：

将Arcmap目录下的vector拖至ModelBuilder中，在Arctoolbo*中，选择分析工具→提取→筛选工具，将其拖至ModelBuilder窗口中，双击筛选框，在出现的对话框中选择输入栅格为vector，输出栅格命名为Vector_Select1.shp，如图6所示，确定，结果如图7：

图6

图7

3〕白水县DEM数据的拼接：

在ModelBuilder中点击

添加数据DEM1和DEM2，选择数据管理工具→栅格→栅格数据集→镶嵌至新栅格工具，将其拖至ModelBuilder窗口，双击镶嵌至新栅格，在出现的对话框中设置如图8，确定，结果如图9所示：

图8图9

4〕利用白水县围对DEM裁切：

在Arctoolbo*中选择spacialanalyst→提取分析→按掩膜提取工具，将其拖至ModelBuilder中，如图10所示。

双击按掩膜提取工具框，在出现的对话框中将DEM作为输入栅格，输入栅格数据或要素数据选为vector-select.shp，输出栅格为E*tract-DEM1,如图11，裁切结果如图12：

图10图11

图12

5〕白水县DEM的投影变换：

选择数据管理工具→投影与变化→栅格→投影栅格工具，将其拖至ModelBuilder窗口中，双击翻开投影栅格对话框，将E*tract-DEM1设为输入栅格，输出栅格命名为E*t-DEM1-Prj，点击输出栅格右边的

按钮，进入空间参考属性对话框，点击选择按钮，浏览坐标系，选择ProjectedCoordinateSystems→GaussKruger→*ian1980→*ian1980GKZone19.prj，重采样技术选择NEAREST,如图13所示，确定。

图13

6〕运行模型:

点击运行按钮运行模型，如图14所示：

图14

三、实验结果

〔1〕最终模型：

图15白水县数据更新变换模型图

〔2〕专题地图：

图16白水县DEM图

实验二：

寻找最正确路径

一、流程图：

图1寻找最正确路径流程图

二、实验步骤：

〔1〕激活工具：

翻开Arctoolbo*，激活spacialanalyst空间分析和3D分析扩展模块〔如图2〕：

自定义→扩展模块

图2

（2）设置环境，加载数据：

1〕翻开Arcmap，添加river.shp、startpot、endpot和DEM数据，如图3；翻开Arctoolbo*，右击选择环境，设置工作空间如图4：

图3图4

2〕点击

按钮翻开模型构建器，执行菜单模型→模型属性→常规，设置模型名称及标签，如图5左所示，确定。

菜单模型→图属性→符号系统→样式2，如图5右，确定。

将数据river.shp、startpot、endpot和DEM数据拖至模型构建器窗口中，如图6：

图5

图6

（3）利用modelbuilder构建模型：

1〕生成坡度本钱数据集：

①选择SpacialAnalyst工具→外表分析→坡度，将坡度工具拖拽至模型构建器窗口中，双击坡度工具框，在出现的对话框中，输入dem，输出栅格命名为slope，如图7左所示，确定。

运行如图7右：

图6

②选择SpacialAnalyst工具→重分类→重分类，将其拖至模型构建器，双击重分类对话框，输入slope，字段为Value，点击分类按钮，分类方法为相等间隔，类别为10，如图8所示，确定。

图8

2〕生成起伏度本钱数据集：

①选择SpacialAnalyst工具→邻域分析→焦点统计，将其拖至模型构建器中，双击焦点统计工具框，输入dem输出栅格命名为QDF，具体设置如图9所示，确定。

右击重分类工具框，选择运行。

图9

②选择SpacialAnalyst工具→重分类→重分类，将其拖至模型构建器中，双击重分类对话框，输入QDF数据层，点击分类按钮，选择相等间隔分类方法，类别选为10，如图10，确定。

分类结果如图11，输出栅格为Reclass-QDF，确定。

图10

图11

3〕生成河流本钱数据集：

选择SpacialAnalyst工具→重分类→重分类，将其拖至模型构建器中，双击重分类工具框，选择river数据层，其他设置如图12，确定。

图12

4〕加权合并单因素本钱数据，生成最终本钱数据集：

选择SpacialAnalyst工具→地图代数→栅格计算器,并将其拖至模型构建器中，双击栅格计算器，计算公式为cost=Recriver〔重分类流域数据〕+〔Reclassslop〔重分类坡度数据〕*0.6+ReclassQFD〔重分类起伏度数据〕*0.4〕，如图13所示，确定。

图13

5〕计算本钱权重距离函数：

选择SpacialAnalyst工具→距离分析→本钱距离，将其拖至模型构建器中，双击本钱距离，输入栅格为startpot，输入本钱栅格为cost，输出距离栅格为CostDis-star1，输出回溯连接栅格为blanklink，如图14所示：

图14

6〕求取最短路径：

选择SpacialAnalyst工具→距离分析→本钱距离，将其拖至模型构建器中，双击本钱路径工具框，在出现的对话框中输入栅格为endpot，目标字段为Id，输入本钱距离栅格为CostDis-star1，输入回溯连接栅格为blanklink，输出栅格为CostPat-endP1，如图15所示，确定。

图15

7〕运行模型：

点击运行按钮运行模型如图16：

图16

〔4〕加载模型运行结果图：

在Arcmap中加载模型运行后的结果，如图18所示：

图17

三、实验结果

〔1〕模型：

图18最正确路径模型图

〔2〕专题图：

图19最正确路径专题图

实验三：

土壤稳定性评估

一、流程图：

DEM

土地利用矢量图

坡向图

坡度图

土地利用栅格图

重分类

重分类

重分类

阴、阳坡等级图阴阳坡等级图

级土地利用等级图

坡度等级图

图1土壤稳定性评估流程图

二、实验步骤：

〔1〕激活工具：

翻开Arctoolbo*，激活spacialanalyst空间分析和3D分析扩展模块〔如图2〕：

自定义→扩展模块

图2

〔2〕设置环境，加载数据：

1）翻开Arcmap，添加DEM和landuse.shp数据，如图3。

翻开Arctoolbo*，右击选择环境，设置工作空间如图4所示：

图3图4

2）点击

按钮翻开模型构建器，执行菜单模型→模型属性→常规，设置模型名称及标签，确定。

菜单模型→图属性→符号系统→样式2，确定，如图5所示；将Arcmap中的DEM及laduse.shp数据拖至模型构建器窗口中，如图6：

图5

图6

（3）构建模型：

1〕提取坡度：

翻开spacialanalyst工具→外表分析→坡度，将坡度工具拖拽至模型构建器窗口中，双击坡度工具框，在出现的对话框中将dem设为输入栅格，输出栅格为slope，如图7，确定。

右击坡度工具框，选择运行，如图8所示：

图7图8

2〕坡度重分类：

选择spacialanalyst工具→重分类→重分类工具并将其拖拽至模型构建器中，双击重分类工具框，在出现的对话框中选择slope作为输入栅格，值字段选择value，如图9，点击分类按钮，类别选为6,方法选为手动，修改中断值如图10，确定。

图9图10

翻开将重分类对话框中的新值一栏依次改为10、8、7、5、3、1，输出栅格命名为Re-slope，如图11所示，确定。

模型如图12。

图11

图12

3〕提取坡向：

翻开spaciaalanlyst工具→外表分析→坡向，将坡向工具拖至模型构建器窗口，双击坡向工具框，在工具对话框中输入l，dem，输出栅格命名为Aspect，如图13，确定。

右击坡向工具框，选择运行，如图14所示：

图13

图14

4〕坡向重分类：

选择spacialanalyst工具→重分类→重分类工具并将其拖拽至模型构建器中，双击重分类工具框，在出现的对话框中选择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，确定。

右击面转栅格，运行，结果如图19：

图18图19

6〕对landusegrid重分类：

选择spacialanalyst工具→重分类→重分类工具并将其拖拽至模型构建器中，双击重分类工具框，在出现的对话框中选择landusegrid作为输入栅格，值字段选为LANDUSE,新值分别改为2、8、6、10，输出栅格命名为Re_land，如图20，确定。

模型如图21：

图20图21

7）空间叠加分析：

选择spacialanalyst工具→叠加分析→加权总和，并将其拖至模型构建器窗口中，双击翻开加权总和对话框，逐个添加Re-slope、Re-Aspect和Re_land，字段均选为VALUE,Re-slope的权重为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

图26

三、实验结果

图27土壤稳定性专题图

实验四：

土壤侵蚀性分析建模

一、流程：

图1土壤侵蚀性分析建模流程图

二、实验步骤：

〔1〕加载数据：

翻开Arcmap并添加数据studyarea.shp、vegetation.shp、soilsgrid.shp,如图2所示：

图2

（2）激活工具：

翻开Arctoolbo*，激活spacialanalyst空间分析和3D分析扩展模块〔如图3〕：

自定义→扩展模块；

图3

（3）数据前期处理：

1〕右击vegetation翻开属性properties→symbology→Uniquevalue，值字段选择vegtype,选择好色带，点击AddAllValues，如图4所示，点击ok；

图4

2〕调整图层顺序为studyarea，vegtation，soilsgrid如图5，保存地图文档；

图5

（4）利用ModelBuilder构建模型：

1〕设置属性与环境：

①点击

翻开模型构建器ModelBuilder如图6所示，执行菜单命令model→modelpropertis设置如图7，选择environment，设置如图8，点击value设置处理围为SameAsLayerStudyArea,如图9，确定；

图6图7

图8图9

②在modelbuiler窗口中执行菜单Model→DiagramPropertis,选择style2，如图10；

图10

2〕加载数据：

将ArcMap中的vegtation，soilsgrid拖至ModelBuilder窗口中。

3〕将DEM数据转为栅格数据：

从Arctoolbo*中将Conversiontools→toraster下的DEMtoRaser拖至ModelBuilder窗口中，如图11所示。

在modelbuilder中双击DEMtoraser,在出现的工具设置对话框中指定输入USGSDEM文件为：

elevation.dem,选择输出路径，确定，如图12：

图11

图12

4〕生成坡度数据：

将Arctoolbo*中3DAnalystTools→RasterSurface→slope拖至modelbuilder窗口中在ModelBuilder窗口中，点击添加连接按钮

将派生数据图框DEMToRa_elev1与工具图框slope连接在一起。

完成后效果如下列图13所示：

图13

※右键点击图框[outputraster]将其改名为:

[坡度图]

5〕坡度数据重分类：

将Arctoolbo*中SpatialAnalystTools→Reclass下的工具Reclassify拖放到ModelBuilder窗口中，点击

按钮将派生数据图框[坡度图]与工具图框Reclassify连接在一起，如图14所示：

图14

在ModelBuilder窗口中，双击工具图框Reclassfy，对坡度图重分类，输入栅格为坡度图，字段为value，点击load装入数据表，将坡度分为10级，具体设置如图15所示：

图15

8〕植被类型矢量数据转栅格：

将Arctoolbo*的ConversionTools→TORaster→featuretoraster拖放到ModelBuilder窗口中，在ModelBuilder窗口中，点击

按钮将数据图框vegetaion与工具图框featuretoraster连接在一起，如图16所示：

图16

双击与Vegetaion相连的工具图框FeaturetoRatser，在出现的对话框中，设置字段为VEGTYPE,选择输入栅格Vegetaion和输出路径，确定，如图16；右击FeaturetoRaster，执行命令Run，如图17；

图16图17

9〕加权叠加分析：

<1>添加工具：

将Arctoolbo*中SpatialAnalystTools→Overlay→WeightedOverly拖至modelbuilder中，如图18；

图18

<2>添加数据：

双击WeightedOverly，在出现的设置框中点击

按钮，在出现的框中输入soilsgrid，字段选为S-value，如图19，确定；同理，添加植被栅格〔字段值为vegtype，如图20〕和Reclass-slop1〔字段值为value，如图21〕，添加结果如图22所示：

图19图20

图21图22

<3>设置因子〔如图23所示〕：

①Soilsgrid：

根据不同土壤类型对土壤侵蚀危险性的影响力，不同的土壤类型给定不同的的数值，数值1表示改天土壤侵蚀危险度较低，9表示较高；

②植被栅格：

根据不同植被类型对土壤侵蚀危险性的影响力，不同的植被类型给定不同的的数值，数值1表示改天土壤侵蚀危险度较低，9表示较高〔water设为restrieted〕；

③Reclass-slop1：

根据不同坡度低的区域发生土壤侵蚀的危险系统较小，坡度较大的区域发生土壤侵蚀的危险系数较大，不同坡度对应不同数值，数值1表示改天土壤侵蚀危险度较低，9表示较高（第10级也设置为9）；

图23

<4>权重配置：

将上述三个因子的权重分别设置为soilsgrid25%，植被栅格25%，Reclass-slop1为50%，如图24所示：

图24

〔10〕最终模型如图25：

图25

三、实验结果：

*****

图26土壤侵蚀危险性专题图

实验五：

水文分析

一、流程图

流域分析

图1水文分析流程图

二、实验结果

******

图2水流方向专题图

图3流水累积量专题图

图4流域专题图

实验六：

找出某种珍贵药材的生长区域

一、

采样点

等高线数据

流程：

数据降雨量分布

温度分布

生成DEM数据

550=<降雨量<=680

10=<温度<=12

洼地填充

坡向数据

水流方向

提取阳坡

水流累积量

沟谷

沟谷缓冲区

数据药材分布区域图〔叠加分析、专题制图、求算面积〕

图1寻找某种珍贵药材的生长区域流程图

二、模型

图2寻找某种珍贵药材的生长区域模型图

三、专题图

图3某种珍贵药材的生长区域专题图

实验七：

地形鞍部的提取

一、

流程图

图1地形鞍部的提取流程图

二、实验结果

图2地形鞍部点专题图

实验八：

沟谷网络的提取

一、

栅格河网

程无洼地流向

流向

分析

沟谷网络

矢量河网

流量分析

洼地计算

DEM数据

实验流程

剔除伪沟谷

洼地填充

图1沟谷网络的提取流程图

二、实验结果

图2沟谷网络专题图

实验九：

TIN及DEM的生成与应用

实验结果：

图1TIN图

图2DEM图

3.TIN的显示：

图3TIN的显示

〔4〕TIN转矢量：

图4TIN转为矢量

〔5〕DEM的应用——坡度栅格图

图5坡度栅格图

〔6〕DEM的应用——坡度剖面曲率

图6坡度剖面曲率图

〔7〕DEM的应用——坡向栅格图

图7坡向栅格图

〔8〕DEM的应用——坡向剖面曲率图

图8坡向剖面曲率图

〔9〕DEM的应用——提取等高线

图9等高线图

〔10〕DEM的应用——地表阴影图

图10地表阴影图

〔11〕DEM的应用——通视性分析

图11通视性分析图

〔12〕DEM的应用——可视区分析

图12可视性分析图

〔13〕DEM的应用——地形剖面分析

图13地形剖面分析图

实验十：

缓冲区分析的应用

〔1〕点要素图层的缓冲区分析

图1点要素图层的缓冲区分析图

〔2〕线要素缓冲区的分析

图2线要素图层的缓冲区分析图

〔3〕面要素缓冲区的分析

图3面要素图层的缓冲区分析图

〔4〕水源污染防治

图4水源污染防治重点区域图

〔5〕受污染地区的分等定级

图5距污染源小于等于0.1的区域图

图6距污染源大于等于0.1小于等于0.15的区域图

图7污染区的分级图

〔距污染源小于等于0.1为1级，小于等于0.15为2级，大于0.15为3级〕

〔6〕城市化的影响围：

图8城镇化影响围专题图

7.实习心得