椭球面上的常用坐标系及其相互关系

1、6§ 6.2椭球面上的常用坐标系及其相互关系6.2.1大地坐标系P点的子午面NPS起始子午面NG啊构成的二面角L,叫做P点的大地经度，由起始子午面起算，向东为正， 叫东经（0°180° ），向西为负，叫西经（0°180° ）。P点的法线Pn与赤道面的夹角B ,叫做P点的大地纬度。 由赤道面起算，向北为正，叫北纬（ 0°90° ）;向南为 负，叫南纬（0。90° ）。大地坐标系是用大地经度 L、大地纬度B和大地高H 表示地面点位的。过地面点P的子午面与起始子午面间的 夹角叫P点的大地经度。由起始子午面起算，向东为正，

2、叫东经（0° 180° ）,向西为负，叫西经（0° -180° ）。 过P点的椭球法线与赤道面的夹角叫P点的大地纬度。由赤道面起算，向北为正，叫北纬（0° 90° ）,向南为 负，叫南纬（0° -90° ）。从地面点P沿椭球法线到椭 球面的距离叫大地高。大地坐标坐标系中，P点的位置用L, B表示。如果点不在椭球面上，表示点的位置除L, B外，还要附加另一参数一一大地高H ,它同正常局H 正常 及正高H正有如下关系H =H正常+匚（高程异常）H =H正+N（大地水准面差距）6.2.2空间直角坐标系以椭球体中心。为原点

3、，起始子午面与赤道面交 线为X轴，在赤道面上与 X轴正交的方向为丫轴，椭 球体的旋转轴为 Z轴，构成右手坐标系 O- XYZ ,在 该坐标系中，P点的位置用X,Y,Z表示。地球空间直角坐标系的坐标原点位于地球质心（地心坐标系）或参考椭球中心（参心坐标系），z轴指向地球北极，x轴指向起始子午面与地球赤道的交 点，y轴垂直于XOZ面并构成右手坐标系。6.2.3子午面直角坐标系设P点的大地经度为L ,在过P点的子午面上，以子 午圈椭圆中心为原点，建立 x,y平面直角坐标系。在该坐 标系中，P点的位置用L , x, y表示。6.2.4大地极坐标系M为椭球体面上任意一点，MN为过M点的子午线，S为连结M

4、P的大地线长，A为大地线在M点的方位角。以 M为极点，MN为极轴，S为极半径，A为极角，这样就构成 大地极坐标系。在该坐标系中 P点的位置用S , A表示。椭球面上点的极坐标 （S , A）与大地坐标（L ,B ）可以互相换算，这种换算叫做大地主 题解算。6.2.5各坐标系间的关系椭球面上的点位可在各种坐标系中表示，由于所用坐标系不同，表现出来的坐标值也不 同。1. 子午面直角坐标系同大地坐标系的关 系a cos B.1 e2sin2 BacosB过P点作法线Pn ,它与X轴之夹角为B , 过P点作子午圈的切线 TP ,它与X轴的夹角 为（90° + B ）。子午面直角坐标 x, y

5、同大地纬 度B的关系式如下：一 2、a(1 -e ) sin B a 2b sin By(1 -e ) sin B =, 1 -e2 sin 2 B WV2. 空间直角坐标系同子午面直角坐标系的关系空间直角坐标系中的 P2P相当于子午平面直角坐标系中的y ,前者的OP2相当于后者的x,并且二者的经度 L相同。X = x cos L |Y =xsin L *Z =y3. 空间直角坐标系同大地坐标系的关系同一地面点在地球空间直角坐标系中的坐标和在大地坐标系中的坐标可用如下两组公 式转换。x = (N + H CosB cosLz = N(1 -e2 *H 】sinB,.yL = arctanx2c

6、 x z + Ne sin B B = arctan/,寸 x2 +y2sin BH = Z -N(1 -e2 )式中：e子午椭圆第一偏心率，可由长短半径按式e2 = a2 一 b2 / a2 算得。N法线长度，可由式 N =a /力e2 sin2 B算得。§ 6.3几种主要的椭球公式过椭球面上任意一点可作一条垂直于椭球面的法线，包含这条法线的平面叫做法截面， 法截面同椭球面交线叫法截线（或法截弧）。包含椭球面一点的法线，可作无数多个法截面,相应有无数多个法截线。 椭球面上的法截线曲率半径不同于球面上的法截线曲率半径都等于 圆球的半径，而是不同方向的法截弧的曲率半径都不相同。6.3.

7、1子午圈曲率半径子午椭圆的一部分上取一微分弧长 DK =ds , 相应地有坐标增量dx ,点n是微分弧dS的曲率中 心，于是线段Dn及Kn便是子午圈曲率半径 M。任意平面曲线的曲率半径的定义公式为：m dSM =-dB子午圈曲率半径公式为：M =虹W3M =与或V3BM说明B =0 '2CMo H(1 q2)1(1 53在赤道上，M小于赤道半径a0 口 <B <90 =2 、.a(1 e ) <.M <c此间M随纬度的增大而增大B =90 'aM 90 _ T& w在极点上，M等于极点曲率半径cM与纬度B有关.它随B的增大而增大，变化规律如下表所

8、示:6.3.2卯酉圈曲率半径过椭球面上一点的法线，可作无限个法截 面，其中一个与该点子午面相垂直的法截面同椭 球面相截形成的闭合的圈称为卯酉圈。在图中 PEE'即为过P点的卯酉圈。卯酉圈的曲率半径 用N表示。为了推导N的表达计算式，过 P点作以O' 为中心的平行圈 PHK的切线PT ,该切线位于垂 直于子午面的平行圈平面内。因卯酉圈也垂直于 子午面，故PT也是卯酉圈在P点处的切线。即 PT垂直于Pn。所以PT是平行圈PHK及卯酉圈PEE '在P点处的公切线。卯酉圈曲率半径可用下列两式表示：NN6.3.3 任意法截弧的曲率半径子午法截弧是南北方向，其方位角为 酉法截弧是东

9、西方向，其方位角为90°或0° 或 180° 。卯270°。现在来讨论方位角为A的任意法截弧的曲率半径 Ra的计算公式。任意方向A的法截弧的曲率半径的计算公式如下：NNRa =22221 一 - cos A 1 - e cos B cos A(7-87)在实际际工程应用中，当半径的球面。取过地面某点的所有方向个球面的半径平均曲率半径R:在一定范围内，把椭球面当成具有适根据测量工作的精度要求，Ra的平均值来作为这个球体的半径是合适的。这c宁 V因此，椭球面上任意一点的平均曲率半径 半径N的几何平均值。a2、=w2 J e )R等于该点子午圈曲率半径 M和卯

10、酉圈曲率6.3.4 平均曲率半径6.3.5子午线弧长计算公式子午椭圆的一半，它的端点与极点相重合；而赤道又把子午线分成对称的两部分。如图所示，取子午线上某微分弧 PP '=dx，令P点纬度为B ,P点纬度为B+dB , P点的子午圈曲率半径为 M ,于是有： dx =MdB从赤道开始到任意纬度B的平行圈之间的弧长可由下列积分求出：BX = MdB0一式中M可用下式表达：M =a0 -a2 cos2B a4 cos4B-a6 cos6B a8 cos8Bm2 3535a。=m。+十一 m4十一m&十m8干 2816128m2m4 157a2 = + + m6 + m8223216

11、甘出.m.37其甲.a4 = + m6 + m881632m6 , m8a6 =+3216m8 a8 128经积分，进行整理后得子午线弧长计算式：a2a.a&agX =a°B -一sin 2B sin4B-sin6B sin8B 2468为求子午线上两个纬度B1及B2间的弧长，只需按上式分别算出相应的X1及X2 ,而后取差：AX =X2 X1 ,该以即为所求的弧长。克拉索夫斯基椭球子午线弧长计算公式：X =111134.861B 16036.480sin 2B 16.828sin 4B0.022sin 6B_ _ .3_.5 _X =111134.861 B _32005.7

12、80sin BcosB _133.929sin3 B cosB-0.697 sin5 BcosB1975年国际椭球子午线弧长计算公式：X =111133.005B 16038.528sin 2B 16.833sin 4B0.022 sin 6B七3 _5 -X =111133.005B -32009.858 sin BcosB-133.960 sin B cosB-0.698 sin BcosB6.3.6底点纬度计算在高斯投影反算时，已知高斯平面直角坐标(X, Y)反求其大地坐标(L, B)。首先X当作中央子午线上弧长，反求其纬度，此时的纬度称为底点纬度或垂直纬度。计算底点纬度的公式可以采用迭

13、代解法和直接解法。(1)迭代法在克拉索夫斯基椭球上计算时，迭代开始时设_ 1B f = X /111134 .8611以后每次迭代按下式计算:Bf1 =(X -F(Bf )/111134.8611F(Bf) -16036.4803sin2Bf 16.8281sin 4Bf 0.0220sin6Bf重复迭代直至B；* -Bf为止。在1975年国际椭球上计算时，也有类似公式。(2)直接解法 1975年国际椭球：:=X / 6367452.133Bf =B (50228976 293697 (2383 22coS -)cos' - cos' 10° sin ： cos -克拉索夫斯基椭球：|.：- : X / 6367588.4969Bf = "(50221746 293622 (2350 22cos2 -) cos2 |： cos2 | -6.3.7 大地线椭球面上两点间的最短程曲线叫做大地线。在微分几何中，大地线 这样的定义："大地线上每点的密切面（无限接近的三个点构成的平面） 法线”，亦即“大地线上各点的主法线与该点的曲面法线重合”。因曲面法线互不相交，故大地线是一条空间曲面曲线。假如在椭球模型表面 A,