1、 本科学生实验报告姓名 尹永义 学号 1443206000113 专业 地理科学 班级 2014B 实验课程名称 地理信息系统概论(实验) 实验名称 几何变换 指导教师及职称 速绍华 (讲师 ) 开课学期 2014 至 2015 学年下 学期云南师范大学旅游与地理科学学院编印一、实验准备实验名称:实验(几何变换 ) 实验时间:2015/04/16实验类型: 验证实验 综合实验 1、 实验目的和要求:1) 了解使用几何变换,扫描文件的仿射变换, 把扫描图像变换到 UTM 坐标系。2) 几何变换后的扫描文件的矢量化过程,3) 卫星图像的仿射变换 2、 实验材料及相关设备:1) 安装有 Arcmap
2、 Arccatalog 软件的 Windows7 电脑一台,2) 数据来源:地理信息系统导论配套光盘中的 GIS 文件中的 datesets_v7 中的第六章3) Hoytmtn.tif,一个包含扫描文件的土壤界限的 TIFF 文件4) Rect_ hoytmtn.tif,习作一矫正后的 TIFF 文件5) Spot-panbil 一景分辨率为十米的 SPOt 全色卫星图像;road.shp,一个带有有 GPS 接收机测得的,并被投影到 UTM 坐标系的 shapefile 文件3、实验理论依据或知识背景:几何变换 几何变换就是利用一系列控制点和转换方程式在投影坐标上配准数字化地图、卫星图像或
3、航空照片的过程。 GIS 中, 几何变换可分为地图到地图的变换和影像到地图的变换。 均方根误差(RMS)是度量几何变换质量的一种定量方法。 变换方法 坐标系统之间进行几何变换有不同的方法。各种方法的区别在于它所能保留的几何特征,以及允许的变化。 仿射变换 仿射变换是假设输入(地图或影像)是均匀失真的,通常被建议用于地图到地图或影像到地图的变换。 如果输入的是已知的具有不均匀分布的失真,如地形位移的航空照片(由于局部地形造成的对象移位) ,建议使用投影变换。 控制点 控制点在确定仿射变换准确度中起关键作用(Bolstad 等 1990) 。控制点的选择因地图到地图的变换和影像到地图的变换而异。
4、地图到地图变换的控制点选择相对比较直接。我们只需要有已知真实世界坐标的点。 影像到地图变换时的控制点通常也叫做地面控制点。地面控制点(GCPs) 是影像坐标(用行和列显示)和真实世界坐标都能识别的点。GCPs 直接从卫星影像选取。因此,其地面控制点的选取并不像在数字化地图上选取 4 个地理控制点那么直接。 均方根 (RMS) 误差 控制点的好坏通常用均方根(RMS)误差来衡量,即对控制点实际位置(真实的)与估算位置(数字化的)之间偏差的估量。 数字化地图上的均方根误差 如果 RMS 误差在可接受的范围之内,通常我们会假设整幅地图的变换也是可以接受的。 如果在数字化控制点或输入控制点的经纬度值时
5、,出现大的差错,该假设就是错误的。 像元值重采样 三种常用的重采样方法是:邻近点插值法、双线性插值法和三次卷积插值法。 邻近点插值法是将原始图像的最邻近像元值填充新图像的每个像元。 双线性插值法把基于三次线性插值得到的 4 个最邻近像元值的平均值赋予新图像的相应像元。 三次卷积插值法则用五次多项式插值法求出 16 个相邻像元值的平均值。二、实验内容、步骤和结果(要求:详细写清楚本次实验的完成的主要内容、具体实施步骤和实验结果。纸张不够可以自行添加。 )1) 连接 ArcCatalog 至第六章数据库,启动 Arcmap,重命名为数据帧 Task1,将 hoytmtn.tif添加到 Task1.
6、单击 Customize ,菜单指向 toolbars,选中 Georeferencing,工具条出现在Arcmap,窗口,2) 放大 hoytmtn.tif,定位四个控制点,用括号圈起控制点,两个在图像的顶部,两个在图像的底部,按顺时针从一到四编号,其中,左上角的点为 1,3) 在第一个控制点周围放大视窗,激活 Georeferencing 工具条里的 Add Control Point 工具。单击于括号中心线相交的交叉点,再单击一次,在控制点上有一个由绿变红的+出现用同样的方法添加其他三个点。4) 在办步骤中,要更新四个控制点的坐标值,单击 Georeferencing 工具条里的 vie
7、w link table。链接表格点部出了四个控制点,Xsource Ysource,Xmap,Ymap 和残差值,其中前两个是扫描图像上的点,的坐标值,后两个是扫描图像上的坐标值。5) 如果控制点被正确的添加到图像,总的均方根误差应该小于 4.0m,如果均方根误差大,选中残差值大的记录并删除,单击 linktable 上的对话框 OK6) 本步骤要纠正便缓缓的土壤文件,选择 Georeferencing,下拉菜单的 Rectify,在出现的对话框中,选择所有参数为默认值,并将纠正后的图像保存在第六章数据库里,命名为rect_hoytmttn.tif7) 本步骤要创建一个新的 shapefil
8、e,用于存储 rect_hoytmttn.tif 的矢量化要素,右击第六章文件夹,指向 NEW,选择 shapefile,在 creat new shapefile 对话框里,输入文件名为hotmtn_trace.shp,要素类型为 polyline,单击 edit 按钮,新文件的坐标系为 Prejected Coordinate systerms /UTM/NAD1297zone11N,单击 OK 退出8) 插入数据帧,命名为 Task2,天假 rect_hoytmttn.tif 和 hotmtn_trace.shp,从前一个中选择Properties,在 Symbology 表里,选择 U
9、nique Values,选择 biuld the attribute table .改变符号颜色为红色,值为零,关闭对话框,右击 rect_hoytmttn.tif 指向 Zoomto Layer .因为 rect_hoytmttn.tif 里面的颜色很浅。9) 在 Arcmap 中,单击 Editor Toolbar,在 Editor 下拉菜单中选择 StartEditing,选择hoytmtn_trace 进行编辑。编辑模式激活了 Arcscan 工具条。10) 本步骤要矢量化参数,从 Vectorization,下拉菜单选择 Vectorization setting,单击Styles
10、,选中 Polygons,并在下一个对话框中单击 OK,单击 Apply,在单击 Apply,关闭Vectorization setting 对话框11) 在 Vectorization 菜单中选择 Generate Features 是添加中心线的图层是 hoytmtn_trace.单击OK 批量式量化的结果保存在 hoytmtn_trac 图层,12) Rect_hoytmtn 左下角有关于土壤的注释,在标准工具上,选择 Nottes,完成图想到地图的转换1) 在 Arcmap 中插入新的数据帧, ;重命名为 Task3,添加 Spot-panbil 和 road。Shp2) 使用快捷菜单
11、的 Zoom 观察 Spot-panbil 和 road。Shp3) 在下拉菜单中的 AutoAdjust 打钩,右击 road,选择 Zoom to layer,单击 Add Controlpoints,4) 在 Georeferencing 中单击 view link table, ,其显示四条记录,每天条记录对应三个链接,5) 在 Georeferencing 菜单中选择 Rectify。在出现的下一个对话框中,进行祥云大小和像云冲采样的选择,同时还要定义输出文件的名称。单击 Save 关闭对话框。6) 添加并查看 rect_spot,经过地理坐标匹配和纠正过的 rect_spot,是一
12、个栅格数据,7) rect_spot 的值域应该是 16100,但如果他的值域是 16-255,则应该使用以下方法纠正;第一,Spatial Analyst 下拉菜单选择 Reclassify。在 Reclassify 对话框中,单击Unique。行的原值为 255,改为 NO data,单击 OK,第二右击热 Reclass of rect _spot,选择prroperties,在 Show 框中,选择 Stretched。改变 Low Lable 为 16,Hign Lable 为 100,单击 OK,右击 Reclass of rect _spot,单击 Data,并选择 Export
13、 Data。然后指定输出名称和保存位置。三、实验小结1、实验中出现过的问题(或错误) 、原因分析(1) 、在做到习作二,矢量化栅格线条时,忘了打开 ARCScan 的扩展模块, ,后面经过同学的帮忙指正,才打开了扩展功能。(2) 、在做几何变换过程中,出现了很多错误,原因是几何方面公式知识理解较少,再加上对程序的不熟悉造成的;(2) 、在很多次操作中,出现脚本错误和在输入输出要素时出现红叉,没能成功完成,原因是细节不注意,有些小细节没做好(3) 、3 次实验课从简单到困难,执行了几何变换的一些重要的基本操作步骤,好多次坐到后面就做不下去了,后面的很多数据都市一前面的为基础。3、 保证实验成功的关键问题(1) 在操作时,遇到做不来的,应该立即向人请教,或者询问老师不可一直纠结于一个问题,否则会导致做不完;(2) 上课认真听讲,保证对每个概念都有所了解,在做实验时才会更加快捷,更加有效率。(3) 在操作时,应该按照步骤认真完成,不能忽略任何一步,否则可能会导致后面的实验失败(4) 多进行练习,熟悉程序和操作方法。指导教师评语和实验得分:实验得分: 签名: 年 月 日