实验三、影像配准及矢量化
一、实验目的
1. 利用影像配准(Georeferencing) 工具进行影像数据的地理配准 2. 编辑器的使用(点要素、线要素、多边形要素的数字化)。 3. 熟悉GRAMIN GPS手持机的基本使用方法。
注意:在基于ArcMap 的操作过程中请注意保存地图文档。
二、实验准备
数据:昆明市西山区 普吉 地形图 1:10000 地形图――70011-1.Tif,昆明市旅游休闲图.jpg (扫描图),请从课程网站下载。
在开始本实验之前,将同学分成3-5个小组,每个小组从实验室借出一部GARMIN GPS手持机,利用课外时间在昆明市范围内采集20个GPS点,并在自行购买和昆明市旅游休闲图上做好标记。数据采集回来后用MapSource 软件通过数据线将数据从GPS手持机中导入到计算机并存储为GPS.txt 文件以备使用。
在借到GPS手持机后,请在老师指导下认真阅读GPS手持机的说明书,掌握基本操作方法。使用完毕后请归还实验室。外出采集数据时请遵守指导老师宣布的安全注意事项。 软件准备:
MapSource (GRAMIN GPS手持机数据导入软件)
ArcGIS Desktop ---ArcMap
三、实验内容及步骤
第1步 地形图的配准-加载数据和影像配准工具
所有图件扫描后都必须经过扫描配准,对扫描后的栅格图进行检查,以确保矢量化工作顺利进行。
打开ArcMap,添加“影像配准”工具栏。
把需要进行配准的影像—70011-1.TIF增加到ArcMap中,会发现“影像配准”工具栏中的工具被激活。
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
1
GIS软件应用-ArcGIS
第2步 输入控制点
在配准中我们需要知道一些特殊点的坐标。通过读图,我们可以得到一些控件点――公里网格的交点,我们可以从图中均匀的取几个点。一般在实际中,这些点应该能够均匀分布。
在”影像配准”工具栏上,点击“添加控制点”按钮。
使用该工具在扫描图上精确到找一个控制点点击,然后鼠标右击输入该点实际的坐标位置,如下图所示:
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
2
GIS软件应用-ArcGIS
用相同的方法,在影像上增加多个控制点(大于7个),输入它们的实际坐标。点击“影像配准”工具栏上的“查看链接表”按钮。
注意:在连接表对话框中点击“保存”按钮,可以将当前的控制点保存为磁盘上的文件,以备使用。
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
3
GIS软件应用-ArcGIS
检查控制点的残差和RMS,删除残差特别大的控制点并重新选取控制点。转换方式设定为“二次多项式”
第3步 设定数据框的属性
增加所有控制点,并检查均方差(RMS)后,在”影像配准”菜单下,点击“更新显示”。执行菜单命令“视图”-“数据框属性”,设定数据框属性
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
4
GIS软件应用-ArcGIS
在“常规”选项页中,将地图显示单位设置为“米”
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
5
GIS软件应用-ArcGIS
在“坐标系统”选项页中,设定数据框的坐标系统为
“Xian_1980_Degree_GK_CM_102E”(西安80投影坐标系,3度分带,东经102
度中央经线),与扫描地图的坐标系一致
更新后,就变成真实的坐标。
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
6
GIS软件应用-ArcGIS
第4步 矫正并重采样栅格生成新的栅格文件
在”影像配准”菜单下,点击“矫正”,对配准的影像根据设定的变换公式重新采样,另存为一个新的影像文件。
加载重新采样后得到的栅格文件,并将原始的栅格文件从数据框中删除。后面我们的数字化工作是对这个配准和重新采样后的影像进行操作的。
通过上面的操作我们的数据已经完成了配准工作,下面我们将使用这些配准后的影像进行分层矢量化。
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
7
GIS软件应用-ArcGIS
第5 步 分层矢量化-在ArcCatalog中创建一个线要素图层
该数据采用的是西安80坐标系统、3度分带
(1) 打开ArcCatalog.在指定目录下,鼠标右击,在“新建”中,选择“个人 Geodatabase”。并修改该Geodatabase数据库的名称(例如test3.mdb)。
(2) 下面将为该Geodatabase创建新的要素类,首先创建一个“等高线”要素类来存储等高线要素。在ArcCatalog中,鼠标右击test3这个个人 Geodatabase,在“新建”中选择“要素类”.
(3) 输入创建的要素类的名称“等高线”,点击下一步。 (4) 点击下一步。
下面将是我们创建新的要素类的关键,为我们的数据定义坐标系统,空间范围,存储要素类型。以及可以在这增加属性字段。
(5) 点击Shape字段。在对话框中将显示详细的选项,我们首先点击“几何类型”,并将要素类型选择为我们需要的类型(我们现在要创建等高线这个要素类,所以应该选择线).
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
8
GIS软件应用-ArcGIS
(6) 点击“空间参考”选项后面的按钮,在“空间参考属性”对话框中的“坐标系”选项页下,将选择合适的坐标系统,点击“选择”按钮。在(Projected Coordinate Systems图所示:
目
录
下
,
选
择
Gauss
Kruger---Xian
1980--
Xian_1980_Degree_GK_CM_102E.prj)。点击增加,现在这些坐标系统信息应该如下
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
9
GIS软件应用-ArcGIS
(7) 再点击“X/Y 域”选项页,在该选项页下为我们的数据定义存储的空间范围。该空间范围需要认真考虑,不仅要考虑你当前的纸制地图的空间范围,还要考虑到将来工作中还会出现的最大的空间范围。
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
10
GIS软件应用-ArcGIS
为确定这个区域XY(最小值),XY(最大值),可以切换到ArcMap中,点击“绘制”工具栏上的“矩形框”按钮,在地图显示区中画一个矩形,使区在更大范围内包含已配准的栅格地图。右键选中这个矩形框,设置“属性”,将填充色设置为“无”,可得到如下的效果:
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
11
GIS软件应用-ArcGIS
在矩形框属性的“大小和位置”选项页中,可获取矩形框左下角和右上角的坐标(X,Y),将这里获取的X,Y值分别填入到上面“空间参考属性”对话框的“X/Y域”选项页 “最小X”、“最小Y”、“最大X”、“最大Y”
输入框中
通过上面的操作我们为创建的要素类定义了正确的坐标系统和空间范围。 (8)
下面我们将为该数据创
建新的属性字段。“高程”,类型设置为“Float”用来存储等高线的高程值。
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
12
GIS软件应用-ArcGIS
(9)
创建了一个线状的要素类。
点击完成这样,我们就
第6步 从已配准的地图上提取等高线并保存到上面创建的要素类中
(1) 切换到ArcMap中, 将新建的线要素图层,加载到包含已配准地形图的数据框中,保存地图文档为Ex3.mxd
(2) 打开“编辑器”工具栏,在“编辑器”下拉菜单中执行“开始编辑命令”,并选择前面创建的“等高线”要素类。确认编辑器中:任务为――新建要素,目标为――等高线,设置图层――等高线的显示符号为红色,并设置为合适的宽度。
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
13
GIS软件应用-ArcGIS
(3) 将地图放大到合适的比例下,从中跟踪一条等高线并根据高程点判读其高程,输入该条等高线的高程。
(4) 进一步练习线要素的其它操作,比如线段的合并、分割、编辑顶点等操作 (5) 可参照以上步骤,从地图中提出多边形要素(比如居民地),并进一步熟悉多边形要素编辑的相关操作。
第7步 根据GPS观测点数据配准影像并矢量化的步骤
数据:扫描地图-昆明市旅游休闲地图(YNKM.JPG)、
Garmin 手持GPS野外采集数据(gpsdata.dbf)-GCS_WGS_1984 地理坐标系
(1). 打开ArcMap,添加扫描地图-YNKM.JPG,打开“影像配准”工具栏(在ArcMap的工具栏的空白区域点击鼠标右键,然后选择“影像配准”)
(2) 参考练习2中最后一步的内容,根据gpsdata.dbf中的内容,将其转换为一个新的图层:GPS.shp,并将其添加到当前数据框中。
添加gpsdata.dbf
执行菜单命令<工具>-<添加X,Y数据>, 在出现的菜单中指定坐标系统为地理坐标:GCS_WGS 1984, X坐标指定为经度(E),Y坐标指定为纬度(N)
在TOC 面板中“显示”视图下,右键选择图层“gpsdata事件”,从右键菜单中
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
14
GIS软件应用-ArcGIS
执行“数据”->“导出数据”,将其导出成为一个新的Shape文件-名称为gps.shp
将GPS.shp 添加到当前的数据框中
这里显示的数据就是第1个小组在野外用GPS获取的控制点,每个控制点在
纸质地图上都有一个唯一的记号。
(3) 在TOC中右键选择图层――YNKM.JPG,在出现的菜单中点击“缩放到图层”,并将其放大到某一尺度下
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
15
GIS软件应用-ArcGIS
(4) 对比第一组同学在进行GPS数据采集时所使用的纸质地图,在地图显示区中找到第1个控制点。
点击“影像配准”工具栏上的控制点选择工具,在扫描地图中,采集第1个GPS
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
16
GIS软件应用-ArcGIS
控制点的位置,点击。
(5)在TOC中右键选择图层――GPS,在出现的菜单中点击“缩放到图层”,并将其放大到某一尺度下,并移动地图,可以方便地找到第1个控制点
(6) 在地图显示区中,我们可以看到第1组同学在野外采集的GPS控制点,找到与扫描地图中对应的那个GPS控制点,点击鼠标。
(7)通过以上操作我们已经完成了第1个控制点的选择。以类似的方法,添加至少4个控制点。
(8)如果操作正确,在完成以上操作后,扫描地图就被配准到了GCS_WGS_1984地理坐标系下。当鼠标在地图显示区移动时,在ArcMap状态栏上就会显示当前位置在GCS_WGS_1984坐标系下的经纬度坐标。
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
17
GIS软件应用-ArcGIS
在添加所有控制点,并检查RMS后,可使用下拉菜单“地理参考”中的“矫正”命令对栅格地图重采样。转换方法使用“二次多项式”。加载重采样后的栅格地图,并将原始的扫描地图从数据模框中删除。
(9) 在ArcCatalog中新建一个要素类,可命名为:Schools,(要素类型为点要素-将用于从扫描图中提取各个学校的位置,添加所需要的字段-用于存储学校的名称,坐标系统设置为WGS_1984地理坐标系-与GPS野外数据相同)
(10) 将Schools图层添加到当前数据框中。
(11)在TOC面板中,右键选择“图层”,在出现的菜单中,点击“属性”,设定数据框的坐标系为GCS_WGS_1984
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
18
GIS软件应用-ArcGIS
(12) 打开“编辑器”工具栏,在“编辑器”下拉菜单中执行“开始编辑”命令,选择要编辑的图层-Schools,确认“编辑器”工具栏中,任务为――新建要素,目标为――Schools。
(13) 将地图放大到合适的比例下,在“编辑器”中选中“草图”工具
,在地图
中找到图例为○文的要素,在其所在位置,点击鼠标,分层提取昆明市的各高校位置数据-(点要素)。打开属性表,修改要素属性,比如:输入学校的名称。
(14)完成编辑后,停止编辑,并将修改结果保存到要素类-Schools中。
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
19
GIS软件应用-ArcGIS
四、实验报告及要求
做出书面报告,包括原理、过程和结果。具体内容如下:
1) 总结屏幕跟踪数字化过程的基本步骤及每一步骤的必要性。 2) 分析数字化过程中误差的来源及减小误差的相关方法。 3) 为什么要对配准后的数据进行重采样?
GSNIAC: http://www.geo-spatial.net/ kecheng.yang@gmail.com
20
因篇幅问题不能全部显示,请点此查看更多更全内容