精选-GPS坐标和国家大地坐标之间的转换

1、GPS标和国家大地坐标之间的转换一、前言WGS-84标系是目前GP新采用的坐标系统，GP新发布的星历参数就是基 于此坐标系统的。 WGS-84标系统的全称是 World Geodical System-84 (世界 大地坐标系-84),它是一个地心地固坐标系统。WGS-84:标系统由美国国防部 制图局建立，于1987年取代了当时GPST采用的坐标系统-WGS-7/标系统而成 为GPS勺所使用的坐标系统。WGS-84标系的坐标原点位于地球的质心，轴指向 BIH1984.0定义的协议地球极方向，轴指向BIH1984.0的启始子午面和赤道的交 点。采用椭球参数为：a=6 378 137m, f= 1

2、/298.257 223 563 。北京54坐标系、西安80坐标系一属于参心坐标系，北京54坐标系采用克 拉索夫斯基椭球参数，长轴a= 6 3 78 2 4 5 米，扁率f=l : 2 98.3;西安80大地系坐标系椭球参数采用国际大 =地测量和地球物理联合19 75后推荐的地 球椭球参数，长轴a= 6 3 7 8 140 米，扁率f1 : 298.257, 大地原点在我西安 市径阳县永乐镇。西安 80坐标系的建立是在54年北京坐标系的基础上完成 的。在实际的工彳中，对于GPS勺测量数据。我们需要将其转换成所需要的54或80坐标系，才能够使用。或是将其转换成相应的地方坐标系。在转换的过程中 需

3、要进行一系列的变换。本文将对其过程做详细的说明。二、转换过程(1)数据测量：在实际操作中，首先进行的是数据的观测。根据实际工作需 要，采用相应的观测方法进行观测，得到合格的测量成果。本文主要是针 对GP毓制网的转换来说明的。(2)平差：在GPSE制网的测量工程中，在进行完基线测量(地面坐标和高程)后，需要对测量结果进行平差，得到相应的平差结果。下面对相应的条件 平差做具体说明：AV-W=0 1L#=L+V 2基础方程和它的解：设有r个平差线性条件方程：/ alLl + a2L2 + an£n + aO = 0bl£l+ b2L2+ .1+ hnLn + bO = 0k rl

4、Ll + r2L2 + .»+ rnLn + rO = 03式中ai, bm（i=1,2,口）为条件方程系数，ao,boro为条件方程常数项。将方 程带入2中，得条件方程为：f alvl + a2v2 + + anvn + wa = 0；blvl + b2v2+ +bnvn + wb = 0 囿I rlvl + r2v2 !F rnvn + wr = 0式中wa,wb-w成为条件方程的闭合差，即fwa = -（alLl + 我 L2 + 皿 + anLn + a0）囱尸直纣金翳L.也5wr = -(rl£l+ r2£2 + .，+ min + r0)A=Lrl r

5、2按求函数极值求拉格朗日乘数法:求导V=PAQA0-ATPV-2K(AV-W)d/dv=2VTP-2KTA=0-1>Jk=qAkK-W=06TTPal a2 - an bl b2 bnN3a=AQA=AP1ATi=1/p(a ika+bikb+mkr)求评差值L #=L+V（3）将得到的平差结果按以下具体步骤进行解算：1 ）将空间直角坐标系转换成大地坐标系，得到大地坐标（ B, L）:L=arctan （X /Y ）B=arctan （Z + Ne2sin B ） / （犬 + Y 2） 0.5H= ( X2+ Y2) 0.5sin - N8用上式采用迭代法求出大地坐标（B, L ）2

6、）将大地坐标系转换成高斯坐标系，得到高斯坐标（X , Y ）按高斯投影 的方法求得高斯坐标。7r皿氏族的-5蜡+广泗3）将高斯坐标系转换成任意独立坐标系，得到独立坐标（ x ' , y '）。在小范围内测量，我们可以将地面当作平面，用简单的旋转、平移便可将高斯坐标换成工程中所采用坐标系的坐标（X、Y '）。以下是具体的不同大地坐标系之间的转换方法一一布尔莎公式坐标系之间的转换：由于不同的大地坐标系之间换算除具有不同空间直角 坐标系统之间换算的七个转换参数，还要增加两个系统采用的地球椭球不同而 产生的两个地球椭球转换参数。不同大地坐标系的换算公式又称大地坐标微分 公式或变

7、换椭球公式。一般采用的解算方法是布尔莎公式：Y,Z.10(N + H) cosB cosL(N + H) cdsB sin LpV(l - e2) + H sin B顾及N=a/B11所以X,Y,Z 各自是B, L, H, a , e2的函数。因此这些自变量的变化 dB、dL、dH da 必引起 X、Y、Z 的变化 dX、dY、dZ。式中 A=121311中可按8右边相对应变量求偏导得到。由于S = -MsinBJ当=McosB114由8得至IJ:ss15/一现在求偏导数Sada 3adz dK=(1 e) sin B=cos B cos La广 a广一 =cos B sin L M J?(1

8、 J) N sin BN=a(1-e2sin 2B)-1/216对其进行进一步的换算，一 (M + H) sin E 8s L A= 一(M + H)sinBsinLC=(M + H) cos B(N/a) cos B cos L (N/a)cosBsinL3N _ Nda a3N _ NEEinBBfffi? 2(l-e2EiD2B)得到具体的表达式：(M - H) cos B sin L(N - H) cos E cos L0:M T 1 一 a i5ii:JEccsLM/(l a)sin2BcosL17cos B cos L cds B sin L sin B18CN/a)(l - a)

9、2sinB (1 - a)smB(Msm3B - 2N)求A1可得sin B cos Li=(M + H) sinB sin L 詈=(M + H) cos B(N + H) cos B fin L今=(M + H) cos B cos L生=0sldxo . T=cosB sin L dHJr口 m t=cos E sin L0H=sin B 日H日N=-cos B cos LdN" cos B cos L19r-sinBrosL/(M + H)- sinBsinL/(M+H) cosB/(M + H)1 A-1-A=-secBsmL/(M + H) secBcosL(M+N)0c

10、osBcosLcosBsinLsinB将以上各式整理，最终得到：20dBdL.dH.sinHcacL 尸 siaE mnLeid£ 旌PP PM4H 尸3£4H rM+H reoeL ffcoeL n 0tM+H)siuB(K+H)cobBcosBcosL cosBsinLsinBsinLcosttanBcosLtanBsinL.NsinBcosBsinL/p” Ne2sinBcosBcosL/pNe2sinBcosBp "/M0N(l-essinaB)Dk+M+He2sinBcosBp0-(1 -e"sin2B') 且M(2-eaBinaB)_ _-smB