实验四、空间数据处理一、实验目的1.掌握空间数据处理(融合、拼接、剪切、交叉、合并)的基本方法,原理。领会其用途。2.掌握地图投影变换的基本原理与方法。3.熟悉ArcGIS中投影的应用及投影变换的方法、技术4.了解地图投影及其变换在实际中的应用。二、实验准备预备知识:ArcToolbox是ArcGISDesktop中的一个软件模块。内嵌在ArcCatalog和ArcMap中。ArcToolbox具有许多复杂的空间处理功能,包括的工具有:数据管理数据转换Coverage的处理矢量分析地理编码统计分析空间间数据处理是基于已有数据派生新数据的一种方法。是通过空间分析方法来实现的。是基于矢量数据进行的,包括如下几种常用的操作:融合,剪切,拼接,合并(并集),相交(交集)。地理坐标系(GeogrpahicCoordinateSystem)地理坐标系使用基于经纬度坐标的坐标系统描述地球上某一点所处的位置。某一个地理坐标系是基于一个基准面来定义的。基准面是利用特定椭球体对特定地区地球表面的逼近,因此每个国家或地区均有各自的基准面。椭球体长半轴a(米)短半轴b(米)Krassovsky(北京54采用)63782456356863.0188IAG75(西安80采用)63781406356755.2882WGS8463781376356752.3142在ArcGIS中基于这三个椭球,建立了我国常用的三个基准面和地理坐标系:GCS_WGS1984(基于WGS84基准面)GCS_BEIJING1954(基于北京1954基准面)GCS_XIAN1980(基于西安1980基准面)投影坐标系(ProjectedCoordinateSystems)投影坐标系使用基于X,Y值的坐标系统来描述地球上某个点所处的位置。这个坐标系是从地球的近似椭球体投影得到的,它对应于某个地理坐标系。投影坐标系由以下参数确定地理坐标系(由基准面确定,比如:北京54、西安80、WGS84)投影方法(比如高斯-克吕格、Lambert投影、Mercator投影)在ArcGIS中提供了几十种常用的投影方法北京1954投影坐标系和西安1980坐标系都是应用高斯-克吕格投影,只是基准面、椭球、大地原点不同。地理变换地理变换是一种在地理坐标系(基准面)间转换数据的方法,当将矢量数据从一个坐标系统变换到另一个坐标系统下时,如果矢量数据的变换涉及基准面的改变时,需要通过地理变换来实现地理变换或基准面平移。主要的地理变换方法有:三参数和七参数法。投影变换当系统所使用的数据是来自不同地图投影的图幅时,需要将一种投影的地理数据转换成另一种投影的地理数据,这就需要进行地图投影变换。实验数据:云南县界.shp;Clip.shp西双版纳森林覆盖.shp西双版纳县界.shp三、实验内容及步骤空间数据处理步骤:将所需要的数据解压到磁盘:如e:\gisdata,设定工作区:在ArcMap中执行菜单命令:工具-选项,在“空间处理”选项页里,点击“环境变量”按钮,在环境变量对话框中的常规设置选项中,设定“临时工作空间”为e:\gisdata第1步裁剪要素在ArcMap中,添数据GISDATA\云南县界.shp,添加数据GISDATA\Clip.shp(Clip中有四个要素)激活Clip图层。选中Clip图层中的一个要素,注意确保不要选中“云南县界”中的要素!点击打开ArcToolbox,指定输出要素类路径及名称,这里请命名为“云南县界_Clip1”指定输入类:云南县界指定剪切要素:Clip(必须是多边形要素)依次选中Clip主题中其它三个要素,重复以上的操作步骤,完成操作后将得到共四个图层(“云南县界_Clip1”,“云南县界_Clip2”,“云南县界_Clip3”,“云南县界_Clip4”)。第2步拼接图层在ArcMap中新建地图文档,加载你在剪切要素操作中得到的四个图层点击打开ArcToolbox在ArcToolbox中执行“追加”命令输出要素:设定为云南县界_Clip1,输入要素:依次添加其它三个图层右键点击图层“云南县界_Clip1”,在出现的右键菜单中执行“数据”-”导出数据”指定导入数据的路径和名称:YNOK.shp通过以上操作我们就完成了将4个图层拼接为一个图层的处理。新建一地图文档,加载数据YNOK.shp,查看图层及打开其属性表看看与“云南县界”中的属性表有何区别。第3步要素融合在拼接图层的基础上继续执行“融合”命令输入要素:指定为YNOK融合字段:选择为“所属州”,将根据这个字段的值对要素进行融合,YNOK图层中“所属州”相同的要素将合并成一个要素以上操作,根据指定字段的值,对现有图层中的要素进行融合,产生新的图层――YNOK_Dissovle,打开并查看其属性表类似地,重复以上过程,并将融合字段指定为:CHINESE,看看结果有何不同第4步图层合并在ArcMap中新建一个地图文档,加载数据GISDATA\西双版纳森林覆盖.shp和GISDATA\西双版纳县界.shp调整图层顺序,将西双版纳县界置于下方打开ArcToolbox,在ArcToolbox执行“联合”命令在联合对话框中输入要素:依次添加“西双版纳森林覆盖”“西双版纳县界”两个图层输出要素类:设置为Union.shp查看输出要素类:Union的的属性表,并检查属性“Type”,其中为“Y”的表示有植被覆盖的区域,右键点击图层Union,修改属性-符号(设置为唯一值图例,字段设置为TYPE)思考题:勐海县的总面积是多少平方公里?其中有森林覆盖的区域面积是多少?没有森林覆盖的区域面积是多少?第5步图层相交在图层合并练习的基础继续在ArcToolbox中,执行“相交”命令在“相交对话框”中输入要素:依次添加“西双版纳森林覆盖”“西双版纳县界”两个图层输出要素类:设置为Intersect.shp查看输出要素类InterSect,并与“西双版纳森林覆盖”及“图层合并”操作所得结果――“Union”进行比较,并进一步思考这类操作适合求解哪一些现实问题。定义地图投影第6步定义投影(1)在ArcMap中新建地图文档,添加第4步成生成的图层:Union.shp(2)在TOC中,右键点击图层“Union”,查看属性,在属性对话框中,点击“源”选项页,查看这图层是什么坐标系(3)打开ArcToolbox,执行命令“定义投影”命令在定义投影对话框中,选择要素类:Union,点击坐标系输入框右边的按钮,在出现的“空间参考属性”对话框中,选择一个地理坐标系,GCS_BEJING_1954(注意:前提是我们已知道图层Union是使用北京1954地理坐标系)点击“选择按钮”,从预定义的坐标系中选择(坐标系统\GeographicCoordinateSystems\Asia\Beijing1954.prj)(4)在TOC中,右键点击图层“Union”,查看属性,在属性对话框中,点击“源”选项页,查看这个图层的坐标系是否已经被指定为“北京1954地理坐标系”第7步投影变换――地理坐标系-北京1954坐标系转换-西安80坐标系(1)在第6步的基础上进行(2)打开ArcToolbox,执行命令“定义投影”命令(3)在“投影”对话框中,依次设定输入要素类为――Union,输出要素类为――Union_PRJ_BJ54.shp,输出坐标系选择为――“BEJING_1954_GK_ZONE_17N”从预定义的坐标系中选择(坐标系统\ProjectedCoordinateSystems\GaussKruger\Beijing1954\Beijing1954GKZone17N.prj)(4)确定后,完成由地理坐标系GCS_BEJING_1954到投影坐标系BEJING_1954_GK_ZONE_17N的变换。(5)请参照以上过程,完成由投影坐标系-BEJING_1954_GK_ZONE_17N到投影坐标系四、实验报告要求回答第4步中的问题―――没有森林覆盖的区域面积是多少?说明原理并给出结果。提示:面积单位是什么?如何转换成平方公里?(1)基于第4步的结果,用图层相减功能求出没有植被的区域,得到一个图层,检查此图层的坐标系是否已是投影坐标系(地图单位是米?),如果不是则将其转换为投影坐标系(2)查找ArcMap中的联机帮助,以了解如何计算多边形的面积提示:用字段计算器