手持GPS坐标系统转换参数求解与设置.docVIP

手持GPS坐标系统转换参数求解与设置 　　【摘 要】手持GPS全球定位系统所提供的是WGS84系统坐标，而我国目前应用的地形图大多属于北京54坐标系或西安80坐标系，不同的坐标系统需要进行转换，求解出DX、DY、DZ、DA、DF这五个转换参数后并输入GPS中，就可在GPS仪器上自动显示该点对应的北京54坐标系（或西安80）的坐标值了。本文给出了转换参数的求解与设置方法。 　　【关键词】GPS；坐标转换参数；求解；设置 　　1.前言 　　导航型手持GPS接收机目前已经成为野外地质工作者经常使用的定位仪器，既操作简单又便于携带，且定位精度可满足野外工作的需要。GPS导航系统所提供的是WGS84系统坐标，而我国目前应用的许多地形图大多属于北京54坐标系或西安80坐标系，因为不同坐标系之间存在着平移和旋转关系，对于同一点的坐标，可能会有几十米到上百米的差异，于是就出现了如何把WGS84坐标系转换为北京54或西安80坐标系的问题。因此，要使手持GPS所测量的数据（WGS84坐标）转换为自己需要的坐标（北京54或西安80），必须求出两个坐标系之间的转换参数。只要计算出五个参数（DX、DY、DZ、DA、DF）并输入GPS中，就可以在GPS仪器上自动显示该点对应的北京54坐标系（或西安80）的坐标值了。 　　在实际应用中，野外地质工作者经常遇到的棘手问题即是上述不同坐标系统转换参数的求解问题。 　　2.转换参数的求解 　　2.1三种常用坐标系统的简要介绍 　　北京54坐标系、西安80坐标系、WGS84坐标系这三个坐标系统是当前国内常用的，它们均采用不同的椭球基准（椭球参数见表1）。北京54坐标系属三心坐标系，大地原点在原苏联的普而科沃，该坐标系采用克拉索夫斯基椭球作为参考椭球，高程系统采用正常高，以1956年黄海平均海水面为基准；西安80坐标系的大地原点在陕西省西安西北的永乐镇，简称西安原点，椭球参数选用1975年国际大地测量与地球物理联合会第16界大会的推荐值，简称IAG75地球椭球；WGS84坐标系的原点位于地球质心，Z轴指向BIH1984.0协议地极（CTP），该坐标系是使用更高精度的VLBL、SLR等成果而建立的，我们常用的手持GPS定位即采用WGS84坐标系。 　　2.2手持GPS五参数的含义 　　通常，我们使用的手持GPS中有五个参数需要进行设置，即DX、DY、DZ、DA、DF。它们代表的含义如下： 　　（DX，DY，DZ）：这三个参数是两坐标系的原点平移参数，即原坐标系的原点在新坐标系中的三维坐标值（这里的“原坐标系”指的是WGS84坐标系，“新坐标系”指的是北京54或西安80坐标系）。 　　DA：原坐标系统参考椭球体长半轴与新坐标系统参考椭球体长半轴之差（a原-a新）。 　　DF：原坐标系统参考椭球体扁率与新坐标系统参考椭球体扁率之差（f原-f新）。 　　两个需要转换的坐标系统一旦确定以后（比如将WGS84坐标转换为北京54坐标），DA、DF即是常数，对于北京54坐标系来说，DA=-108，DF=+0对于西安80坐标系来说，DA=-3，DF=-0.0000000025（单位为m）。 　　DX，DY，DZ三个参数在不同地区，其值是不同的。因此，需要根据所处地区的若干个控制点坐标来求解出这三个参数的值。 　　下面以北京54坐标系为例，介绍求解坐标转换参数的方法。 　　2.3收集已知控制点资料 　　（1）首先在所处区域内或附近收集已知点（等级点或GPS“B”级网网点）的坐标。数量最好3个以上，分布要均匀，点位要求最好是周围无电磁波干扰、视野开阔、卫星信号要强。 　　到当地测绘管理部门收集或购买这些点的北京54坐标系的大地经纬度（B、L）或高斯平面直角坐标（X、Y），海拔高程h，高程异常值ξ（即大地水准面相对椭球面的高度）以及相应的WGS84坐标系的大地经纬度（B、L），大地高程H。 　　（2）如果只收集到了北京54坐标，其相应的WGS84坐标可以通过GPS测量来获得。 　　方法是：按照后文第3章“GPS坐标系统转换参数的设置”中第3步叙述的“坐标系统”选项设置为“WGS84”后，对各已知控制点进行测定WGS84坐标系下的B、L、H值，以供下面计算转换参数时使用。 　　2.4各已知点的三维直角坐标的计算 　　WGS84坐标转换为北京54坐标，实际上是两个坐标系的空间三维直角坐标的相互转换。因此，首先要计算出各自的三维直角坐标。 　　计算公式： 　　对于同一空间点，大地坐标系与空间直角坐标系有下列转换关系式： 　　X=（N+H）cosBcosL 　　Y=（N+H）cosBsinL 　　Z=[N（1-e2）+H]sinB 　　其中： 　　N=a/（1-e2sin2B）0.5（椭球的卯酉圈曲率半径） 　　H=相对椭球面