探讨不同坐标系的坐标转换方法.docxVIP

探讨不同坐标系的坐标转换方法 摘要:本文主要讲述了应用cass成图软件，用图解的方法进行转 换参数的计算及图形、包括控制点在内的数据转换，并举例进行了 说明。 关键词 地籍调查坐标系转换 目前，城镇地籍调查成果还是把点的平面坐标与高程分开处理，两 种坐标系下的成果转换也是把平面坐标与高程分别转换。 本文就城镇地籍调查中不同坐标系的成果，尤其是控制点成果，如 何转换进行探讨。 不同坐标系下坐标转换的理论基础 不同坐标系的坐标转换与换带计算是有区别的。换带计算是在同一 个参考椭球基础上进行坐标变换，本文所指的不同坐标系的坐标转 换是在不同的参考椭球基础上的坐标转换。 1980年国家大地坐标系与1954年北京坐标系。在参考椭球的大小、 定位及起始方向的确定等方面均存在差异。因此两种坐标系下的测 量成果转换，不但需要进行平移、旋转，通常还需要进行尺度的缩 放。 平面坐标是釆用高斯一一克吕格投影的高斯平面直角坐标。坐标转 换参数有四个。其中二个坐标平移参数Ax、Ay, 一个旋转参数△ a, 一个尺度参数k。要求得这四个转换参数，至少需要两个同时具 有两套坐标系统坐标的控制点（公共点），按式（1）计算转换参数。 如果公共点超过二个，则由（1）式按最小二乘原理计算。 x 二x +△ + k(x- x )cosAa-k(y-y ) sinAa y = y +△, +k (y-y ) cosAa+k (x- x ) sinZXaC) 式中x、y为新坐标系下的坐标，x, y是旧坐标系下的坐标，x , yo是计算坐标原点(旧坐标系)o 一般的成图软件都具有计算坐标转换参数的功能。但有的成图软件 显示的转换参数取位不够，如cass成图软件，其坐标平移参数和 尺度参数仅显示3位小数，旋转参数取位到0. 这对于一般的 图形转换来说，其精度可以满足。但对于控制点的转换而言，其精 度就显然不够。如果要转换四等控制点，则尺度参数要达到6位小 数，旋转参数取位到0. 01”。也有的文献给出的坐标转换计算公式 不正确。 用cass成图软件进行不同坐标系成果的转换研究 在求出转换参数后。先对图形进行旋转、缩放，然后再进行平移。 通过转换参数不仅可以转换图形，而且可以转换数据。如果按式(1) 直接进行计算转换，工作量比较大。即使是利用程序进行计算，也 必须先把原坐标系下的控制点数据逐个输入，手工输入难免有错误 出现；另外，如果控制点数量较多的话。数据输入也是一件不容易 的事。对于已有控制点的转换。如何更快捷、方便和准确呢？这就 是本文主要讨论的内容。 这里通过在全国应用比较广泛的cass成图软件来讨论控制点数据 转换问题，其基本思路是：把原控制点数据展绘到图上，先对图形 进行缩放、旋转和平移，然后在新坐标系下的图形中提取控制点数 据。这些工作在cass成图软件中都有相应的菜单或命令，通过图 形处理的方式来达到成果转换的目的。简而言之，就是通过图解的 方法来实现不同坐标系成果的转换。下面介绍转换步骤和方法。 1求转换参数 式(1)的计算在cass成图软件可以全部实现。 分别展绘同时具有两套坐标系统坐标的控制点(公共点)，量取 任意两个点之间的边长和坐标方位角： 计算转换参数 同一条边的两个坐标方位角之差即为旋转参数，两个长度之比即为 尺度参数。如果有n个公共点，则可求得n=n(n-l)/2组转换参数, 若这n组参数之间最大的差异小于同等级控制网起算数据的要求， 则求其相应的平均值为两套坐标系统的转换参数。 2. 2对原有控制点数据进行缩放和旋转 先执行“新建图形文件”命令：然后用“插入图块”命令，把展绘 了原坐标系统控制点的图形插入到新建图形文件中。在图1所示的 “插入”对话框中，浏览到待插入的“控制点图”.然后在“缩放 比例”栏选择“统一比例”复选框.在X栏内输入缩放参数；在“旋 转”栏输入旋转角度，注意若是逆时针旋转，旋转角度为负，而顺 时针旋转时旋转角度为正。 图1 “插入图块”命令对话框 2. 3对图形进行平移 选择一个公共点(最好是测区中部位置)，把经过缩放和旋转的图形 平移到新的坐标系下。 4提取控制点数据 cass成图软件“工程应用”下拉菜单中有“控制点生成数据文件” 子菜单。 执行“控制点生成数据文件”子菜单，在屏幕弹出的“输入坐标数 据文件名”对话框中输人新坐标系统下控制点数据文件的保存路径 即可。 该命令能自动识别控制点类型，比如说，gps控制点，导线点，图 根点，水准点等。下面的数据是从某幅地形图中提取的部分控制点 数据文件： c02 —b25, 530702. 150, 2839228. Ill, 301? 280 2 c08—汝崇一 1, 530697. 040, 2839257. 985, 285360 c09 —a6, 530606. 8