WGS84与北京54坐标系之间的转换

WGS84 与北京 54 坐标系之间的转换第一章绪论l.i 概述坐标转化并不是一个新的课题,随着测绘事业的发展，全球一体化的形成，越来越要求全球测绘资料的统一。尤其是在坐标系统的统一方面.原始的人地测量工作主要是依靠光学仪器进行,这样不免受到近地面人气的影响,同时受地球曲率的影响很人，在通视条件上受到很人的限制,从而对全球测绘资料的一体化产生巨大的约束性。另外由于每一个国家的大地坐标系的建立和发展具有一定的历史特性，仅常用的大地坐标系就有150 余个。在同一个国家,在不同的历史时期由于习惯的改变或经济的发展变化也会采用不同的坐标系统。例如:在我国建国之后，为了尽快搞好基础建设,我国采用了应用克氏椭球与我国实际相结合的北京 54 坐标系;随着经济的发展北京 54 坐标系的缺陷也随之被表露的越来越明显，特别是对我国经济较发达的东南沿海地区的影响表现得更为明显，进而我国开始研究并使用国家 80坐标系。在实际生活中，在一些地区由于国家建设的急需，来不及布设国家统一的大地控制网，而建立局部的独立坐标系。而后，再将其转换到国家统一的人地控制网中，这些坐标系的变换都离不开坐标值的转化.在国际上，随着 1964 年美国海军武器实验室对第一代卫星导航系统一 NNSS 的研制成功，为测绘资料的全球一体化提供了可能。到 1972 年，经过美国国防部的批准，开始了第二代卫星导航系统的开发研究工作，即为现在所说的 GPS。此套卫星导航系统满足了全球范围、全天候、连续实时以及三维导航和定位的要求.正是由于 GPS 卫星的这些特性，这种技术就很快被广人测绘工作者接受。是由于坐标系统的不同，对 GPS 技术的推广使用造成了一定的障碍。这样坐标转换的问题再一次被提到了重要的位置。为了描述卫星运动，处理观测数据和表示测站位置，需要建立与之相应的坐标系统。在 GPS 测量中，通常采用两种坐标系统，即协议天球坐标系和协议地球坐标系。其中协议地球坐标系采用的是 1984 年 世 界 大 地 坐 标 系(WordGeodeticSystem1984WGS-84)其 主 要 参 数 为：长 半 轴a=6378137;扁率 21:298.257223563.而我国采用的坐标系并不是 WGS-84 坐标系而是 EJ-54 坐标系，这个坐标系是与前苏联的 1942年普耳科沃坐标系有关的，其主要参数为：长半轴 a=6378245;扁率 21:298.3.这就使得同一点在不同的坐标系卞有不同的坐标值，这样使测绘资料的使用范鬧受到很人的限制,并且对 GPS 系统在我国的广泛使用造成了一定的约束性,对我国的测绘事业的发展不利。为了解决这个问题，我国的测绘工作者做了大量的工作,并且已经有许多使用的商品化软件投入市场，其中以七参数法用的最为普遍。本次毕业设计也正是在这个理论基础之上进行的,但是由于此方法中参数较多，对于一部分 GPS 用户来说不容易使用足够的已知点坐标求解。并且现在国内所使用的坐标转换程序人多使用国外软件，除了成本较高以外，也不利于推广使用，同时人多软件中没有对参数进行检验的过程，为了能够较好的解决这个问题,为了使我国的 GPS 用户没有后顾之忧，提出了这个课题。1.2本课题的研究任务研究本课题的任务是基于七参数法原理，利用相应的数学知识（主要使用到了最小二乘法）求解出由 WGS-84向 BJ-54坐标系转化的七参数值，再计算相应的 BJ-65 坐标系的坐标，最后通过高斯正算公式完成高斯投影坐标的计算。通过本次毕业设计基本完成了七参数的求解工作，所得出的结果可用于不同区域内，精度要求不太苛刻的、不同空间直角坐标系之间的坐标值的转化，可以减少已知点的数目，降低使用条件和计算成本，以便于在国内推广应用。第二章 GPS 概述2.1全球定位系统简介GPS 系统，即全球定位系统，是由美国陆海空三军联合研制的一种新的卫星导航系统。它的结构包括：1.GPS 卫星星座。由分布在 6个轨道平面内的 24颗 GPS 卫星组成。2.地面监控系统。包扌舌一个主控站、三个注入站和五个监控站及其它的通信系统和辅助系统,其主要功能是收集数据，编算导航电文，向用户发送广播星历及精密星历。3.GPS 信号接收机。它是用来接收、跟踪、变换和测量来自 GPS 卫星的信号，提供用户状态参数的一种电子设备。2.2GPS 卫星测量原理定位应用基本原理：GPS 接收机产生与卫星发射相同的伪随机码，由机内延迟锁相环使复制码与卫星输入码对准，求出信号传播时间，推算出卫星至接收机的距离。通过接收来自 4 颗或 4 颗以上卫星的信号，使用空间后方交会原理，求出接收机安置点的三维坐标和参数，达到定位的目的。测量应用是使用两台或两台以上的接收机在待测点上同步观测卫星。通过测定载波相位差，依据相位的线性组合（单差，双差，三差）和求出的整周模糊度，来测定两点或多点空间坐标的三维坐标增量。其精度可达 12PPM（相对定位精度）。2.3GPS 在工程上的应用GPS 是多功能、高效、快速和高精度的定位系统。该系统应用于工程测量是一种崭新的技术方法和手段。它具有测量精度高，观测时间短，数据管理快，成本低等优点。作为 GPST 程控制网，要求点位精度都要达到施工所需的精度，由于 WGS-84 椭球人地高转换为正高或正常高还没有达到工程要求，故 GPS 控制网主要用于平面工程控制网。第三章坐标系统3.1WGS-84坐标系WGS-84 坐标系是美国国防部研制确定的人地坐标系，是一种协议地球坐标系。WGS-84 坐标系的几何定义是：原点是地球的质心，空间直角坐标系的 Z 轴指向 EIH（1984o0）定义的地极（CTP）方向，即国际协议原点 CIO,它由 IAU 和 IUGG 共同推荐。X 轴指向 BIH 定义的零度子午面和 CTP 赤道的交点，Y 轴和乙 X 轴构成右手坐标系。WGS-84 坐标系是对 NSWC9Z-2（NNSS 卫星多普勒定位系统的一个参考坐标系）的修正。WGS-84 椭球采用国际人地测量与地球物理联合会第 17届大会测量常数推荐值，釆用的两个常用基本几何参数：A=6378137m;fM:298.257223563.3.2国家大地坐标系我国目前常用两个国家大地坐标系，即 1980年国家大地坐标系（C80）和 1954年北京坐标系（旧 EJ-54）。两者都是参心大地坐标系统3.2.11954年北京坐标系旧 BJ-54系是建国初期我国根据苏联 1942 年普尔科夫坐标系的起算数据传算过来的坐标系，它是普尔科夫系的延伸。归结其要点为：1.北京坐标系的参考椭球是属于克拉索夫斯基椭球常用基本参数如下：长半轴 a=6378245m2m扁率 f=l:298.32.lidBJ-54坐标系的人地原点在前苏联的普尔科夫；3.采用多点定位进行了椭球定位；（1）.高程基准为 1956 年青岛验潮站求出的黄海平均海水面；（2）.高程异常以前苏联 1955 年人地水准面重新平差结果为起算数据，按我国天文水准路线推算而得。3.2.21980年西安大地坐标系1978年在西安召开的全国天文人地网整体平差会议上，专家们建议建立我国新的人地坐标系，以改正旧 BJ-54坐标系的诸多缺点，1980 年完成天文人地网平差，故称为 1980年国家大地坐标系（简称 C80）。归结 1980年大地坐标系的要点为：1.该坐标系人地原点定在我国中部，即陕西省泾阳县永乐镇。2.采用 IUA 和 IUGG1975年推荐的地球椭球参数：长半轴 a=6378140m 扁率 fH:298.2573.定向明确：1980年国家人地坐标系的椭球短轴平行于地球质心指向地极原点 JYD（1968o0）方向，起始大地子午面平行于格林尼治平均天文台的子午面：3.椭球定位参数以我国范閑内高程异常值平方和等于最小为条件求定；4.大地点高程以 1956年青岛验潮站求出的黄海平均海水面为基准；5.建立两套 1980 年国家人地坐标系和地心坐标系。前者是在后者的基础上通过精确求定位坐标变化参数，换算成地心坐标。3.2.3新 1954年北京坐标系（新 BJ-54）新 1954 年北京坐标系（整体平差转换值）是由 1980新国家人地坐标系派生得来的，是作为老54系到 80系的过渡而存在的。将其与 80系和旧 54系比较，归结其要点为：1.采用克拉索夫斯基椭球参数：长半轴 a=6378245m扁率 f=l:298.32.多点定位。参心虽和老 54系参心不一致，但十分接近。3定向明确。坐标轴方向和起始人地子午面与 80系相同。坐标轴的旋参数等于零。4.人地原点位于陕西省泾阳县永乐镇，但与 80系人地原点大地起算数不同。5.大地点高程以 1956年青岛验潮站求出的黄海平均海水面为基准。6.提供的坐标是 1980年国家人地坐标系整体平差转换值，坐标精度和 80系的坐标精度完全一样。新老 54系不存在椭球差异和定位差异，两系统同一点坐标的不同主要是由于一个是全国统一平差的结果，另一个是局部平差结果的缘故产生的。3.3不同空间大地直角坐标系的转换一般包括地心直角坐标系与参心直角坐标系转换及不同参心直角坐标系的转换。根据转换参数的不同转换，模型有很多种，如三参数法，七参数法，九参数法，十参数法等。下面简介三参数法和七参数法。1三参数法两个坐标系的坐标轴平行，仅有原点不同，也就说只有平移参数，采用三参数法。模型如下：1.七参数坐标转换模型两个坐标系除平移参数外，坐标轴还存在旋转参数及尺度比参数，共七个参数，转换模型有布尔莎、莫洛金斯基及范士公式。本文简介布尔莎公式，其余两种参见有关文献。模型如下:其中，m 为尺度比参数，为旋转参数，为平移参数。第四章 GPS 工程网数据后处理的过程和内容4.1后序数据处理简介 GPSGPS 相对定位前序数据处理方法解求出两个端点的人地坐标差（dE,dL.dH）或属于空间直角坐标系的坐标差（dX.dY,dZ）,般定义被观测的两端点的边线为 GPS 的观测基线，所得到的坐标差即为相应基线的基线解。同时，GPS 绝对定位的精度很低，不能在工程测量中加以应用，GPS 观测网是一个局部独立的自由网，其网点的坐标基准和方位基准相对于实际存在的各种国家统一的坐标系或局部坐标系都是未知的。GPS 相对定位的基线解由 GPS 接收机的随机软件或第二家软件开发商提供的基线解算软件解求。这以后的全部工作便属 GPS 观测后序数据处理的领域，即 GPS 测量后序数据处理就是对 GPS 基线结果具体应用。GPS 测量后序数据处理的目的即为：1.提取 GPS 网的基线结果；2.GPS 基线网的独立平差：3.GPS 基线网测量成果到地面坐标系的转换。卜面将按照计算步骤的先后，对 GPS 后序数据处理每个部分的内容予以简述。4.2GPS 后序数据处理过程1.挑选并汇总独立基线GPS 测量时，n台接收机的同步观测值可以两两结合，通过计算将得到多于(ml)条的基线。但是，其中只有(ml)条基线是独立的，进行平差计算的时候只能让(ml)条独立基线参与平差计算。2.提取三维基线解并组成三维基线结果文件在基线解算软件提供的基线解文件中，除了包含有两个点的点名、坐标差、测量的时段号，还记录了许多用作他途的其它信息。应该把分散在各个解算结果文件中的真正有用的信息提取出来，并按照较为合理的格式组织汇总到一个三维基线结果文件中去，作为输入数据提供给后序有关的计算，提取三维基线解就是提取后续各种数据处理需要的四个基本要素：a.基线两端点的点名；b.基线端点的三维近似坐标；c.基线向量及其方差、协方差；d.测量的时段号。3.GPS 网点三维坐标的概算对 GPS 网点进行三维坐标的概算，以某一网点的单点定位坐标为基准，根据解算出来的坐标差，推倒出其它两点的坐标，从而为 GPS 网三维平差提供相对关系较为精确的 GPS 网点三维近似坐标。4.三维基线投影到高斯平面上将三维基线解的投影分解为基线两个端点的投影。利用高斯投影变换的公式将一个端点投影到高斯平面上，另一个端点用第一个端点坐标加上基线解向量作为三维坐标进行投影。然后重组基线解。此外还存在投影高度的问题，简易作法就是将标准椭球的长半径加上投影面高度作为与右影高程面相吻合的投影椭球体的长半径。5.GPS 网点二维网坐标概算经过坐标概算之后的 GPS 网点三维近似坐标直接投影到高斯平面上，所得到的 GPS 网点的二维近似坐标之间相对关系同样较为精确。但是为了与地面局部坐标系相联系，还应该利用一个或一个以上的地面已知点将 GPS 网平移，旋转到与地面已知点较为接近的位置。6.GPS 三维及二维网平差相对与实际存在的各种地面坐标系，GPS 测量得到的基线网的位置和方位基准都是未知的，该网与地面局部坐标系之间存在平移，旋转和比例三种系统转换参数；其次，观测网中总不可避免存在着测量误差。因此，必须利用足够的地面已知点对 GPS 网进行平差或拟合，把 GPS 测量结果转换到地面局部坐标系中。4.3GPS 数据后序处理的结果经过 GPS 数据后序处理，将 GPS 测量成果(三维基线解)转换到局部坐标系中，得到我们需要的点的地面局部坐标系的坐标及所需边长。第五章 WGS84BJ54 两种空间直角坐标系的转换目前，GPS 被大量地应用到旧城市控制网的改造和扩展中，工程控制网的建立，测屋成果一般都要同时提供 54系和80系高斯平面坐标以及独立坐标，而 GPS 测量成果属于 WGS-84 系的，因此将 WGS-84 系转化为54系或 80系就显得非常必要和迫切，下面就试着讲座一下由 WGS-84 系到 54系高斯平面坐标的转换过程与方法；并试着进行了计算程序的编制。5.1数学模型的建立WGS84椭球与 BJ54坐标系所属的克拉索夫斯基椭球有差异，因此要将 WGS84系空间直角坐标系转化到 54系高斯平面坐标，首先得完成 WGS-84 椭球到克拉索夫斯基椭球的转代，数学模型采用布尔莎公式：对该公式精加变换解算这七个参数，至少要用到三个已知点，采用间接平差模型进行解算：其中：V 为残差矩阵；X 为未知七参数；A 为系数矩阵；L 为闭合差解之：解得七参数，每输入一坐标值，就能求出它在新坐标系中的坐标。有时在转换精度要求不太高时，也可采用五参数布尔模型（去掉旋转参数）这样，模型就为:解算时只需两个已知点即可，解算方法与七参数法完全相同。5.2由空间直角坐标到大地坐标的转换5.1节只求解出点在克拉索夫斯基椭球的坐标，它还不能作为工程网的依据：要将它转化为高斯平面坐标必须还要将其先从空间直角坐标转化为大地坐标：由控制测量学我们知道：显然很容易知道：但人地纬度 B的计算比较复杂，通常采用迭代法：如图所示：9由图可知:上式右端有待定量 E,需迭代计算，迭代时可取用 E 的初值 B1 计算 N1 和 smB1,将上式进行第二次迭代直至最后两次 B 值之差小于允许误差为止。计算出了 E 值，大地高也可得出5.3高斯平面坐标的计算得到了点的人地坐标，就可以将其转化为某投影带的高斯坐标：因为BT-54 属于克拉索夫斯基椭球，我们将克拉索夫斯基椭球参数代入高斯投影正算公式：得到更适用于电算的高斯坐标计算的实用公式：式中它们的计算精度，即平面坐标可达 0.001111,注意：椭球不同的不能应用这个实用电算公式。5.4程序的编制程序框图：1坐标转换总程序图：2坐标转换框图:3组成误差方程框图:第六章数据处理结果的比较本章对提供的 GPS 三维无约束平差结果进行了两种方法的数据处理：第一种是布尔莎七参数坐标转换；第二种是布尔莎五参数坐标转换。它们均是从 WGS-84人地坐标系转化到 BJ-54高斯投影坐标。（一）已知数据1.已知 WGS-84人地坐标系坐标2.已知公共点坐标3.已知平差后高斯投影坐标（二）坐标转换1.转换参数比较七参数五参数X-79.012056-76.165648Y104.541645110.229246Z19.39996913.012207111-0.000001-0.000001ex0.319096ey-0.081277ez0.3194050.2849462.转换结果比较转换结果的比较可参考实例的附三。（三）总结由通过实例的计算町知，布尔莎七参数和五参数公式求出来的转化结呆在小范 I 制面积内近似相同，七参数公式适用范围人，而五参数公式虽然适用范围小，但由于公式更简洁，理所当然地更适合用于小范围的坐标转化。