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

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

6.2.1大地坐标系

点的子午面NPS与起始子午面NGS所构成的二面角，叫做点的大地经度，由起始子午面起算，向东为正，叫东经（0°～180°），向西为负，叫西经（0o～180°）。点的法线与赤道面的夹角，叫做点的大地纬度。由赤道面起算，向北为正，叫北纬（0°～90°）；向南为负，叫南纬(0°～90°)。

大地坐标系是用大地经度L、大地纬度B和大地高H表示地面点位的。过地面点P的子午面与起始子午面间的夹角叫P点的大地经度。由起始子午面起算，向东为正，叫东经（0°~180°），向西为负，叫西经（0°~-180°）。过P点的椭球法线与赤道面的夹角叫P点的大地纬度。由赤道面起算，向北为正，叫北纬（0°~90°），向南为负，叫南纬（0°~-90°）。从地面点P沿椭球法线到椭球面的距离叫大地高。大地坐标坐标系中，点的位置用,表示。如果点不在椭球面上，表示点的位置除,外，还要附加另一参数——大地高，它同正常高及正高有如下关系

6.2.2空间直角坐标系

以椭球体中心为原点，起始子午面与赤道面交线为轴，在赤道面上与轴正交的方向为轴，椭球体的旋转轴为轴，构成右手坐标系-，在该坐标系中，点的位置用表示。

地球空间直角坐标系的坐标原点位于地球质心（地心坐标系）或参考椭球中心（参心坐标系），z轴指向地球北极，x轴指向起始子午面与地球赤道的交点，y轴垂直于XOZ面并构成右手坐标系。

6.2.3子午面直角坐标系

设点的大地经度为，在过点的子午面上，以子午圈椭圆中心为原点，建立平面直角坐标系。在该坐标系中，点的位置用,表示。

6.2.4大地极坐标系

为椭球体面上任意一点，为过点的子午线，为连结的大地线长，为大地线在点的方位角。以为极点，为极轴，为极半径，为极角，这样就构成大地极坐标系。在该坐标系中点的位置用,表示。

椭球面上点的极坐标（,）与大地坐标(,）可以互相换算，这种换算叫做大地主题解算。

6.2.5各坐标系间的关系

椭球面上的点位可在各种坐标系中表示，由于所用坐标系不同，表现出来的坐标值也不同。

1.子午面直角坐标系同大地坐标系的关系

过点作法线，它与轴之夹角为，过点作子午圈的切线，它与轴的夹角为（90°+）。子午面直角坐标同大地纬度的关系式如下：

2.空间直角坐标系同子午面直角坐标系的关系

空间直角坐标系中的相当于子午平面直角坐标系中的，前者的相当于后者的，并且二者的经度相同。

3.空间直角坐标系同大地坐标系的关系

同一地面点在地球空间直角坐标系中的坐标和在大地坐标系中的坐标可用如下两组公式转换。

式中：e——子午椭圆第一偏心率，可由长短半径按式算得。

N——法线长度，可由式算得。

§6.3 几种主要的椭球公式
过椭球面上任意一点可作一条垂直于椭球面的法线，包含这条法线的平面叫做法截面，法截面同椭球面交线叫法截线（或法截弧）。包含椭球面一点的法线，可作无数多个法截面，相应有无数多个法截线。椭球面上的法截线曲率半径不同于球面上的法截线曲率半径都等于圆球的半径，而是不同方向的法截弧的曲率半径都不相同。

6.3.1子午圈曲率半径

子午椭圆的一部分上取一微分弧长，相应地有坐标增量，点是微分弧的曲率中心，于是线段及便是子午圈曲率半径。

任意平面曲线的曲率半径的定义公式为：

子午圈曲率半径公式为：

或

与纬度有关．它随的增大而增大，变化规律如下表所示：

		说 明
		在赤道上，小于赤道半径 此间随纬度的增大而增大 在极点上, 等于极点曲率半径

6.3.2卯酉圈曲率半径

过椭球面上一点的法线，可作无限个法截面，其中一个与该点子午面相垂直的法截面同椭球面相截形成的闭合的圈称为卯酉圈。在图中即为过点的卯酉圈。卯酉圈的曲率半径用表示。

为了推导的表达计算式，过点作以为中心的平行圈的切线，该切线位于垂直于子午面的平行圈平面内。因卯酉圈也垂直于子午面，故也是卯酉圈在点处的切线。即垂直于。所以是平行圈及卯酉圈在点处的公切线。

卯酉圈曲率半径可用下列两式表示：

6.3.3 任意法截弧的曲率半径

子午法截弧是南北方向，其方位角为0°或180°。卯酉法截弧是东西方向，其方位角为90°或270°。现在来讨论方位角为的任意法截弧的曲率半径的计算公式。

任意方向的法截弧的曲率半径的计算公式如下：

（7-87）

6.3.4 平均曲率半径

在实际际工程应用中，根据测量工作的精度要求，在一定范围内，把椭球面当成具有适当半径的球面。取过地面某点的所有方向的平均值来作为这个球体的半径是合适的。这个球面的半径——平均曲率半径R：

或

因此，椭球面上任意一点的平均曲率半径等于该点子午圈曲率半径和卯酉圈曲率半径的几何平均值。

6.3.5 子午线弧长计算公式

子午椭圆的一半，它的端点与极点相重合；而赤道又把子午线分成对称的两部分。

如图所示，取子午线上某微分弧，令点纬度为，点纬度为，点的子午圈曲率半径为，于是有：

从赤道开始到任意纬度的平行圈之间的弧长可由下列积分求出：

式中M可用下式表达：

其中：

经积分，进行整理后得子午线弧长计算式：

为求子午线上两个纬度及间的弧长，只需按上式分别算出相应的及，而后取差：，该即为所求的弧长。

克拉索夫斯基椭球子午线弧长计算公式：

1975年国际椭球子午线弧长计算公式：

6.3.6 底点纬度计算

在高斯投影反算时，已知高斯平面直角坐标（X，Y）反求其大地坐标（L，B）。首先X当作中央子午线上弧长，反求其纬度，此时的纬度称为底点纬度或垂直纬度。计算底点纬度的公式可以采用迭代解法和直接解法。

（1）迭代法

在克拉索夫斯基椭球上计算时，迭代开始时设

以后每次迭代按下式计算：

重复迭代直至为止。

在1975年国际椭球上计算时，也有类似公式。

（2）直接解法

1975年国际椭球：

克拉索夫斯基椭球：

6.3.7 大地线

椭球面上两点间的最短程曲线叫做大地线。在微分几何中，大地线（又称测地线）另有这样的定义：“大地线上每点的密切面（无限接近的三个点构成的平面）都包含该点的曲面法线”，亦即“大地线上各点的主法线与该点的曲面法线重合”。因曲面法线互不相交，故大地线是一条空间曲面曲线。

假如在椭球模型表面,两点之间，画出相对法截线如图所示，然后在,两点上各插定一个大头针，并紧贴着椭球面在大头针中间拉紧一条细橡皮筋，并设橡皮筋和椭球面之间没有摩擦力，则橡皮筋形成一条曲线，恰好位于相对法截线之间，这就是一条大地线。由于橡皮筋处于拉力之下，所以它实际上是两点间的最短线。

在椭球面上进行测量计算时，应当以两点间的大地线为依据。在地面上测得的方向、距离等，应当归算成相应大地线的方向、距离。