2000国家大地坐标系转换的研究.docVIP

2000国家大地坐标系转换的研究 　　【摘要】如何有效的利用原有的测绘成果和测绘资料，建立新旧坐标系间的联系，研究它们间的影响因素，将是一项紧迫而且富有实际意义的课题。通过对每个坐标系之间相互转换的方法进行分析，对二维参心坐标系转换到三维地心坐标框架进行了一些研究，以空间直角坐标为过渡，采用一种以正常高为初值，迭代方法逼近参心坐标基准大地高，进而解算地心坐标系与参心坐标系转换七参数的方法。 　　【关键词】坐标转换；参心坐标系；地心坐标系； 　　1引言 　　由于1954北京坐标系和1980西安坐标系在我国存在了较长的时间，以前的很多测量成果都是以这两个坐标系为基础获得的，但随着卫星定位技术的快速发展，卫星定位技术赢得了广大测绘工作者的信赖，并广泛地应用于大地测量、工程测量、航空摄影测量、运载工具导航和管制、地壳运动监测、工程变形监测、资源勘察、地球动力学等多种学科，我国大地坐标系统建设进入了一个新的阶段。这也就存在着直接观测获得的成果（2000 国家大地坐标系成果）与落后的坐标系统之间的衔接问题。 　　2 1954北京坐标系、1980西安坐标系的定义 　　1）1954年北京坐标系为参心大地坐标系，大地上的一点可用经度L、纬度B和大地高H定位，它是以克拉索夫斯基椭球为基础，经局部平差后产生的坐标系，1954年北京坐标系可以认为是前苏联1942年坐标系的延伸。它的原点不在北京而是在原苏联的普尔科沃。 　　2）1980西安坐标系是为进行全国天文大地网整体平差而建立的。主要有以下特点： 　　（1）大地原点在陕西省泾阳县永乐镇； 　　（2）它是参心坐标系，椭球短轴Z轴平行于地球质心指向地极原点方向，大地起始子午面平行于格林尼治平均天文台子午面； 　　（3）X轴在大地起始子午面内与Z轴垂直指向经度0方向；Y轴与Z、X成右手坐标系；椭球参数采用IUG1975年大会推荐的参数为： 　　长轴： ； 　　扁率： 　　地心引力常数： 　　地球自转角速度： 　　椭球定位时按我国范围内高程异常值平方和最小为原则求解参数； 　　（4）多点定位； 　　（5）大地高程以1956年青岛验潮站求出的黄海平均海面为基准。 　　3CGCS 2000坐标系的定义方法 　　CGCS2000国家大地坐标系的定义包括坐标系的原点、3个坐标轴的指向、尺度比以及地球椭球的4个基本常数的定义。CGCS2000国家大地坐标系的原点为包括海洋和大气的整个地球质量中心；CGCS2000国家大地坐标系的Z轴由原点指向历元2000.0的地球参考极方向，该历元的指向由国际时间局给定的历元为1984.0的初始指向推算，定向的时间演化保证相对于地壳不产生残余的全球旋转，X轴由原点指向格林尼治参考子午线与地球赤道面（历元2000.0）的交点，Y轴与Z轴、X轴构成右手正交坐标系。采用广义相对论意义下的尺度。值得指出的是，这里a、f采用的是GRS80值，ω、GM采用的是IERS推荐值。 　　4 坐标转换的基本方法 　　基于旧测绘成果的继续利用，新国家大地坐标系统建立起来后，需将原有其它坐标系成果转换到新坐标系系统下，这就需要选取合理的转换模型来进行不同坐标系统之间的转换。通过选区分布均匀、点位稳定且具有一定密度的两种坐标系的重合点，计算坐标转换参数，实现坐标转换，并通过坐标转换精度评定对坐标转换质量进行评价。 　　1）高斯平面下转换 　　（1）平面四参数转换 　　由至少两个同名点求解两个平移参数，一个旋转参数，一个尺度参数。 　　其中： 为平移参数， 为旋转参数，k为尺度参数，x2，y2为输出平面直角坐标。 　　根据重合点位分布情况，选择有代表性的高精度控制点作为重合点和检验点。用最小二乘法求解相似变换参数。 　　（2）二次多项式拟合 　　由最小二乘参数求解法。 　　（3）比例缩放法 　　由参心（或地心）坐标变换为独立坐标（抵偿坐标）公式 　　公式逆转换： 　　式中， 为缩放系数；H为投影面高度，起算面为愿投影面； 为测区中心的平均曲率半径； 、 为参心（地心）3度带坐标； 、 为独立坐标系坐标； 、 为测区中心点3度带坐标。 　　2）在空间坐标系下转换 　　Bursa模型一般含3个平移参数、3个旋转参数和一个尺度参数。由若干个（一般多于3个）公共点坐标求得上述7个转换参数后即可进行不同坐标系的转换。相应的坐标变换公式为： 　　（1） 　　其中m为尺度变化参数； 　　、 、 为三个旋转参数； 　　、 、 为三个平移参数。 　　为了得到7个转换参数，至少需要3个公共点，当多于3个公共点时，可按最小二乘法7个参数的最或是指。采用Bursa七参数转换模型。坐标转换中误差应小于0.05m。 　　3）大地高及坐标转换七参数计算