第四章3椭球面上的弧长大地线计算.pptxVIP

上一讲应掌握的内容1、椭球面上法截线有关概念 法截线 子午圈(线)卯酉圈 (线)曲率 曲率半径2、子午线曲率半径 在赤道上M小于a，M随纬度增加而增加3、卯酉线曲率半径 在赤道上N等于a，N随纬度增加而增加4、任意法截弧的曲率半径 主曲率半径M及N分别是 RA的极小值和极大值。R0＝M，R90＝N5、平均曲率半径 （不是地球平均半径）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）系数对于75国际椭球，子午线弧长的具体计算公式为：二、由子午弧长求大地纬度由子午弧长求大地纬度可以采用迭代解法和直接解法迭代解法:直接解法：对于克氏椭球三、平行圈弧长公式任一平行圈都是半径相等的圆，纬度为B的平行圈的半径为：子午线弧长和平行圈弧长变化的比较B+dBd?MdBBNcosBdLLL+dL四、椭球面梯形图幅面积的计算只要求出经纬格网的面积，就可度量出椭球面的整个面积或局部区域的面积。上式利用二项式展开并积分，得：四、椭球面梯形图幅面积的计算（续）即：椭球面积若取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三个点处所测得的角度(各点上正法截线之夹角)将不能构成闭合三角形。为了克服这个矛盾，在两点间另选一条单一的大地线代替相对法截线，从而得到由大地线构成的单一的三角形。 二、大地线的定义与性质定义1：椭球面上两点间的最短程曲线叫大地线（测地线）。定义2：大地线是主法线与曲面法线处处重合的曲线。 大地线的性质：1、曲面上连接任何两点的最短弧线必为大地线；2、大地线上任何点的密切平面就是该点的法截面；3、大地线的测地曲率等于0；4、大地线位于相对法截线之间，并靠近正法截线，它与正法截线间的夹角 5、在椭球面上进行测量计算时，应当以两点间的大地线为依据。在地面上测得的方向、距离等，应当归算成相应大地线的方向、距离。 长度差异可忽略（40km以内）， 但方向差异需改化。三、大地线的微分方程（推导）由球面余弦定理得：三、大地线的微分方程四、大地线的克莱劳方程 在旋转椭球面上，大地线各点的平行圈半径与大地线在该点的大地方位角的正弦的乘积等于常数。式中常数C也叫大地线常数  克莱劳方程的应用 对于同一条大地线上的各点因要保持同一个C值，而使点上的大地线方位角与平行圈半径两者间的变动互为制约。当大地线穿越赤道时 赤道是大地线 C = a当大地线达极小平行圈时 子午圈是大地线 C = 0由克莱劳方程可以写出 用以检查纬度和方位角计算的正确性 习题在克拉索夫斯基椭球面上，某一大地线穿越赤道时的大地