GIS矢量数据分析及栅格数据分析实验

1、YU NNAZ NORMAL UNIVERSITY本科学生实验报告姓名 尹永义 学号1443206000113专业地理科学班级 2014B _实验课程名称 地理信息系统 概论（实验）实验名称矢量数据分析与栅格数据分析指导教师及职称开课学期2014 至一 2015 学年学期云南师范大学 旅游与地理科学学院编印一、实验准备实验名称：实验（矢量数据分析与栅格数据分析）实验时间：2015/05/14实验类型：验证实验 综合实验1、实验目的与要求：1）了解使用矢量数据分析，栅格数据分析的基本工具。包括建立Buffer、Overlay，Select与 Clip 等。2）用Calculate Geometr

2、y计算面积与周长。以及多组多边形的地图叠置操作， 点线距离测量的两种方法，与空间自相关等。3）同时还有局域运 算，用Combine进行的局域运算，邻域运算，分区运算，数 据查询中的自然距离量测运算，栅格数据的提取操作等。2、实验材料及相关设备：1）安装有 Arcmap Arccatalog 软件的 Windows7 电脑一台，2）数据来源：地 理信息系统导论配套光盘中的 GIS文件中的datesets_v7 中的第十一章与十二章3）数据 Landuse, Soil,与 Sewers的 Shp文件，4）记录有博伊西国 家森林19086一 1996年的森林火情况的Boise_Fire文件， 198

3、6年的林火记录 Fire1986,Fire1992等，Deer.Shp,Edge.Shp,Adabg00.Shp5）一个象元大小30m的高程栅格Emidalat,6）一个包含四个坡度等级的坡度栅格Slope_Gd,7）一个有平地与四个主要方向的坡向栅格 Aspect_Gd,8）一个表示爱达荷州年平均降水量的栅格，Precipgd一个流域栅格，9）一个表示河流的栅将Strmgd,一个高度分带的栅格Elevgd3、实验理论依据或知识背景：矢量数据分析 矢量数据以点、线与面空间要素为输入数据。分析结果的准确性取决于空间特征的位置及形状的准确性。拓扑关系是一些矢量数据分析(如建立缓冲区与叠置分析)的一

4、个因素。基于邻近(Proximity)概念，建立缓冲区可把地图分为两个区域：一个 区域位于所选地图要素的指定距离之内，另一个区域在指定距离之外。 在指定距离之内的区域 称为缓冲区。围绕点建立缓冲 区产生圆形缓冲区。围绕线建立缓冲区形成一系列围绕每条线段的长条形缓冲带。围绕多边形建立缓冲区则生成由该多边形边界向外延伸的缓冲区。对线要素建立缓冲区未必在线两侧都有缓冲区，可以只在线的左侧或右侧建立缓冲区。缓冲距离(又叫缓冲大小)未必为常数，可以根据给定字段取值而变化。缓冲区边界也可以被融合掉，使得缓冲区之间没有叠置区。地图叠置操作是将两个要素图层的几何形状与属性组合在一起，生成新 的输出图层。输出图

5、层的几何形状代表来自各输入图层的要素的几何交集。输出图层的每个要素包含所有输入图层的属性组合，而这种组合不同于 其邻域。所有叠置方法都是基于布尔连接符的运算，即 AND、OR与XOR。 若使用AND连接符，则此叠置操作为求交(In tersect)。若使用OR连接符，则此叠置操作称为联合(Union)。若使用 XOR 连接符，则此叠置操作称为对称差异(Symmetrical Differenee)或差异(Differenee )。若使用以下表达 式(Input Layer)AND ( Identity Layer) OR( Input Layer)，贝该叠置操作称为识别(Ide ntity)或

6、减去(Minus )。模式分析是关于二维空间点要素空间分配的研究。在整体水平上，模式分析可以揭示某分布模式是随机、离散还是集聚的。在局部水平上，模式分析可以检测出分布模式中是否含有高值或低值的局部集聚。模式分析包括点模式分析、量测空间自相关的莫兰指数(Moran s I)与 量测高/低聚集度的G统计量。栅格数据分析栅格数据分析是基于栅格像元与栅格的。栅格数据分析能在独立像元、像元组或整个栅格全部像元的不同层次上进行。一些栅格数据运算使用单一栅格，而另一些则使用两个或更多栅格数据。栅格数据分析也应考虑像元数值类型(数字型数值，类别型数值)。 栅格数据分析环境包括分析的区域范围与输出像元大小。局域

7、运算局域运算是一个像元接一个像元运算，建立栅格数据分析的核心。 局域运算由单个或多个输入栅格生成一个新的栅格。格局域运算:单一栅格假定以单一栅格为源数据，基于输入栅格的像元值，局域运算通过空间数学函 数计算输出栅格的每个像元值。由于可以用多个 栅格图层进行运算，所以局域运算相当于基于矢量的地 图叠置操作。除了可用于独立 栅格的数学公式外，其他的基于输入栅格的像元值或其频率的度量也都可存储于输出栅格。邻域运算邻域运算，涉及一个焦点像元与一组环绕像元。环绕像元是按其相对于焦点像元的距离与（或）方向性关系来选定的。邻域运算得到的既可以是最小值、最大值、值域、总与、平均值、中值、 标 准差等统计值，也

8、可以是众数、少数与种类数等测量值列表。常见的邻域类型 有矩形、圆形、环形与楔形。分区运算分区运算用于处理相同值或相似要素的像元分组。这些组称为分区。分区可以是连续的 或不连续的。分区运算可对一个或两个栅格进行处理。若为单个输入栅格，分区运算量测每个分区的几何特征，如面积、周长、 厚 度（Thickness）与重心。给定两个栅格（一个输入栅格与一个分区栅格），要求以分区栅格的区 域 为范围对输入栅格进行分区运算生成输出栅格，输出栅格对分区栅格的每个分区概括了输入栅格的像元值。自然距离量测运算距离可以表达为自然距离与耗费距离。自然距离量测运算是计算与源像元的直线距离。配置与方向配置栅格中的像 元值

9、对应于距该像元最近的源像元。方向栅格中的像元值对应于距它最近的源像元的方向值。其他的栅格数据运算1. 栅格数据管理的操作包括剪取（Clip）与镶嵌（Mosaic）。2. 栅格数据提取是指从一个现有栅格提取数据生成一个新的栅格。提取栅格数据的工具可以是一个数据集、图形对象或查询表达式。3. 栅格数据的综合 归纳包括聚合（Aggregate）与 区组（Regiongroup）。基于矢量与基于栅格的数据分析的比较矢量数据分析与栅格数据分析是 GIS分析的两种基本类型。GIS软件包 不能在相同操作 中同时进行这两种分析，因此被分开处理。一般原则是，对于项目，选择有效的与适当的数据分析类型。、实验内容、

10、步骤与结果（要求：详细写清楚本次实验的完成的主要内容、具体实施步骤 与实验结果。纸张不够可以自行添加。）1）启示 Arccatalog,连接到 Chapll 数据库。启动 Arcmap,添加 Sewers.Shp, .Shp, Soil.Shp, 到图层中，将图层改名为Taskl,2）首先建立的缓冲区。单击并打开 Arctoolbox窗口。从快捷菜单中设置Environments，将数据库Chap11设置为当前工作空间，。在工具箱Analysistools/Proximity内双击Buffer工具。再出现的 Buffer对话框中，选择 Sewers为输入要素，.Shp作为输出要素集，输入300

11、米作为距离，选择All为 Dissolved Type，打开的 Sewerbuf 属性表。3）接着进行Soils丄anduse与Sewerbuf地图叠置，在工具箱An alysistools/Overlay 内双击 In tersect 工具。选择 Soils ,La nduse 与 Sewerbuf 为输入要素，输入Final.Shp，作为输出要素类，然后单击 Ok。4）将输出要素类命名为Sites.Shp,并单击用于输入表达式的SQL，输入表达式。5）打开Sites属性表，把Sites.Shp转换为Geodatabase要素类，6）字段类型，输入11为字段精度，3为字段尺度，单击Ok。以同

12、样的方法， 将Leng为新字段添加到中。7）在Sites属性表中，右击 Shape_Area并选择Caculate Geometry。在对话框，选中Area为性质，平方米为单位，单击Ok。8）右击Shape_Length并选择Caculate Geometry。在对话框，选中Length为性质，平方米为单位，单击Ok。2多组分多边形的地图叠置1）在Arcmap中插入新的数据帧，重命名为Tsk2，添加Regio ns，打开Boise_Fire属性表，2）首先，进行与的联合，在An alysistools/Overlay工具箱内双击 Un io n工具。选择Fire1986,与Fire1992为输

13、入要素，输入 Regoins作为要素 数据集的输出要 素 类，单击Ok3）在 An alysistools/Overlay 工具箱内双击 In tersect 工具。选择 1986,与 Fire1992 为输入要素，输入Fire_Intersect作为输出要素类，单击 Ok3点与线之间的距离量测1）在Arcmap中插入新的数据帧，重命名为Tsk3，添加Deer.Shp,.Shp,2）右击，指向Jion And Relates并选择Jion在Jion Data对话框中单击第一个下拉箭头，并 选择。指定Deer_Edge Shp为输出文件，单击 Ok3）右击Deer_Edge并打开属性表，4计算整

14、体与局部G统计量1）在Arcmap中插入新的数据帧，重命名为Tsk4，添加Adabg00.Shp,2）右击 Adabg00,并选择 Properties，在 Symbology 栏中，选择 Quantities/Graduated Colors来显示Latino字段值，3）打开 Arctoolbox， 首先计算整体 G 统计量，在 Spatial Statistics Tools/Analyzing Patterns工具箱内双击 High/Low Clustering工具。选择 Adabg00为输入要素， 选择 Latino 为输入字段，同 时 GeneralReport勾选前的复选 框，单击

15、 OK4）从 Geoprocessing菜单选择 Result,在 CurrentSession下，扩展 High/LowClustering,然后双击 HTML Report File5）双击 Spatial Statistics Tools/Maping Clusters 工具，选择 Latino 为输入字段， 输入Local_G.Shp作为输出要素类，6）在Arcmap中添加，Local_G,打开Local_G属性表，5 执行 Select 与 Clip1）在Arcmap中插入新的数据帧，重命名为Tsk5，添加JeffersonANDAmscmtype_PUB_24K_POINT ，2）

16、在Analysis Tools/Extract工具箱内双击 Select工具。在Select对话框中， 输入Amscmtype_PUB_24K_POINT 为输入要素，命名输出要素类为 action_completed，键入 SQL 表达式3）在 Analysis Tools/Extract 工具箱内双击 Clip 工具。选择 action_completed 为输入要素，Jefferson为剪取要素，ac_Jefferson为输出要素类，单击Ok栅格数据分析10）启动arccatalog,连接到Chap12数据库，在目录树中，从的快捷菜单中选择properties,raster datase

17、t properties 对话显示出 emidalat 属性有 186 列，214 行，像元大小 为 30m，11）启动 ArcMap添加Emidalat,至图层，图层重命名为 TASK1 &3 打开AtcToolbox右击Arctoolbox，选择Evironment,设定第十二章数据为当前与过期工作空间。双击SpatialAmalyst Tools/Math工具及下的TIMEs，在出现的对话框中，选择Emidalat，为12）输入栅格或常量值为1,;输入3。28为输入栅格或常量值为2,在当前工作 空间保存输出栅格为Emidaft。除了用ft度量单位外，Emidalat,与Emidaft是

18、一样的，单击ok。13）在 arcmap的 insert 中选择 Data Frame.并重命名为 Task2把 slope_gd 与aspect_gd添加到其中14）双击Spatial Amalyst Tools /local工具集中的combine工具，在出现的对话框中，选择slope_gdAndaspect_gdAnd作为输出栅格输入slp_asp为输出栅格，单击 ok，运行操作，Slp_asp显示出对输入值的每个独特组合独有一个独特输出冷打开其属性表查看独特组合及其计数，15）在 Spatial Amalyst Tools/Neighborhood 工具集中双击 Focal stati

19、stic，随后的对 话框中，选择mean为统计类型，单击Ok16）在 arcmap 的 insert 中选择 Data Frame 并重命名为 Task4 把 precipgd 与 hucgd 添加到其中17）双击 Spatial Amalyst Tools/zonal 工具集中的 zonal Statistic 工具，在出现的对 话框中，选择hucgd为输入栅格，选择precipgd为数值栅格，保持 huc_precip为输出栅格，选 择mean为统计类型，单击0K，18）为显示每个流域的平均降水量，可通过如下操作实现；右击huc_precip,选择 properties 在 Show 栏中

20、，单击 Unique Values,然后单击 OK,19）在 arcmap 的 insert 中选择 Data Frame .并重命名为 Task5把 strmgd 与 elevgd 添加到其中20）双击 Spatial Amalyst Tools/Distanee 工具集中的 Euclidean Distanee 工具，再出现的对话框中，选择 strmgd为输入栅格，将输出保存为 strmdist,然后单击OK, strmdist显示了在 strmgd中对河流的连续距离分区，21）下一步是生成一个距离河流 200m的区域的新栅格，双击Spatial Amalyst Tools/Reelass工具集中的Reclassify共聚焦，在其对话框中，选择strmgdist做 为输入栅格，并单击Classify按钮，在其后的对话框中，更改分类数为2,输 入200作为第一个断点值，单击 OK。22）双击Spatial Amalyst Tools /local工具集中的combine工具，再出现的对话框中，选择rec_gd与Elev_gd为输入栅格，将habitatl并保存为输出栅格，单击 OK 提取栖息地，双击 Spatial AmalystTools/Extractio