坐标转换说明
GPS接收机接收到GPS(大地坐标:经度、纬度和高度值)信号后,并不利于显示,需要将大地坐标进行转换,现选用东北天坐标系(也叫站心坐标系)作为显示的依据。
GPS接收机接收到的第一个信号L(经度)、B(纬度)和H(高度),作为东北天坐标系的原点。当接收到第二个信号时L
1、B1和H1,应用坐标转换公式,转换到东北天坐标系下进行显示。依次类推,凡是接收到的GPS信号都转换到东北天坐标系下进行显示,在东北天坐标系下预测出来的坐标值通过坐标转换公式在显示屏上显示大地坐标(经度、纬度和高度)。
1.大地坐标与直角坐标的相互转化
对空间某一点,大地坐标系(L,B,H)到直角坐标系(X,Y,Z)的转换关系如下:
X(NH)cosBcosLY(NH)cosBsinL
(1) Z[N(1e2)H]sinB由直角坐标系(X,Y,Z)转化到大地坐标系(L,B,H)的公式如下:
Larctan(Y/X)Barctan{Z(NH)/[X2Y2(N(1e2)H]}
(2)
HZ/sinBN(1e2)式中:Na/1e2sin2B,N为该点的卯酉圈曲率半径;e2(a2b2)/a2,a、b、e分别为该大地坐标系对应参考椭球的长半轴、短半轴和第一偏心率。长半
9013。 轴a=63781372m,短半轴b=6356.7523142km,e20.0066943799从公式(2)看出,经度比较容易求得,纬度和高度必须通过迭代计算获直接计算得到。迭代计算的次序为:NHB,通常迭代四次可以达到H优于0.001m,B优于0.00001的计算精度;教科书中给出的直接法计算公式比较繁琐,有的计算公式的应用条件受到一定限制,例如要求大地高度小于10000m时,才能使B、H达到上述计算精度,有的直接计算公式精度较低。
根据[张华海]提供的方法,本文建议采用该方法将直角坐标(X,Y,Z)转变成大地坐标(L,B,H)。该方法的公式形式比较简便,B、H的计算精度高;用计算出的具有一定精度的B0,直接求出H,一次性计算出满足精度要求的H;再将H值代入公式(2)中,求出B值。
令uarctan(Za/(X2Y2b)),a、b分别为长半轴和短半轴。将u代入下式,求出B0:
B0arctan((Zbe2sin3u)/(X2Y2ae2cos3u))
Na/1e2sin2B0
Larctan(Y/X);HX2Y2(ZNe2sinB0)2N;Barctan{Z(NH)/[X2Y2(N(1e2)H)]}其中:e为第二偏心率,e20.00673949674227。通过上市就可以得到精度较高的大地坐标(LBH)。
2.直角坐标与东北天坐标的相互转化
以GPS接收到的第一点作为东北天坐标系的原点,以通过坐标原点且指向天顶的法线为z轴(指向天顶为正),以子午线方向为y轴(向北为正),x轴指向东,且与y、z轴垂直(向东为正)。
(NH)cosBcosLX1sinBcosLsinLcosBcosLyY=sinBsinLcosLcosBsinLx+(NH)cosBsinL 120sinBZ1直角cosBz东北天[N(1e)H]sinB
(3)
设坐标换矩阵表示为:
sinBcosLsinLcosBcosL
(4) RsinBsinLcosLcosBsinL0sinBcosB(3)式可化简为:
(NH)cosBcosLX1yY=Rx+(NH)cosBsinL1
(5) 2Z1直角z东北天[N(1e)H]sinBˆXX(NH)cosBcosL1ˆ令YY1(NH)cosBsinL,则可得到东北天坐标:
2ˆZZ1直角[N(1e)H]sinBˆXyx1ˆ
(6) =RYˆz东北天Z其中: R1sinBcosLsinBsinLcosB
(7) sinLcosL0cosBsinLsinBcosBcosL将公式(6)展开得到东北天坐标:
ˆcosLYˆxsinLXˆsinLsinBYˆcosBZˆ
(8) ysinBcosLXˆˆˆzcosBcosLXcosBsinLYsinBZ备注:
纬度1度合:110.94km 纬度1分合:1.84km 纬度1秒合:30.8m 某一纬度下的经度1度合:纬度1度×cos(纬度)
参考文献:
1 董绪荣,张守信,华仲春.GPS/INS组合导航定位及其应用.北京:国防科技大学出版社,1998 2 胡伍生,高成发.GPS测量原理及其应用.北京:人民交通出版社,2004 3 张华海,郑南山,王军,李景芝.由空间直角坐标计算大地坐标的简便公式.全球定位系统.2002,4:9-12