坐标转换程序,可以进行地方坐标与WGS-84坐标的相互转换,同时具备参数求解功能。
下面对这个工具进行介绍:
一般而言椭球之间的转换,比较严密的是用七参数法(或称布尔莎模型),即X平移,Y平移,Z平移,X旋转,Y旋转,Z旋转,尺度变化K。要求得七参数就需要在一个地区需要3个以上的已知点(7个参数至少7个方程可解,所以需要三个点列出9个方程),如果区域范围不大、最远点间的距离不大于30Km(经验值)的情况可以用三参数,即X平移,Y平移,Z平移,而将X旋转,Y旋转,Z旋转,尺度变化K视为0,所以三参数只是七参数的一种特例。
七参数模型的实质是用一个局部椭球去拟合地方坐标系的形态;所以转换后获得的地方椭球高就是水准高。当然我们也可以把平面和高程两个方向分别进行拟合。例如平面用四参数模型拟合,高程方向则用二次曲面等模型来拟合。这样分开处理的模式相对七参数模型自由度更高。但是由于四参数模型参数较少,表达能力较弱,通常只用于小区域坐标转换。
综上所述,从实用的角度出发,坐标转换程序提供了两种转换策略供给客户选择使用:
1.七参数模型,一步得到地方平面和水准数据。
2.四参数加高程拟合模型,分两步得到地方平面和水准数据。
二、使用软件进行坐标转换
本软件将坐标转换参数以文件管理,用户可以将定义好的一组转换参数作为一个文件(*.dam)保存下来,下次使用时从文件菜单中选择打开这个文件来调用所有已有的转换参数。
涉及到坐标转换参数一般指的是 椭球参数,投影参数,七参数,四参数,高程拟合参数,水准格网文件。所有这些参数集成到如下界面进行输入。输入完成后,取一个文件名称,点击"保存",即可在程序目录下的“GeoPath”目录下生成*.dam参数文件。
图1.参数输入与设置界面
在该界面完成输入参数,或者点击 "V"下拉按钮选中调入已有的坐标转换文件后,点击确定,即可回到主界面进行坐标的正反算:
图2.坐标转换主界面
程序还提供了简单的批量转换功能,在主界面切换到"文件"选项卡,即可调入预先准备好的文件,进行批量正算或者反算:
图3.批量转换功能界面
三、参数计算
当用户有一组控制点(该组点具有WGS84坐标和当地坐标),即可使用本软件进行参数求解。如前面所述,本软件提供七参数模型参数和四参数加高程拟合模型的参数求解,两种模型的求解均在同一个界面完成,方便用户对比选用不同模型的精度。主界面点击"工具" ->"参数计算器" 即可打开参数计算界面(如果未输入椭球和投影参数,将提示并弹出参数设置界面)
图4.参数计算界面
参数计算的流程是:
1.输入基本参数:即先输入当地椭球和投影参数。
2.导入数据:逐个添加坐标点对数据或者先准备好文本格式的数据文件再点击“导入文件”(备注:文件格式为Name,B,L,H,x,y,h)
3.计算参数:本软件支持两种模式的坐标转换点击"七参数"或"四参数+高程拟合"若用后者,请先选定高程拟合模型。
4.查看结果:在结果栏会输出计算出的参数,用户可以复制出以保存。
5.运用参数:点击"生成坐标系统",查看转换参数,椭球参数及投影参数。确定无误后,可以输入名称并保存为*.dam文件,用于中海达其他软件。
注意:请切换到"图形"界面查看起算点的几何分布情况。
四、运用实例
例1(用户使用的当地坐标是国家坐标系)
用户在一个A测区内使用RTK GPS接收机测量了一组WGS-84的坐标,现在希望将其转换为北京54。用户有A测区的四个控制点,这些控制点既有WGS-84坐标,也有北京-54坐标。同时已知地方坐标的投影参数。
分析:
用户的控制点数量超过3个, 且坐标为国家坐标,可以直接使用七参数
步骤:
1、进入参数设置界面,新建坐标转换文件" A市54坐标系",便于下次使用转换是不用重新输入,并首先选择椭球为BJ54并输入投影参数:
图5.参数输入
2、进入参数求解工具界面,录入四个控制点坐标数据,也可以实现按照格式准备好txt文件进行导入(推荐方法):
图6.控制点数据输入
3、切换到图形界面查看控制点分布形态:
图7. 查看控制点分布
如上图显示,该组控制点分布呈直线形态,形态不是很好,我们可以首先尝试用七参数模型,同时对比一下用四参数加高程拟合模型的方法。
4、点击“计算七参数”,即可完成七参数的计算,界面将输出参数的残差及具体参数,如果残差超过2个cm,表格将以黄色提醒。
图8. 参数结果及残差信息
在这一步,可以点击“计算四参数与高程拟合”,对比两种不同的转换流程的精度。根据残差信息进行选用。
5、确定参数计算无误后,点击"生成坐标系统",即可弹出自动形成的坐标转换参数集,检查无误,即可点击保存,生成*.dam文件使用。
图8.参数结果自动录入
例2 用户要将WGS-84转换到地方坐标(中央子午线已知121:30)
分析:首先我们需要七参数转换的特点要求坐标之间不能相差太大。(这点很重要,可以查看相关大地测量的相关知识,有关椭球转换时的大角度与小角度转换问题)
WGS-84的坐标:31:25:06.30317,121:11:43.96187,14.0045
上海城市坐标为:20317.9891,-25826.1791,5.556
处理:使用WGS-84坐标不加任何参数转换为平面坐标为
3477405.6555,471050.9167,14.0045
可以看到X的大数差-3457000 Y相差:-500000 左右
可得投影中的东和北方向加常数(默认值加上上面的差值)
北方向:0-3457000=-3457000 东方向:500000-500000=0
可知投影参数如下设置:
求得结果如下: