第四篇 3椭球面上的弧长(大地线)计算.ppt

第四章 Ⅲ椭球面上的弧长计算 ——子午线弧长计算公式 ——由子午线弧长求大地纬度 ——平行圈弧长计算公式 ——椭球面梯形图幅面积的计算 ——大地线的概念与克莱劳方程 上一讲应掌握的内容 1、椭球面上法截线有关概念 法截线 子午圈(线) 卯酉圈 (线) 曲率 曲率半径 2、子午线曲率半径 在赤道上M小于a，M随纬度增加而增加 3、卯酉线曲率半径 在赤道上N等于a，N随纬度增加而增加 4、任意法截弧的曲率半径 主曲率半径M及N分别是 RA的极小值和极大值。R0＝M，R90＝N 5、平均曲率半径 （不是地球平均半径） 6、椭球面上几种曲率半径的关系 子午圈曲率半径M及卯酉圈曲率半径N，是两个互相垂直的法截弧的曲率半径，这在微分几何中统称为主曲率半径。 §4.4 椭球面上的弧长计算 在高斯投影计算和弧度测量计算中， 往往用到子午线弧长和平行圈弧长 一、子午线弧长计算公式 子午线弧长计算公式推导（方法1） 子午线弧长计算 对于75国际椭球，子午线弧长的具体计算公式为： 如果以B＝90°代入，则得子午椭圆在一个象限内的弧长约为10 002 137m。旋转椭球的子午圈的整个弧长约为40 008 549.995m。即一象限子午线弧长约为10 000km，地球周长约为40 000km。 为求子午线上两个纬度B１及Ｂ２间的弧长，只需按上式分别算出相应的X１及X２，而后取差： Ｘ1 2＝Ｘ２－Ｘ１，该ΔＸ即为所求的弧长。 子午线弧长计算（续） 对于400km以下的子午线弧长计算公式 （精确到：0.001m ） 若弧长＜45km,仅取△B一次项，即可精确到：0.001m 当弧长甚短(例如X≤45km，计算精度到0.001m)，可视子午弧为圆弧，而圆的半径为该圆弧上平均纬度点的子午圈的曲率半径Mm 子午线弧长计算公式推导(方法2) 将被积函数按二项式作级数展开，有： 代入上式，积分得： 子午线弧长计算公式（2）系数 由子午弧长求大地纬度可以采用迭代解法和直接解法 迭代解法: 直接解法： 三、平行圈弧长公式 任一平行圈都是半径相等的圆， 纬度为B的平行圈的半径为： 四、椭球面梯形图幅面积的计算 只要求出经纬格网的面积，就可度量出椭球面的整个面积或局部区域的面积。 四、椭球面梯形图幅面积的计算（续） 椭球面积 若取L2-L1=2π，B2=π/2，B1=0 可得出北半球的椭球面面积； 再乘以2得到椭球面的总面积。 克氏椭球总面积约为5.1亿km2 梯形图幅的面积，直接代入经纬度计算。 对于形状不规则的区域面积可将该区域分割成若干个规则格网的子块，然后求和。 §4.5 大地线的概念与计算 两点间的最短距离，在平面上是两点间的直线，在球面上是两点间的大圆弧，那么在椭球面上又是怎样的一条线呢? 它应是大地线。 一、相对法截弧 法截弧：由椭球面上A点的法线与B点所确定的法截面与椭球面相割得到的曲线称为A到B的法截弧。 相对法截弧： A到B的法截弧与B到A的法截弧。 互为正反法截弧 当A，B两点位于同一子午圈或同一平行圈上时，正反法截线则合二为一，但不是平行圈。 在通常情况下，正反法截线是不重合的。因此在椭球面上A，B，C三个点处所测得的角度(各点上正法截线之夹角)将不能构成闭合三角形。为了克服这个矛盾，在两点间另选一条单一的大地线代替相对法截线，从而得到由大地线构成的单一的三角形。 二、大地线的定义与性质 定义2：大地线是主法线与曲面法线处处重合的曲线。 三、大地线的微分方程（推导） 四、大地线的克莱劳方程 克莱劳方程的应用 对于同一条大地线上的各点因要保持同一个C值，而使点上的大地线方位角与平行圈半径两者间的变动互为制约。 当大地线穿越赤道时 赤道是大地线 C = a 当大地线达极小平行圈时 子午圈是大地线 C = 0 由克莱劳方程可以写出 用以检查纬度和方位角计算的正确性 习题 在克拉索夫斯基椭球面上，某一大地线穿越赤道时的大地方位角A=30o，试求该大地线所能达到的最小平行圈上，经差l=6′52. 53″的平行圈弧长（精确至0.01m） 结束 再见! * * 主曲率半径的计算公式 不同的椭球元素对应不同的系数 写在黑板 主曲率半径的计算公式中的系数 为了便于积分通常将正弦的幂函数展开为余弦的倍数函数。如： 子午圈曲率半径公式可变为： 积分，整理得： 板书 写 写 对于75国际椭球，子午线弧长的具体计算公式为： 二、由子午弧长求大地纬度 对于克氏椭球 子午线弧长和平行圈弧长变化的比较 L L+dL B B+dB MdB NcosBdL d? 上式利用二项式展开并积分