ArcGIS操作步骤
1.问题提出
纸质地图是GIS的重要数据源之一,在目前大多数应用的数据获取阶段都是采用先将纸质地图扫描,然后再进行屏幕矢量化得到矢量数据,对数据进行编辑,如图形编辑,建立拓扑,属性编辑等,以备后续的分析应用。在本节中将通过一个实例来学习在ArcGIS中如何进行屏幕矢量化及拓扑建库。
2.数据准备
现有一张扫描的地图,如图3.4 1所示,每个多边形表示不同的地块,为了对这些地块进行一些复杂的分析,需要把这幅扫描地图转化为矢量形式,并建立多边形间的拓扑关系。扫描地图文件名为cropmap.jpg。
图像配准
图像配准实际上是对扫描的地图建立空间坐标系,这在ArcGIS中叫做Georeference(空间参考)。当我们在ArcMap中添加cropmap.jpg图层时,系统会提示该图层没有空间参考信息,需要通过图像配准为该图层添加空间参考信息。
1.连接数据
2.打开数据
启动ArcMap,添加cropmap.jpg数据
3.添加控制点
选取图像上线划的交点为控制点,如图3.4 3所示,上排的点从左到右依次被编号为:21,22,23;下排的点从左到右依次被编号为:31,32,33。由于ArcMap是动态配准的,所以首先必须根据左上的21和右下的33两个控制点确定图象的范围。
4.为控制点赋值
这六个控制点的坐标是利用GPS测量得到的,采用的是Geocentric Datum of Australia1994 (GDA94)大地坐标系,X为经度,Y为纬度,具体坐标值见表3.4 1,DMS为以度分秒表示,DD为以十进制度数表示。首先添加控制点21和33以确定地图的范围,然后添加控制点23,31,22和32,在添加并数字化确定地图范围的控制点21和33后,其他几个点的添加顺序就无所谓了。由于地图是动态的配准,即随着控制点的添加逐步纠正,地图会进行诸如旋转等的校正。在ArcMap的主菜单中添加Georeferencing(地理配准)工具条。
5.输入控制点坐标
在Georeferencing工具条上点击(Add Control Points) 图标,鼠标会变成十字丝形状。先左键点击21点,然后右键点击21,在弹出菜单中选择“Input X and Y”,在弹出的“Enter Coordinates”对话框中参照表3.4 1中的DD形式的值输入控制点21的空间坐标,如图3.4 4所示,点击“OK”按钮结束控制点坐标输入。
提示:如果无法精确点击21点的话,可以先将原图放大。
6.重复步骤5
依次添加控制点33,22,23,31,32的坐标。
注意:如果在控制点添加过程中出现错误,只能在Georeferencing工具栏的下拉菜单中选择Clear All Control Point删除已添加坐标值的控制点,全部重新添加。
思考:为什么要先添加控制点21和33的坐标值。
7.图像配准
当6个控制点的坐标都添加完毕后,点击Georeferencing/Recftfy进行图像配准的操作,在弹出的Save As对话框中采用缺省设置,将配准好的图象以rectifcropmap.tif文件名存放好。
思考:当前的栅格图像的地理参考是什么?
8.新建个人地理数据库
在ArcCatalog新建个人地理数据库命名为CropDB。
9.新建要素集
新建要素集以存储数字化地图,将新要素集命名为CropsGDA。
10.选择坐标系统
GPS采集的控制点是GDA94坐标系,因此将CropsGDA要素集的坐标设定为GDA94坐标系。选择Geographic Coordinate Systems/Australia and New Zealand/Geocentric Datum of Australia 1994,如图3.4 7所示。点击“下一步”,后面的选项都采用缺省设置,最后点击Finish。
11.在CropsGDA要素集里新建CropsLine线要素类
12.数字化区域边界
13.多边形追踪建立拓扑
利用ArcToolBox中的“数据管理工具箱”中的“要素工具集”中的“要素转面”工具将CropLine线要素转为名为CropPoly的面要素。转换结果如图3.4-14所示。
14.添加属性
在ArcCatalog打开为CropPoly要素的属性表添加一个名为TypeID的字段,字段类型为短整型(Short Integer)。再对照cropmap.jpg文件中每个地块的编号标准为每一个地块的TypeID字段赋值。
15.新建属性表
在ArcCatalog中为地理数据库CropDB新建一个表,命名为mgmt,表中内容如表3.6-1所示:
15.为属性表mgmt添加内容
在ArcMap中添加mgmt,将该图层设置为编辑状态,添加记录。
16.连接属性
在ArcMap中将右键点击CropPoly图层,选择“连接和关联”-“连接”,将表mgmt和面要素图层CropPoly连接起来,使得该图层的属性表中能够显示农作物类型名称。
【近期GIS制图三天两晚课程预告】
2021年10月28日-31日 济 南