常用坐标系之间的关系与转换VIP

1、sBB 7-23Z=N (I护sinA7.5常用坐标系之间的关系与转换一、大地坐标系和空间大地直角坐标系及其关系大地坐标系用大地纬度企丈地经度L和丈地高H来表示点的位置。这种坐标系是经典大地 测量甬:画用座标紊7展据地图投影的理论，大地坐标系可以通过一定的投影转化为投影平面上 的直角坐标系，为地形测图和工程测量提供控制基础。同时，这种坐标系还是研究地球形状和 大小的 种有用坐标系。所以大地坐标系在大地测量中始终有着重要 的作用.空间大地直角坐标系是一种以地球质心为原点购亘跑®坐标系，一般用X、化Z表 示点 BSSTSTT逐碇SS靶菇飞雨H绕祺扭转冻其轨道平面随时通过 地球质心。对它们

2、的跟踪观测 也以地球质心为坐标原点，所以空间大地直角坐标系是卫星大地测量中一种常用的基本坐标系。 现今，利用卫星大地测量的手段*可以迅速地测定点 的空间大地直角坐拯，广泛应用于导航定位等空间技术。同时经过数学变换，还可求出点的 大地坐标I用以加强和扩展地面大地网，进行岛屿和洲际联测，使传统的大地测量方法发生了 深刻的变化，所以空间大地宜角坐标系对现今大地测量的发展具有重要的意义。如图7-23所示尸点的位置用空间大地直角坐标X,匕Z)表示，其相应的大地坐 标为(E,勿将该图与图？ 一 5 加以比较可见，图7-5中的子午椭圆平面相 当于图7-23中的OJVP平面.其中PPz=Z.相 当于图7-5中

3、的j7； OP3相当丫于图7-5中的 小两平面的经度乙可视为相同，等于叽于是可以直接写出X=jrcQsif Y=jrsinL, Z=y将式(7-21).式(7-20)分别代入上式，井考虑式(7-26)得X=NcosAcosZr ”Y=NcQsBsinL >(778)上式表明了 2种基本坐标系之间的关系。、大地坐标系和空间大地直角坐标系的转换1 .由大地坐标求空间大地直角坐标当已知椭球面上任一点P的大地坐标(民勿时，可以按式(7-78)直接求该点的空间大地 直角坐标(X, Y, Z)o如果P点不恰好位于椭球面，例如位于大地高为H的H点处，此时由大地坐标求空间 大地直角坐标的公式则为(N+H

4、) cosBcosL (N 十 H) cosBsinL 、(7-79)N(le2)+刃sinBj2 .由空间大地直角坐标求大地坐标当已知X、y、Z反求B、L、H时，可以采用直接解法或迭代解法。由公式(7-79)第、二两式得X 二y 二L 二 arc tan §(7-80)利用上式可直接由空间大地直角坐标X、Y求出大地经度公。为了求出B和H,还应对公式作些变化，以适应迭代计算的需要。由公式(7-79)第一、 三式得(N 十 H) ZcosL*a” =(l_e2)式中，cos ZA仍由式(7-79)得出cos, 一 (n+h)cosB yy2_l_y2代入前式B 二 arctan2(7-

5、81)V I又由式(7一79)得H =7、一 N(7-82)cosB式(7-81).式(7-82)就是求3、H的迭代公式。迭代开始时设Nq = QHo=+ W+Z2_ 应Bo= arc tan2VX2aY2随后，每次迭代按下列公式进行V 1, sin 叨-cosB 宜至b-5和Hi-H小于要求的限值为止。一般，在要求H精确至0.001m、占精确 至 0.000。/时，需要迭代4次。三、不同空间大地直角坐标系的换算利用“GPS”定位所获取的点位属于空间大地直弟坐标系。可是由于各国所采用的参考椭球及其定位不同，参考椭球中心也不和地球质心重合，所以世界上存在着各不相同的空间大地直角坐标系。为了将“G

6、PS”定位成果转换成各自需用的成果，就出现了不同空间大地直角坐标系的换算。这在“GPS”定位的数据处理中，应用十分广泛.在高等数学的解析几何里，曾经论证了二维直角坐标系中，当坐标轴旋转角度。时(图7-24),用旧系坐标表示新系坐标的公式为.>(7-33)弄二一日 sma-FAia costfj在三维空间直角坐标系中，新、旧两坐标系的变 换需要在3个坐标平面上，分别通过3次转触才能完 成。如图7-25所示.2个空间大地直角坐标系和0 一心Yin Ze，它们的原点一致，但相应的坐标轴 互不平行，存在微小差异。按以下步骤进行转轴可以将。一心丫旧乙日转换成o-x新第一*保持0Z时轴不动，绕其将O

7、X” ox日轴金定转微小角度殳，旋转后的坐标轴设为图7OX OY OZ 则有X* XA cose.+Yjg sinefY1 = 一XIH sin®十 F 旧 cosesZTB第二，保持OF轴不动，绕其将OZJOX,轴旋转微小角度旋转后的坐标轴设为OX”、 OYOZ则有X”二 X, cos®Z'sin八Yr=r>(7-85)Z" 二X'sin£Y+Z'cos&Y,第三，保持OX”轴不动，绕其将。严、OZ”轴旋转微小角度匕，旋转后的坐标轴设为OX*、 0Y祈、OZ觥则有X八=X *=y”coss+Z”sin£x

8、 (7-86)二 一 y"sinEx+Z"cos&x这样，将O X旧Y旧Z旧分别绕3个坐标轴旋转了 3个微小角度£z、弘5,使其和 O天新丫新Z新重合。£x、£丫、£z称为欧勒角。将式(784)代入式(785),再代入式(7 86),由于取、£丫、£乙是秒级微小量，略去其 正弦、余弦函数展开式中2次及以上各项，得Xgj 二Xg+£疔旧 r£yZ 日丫新=丫旧旧+e：xZ日2新二2旧十£丫 乂旧一 5 丫旧当新、旧2个坐标系的原点不相一致时，还需根据坐标轴的平移原理，将旧系原点移

9、 至新系原点，其变化公式为X新=Xo+XI日+£疔口 £ 丫乙日(7-87)Z=Zo + Z用+£yX用一包丫旧 J式中，X。、y。、Z。称为3个平移参数，是旧坐标系原点在新坐标系中的3个坐标分若再考虑两个坐标系的尺度比例也不一致，即存在有尺度变化的参数，设为虹则有 Xgf=Xo+ (1+上)Xia+®Y 旧一£yZ(,日丫新二丫。+ (1+”)丫旧一叨八十谄旧(7-88)Zjfi=Zo+ (1+&) Zig+Cy"S-ExYlo .上式即为布尔莎公式。公式中存在7个参数：3个平移参数X。、丫。和Z。，3个旋 转参数昭为、切1

10、个尺度变化参数虹习惯上称这种换算法为七参数法。七参数法除布尔莎 公式外，还有莫洛琴斯基公式和范士公式等。由公式(7-88)可知，由一个坐标系换算成另一个坐标系，必须知道其转换参数。转换 参数可以通过联测一些公共点获得，因为通过公共点联测，可以得到这些公共点在新、旧 2个坐标系中的坐标值，于是就可以利用公式(7-88)求出转换参数。当公共点数较多时，观 测方程式个数就大丁所求参数个数，这时还可根据测董平差原理列立观测值的误差方程式， 组成并解算法方程，求得转换参数。四、不同大地坐标系的换算地面点在椭球面上的位置，是由一定元素和定位的椭球所规定的.如果选择的椭球元 素和 定位发生变化，地面点在椭球

11、面上的大地坐标必将随之变化根据椭球元素和定位的变化推求 点的大地经纬度和大地高变化的公式，叫做大地坐标微分公式，它是不同大地坐标换算的基 础，下面首先来推导大地坐标微分公式。由公式(7-79)可以看出，点的空间大地直角坐标是椭球几何元素(用长半径和扁率产表 示)和椭球定位元素(品、L、H)的函数口当椭球元素和定位结果发生了变化时，点的空间大地 直角坐标必然发生变化*取式(7-79)的全微分，即dX器也+就|/ +辍盼超舄心=和+副，+豹B得亚+各1H >(7-89)超=融+一+用曲+知+豹”. 考虑到 (JN d r,J> NUN 37Vde 3 f找、_L ".石二厉石

12、二厉卫(1 一代曲一刃石 -)2匚严历需二磊彳复 tl-ehin2B)3 =sinBcosB则根据式(7-79)可以求出3X二 AcosBcL= -cosBcosLoUzcosBcosL= yAjcosBcosLsin2BcosBcosL (N+H) sin£cosZ= (M+H) sinBcosL3XdL = (N+H) cosBsinLdXdH=cosBcosL将以上5式代入式(7-89)第1式得(M+彳/) ynBcosLdBdX二NcosBcos Za虫+McosBcos A sin?H 7A7一一(N+H) cosBsinLd£+cosBcosLd/l同理dy=N

13、cosBsinL 4-A/cosBsinLsin2B 7A7一 (M+H) sinBsin A dB>(7-90)a1-j+ (M+H) cosBcosLdL+cosBsinLdHdZ=N (1e2) sinB M (l+coszBA2sin2B) sinB 县：a- 1J+ (M+H) cosBdB+sinBdH若以dH、45、d Za为未知数解算以上3式，则得dH=cosBcosLdX+cosBsinLdY+sinBdZ一N (1一八2sin2B)乎+M (l-e2sin2B) sin2B曲二肚*芳一 sinBcosLdX一sinBsinLdy-4-cosBdZ> (7-91)

14、+WsinBcosB 乎+M (2-? siiM) sinBcosB 管dL- xt rj (secBsin£dX+secBcosLdy) h+H式中，血、”表示椭球元素(长半径、扁率)的变化；dx、1|又、dz表示楠球中心的变化，即椭球 定位的变化。因此，式(7-91)就是由于椭球元素和定位变化引起点的大地坐标变化的公式，亦 即大地坐标微分公式。将上式代入下式，即得不同大地坐标系的换算公式力新二A旧+d ABJBZdB (7-92)日新二H旧+dH：当考虑欧勒角和尺度变化参数时，可将式(7-88)写成如下形式dX=X新一/=&+宓旧+莎旧一翻旧dK=yM 一 丫旧二人+羽旧

15、rzXm +翻旧dZ=Zfr %=Zo+ Z 旧 +&X 旧一“Y旧上式等号右端的X旧、Z13用式(7-79)等号右端的函数代入后，再将上式代入式(791),经过整理可得广义大地坐标的微分公式，da>(7-93)(i/=cosBcosLXo+cosBsinjLya+sinBZoVe2sinBcosBsinZxx+N/sinBcosBcosLy+N (1一e2sin2B) kN (1一e2sin2B) +M (1 /sirfB) sin2BdB二(一sinBcos A X°一sinBsinZyo+cosBZo) 一sinLex M 十 H 4-cosLeye2sinBcosBA+/：N/sinBcosB 屯+M (2护 sin