卯酉圈曲率半径

1大地测量学基础第四章地面观测值归算至椭球面2第四章地面观测值归算至椭球面4.1地球椭球的基本性质4.2将地面观测值归算至椭球面4.3大地测量主题解算概述本章重点：地球椭球几何性质、地面观测值归算本章难点：地面观测值归算至椭球面、大地主题解算34.1地球椭球的基本性质4.1.1椭球面上的几种法截线的曲率半径法截面和法截线：过椭球面上任意一点可作一条垂直于椭球面的法线，包含这条法线的平面叫作法截面；法截面与椭球面的交线叫法截线。4一、子午圈曲率半径4.1.1椭球面上的几种法截线的曲率半径子午椭圆的一部分上取一微分弧长DK＝dS，相应地有坐标增量dx，点n是微分弧dS的曲率中心，于是线段Dn及Kn便是子午圈曲率半径Ｍ。由任意平面曲线的曲率半径的定义公式，易知5一、子午圈曲率半径6一、子午圈曲率半径7二、卯酉圈曲率半径过椭球面上一点的法线，可作无限个法截面，其中一个与该点于午面相垂直的法截面同椭球面相截形成的闭合的圈称为卯酉圈。卯酉圈的曲率半径用N表示。卯酉圈曲率半径8卯西圈曲率半径恰好等于法线介于椭球面和短轴之间的长度，亦即卯酉圈的曲率中心位于椭球的旋转轴上。二、卯酉圈曲率半径9二、卯酉圈曲率半径10三、主曲率半径的计算主曲率半径:子午圈曲率半径Ｍ及卯酉圈曲率半径N，是两个互相垂直的法截弧的曲率半径，称为主曲率半径。2322211)sin)((BeeaM21221)sin(BeaN11四、任意法截线的曲率半径AMANMNRA22sincos按泰勒级数展开：)coscos(AANRA4422112四、任意法截线的曲率半径引入平均曲率半径Ｒ代入上式得注：子午圈曲率半径和卯酉圈曲率半径的变化规律13五、平均曲率半径在一定范围内，把椭球面当做球面来处理，推出这个球面的曲率半径－平均曲率半径：MNR或)(22221eWaVNVcWbR14六、M、N、R之间的关系一般情况下NRM在极点上都等于极曲率半径ｃ曲率半径NRM公式1Vc2Vc3Vc121Wea221Wea321WeacMRN909090154.1.2椭球面上的弧长计算一、子午线弧长计算公式取子午线上微分弧PP’=dxP点的子午圈半径为Ｍ赤道到任意纬度Ｂ的平行圈之间的弧长MdBdxBMdBX016二、由子午弧长求大地纬度1.迭代解法2.直接解法17三、平行圈弧长公式旋转椭球体的平行圈是一个圆，其短半轴ｒ就是圆上任意一点的子午面直角坐标ｘ如果平行圈上有两点，它们的经度差平行圈弧长公式：18三、平行圈弧长公式平行圈弧长随纬度变化的微分公式可近似地写为由于式中194.1.3大地线一、相对法截弧设在椭球上任取两点A和B，纬度分别为B1和B2，且二者不等，过A、B两点分别做法线与短轴和赤道面相交。2211BnQOnBnQOnbbaasinsin如图又则2NeQn222121BeNOnBeNOnbasinsin20一、相对法截弧正、反法截线：A点照准B，照准面同椭球面的交线AaB，叫做A点的正法截线，或者B点的反法截线；同样，B点照准A，照准面同椭球面的交线BbA，叫做B点的正法截线，或者A点的反法截线。21一、相对法截弧某点的纬度愈高，其法线与短轴的交点愈低，当A、B两点位于同一子午圈或同一平行圈上时，正反法截线则合二为一。在通常情况下，正反法截线是不重合的。AB方向在不同的象限时，正反法截线关系如图：22二、大地线的定义和性质大地线:大地线是一条空间曲面曲线，是椭球面上两点间的最短线。大地线上每点的密切面(无限接近的三个点构成的平面)都包含该点的曲面法线，大地线上各点的主法线与该点的曲面法线重合。不在同一子午圈或同一平行圈上的两点的正反法裁线是不重合的，它们之间的夹角△；大地线是两点间惟一最短线，而且位于相对法截线之间，并靠近正法截线。23二、大地线的定义和性质在一等三角测量中，数值可达干分之一二秒，可见在一等或相当于一等三角测量精度的工程三角测量中是不容忽略的。大地线与法截线长度之差只有百万分之一毫米，所以在实际计算中，这种长度差异总是可忽略不计的。24三、大地线的微分方程和克莱劳方程设P为大地线上任意一点，其经度为Ｌ，纬度为Ｂ，大地线方位角为Ａ。当大地线增加dＳ到P1点时，则上述各量相应变化dL，dB及dA。所谓大地线微分方程，即表示dL、dB和dA与dS的关系。dS在子午圈上的分量dS在平行圈上的分量MdBpp12BdLNrdLppcos225三、大地线的微分方程和克莱劳方程三角形PP2P1是一微分直角三角形dSBNAdLAdsBdLNdSMAdBAdSMdBcossinsincoscoscos26三、大地线的微分方程和克莱劳方程大地线微分方程BdSNAdAdSBNAdLdSMAdBtansincossincos27三、大地线的微分方程和克莱劳方程由和得由于两边积分，易得克莱劳定理克莱劳定理表明：在旋转椭球面上，大地线各点的平行圈半径与大地线在该点的大地方位角的正弦的乘积等于常数。284.2将地面观测值归算至椭球面将地面观测元素（包括方向和距离等)归算至椭球面。在归算中有两条基本要求：①以椭球面的法线为基准；②将地面观测元素化为椭球面上大地线的相应元素；294.2.1将地面观测的水平方向归算至椭球面将地面观测的水平方向归算至椭球面，进行三差改正：垂线偏差改正标高差改正截面差改正30一、垂线偏差改正在椭球面上则要求以该点的法线为依据。在每一三角点上，把以垂线为依据的地面观测的水平方向值归算到以法线为依据的方向值而应加的改正定义为垂线偏差改正——uu31一、垂线偏差改正以AO方向作为参考方向。以垂线AZ1为准，照准M点得OR1；以法线AZ为准，则得OR。由此可见，垂线偏差对水平方向的影响是(R一R1)，这个量就是垂线偏差。u32二、标高差改正标高差改正:又称由照准点高度而引起的改正。不在同一子午面或同一平行圈上的两点的法线是不共面的。这样，当进行水平方向观测时，如果照准点高出椭球面某一高度，则照准面就不能通过照准点的法线同椭球面的交点，由此引起的方向偏差的改正叫做标高差改正。h33二、标高差改正标高差改正的计算公式h34三、截面差改正在椭球面上，纬度不同的两点由于其法线不共面，所以在对向观测时相对法截弧不重合，应当用两点间的大地线代替相对法裁弧。这样将法裁弧方向化为大地线方向应加的改正叫截面差改正。s35三、截面差改正s364.2.2将地面观测的长度归算至椭球面一、基线尺量距的归算1.垂线偏差对长度归算的影响37一、基线尺量距的归算2．高程对长度归算的影响如果将上式展开级数，取至二次项，有38一、基线尺量距的归算地面基线长度归算到椭球面上长度的公式)()(12211021HHuuRHSSm39二、电磁波测距的归算电磁波测距仪测得的长度是连接地面两点问的直线斜距，也应将它归算到参考椭球面上40二、电磁波测距的归算将上式按反正弦函数展开级数，舍去五次项，得414.3大地测量主题解算概述4.3.1大地主题解算的一般说明大地元素:大地经度L、大地纬度B、两点间的大地线长度S及其正反大地方位角A12、A21。大地主题解算：如果知道某些大地元素推求另一些大地元素，这样的计算问题就叫大地主题解算，大地主题解算有正解和反解。424.3.1大地主题解算的一般说明大地主题正解：已知Pl点的大地坐标(L1，B1)，P1至P2的大地线长S及其大地方位角A12，计算P2点的大地坐标(L2，B2)和大地线S在P2点的反方位角A21，这类问题叫做大地主题正解。434.3.1大地主题解算的一般说明大地主题反解：如果已知P1和P2点的大地坐标(L1，B1)和(L2，B2)，计算P1至P2的大地线长S及其正、反方位角A12和A21这类问题叫做大地主题反解。444.3.1大地主题解算的一般说明大地主题解算根据不同理论基础的分五类：1．以大地线在大地坐标系中的微分方程为基础，直接在地球椭球面上进行积分运算。2．以白塞尔大地投影为基础。3．利用地图投影理论解算大地问题。4．对大地线微分方程进行数值积分的解法。5．依据大地线外的其他线为基础。454.3.1大地主题解算的一般说明1．以大地线在大地坐标系中的微分方程为基础，直接在地球椭球面上进行积分运算。大地线微分方程BNBdSdABNAdSdLMAdSdBsincoscossincos462．以白塞尔大地投影为基础白塞尔大地主题解算的步骤：1)按椭球面上的已知值计算球面相应值，即实现椭球面向球面的过渡；2)在球面上解算大地问题；3)按球面上得到的数值计算椭球面上的相应数值，即实现从圆球向椭球的过渡473．利用地图投影理论解算大地问题。如在地图投影中，采用椭球面对球面的正形投影和等距离投影以及椭球面对平面的正形投影(如高斯投影)，它们都可以用于解算大地主题，这类解法受距离的限制，只在某此特定情况下才比较有利。4．对大地线微分方程进行数值积分的解法。这种解法直接进行数值积分计算以解决大地主题的解算。常用的数值积分算法有高斯法，龙格－库塔法．牛顿法以及契巴雪夫法等。这种算法易于编写程序，适用任意长度距离。其缺点是随着距离的增长，计算工作量大，且精度降低，而在近极地区，这种方法无能为力。485．依据大地线外的其他线为基础。连接椭球面两点的媒介除大地线之外，当然还有其他一些有意义的线，比如弦线、法截线等。利用弦线解决大地主题实质是三绝大地切量问题，由电磁波测距得到法截线弧长。所以对三边测量的大地主题而言，运用法截弧进行解法有其优点。当然，这些解算结果还应加上归化至大地线的改正。494.3.2勒让德级数式在过已知点P1(L1，B1)且在该点处大地方位角为A12的大地线S上任意一点P2的大地坐标(L2，B2)及其方位角A21必是大地线长度S的函数S＝0时，这些函数值等于P1点的相应数值因此，可在已知点P1点(S＝0)上，按马克劳林公式将Pl和P2点的纬度差、经度差及方位角之差展开为大地线长度S的幂级数。504.3.3高斯平均引数正算首先把勒让德级数在P1点展开改在大地线长度中点M展开，以便级数公式项数减少，收敛快，精度高；其次，考虑到求定中点M的复杂性，将M点用大地线两端点平均纬度及平均方位角相对应的m点来代替，并借助迭代计算，便可顺利地实现大地主题正解。514.3.4高斯平均引数反算大地主题反算是已知两端点的经、纬度L1，B1及L2，B2，反求两点间的大地线长度S及正、反大地方位角A12和A21。这时，由于经差△L、纬差△B及平均纬度Bm均为已知，故可依正算公式很容易地导出反算公式。524.3.5白塞尔大地主题解算方法将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上，继而在球面上进行大地主题解算，最后再将球面上的计算结果换算到椭球面上。由此可见，这种方法的关键问题是找出椭球面上的大地元素与球面上相应元素之间的关系式。同时也要解决在球面上进行大地主题解算的方法。53本章小结1.地球椭球的几何性质。几个重要概念：法截线、子午圈曲率半径、卯酉圈曲率半径、大地线、相对法截弧。2.地面观测值归算到椭球面的原理及过程：方向归算、长度归算。3.大地测量主题解算方法。54思考题：1、分别解释大地测量主题解算正解及反解?2、说明以大地线微分方程为基础的大地测量主题解算的步骤。3、说明白塞尔大地主题解算方法的基本思想。4、白塞尔大地测量主题解算的三个投影条件是什么？55下课了。。。休息一会儿。。。