手持GPS坐标转换参数在EXCEL中的计算方法.doc

手持GPS坐标转换参数求解方法及在EXCEL中的实现 摘要：介绍了手持GPS坐标转换参数求解的方法步骤，并通过具体实例详述了在EXCEL中的实现。 关键词：手持GPS 坐标转换参数 方法步骤 EXCEL 实现 全球定位技术以其精度高、速度快、费用省、全天候、操作简便等优良特性在矿产资源调查、林业管理、农业、气象、测绘、航海航空等领域有着广泛应用，尤其是手持式GPS接收机作为新型定位、导航的现代化工具，在矿产资源调查中发挥着重要作用。但GPS定位技术所采用的坐标系统为WGS－84坐标系统，空间点位描述用大地坐标B，L，H表示，即纬度、经度和大地高，大地高基准面为WGS－84椭球面，而我们国家所采用的平面坐标系统为1954北京坐标系统（简称北京54）或1980年西安坐标系统（简称国家80），投影方式采用6o带或3o带高斯投影，高程系统为正常高系统。国家大地控制成果及1：5万、1：1万等基本地形图成果均采用国家坐标系统，点位描述为x、y、h。因些，在实际应用中，如何将WGS－84坐标系统下的B，L，H转换为国家坐标系统下的x，y，h是非常重要的，由于不同地区高程异常ξ大小的差异，其转换参数也所有不同。可以将求解好的参数输入手持式GPS接收机相应项目中自动完成转换。下面就转换参数求解方法步骤和在EXCEL中的实现进行了详细介绍。 1 GPS坐标转换参数求解的方法步骤 将WGS－84坐标系统下的B、L、H转换为国家坐标系统下相应的x、y、h，需要求解的参数有DA、DF、DX、DY、DZ等五个。 1.1　DA求解 DA＝aWGS-84－a北京54　或　DA＝aWGS-84－a国家80 （1） 即WGS-84椭圆长半轴与克拉索夫斯基椭圆长半轴（北京54）之差或WGS-84椭圆长半轴与1975年国际椭圆长半轴（国家80）之差。 1.2　DF求解 DF＝αWGS-84－α北京54　或　DA＝αWGS-84－α国家80 （2） 即WGS-84椭圆扁率与克拉索夫斯基椭圆扁率（北京54）之差或WGS-84椭圆扁率与1975年国际椭圆扁率（国家80）之差。 1.3　DX、DY、DZ求解 根据公共点在WGS－84坐标系统下的大地坐标BWGS-84、L WGS-84、H WGS-84和国家坐标系统下的高斯平面坐标x、y、h求解，具体步骤是，首先将x、y通过高斯反算转换为大地坐标纬度、经度，高程h不变，再将两套大地坐标值分别换算为空间直角坐标X、Y、Z，最后求解出DX、DY、DZ。 1.3.1高斯坐标反算 将国家平面坐标通过高斯反算转换为大地坐标纬度和经度，高程值不变，即B国、L国、H国。中央午线经度按公式（3）或（4）求解，其中L值可采用WGS－84坐标系统的经度值，高斯反算可按公式（5）进行，不同坐标系统相应椭球参数值见表1。 6°带L0的计算： N6＝[L/6]+1 L0=N6×6°-3° (3) 3°带L0的计算： N3＝[（L－1.5°） 表1 不同坐标系统相应椭球参数值 符号 1954年北京坐标系 1980年国家坐标系 WGS-84坐标系 备注 ɑ 6 378 245m 6 378 140m 6 378 137m 椭圆长半径 α 1/298.3 1/298.257 1/298.257 223 563 椭圆扁率 e2 6.693 421 622 965 949×10-3 6.694 384 999 587 952×10-3 0.006 694 379 901 3 第一偏心率 e′2 6.738 525 414 683 497×10-3 6.739 501 819 472 927×10-3 0.006 739 496 742 27 第二偏心率 C 6 399 698.901 782 711 1m 6399 596.651 988 010 4m k0 1.570 460 641 22×10-7m-1 1.570 486 874 73×10-7m-1 k1 5.051 773 902×10-3 5.052 505 593×10-3 k2 2.983 868×10-5 2.984 733 5×10-5 k3 2.415×10-7 2.416×10-7 k4 2.2×10-9 2.2×10-9 1.3.2 大地坐标转换为空间直角坐标 大地坐标转换为为空间直角坐标的计算公式如（6）。 （6） 式中：N——法线长度，可由式算得。 将WGS-84坐标系统中的BWGS-84、L WGS-84、H WGS-84和通过高斯反算后国家坐标下的大地坐标B国、L国、H国分别代入公式（6）求出XWGS-84、YWGS-84、ZWGS-84和X国、X国、X国。 1.3.3 DX、DY、DZ求解 2