浅谈WGS84与北京54之间的转换.docVIP

浅谈WGS84与北京54之间的转换 葛岚 （德清县城乡规划所 湖州313200） 摘 要：GPS测量得到的是WGS84坐标系下的坐标，而实际应用中较多使用的是北京54坐标，如何实现WGS84坐标与北京54坐标系的转换，一直是GPS应用中关心的热点。本文详细介绍了GPS定位结果转换至北京54平面坐标系的两种坐标转换模型，并对实验结果进行了分析比较。 关键词：GPS；平面转换；空间转换 1前言测绘信息网 随着GPS技术的快速发展，GPS定位精度的不断提高，使得GPS技术在测量中的应用也越来越广泛。由于GPS卫星星历表示在WGS-84坐标系中，因此算得的GPS定位结果也直接表示在WGS84全球坐标系中。而我国测绘成果大多表示在北京54坐标系中，它以克拉索夫斯基椭球为参考椭球，投影方式为高斯克吕格投影，以3度或者6度带划分整个中国区域。本文将从平面转换模型和空间转换模型两种方法中详细介绍这两种坐标系之间的转换实现的方法。 2平面转换模型测绘信息网 平面转换可以采用两种方法：一是将WGS84的大地坐标直接采用高斯投影的方法获得平面坐标，然后将该平面坐标经过平移、旋转和缩放转换至北京54坐标；二是假设北京54坐标系中的克拉索夫斯基椭球的中心与坐标轴方向与WGS84椭球相一致，将GPS测得的空间坐标利用克拉索夫斯基椭球参数换算至北京54椭球下的大地坐标，然后再将该大地坐标经高斯投影获得其在北京54椭球下的投影坐标，最后将该投影坐标经平移、旋转和缩放与实际的北京54坐标相拟合。本文主要对第二种方法进行了论述。 （1）由克拉索夫斯基椭球参数，将GPS测得的WGS84空间直角坐标（X,Y,Z）84换算至北京54椭球下的大地坐标形式（B,L,H）54： （1式） 式中B的计算需要采用迭代法，迭代时取B的初值,用B的初值B0代入计算B1，反复迭代直到最后两次B值之差小于允许误差为止。 （2）将（1式）中所求得的北京54椭球下的大地坐标（B,L）54投影为高斯坐标（x,y）‘54测绘信息网 （2式） 以上步骤是在假定54椭球与WGS84椭球的中心与坐标轴相同的前提下进行的，因此投影后获得的平面坐标还应考虑旋转平移缩放问题，以与实际相吻合。若GPS测定的大量点中，已知部分点在北京54坐标系下的坐标为,这些点通过以上方法计算所得的平面坐标为，则通过平面转换模型，两坐标系之间有如下关系： （3式） 其中为坐标平移量，m为尺度因子，为两平面坐标系的旋转角度。 为求出（3式）中的平移量、尺度因子和旋转角度，至少需要已知两个GPS点的北京54坐标，如多于两个点可利用最小二乘平差求出参数。 对所有的GPS点经过以上两个步骤及（3式），即可求出其在北京54坐标系下的平面坐标。水准高程可由大地高h84减去大地水准面差距求得。 平面模型原理简单，数值稳定可靠，但仅适用于范围较小的工程使用，对于大范围的GPS测量应使用空间转换模型。测绘信息网 3空间转换模型 由于WGS84椭球与北京54坐标系所属的克拉索夫斯基椭球有差异，因为要将WGS84空间坐标系转换为北京54直角坐标系，首先需完成WGS84椭球到克拉索夫斯基椭球之间的转换，两不同参考椭球之间的空间转换可采用布尔萨公式： （4式） 其中为两空间坐标系之间的三个平移参数，为三个旋转参数，m为尺度因子。测绘信息网 当已知至少三个GPS点在北京54坐标系下的坐标时，由（4式）利用最小二乘法可求出两不同空间坐标系之间的转换七参数。求出七参数后，对于每一个GPS点的空间直角坐标（X,Y,Z）84可通过以下几步转换为北京54坐标系下的平面坐标。 （1）利用（4式）将GPS测定的（X,Y,Z）84换算为（X,Y,Z）54； （2）根据54参考椭球，利用（1式）将（X,Y,Z）54换算为大地坐标（B,L,H）54， （3）由高斯投影正算（2式）将（B,L）54投影为高斯坐标（x，y）54 空间转换模型适用于大范围GPS测量使用，但在实际施工过程中，根据施工经度的要求，又分为三种情况：利用（4式）转换的，需要解求七个参数，故称为七参数转换模型；若其中的缩放比例不变，不需求尺度因子，则称为六参数转换；如尺度因子和旋转参数均不考虑，则称为三参数转换。测绘信息网 对于七参数模型的求解，至少需要三个公共点；六参数模型的求解也至少需要3个公共点，因为尽管两个公共点有6个坐标分量，可列出6个观测方程，但这6个坐标分量种，只有5个是独立的；三参数模型至少需要一个公共点。 4解算流程测绘信息网 平面转