测量中几种主要椭球公式.doc

几种主要的椭球公式 过椭球面上任意一点可作一条垂直于椭球面的法线，包含这条法线的平面叫做法截面，法截面同椭球面交线叫法截线（或法截弧）。包含椭球面一点的法线，可作无数多个法截面，相应有无数多个法截线。椭球面上的法截线曲率半径不同于球面上的法截线曲率半径都等于圆球的半径，而是不同方向的法截弧的曲率半径都不相同。 子午圈曲率半径 子午椭圆的一部分上取一微分弧长，相应地有坐标增量，点是微分弧的曲率中心，于是线段及便是子午圈曲率半径。 任意平面曲线的曲率半径的定义公式为： 子午圈曲率半径公式为： 或 与纬度有关．它随的增大而增大，变化规律如下表所示： 说 明 在赤道上，小于赤道半径 此间随纬度的增大而增大 在极点上, 等于极点曲率半径 卯酉圈曲率半径 过椭球面上一点的法线，可作无限个法截面，其中一个与该点子午面相垂直的法截面同椭球面相截形成的闭合的圈称为卯酉圈。在图中即为过点的卯酉圈。卯酉圈的曲率半径用表示。 为了推导的表达计算式，过点作以为中心的平行圈的切线，该切线位于垂直于子午面的平行圈平面内。因卯酉圈也垂直于子午面，故也是卯酉圈在点处的切线。即垂直于。所以是平行圈及卯酉圈在点处的公切线。 卯酉圈曲率半径可用下列两式表示： 任意法截弧的曲率半径 子午法截弧是南北方向，其方位角为0°或180°。卯酉法截弧是东西方向，其方位角为90°或270°。现在来讨论方位角为的任意法截弧的曲率半径的计算公式。 任意方向的法截弧的曲率半径的计算公式如下： （7-87） 平均曲率半径 在实际际工程应用中，根据测量工作的精度要求，在一定范围内，把椭球面当成具有适当半径的球面。取过地面某点的所有方向的平均值来作为这个球体的半径是合适的。这个球面的半径——平均曲率半径R： 或 因此，椭球面上任意一点的平均曲率半径等于该点子午圈曲率半径和卯酉圈曲率半径的几何平均值。 子午线弧长计算公式 子午椭圆的一半，它的端点与极点相重合；而赤道又把子午线分成对称的两部分。 如图所示，取子午线上某微分弧，令点纬度为，点纬度为，点的子午圈曲率半径为，于是有： 从赤道开始到任意纬度的平行圈之间的弧长可由下列积分求出： 式中M可用下式表达： 其中： 经积分，进行整理后得子午线弧长计算式： 为求子午线上两个纬度及间的弧长，只需按上式分别算出相应的及，而后取差：，该即为所求的弧长。 克拉索夫斯基椭球子午线弧长计算公式： 1975年国际椭球子午线弧长计算公式： 底点纬度计算 在高斯投影反算时，已知高斯平面直角坐标（X，Y）反求其大地坐标（L，B）。首先X当作中央子午线上弧长，反求其纬度，此时的纬度称为底点纬度或垂直纬度。计算底点纬度的公式可以采用迭代解法和直接解法。 （1）迭代法 在克拉索夫斯基椭球上计算时，迭代开始时设 以后每次迭代按下式计算： 重复迭代直至为止。 在1975年国际椭球上计算时，也有类似公式。 （2）直接解法 1975年国际椭球： 克拉索夫斯基椭球： 大地线 椭球面上两点间的最短程曲线叫做大地线。在微分几何中，大地线（又称测地线）另有这样的定义：“大地线上每点的密切面（无限接近的三个点构成的平面）都包含该点的曲面法线”，亦即“大地线上各点的主法线与该点的曲面法线重合”。因曲面法线互不相交，故大地线是一条空间曲面曲线。 假如在椭球模型表面,两点之间，画出相对法截线如图所示，然后在,两点上各插定一个大头针，并紧贴着椭球面在大头针中间拉紧一条细橡皮筋，并设橡皮筋和椭球面之间没有摩擦力，则橡皮筋形成一条曲线，恰好位于相对法截线之间，这就是一条大地线。由于橡皮筋处于拉力之下，所以它实际上是两点间的最短线。 在椭球面上进行测量计算时，应当以两点间的大地线为依据。在地面上测得的方向、距离等，应当归算成相应大地线的方向、距离。