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日 济 南
联系客服
