基于VC++的常用大地坐标转换程序实现.doc

基于VC++的常用大地坐标转换程序实现 赵 丹 （中铁第四勘察设计院集团有限公司，湖北武汉，430063） Geodetic Coordinates Transformation Programming Based on VC++ Zhao Dan 摘要：介绍了工程应用中常见的大地坐标系统，并介绍了三参数、四参数、七参数坐标转换模型，以及大地坐标系与空间直角坐标系之间的转换模型，同时利用VC++开发工具设计了实用的坐标转换程序，实现了相应的坐标转换功能。 关键字： 坐标转换；WGS-84；BJ54；GZD80；VC++； 1 引言 常用的坐标转换一般包括各种空间直角坐标系与大地坐标系，地心空间直角坐标系与参心空间直角坐标系，以及不同参心空间直角坐标系之间的相互转换。当不同坐标系之间存在严密的数学转换模型时，可以采用相应的模型之间进行坐标转换。当不同坐标系之间不存在明确的函数关系时，常采用三参数、四参数、七参数坐标转换模型。本文介绍常用的坐标转换模型，同时利用VC++开发工具设计了实用的坐标转换程序，实现了相应的坐标转换功能。 2 坐标转换模型 2.1 七参数法 工程应用中常常需要将不同参考椭球下的空间直角坐标转换。对于既有平移、旋转，又有缩放的两个空间直角坐标系的坐标转换，相应地存在3个平移参数、3个旋转参数以及1个尺度参数，共计7个参数，称为七参数模型。相应的坐标转换模型为： 式中，为3个平移参数；为3个旋转参数，为尺度参数。为了求得这7个转换参数，至少需要3个公共点，当多余3个公共点时，按最小二乘法求得7个参数的最或然值。 2.2 三参数法 三参数法是七参数法的特殊情况，当测区所覆盖的范围不大，精度要求不高时，可以只考虑坐标轴的平移，而不考虑旋转和缩放，相应地只存在3个平移参数，即三参数模型。相应的坐标转换模型为： 2.3 四参数法 由于高斯投影带的划分不同，会造成常用的大地坐标系（如BJ54,GZD80,WGS-84等）在某些地区的投影不能满足工程要求的精度时，此时就需要建立地方独立坐标系。此时面临的问题是如何将常用的大地坐标系转换成地方独立坐标系。这里就要用到四参数模型进行转换，相应的坐标转换模型为： 2.2 空间直角坐标系与大地坐标系的相互转换[1] 同一地面点在空间直角坐标系中的坐标和在大地坐标系中的坐标可用如下两组公式进行转换。 从大地坐标转换到空间直角坐标的转换模型如下式： 利用点的空间直角坐标到大地坐标的转换模型如下式： 式中：e为参考椭球第一偏心率，可由长短半径按式算得。N为法线长度，可由式算得。大地纬度的计算比较复杂，通常采用迭代法。迭代时可取，用的初值计算和，按进行第二次迭代，直至最后两次值之差小于限差为止。 3 常用大地坐标系 3.1参心坐标系 参心大地坐标的应用十分广泛，它是经典大地测量的一种通用坐标系。根据地图投影理论，参心大地坐标系可以通过高斯投影计算转化为平面直角坐标系，为地形测量和工程测量提供控制基础。由于不同时期采用的地球椭球不同或其定位与定向不同，在新中国历史上先后建立了北京54坐标系（旧）（简记BJ54旧）、1980年国家大地坐标系（也称西安80坐标系，简记GZD80）和1954年北京坐标系（新）（简记BJ54新）等三种参心大地坐标系[2]。 1．1954年北京坐标系（BJ54(旧)） 1954年总参谋部测绘局在有关方面的建议和支持下，采取中国和前苏联远东一等锁相连接、平差，这样传过来的坐标定名为1954年北京坐标系（旧）。几十年来，中国在该坐标系上完成了大量的测绘工作，实施了天文大地网局部平差，通过高斯-克吕格投影，得到点的平面坐标，测绘了各种比例尺的地形图。但是随着测绘新理论、新技术的不断发展，人们发现该坐标系存在诸多缺点。另外，该坐标系是按分区平差逐步提供大地点成果的，在分区的结合部产生了较大的不符值。但该坐标系确实在测绘生产中发挥了巨大的作用，至今仍在一些部门使用。1954 年北京坐标系（旧）的缺陷是: ①椭球参数有较大误差; ②参考椭球面与我国大地水准面存在着自西向东明显的系统性的倾斜; ③几何和物理大地测量应用的参考面不统一; ④定向不明确; ⑤当时使用的仪器、测量方法的落后,致使在大面积长距离传递中误差累计较大,且系统只是进行了局部的平差。因而不可避免地出现一些矛盾和不够合理的地方。后者优点是: ①椭球短轴平行于地球地轴; ②起始大地子午面平行于格林尼治天文台起始子午面; ③椭球面同似大地水准面在我国境内最为密合; ④系统经过了整体平差。所以向1980年国家大地坐标系过渡及1980年国家大地坐标系的全面使用就成了必然。 2．1980年国家大地坐标系（GDZ80） 为了解决1954年北京坐标系中存在的椭球参数不够精确、参考椭球与