基于Google Maps 与CASS9.0 的小区域地形图测绘方法应用分析

2019-04-23郝长春郑攀攀

陕西水利 2019年3期

郝长春，郑攀攀

（安徽省水利水电勘测设计院勘测分院，安徽蚌埠 233000）

0 引言

近年来，随着基础设施建设的快速发展，如何快速的测绘满足工程要求的地形图一直以来都是测绘行业研究的重要课题。地形图测绘技术的发展使测绘手段不断更新，从全站仪配合电子平板模式、GPS-RTK 配合草图模式、大型航空摄影模式，到低空无人机模式。针对使用频率高、更新速度快、植被覆盖率低、资金投入少的小区域地形测绘，上述方法均有不足之处。探寻采用高清高分辨率的Google Maps 卫星影像结合GPS-RTK 及CASS9.0 技术解决小区域大比例尺地形测绘的方法，能有效的降低生产成本，提高作业效率，加快地形图的更新速度。

1 Google Maps 和南方 CASS9.0 技术

1.1 Google Maps 技术

卫星地图又称卫星遥感地图，是通过卫星在太空中探测地表物体的电磁波，从中提取物体的位置、几何形状及部分属性信息，再将这些信息通过转换、识别，得到对应物体的图像，即卫星地图[1]。Google Maps 是 Google 公司在 2005 年推出的电子地图产品，能提供不同分辨率的卫星影像[2]。与传统地图不同，卫星地图是以卫星为媒介，向用户及时反馈地球表面真实面貌的图像，具有更强的时效性、直观性和更高的精确度。Google 卫星影像地图为WGS84 大地坐标系，分辨率0.5 m～10 m，可以下载1～19 级的各种影像，高清影像的精度达到0.5 m[3]。卫星地图还具有更新速度快的特点，Google Maps 区域地图每一个季度更新一次，全球地图每18 个月跟新一次，在线卫星地图每天更新一次，每季度全部跟新一次，非常适合于地形图的调绘和更新。

1.2 南方CASS9.0 技术

CASS9.0 是广东南方数码科技股份公司2018 年全新发布的地籍地形成图软件，是基于全球公认的最优秀的图形与设计平台AutoCAD 软件的基础上通过二次开发完成的[4]。针对栅格图像，南方CASS 提供了图像的插入、纠正、裁剪等诸多功能。本文将基于南方CASS 在栅格图像处理方面的技术支持，讨论Google 影像在地形测绘中的应用。

2 Google Maps 影像与南方CASS9.0 技术的结合

2.1 Google Maps 影像的获取

根据测绘工程项目的特点和要求从Google 卫星影像客服端下载相应级别和范围的影像。以成都水经注信息科技有限公司推出的水经注万能地图下载器为例，该软件是基于谷歌、天地图等数十种地图服务器的多线程卫星地图下载器。通过该下载器可以实现测区规划(在线规划或者导入下载范围)、中央经线设置、坐标转换参数计算(WGS-84 与国家坐标系之间的七参数)、影像级别选择、自动下载拼接等功能完成影像图的下载工作[5]。影像在拼接导出的过程中有多中影像格式可供选择，如GeoTIF、GeoJPG、PNG 等。随同影像一起下载的还有.TXT 格式的WGS-84 大地坐标文件，如图1，它是影像所在矩形边框4个拐角的坐标(通常为十进制格式)。通过七参数投影转换，可将WGS-84 坐标转换为国家标准的高斯平面直角坐标，如图2，该成果可以用于后期图像的纠正。

图1

图2

2.2 在CASS9.0 中加载卫星影像

南方CASS9.0 不但具有地形成图、地籍成图、工程测量等强大功能，而且还提供了栅格影像的处理功能。影像加载的基本步骤如下[6]，如图3。

（1）打开CASS9.0 的绘图桌面，按照工程要求在“绘图处理”中设置作业区域的比例尺。

（2）通过“绘图处理”中的“展野外测点点号”将影像4 个拐角的高斯平面直角坐标和采集的野外地物平面特征点展绘在绘图区，作为影像纠正的依据和检验点，如图4。

图3

图4

（3）在“工具”中选择“光栅图像”，点击“插入图像”，打开“图像管理器”，单击“附着”，找到影像所在的位置，加载“.tif”格式的影像，进入“图像”对话框，勾选“插入点”和“缩放比例”，点击“确定”进入绘图区域，在影像拐点和野外采集的特征点附近任意选择影像的插入点和缩放比例。

2.3 基于CASS9.0 的卫星图像纠正

根据影像拐点的高斯平面坐标(GJ+序号)和野外采集的平面特征点坐标(JZ+序号)，对插入的Google 卫星影像进行纠正，使影像上地物的位置和几何形状与实际趋于一致，使平面精度达到地形测绘的要求。

（1）进入“工具”，“光栅图像”，选择“图像纠正”，根据提示选取要纠正的图像，进入“图像纠正”对话框。

（2）在“图像纠正”对话框的“图面”项，点击“拾取”，选取影像边框左上角拐点的坐标，再在“图像纠正”对话框的“实际”项，点击“拾取”，选取影像边框左上角对应的GJ1 点的高斯平面直角坐标，最后点击“添加”，将该组点添加到“已采集控制点”的列表里。以此类推，将其余组的拐角点加入到“已采集控制点”列表里面。

（3）针对野外采集的平面特征点，拾取图面坐标时，只需将卫星影像适当放大，找到野外采集时特征点的平面位置，拾取该点图面坐标。“实际”坐标即是野外实地采集的特征点的坐标，如JZ1 等，以此类推，加载剩余野外控制点的坐标。

（4）纠正方法包括赫尔默特法、仿射变换、线性变换、二次变换、三次变换，其中赫尔默特要求至少3 个控制点，仿射变换至少4 个控制点，线性变换至少5 个控制点，二次变换至少7 个控制点，三次变换至少11 个控制点，这些方法均包括1 个多余控制点。在实际工作中可根据控制点情况进行选择。

（5）“误差”主要是给出图像纠正的精度。

3 实例分析

3.1 工程概况

为了保障安徽段淮河在汛期的行洪安全，现对该段的保庄圩堤防的达标进行勘测。胡姚保庄圩是其中的一个圩子，因现有基础资料缺乏，相关地形、地貌资料不足，为满足本工程科研报告设计需要，要对该圩的堤防进行测量，提供工程需要的1∶2000 地形图。目前收集到的资料主要是该圩2005 年的1∶5000 竣工数字化地形图，坐标系为1954 北京坐标系。由于该资料年代较久，现时性差，成图不够规范，根据工程需要，拟采用卫星影像为底图，对原有地形图进行修测。

3.2 影像加载精度分析

根据测区范围应用谷歌卫星地图下载器X2 下载测区卫星影像及影像4 个拐角WGS84 坐标，通过测区的布尔莎七参数模型，将其转换为测区1954 北京坐标系。在野外，应用GPS-RTK 连接省CORS 系统，设置好测区七参数，采集测区影像范围内清晰明显的平面特征点，如低矮的房角，围墙拐角，水泥地拐角等。从中选取8 个控制点用于卫星影像的纠正，其中包括4 个影像拐角点和4 个实测平面特征点，采用二次变换的方法进行纠正，纠正误差均为0，满足纠正要求，如表1，图5 和图6。