第四章 地球椭球及其数学计算

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,上一堂课内容回顾,地球,形状,大地水准面与地球椭球,大地水准面（一次逼近）,地球椭球（二次逼近）,参考椭球,地球重力场,地球引力（位）,离心力（位）,重力（位）,垂线偏差,地球磁场,地磁七要素,上一堂课内容回顾,重力场与地磁场的异同,相同点：,均为人为不可控天然稳定场（对导航有利）,均随空间位置变化，具有一定的空间分布规律,都随时间发生微小变化,不同点：,重力场和地磁场成因不同、性质不同；,重力场为,单极场,，地磁场近似为,偶极场,；,重力场为,强场,，地磁场为,弱场,；,地磁场为,强时变，,重力场,弱时变；,第四章 地球椭球及其数学计算,张小红,武汉大学测绘学院,导航学,第四章 地球椭球及其数学计算,4.1,地球椭球的几何参数及其相互关系,4.2,大地坐标系、空间直角坐标系及其相互关系,4.3,地心纬度、归化纬度及其与大地纬度间的关系,4.4,地球椭球上的曲率半径,4.5,椭球面上的弧长计算,4.6,法截线与大地线,4.7,大地主题解算,4.8,导航中大地线长度的计算方法,4.9,把地面观测值归算至椭球面,第四讲,第四章 地球椭球及其数学计算,第四章 地球椭球及其数学计算,第一节,地球椭球的几何参数,及其相互关系,4.1,地球椭球的几何参数及其相互关系,椭球上的点和线,地球椭球是一个具有合适的形状和大小的椭圆绕短轴旋转一周后所形成的一个,旋转椭球,北极,N,和南极,S,椭球中心,O,赤道平面（赤道圈）,子午面（子午圈）,平行圈或纬圈,旋转椭球体的特点,对称性,过任意一点的子午圈的形状和大小相同,平行圈（纬圈）和赤道圈都是正圆,子午圈的形状和大小 决定了地球椭球的形状和大小,4.1,地球椭球的几何参数及其相互关系,椭球的基本几何参数,椭球长半径,椭球短半径,椭球的扁率,椭球的第一偏心率,椭球的第二偏心率,上述,5,个参数中任选两个参数就能表示椭球的形状和大小，但其中至少有一个长度参数 ，通常选 和,其中,a,，,b,称为长度元素，扁率反映了椭球体的扁平程度。,偏心率是子午椭圆的焦点离开中心的距离与椭圆半径之比，它们也反映椭球的扁平程度，偏心率愈越大，椭球越扁,4.1,地球椭球的几何参数及其相互关系,椭球几何参数间的相互关系,4.1,地球椭球的几何参数及其相互关系,辅助参数（为简化后续公式推导）,极点处的子午曲率半径,第四章 地球椭球及其数学计算,第二节,大地坐标系、空间直角坐标系,及其相互关系,4.2,大地坐标系、空间直角坐标系及其相互关系,大地坐标系,大地坐标系是,大地测量学与导航学,中常用的一种坐标系，亦称,地理坐标系,或,椭球坐标系,它是以经过椭球定位后的地球椭球上所定义的,点线面,为参考的一种坐标系。,地面一点的大地坐标（,B,，,L,，,H),大地纬度,B (N/S 090),大地经度,L (E/W 0180),大地高,H,4.2,大地坐标系、空间直角坐标系及其相互关系,空间直角坐标系是大地测量与导航计算常用的坐标系,空间直角坐标系定义,坐标原点,O,：,位于总地球椭球(或参考椭球)中心,Z,轴：,与地球平均自转轴相重合，指向某段时间的平均北极点；,X,轴：,指向由平均格林尼治天文台和平均自转轴所确定的子午面与赤道面的交点,Ge,；,Y,轴：,垂直于,X,轴和,Z,轴构成,右手系,4.2,大地坐标系、空间直角坐标系及其相互关系,子午面直角坐标系,（大地坐标与空间直角坐标系转换所需的,中间坐标系,）,P,点为空间某点,P,沿法线方向,在地球椭球上的投影点,，以过,P,点的子午椭圆中心为原点，建立一个平面直角坐标系，,x,轴与子午椭圆的长轴重合，,y,轴与椭圆的短轴重合,。,在该坐标系中，,P,点的位置用,( x, y),表示,过,P,点作子午椭圆的切线,TP,，切线的斜率为,4.2,大地坐标系、空间直角坐标系及其相互关系,子午面直角坐标系,引入辅助参数,代入,x,4.2,大地坐标系、空间直角坐标系及其相互关系,大地坐标转空间直角坐标,在椭球面上的点,不在椭球面上的点（推导）,4.2,大地坐标系、空间直角坐标系及其相互关系,大地坐标系到空间直角坐标系的转换推导思路,建立空间直角坐标系,建立子午面直角坐标系（中间过渡）,推导子午面直角坐标和大地纬度与椭球有关参数之间的关系,找到空间直角坐标和子午面直角坐标之间的相互关系,建立空间直角坐标和大地坐标之间的关系,4.2,大地坐标系、空间直角坐标系及其相互关系,空间直角坐标转大地坐标,迭代公式,迭代初值：,或,4.2,大地坐标系、空间直角坐标系及其相互关系,空间直角坐标系转大地坐标系,直接公式,低精度直接公式,高精度直接公式,纬度的精度可达,大地高的误差小于,第四章 地球椭球及其数学计算,第三节,地心纬度、归化纬度及其,与大地纬度间的关系,4.3,地心纬度、归化纬度及其与大地纬度间的关系,4.3,地心纬度、归化纬度及其与大地纬度间的关系,空间直角坐标与归化纬度间的关系,P,点在子午面直角坐标系中的坐标,4.3,地心纬度、归化纬度及其与大地纬度间的关系,4.3,地心纬度、归化纬度及其与大地纬度间的关系,大地纬度、地心纬度、归化纬度之间的差异很小，经过计算，当,B=45,时：,第四章 地球椭球及其数学计算,第四节 地球椭球上的曲率半径,4.5,椭球,面,上的,弧长计算,基本知识,三角函数级数展开,4.5,椭球,面,上的,弧长计算,基本知识,弧度和度的定义,角度是表示角的大小的量，通常用度或弧度来表示,角度制：,规定周角的,360,分之一为,1,度的角,弧度制：,规定长度等于半径的弧长所对的圆心角为,1,弧度,4.4,地球椭球上的曲率,半径,子午圈曲率半径,M,对于一条平面曲线,其曲率半径可用下式计算,4.4,地球椭球上的曲率,半径,子午圈曲率半径,M,4.4,地球椭球上的曲率,半径,子午圈曲率半径,M,B,M,说明,极点处的子午曲率半径,4.4,地球椭球上的曲率,半径,卯酉圈,过椭球面上任意一点,P,可作一条垂直于椭球面的法线,PF,，包含这条法线的平面叫作,法截面,，法截面与椭球面的交线叫,法截线,过椭球面上一点的法线，可作无限个法截面，其中,与子午面垂直,的,法截面称为,卯酉面,，卯酉面与椭球面的交线称为,卯酉圈,卯酉圈的曲率半径通常用,符号,N,表示,4.4,地球椭球上的曲率,半径,卯酉圈曲率半径,N,麦尼尔定理：,假设通过曲面上一点引两条截弧，一条为法截弧（卯酉圈），一条为斜截弧（平行圈），且,在该点上这两条截弧具有公共切线,，这时斜截弧在该点处的曲率半径,r,等于,法截弧的曲率半径,N,乘以两截弧平面夹角,B,的余弦,卯酉圈曲率半径,N,就是从,P,点至法线与椭球短轴的交点,F,间的距离,(,课后推导,),4.4,地球椭球上的曲率,半径,卯酉圈曲率半径,N,B,M,说明,4.4,地球椭球上的曲率,半径,子午圈和卯酉圈圈曲率半径级数展开（实际计算）,或,顾及,8,次项,的计算公式一般已能保证,mm,级,的计算精度,（参见,4-47,式）,4.4,地球椭球上的曲率,半径,任意方向法截弧的曲率半径,子午法截弧,是南北向，方位角为,0,或,180,卯酉法截弧,是东西向，其方位角为,90,或,270,子午法截弧和卯酉法截弧在,P,点处正交,过,P,点的,子午曲圈率半径,M,和,卯酉圈曲率半径,N,称为曲面在该点的两个主曲率半径,4.4,地球椭球上的曲率,半径,任意方向法截弧的曲率半径,尤拉公式,级数展开,A,为任意法截弧的大地方位角,4.4,地球椭球上的曲率,半径,4.4,地球椭球上的曲率,半径,平均曲率半径,平均曲率半径就是过该点的所有的法截弧的曲率半径的,算术平均值,积分,椭球面上任一点处的平均曲率半径就等于该处的子午圈曲率半径与卯酉圈曲率半径的几何平均值,4.4,地球椭球上的曲率,半径,M,、,N,、,R,的关系,曲率半径,N,R,M,公式,4.4,地球椭球上的曲率,半径,M,、,N,、,R,的数值表,B,N,(m),R,(m),M,(m),6,378 245,6 356 863,6 335 553,6 379 675,6 359 714,6 339 816,6 383 588,6 367 518,6 351 488,6 388 945,6 378 209,6 367 491,6 394 315,6 388 936,6 383 561,6 398 255,6 396 811,6 395 368,6 399 699,6 399 699,6 399 699,参数：克拉索夫斯基椭球体,第 五 讲,第四章 地球椭球及其数学计算,上一堂课内容回顾,地球椭球的几何参数及其相互关系,大地坐标系、空间直角坐标系及其相互关系,地心纬度、归化纬度及其与大地纬度间的关系,椭球上的曲率半径,子午圈曲率半径,M,法截弧,、,卯酉圈的定义，卯酉圈,曲率半径,N,任意方向上的法截弧曲率半径,RA,平均曲率半径,R,N,R,M,均随纬度增大而增大,B=90,时，,N=R=M=c,第四章 地球椭球及其数学计算,第五节,椭球面上的弧长计算,4.5,椭球,面,上的,弧长计算,子午线弧长,子午椭圆的一半，它的端点与极点重合,赤道又把子午线分成对称的两部分,椭球上每个子午圈的形状和大小均相同,计算子午线弧长与该子午圈的经度无关,推导从赤道沿子午线至任一纬度,B,的子午线弧长计算公式就可以计算出所需要的子午弧长,4.5,椭球,面,上的,弧长计算,子午线弧长计算公式,其中,4.5,椭球,面,上的,弧长计算,子午线弧长计算公式,小于,0.1mm,，可以忽略,4.5,椭球,面,上的,弧长计算,子午线计算公式,克拉索夫斯基椭球,：,1975,年国际椭球,：,CGCS2000,地球椭球,：,WGS 84,椭球,：,4.5,椭球,面,上的,弧长计算,当子午线很短时，例如子午线两端的纬差 时，可将子午线视为圆弧。其曲率半径采用两端的平均纬度处的子午曲率半径,子午线弧长公式可简化为：,1,弧度对应的度，分，秒值,4.5,椭球,面,上的,弧长计算,不同地球椭球上的子午线弧长,椭球,克拉索夫斯基,椭球,5540944.463,2212405.723,3328538.740,1975,年国际椭球,5540849.645,2212367.296,3328482.349,CGCS2000,椭球,5540847.039,2212366.253,3328480.786,WGS 84,椭球,5540847.056,2212366.256,3328480.800,北纬 至北纬 的子午线弧长,4.5,椭球,面,上的,弧长计算,平行圈（圆）弧长计算,经度差,4.5,椭球,面,上的,弧长计算,由子午线弧长求大地纬度,迭代解法,直接解法,4.5,椭球,面,上的,弧长计算,曲率半径,M,、,N,及弧长,S,、,S,随纬度,B,的变化,6335439.327,6378137.000,1842.9046,30.7151,1855.3248,30.9221,6337358.121,6378780.844,1843.4628,30.7244,1827.3227,30.4554,6342888.483,6380635.807,1845.0715,30.7512,1744.1181,29.0686,6351377.104,6383480.918,1847.5407,30.7923,1608.1047,26.8017,6361815.827,6386976.166,850.5772,30.8430,1423.2369,23.7205,6372955.926,6390702.844,1853.9177,30.8970,1194.9292,19.9155,6383453.858,6394209.174,1856.8715,30.9479,930.0000,15.5000,6392033.193,6397072.488,1859.3671,30.9895,636.4424,10.6074,6397643.327,6398943.460,1860.9990,31.0167,323.2248,5.3871,6399593.626,6399593.626,1861.5663,31.0261,0.0000,0.0000,1850.5772,4.5,椭球,面,上的,弧长计算,第四章 地球椭球及其数学计算,第六节,法截线与大地线,4.6,法截线与大地,线,法截线面与法截线,过椭球面上任意一点,P,可作一条垂直于椭球面的法线,PF,，包含这条法线的平面叫作,法截面,，法截面与椭球面的交线叫,法截线,4.6,法截线与大地,线,相对法截线,地球椭球是一个旋转椭球，过椭球面上任意一点,D,作椭球面的法线时，该法线必定位于过,D,点的子午面上，与椭球短轴交于,D,1,D,点在子午平面直角坐标系中的,y,坐标为,4.6,法截线与大地,线,相对法截线,如果椭球上另有一点,E,，（ ， ），过,E,点作椭球的法线（该法线位于过,E,点的子午面内），与椭球短轴交于,E,1,4.6,法截线与大地,线,相对法截线,当椭球面两个点不在同一子午圈上，也不在同一平行圈上时，两点间就有两条法截线存在,D,点,(D,到,E),正法截线,D,点,(E,到,D),反法截线,称为相对法截线,不在同一子午圈或同一平行圈上的两点的正反法截线不重合,正法截线：,过某一点的法线和椭球面上的另一点作一法截面，与椭球相交的线为该点的正法截线,反法截线：,过另一点的法线和该点作一法截面，与椭球相交的线为该点的反法截线,4.6,法截线与大地,线,相对法截线,不在同一子午圈或同一平行圈上的两点的正反法截线,不重合,某点的纬度越高，其法线与短轴的交点离开椭球中心就越远,法截线,BbA,偏上（北），法截线,AaB,偏下（南）,当,A,，,B,两点位于同一子午圈或,同一平行圈时，正反法截线重合,4.6,法截线与大地,线,相对法截线,北半球，,A,点在,B,点的南面，则,A,的正法截线也在南面，而反法截线在北面,南半球则相反,北向,东向,正法截线,反法截线,AB,方向在不同象限时，正反法截线的关系,4.6,法截线与大地,线,大地线的提出,正反法截线不重合,椭球面上,A,，,B,，,C,三个点处所测得角度,(,各点上正法截线之夹角,),不能构成闭合三角形（图中实线部分的观测角）,为了克服这个矛盾，在两点间另选一条,唯一的,大地线,代替相对法截线，从而得,到由大地线构成的封闭椭球面三角形,内角和不等于,180,4.6,法截线与大地,线,两点之间的最短距离,平面上是连接两点的直线端,球面上是连接两点的大圆弧段（圆心和球心重合的圆就是大圆）,在椭球面上呢？,4.6,法截线与大地,线,大地线定义,定义,1,：大地线上每点的,密切平面,（无限接近的三个点构成的平面）都包含该点的曲面法线，即大地线上各点的主法线与该点的曲面法线重合。大地线又称为,测地线。（微分几何中的定义）,定义,2,：,椭球面上两点间距离最短的曲线叫,大地线,在平面上大地线就是一条直线,在球面上大地线是一段大圆弧,（圆心和球心重合的圆就是大圆）,4.6,法截线与大地,线,大地线的性质,大地线是两点间唯一最短线,大地线位于相对法截线之间,大地线更靠近正法截线,大地线上任何点的密切平面就是该点的法截面,两点的正反法截线不重合，它们之间的夹角,，在一等三角测量中,(,平均点间距,40km),可达千分之四秒,大地线与正法截线的夹角,大地线与法截线长度之差非常小，可忽略不计,在椭球面上进行测量和计算时，应当以两点间的,大地线为依据,地面上测得的方向，距离等，应归算到相应,大地线的方向和距离,4.6,法截线与大地,线,大地线的微分方程,指的是大地线的弧长,S,与所定义坐标系中的自变量间的微分关系,在空间直角坐标系中,dS,与,dx,，,dy,，,dz,间的关系,在大地坐标系统中,dS,与,dL,，,dB,，,dA,间的关系,4.6,法截线与大地,线,球面三角形余弦定理,若已知球面三角形,ABC,的三頂点是,A,、,B,、,C,，其所对应的三边分別是,a, b, c,，則有球面上的角余弦定理,:,cos,A,= -,cosBcosC,+,sinBsinCcosa,4.6,法截线与大地,线,大地线的微分方程,球面三角形,(,p,1,p,3,N),余弦定理,:,微小量,4.6,法截线与大地,线,大地线微分方程,4.6,法截线与大地,线,大地线的克莱劳方程,积分,克莱劳定理,：,4.6,法截线与大地,线,大地线的克莱劳方程,克莱劳定理,：,当大地线穿越赤道时，,B=0,，,r=a,，,A=A,0,当大地线达极小平行圈时，,A=90,，设,r=r,0,，,B=B,0,某一大地线常数等于椭球半径与该大地线穿越赤道时的大地方位角正弦的乘积，或者等于该大地线上具有最大纬度的那一点的平行圈半径,可用于检测纬度和大地方位角计算的正确性,思考：赤道是不是大地线，子午圈是不是大地线？,4.6,法截线与大地,线,大地线的克莱劳方程,图中有一条大地线，该大地线在与赤道相交的,E,点处的大地方位角为,A,，则该大地线将与半径 的平行圈相切于,F,点，然后调头向南，而不会跑到纬度更高的地区去,4.6,法截线与大地,线,相对法截线夹角,4.6,法截线与大地,线,截面差,A,点的大地线与正法截线,的,夹角（,截面差,）：,可忽略不计,在一、二等,大地测量中,，截面差具有系统性，,为防止误差的不断积累,，,进行方向观测值的归算时均需施加截面差改正,导航中很少,在椭球面上进行精密方向观测值,的,计算,4.6,法截线与大地,线,大地线与法截线的长度之差,第四章 地球椭球及其数学计算,第七节,大地,主题解算,4.7,大地主题解算,大地主题解算分为:,短距离(400,km,),中距离(1000,km),长距离(1000,km,以上),4.7,大地主题解算,大地主题解算,在涉及描述地球表面或近地空间相关学科，如,大地测量,、,导航制导、远程导弹的弹道解算、远洋航海、航天测控和地球物理,等领域的问题时，经常会遇到计算椭球面上两点之间距离、方向等问题,是弹道计算的核心算法,在航路规划中，大地主题解算是最基本的运算，其运算的速度，运算结果的精度将直接关系到航路规划的效果,4.7,大地主题解算,大地主题解算,解算时涉及复杂的数学运算与推导，以前没有电子计算机，处理非常繁琐，先后有一大批学者致力于大地主题解算的研究,1806,年法国著名数学家勒让德首先采用泰勒级数将大地线微分方程展开为大地线长度的升幂级数，用以解算大地主题,历史上研究大地主题解算的算法有许多种，据不完全统计，有,100,多种，如,Gauss,法、,Helmer,t,法、,Bessel,法和,Boring,法等,手工推导得到的直接解算公式，形式较复杂，使用不甚方便,随着计算机技术的发展，大地主题解算更简单，更准确,4.7,大地主题解算,大地线在大地坐标系中的微分方程,其精确值不能直接计算，必须进行逼近求解,4.7,大地主题解算,4.7,大地主题解算,4.7,大地主题解算,在起点,P,1,处展开的级数算法,一阶导数：,二阶导数：,再次求导,三阶导数：,四阶导数：,五阶导数：,将所有求得的导数代入,麦克劳林级数,展开式,4.7,大地主题解算,在起点,P,1,处展开的级数算法,纬度计算公式：,4.7,大地主题解算,在起点,P,1,处展开的级数算法,经度计算公式：,4.7,大地主题解算,在起点,P,1,处展开的级数算法,方位角计算公式：,4.7,大地主题解算,提出高斯平均引数法,4.7,大地主题解算,高斯平均引数法,高斯平均引数法在大地线的,中点处,M,同时向前和向后用级数展开，然后再合并在一起,推导思路：,首先把勒让德级数在,P,1,展开改为在大地线中点,M,处展开，使公式项数减少，加快收敛，提高精度,将,M,点用大地线两端点平均纬度及平均方位角相对应的,m,点来代替，并迭代计算,优点,可以使大地线的长度减半,可消除公式中的偶阶项,同时适用于大地主题正、反算,4.7,大地主题解算,高斯平均引数算法正算公式,相减,4.7,大地主题解算,高斯平均引数算法正算公式,4.7,大地主题解算,高斯平均引数算法正算公式,怎么解决？,引入平均纬度和方位角,未知，需迭代,4.7,大地主题解算,高斯平均引数算法正算公式,相加除,2,=,同理得：,4.7,大地主题解算,高斯平均引数算法正算公式,前面已得到,4.7,大地主题解算,高斯平均引数算法正算公式,4.7,大地主题解算,高斯平均引数算法正算公式,代入,略去高阶项,4.7,大地主题解算,高斯平均引数算法正算公式（,适用于,120km,以内,）,4.7,大地主题解算,高斯平均引数算法正算公式,4.7,大地主题解算,高斯平均引数算法正算公式（,适用于,70km,以内,）,其中主项为：,主项迭代,3,次,改正项迭代,12,次,4.7,大地主题解算,4.7,大地主题解算,高斯平均引数算法反算公式,将正算公式左右移项后得：,4.7,大地主题解算,高斯平均引数算法反算公式,将正算公式左右移项后得：,省去高阶项，右式中用主项代替，主项为：,主项,用主项代入,4.7,大地主题解算,高斯平均引数算法反算公式（,适用于,200km,以内,）,令：,得：,无需迭代,直接求得,4.7,大地主题解算,高斯平均引数算法反算公式（,适用于,200km,以内,）,代入至下式,4.7,大地主题解算,高斯平均引数算法反算公式（,适用于,70km,以内,）,公式简化得：,其中有：,其余算法和,200km,距离的反算一样,4.7,大地主题解算,高斯平均引数算法反算公式总结和要点,反算公式由正算公式左右移项得到,利用主项替换公式右边的未知项,无需迭代，直接求得,对于短距离大地线，通过省去高次项，可以简化计算公式,4.7,大地主题解算,白塞尔大地主题解算法,按照一定的法则将椭球面上的大地元素投影至一个辅助球上,，,实现椭球面向球面的过渡,在圆球上采用简单的球面三角公式进行计算，推求大地元素,将计算结果转换至椭球面上,，,实现球面向椭球面的过渡,特点：,计算精度不会随着距离的增加而迅,速下降,，,适用于,短距离和,长距离大地主题解算,算法：,白塞尔大地主题算法,4.7,大地主题解算,白塞尔大地主题解算方法,1825,年德国大地测量学家、天文学家白塞尔（,F.W.Bessel,）提出了一种长距离大地主题解算公式。,基本思想是,：,首先,按照一定法则将椭球面上的已知大地元素“投影”至一个单位球面上，然后再在这个辅助球面上用球面三角公式进行严密的坐标正反算，最后再将计算结果再转换至椭球面上来,。由于地球椭球的扁率仅为,1/297,，与圆球非常接近，因而从椭球面至圆球（或从圆球至椭球）的投影的“变形”就很小，即使采用一个相对较为简单的公式就能适用于长距离的大地主题解算。,此后又有不少学者对该公式进行了改进。,特点：,计算精度不会随着距离的增加而迅速下降,，,适用于,短距离和,长距离大地主题解算,4.7,大地主题解算,直接对大地线微分方程进行数值积,分,采用数值计算的方法以适当的步长对大地线微分方程进行数值积分,优点：易于编程实现，适用于不同长度的大地主题解算,缺点：,计算工作量大，随着距离的增加精,度可能会有所下降，此外在两极地区效果也不好,算法：,龙格,库塔法，阿达姆斯法,4.7,大地主题解算,其它方法,转换至高斯平面进行计算,直接在空间中计算点的几何关系（卫星大地测量）,第四章 地球椭球及其数学计算,第八节,导航中大地线长度的计算方法,4.8,导航中大地线长度的计算方法,导航中大地线长度的计算,需求：由两点大地坐标求大地线,特点：精度要求不高，公式要求简单方便，计算速度快,方法：,Andoyer,-Lambert,法,，,大椭圆法,4.8,导航中大地线长度的计算方法,Andoyer,-Lambert,法,该法是,以适用于中长距离的白塞尔大地主题解算公式为基础进行简化和近似处理,得到的,大地线微分方程：,替换,4.8,导航中大地线长度的计算方法,Andoyer,-Lambert,法,大地线微分方程：,泰勒展开：,4.8,导航中大地线长度的计算方法,Andoyer,-Lambert,法,积分,未知量，需要求解,4.8,导航中大地线长度的计算方法,Andoyer,-Lambert,法,未知，先放着，稍后做处理,4.8,导航中大地线长度的计算方法,Andoyer,-Lambert,法,4.8,导航中大地线长度的计算方法,Andoyer,-Lambert,法,4.8,导航中大地线长度的计算方法,Andoyer,-Lambert,法,4.8,导航中大地线长度的计算方法,Andoyer,-Lambert,法,4.8,导航中大地线长度的计算方法,Bowring,公式,英国的,Bowring,与,1981,年按椭球面对球面的正形投影,导出了一个崭新地椭球面上主题解算公式。,易于编程,特别适用于可编程计算器上解算及实时应用中。,在美国导航学会第三十九届年会上,该法被建议应用于导航领域。推荐此法作为短距离大地主题解算的最好公式。,4.8,导航中大地线长度的计算方法,Bowring,正解公式,4.8,导航中大地线长度的计算方法,Bowring,反解公式,4.8,导航中大地线长度的计算方法,大椭圆法解算大地线,请同学们课后自学！,4.8,导航中大地线长度的计算方法,4.8,导航中大地线长度的计算方法,顾及高程影响的大地线长度计算方法,扁率不变,4.8,导航中大地线长度的计算方法,顾及高程影响的大地线长度计算方法,在辅助椭球上进行大地线的计算,第四章 地球椭球及其数学计算,第九节,把地面观测值归算至椭球面,4.9,把地面观测值归算至椭球,面,为什么要把,地面观测值归算至椭球,面,在传统大地测量中，需要先把地面测站上所测得的方向观测值和距离观测值等归算至椭球面上，然后才能在椭球面这个规则的数学曲面上进行数学计算。,在导航学中，也可以沿用这种方法。当然在导航学中，也可采用另一种方法：完全脱离地球椭球的概念，直接在空间直角坐标系中按照点与点之间的几何关系来进行数据处理。如有必要再将这些空间直角坐标系化算为大地坐标,4.9,把地面观测值归算至椭球,面,距离归算值，具体推导见课本,4.9,把地面观测值归算至椭球,面,偏差改正值,4.9,把地面观测值归算至椭球,面,标高差改正值,第四章 地球椭球及其数学计算,4.1,地球椭球的几何参数及其相互关系,4.2,大地坐标系、空间直角坐标系及其相互关系,4.3,地心纬度、归化纬度及其与大地纬度间的关系,4.4,地球椭球上的曲率半径,4.5,椭球面上的弧长计算,4.6,法截线与大地线,4.7,大地主题解算,4.8,导航中大地线长度的计算方法,4.9,把地面观测值归算至椭球面,