椭球面上的测量计算(IV).ppt

《椭球面上的测量计算(IV).ppt》由会员分享，可在线阅读，更多相关《椭球面上的测量计算(IV).ppt（67页珍藏版）》请在三一办公上搜索。

1、7-1 地球椭球的几何参数及其相互关系7-2 椭球面上的常用坐标系及其相互关系7-3 椭球面上的几种曲率半径7-4 椭球面上的弧长计算7-5 大地线7-6 将地面观测值归算至椭球面7-7 大地测量主题解算简介（*）7-8 椭球面上三角形的解算（增加）,第七章 椭球面上的测量计算（教材第六章“GPS卫星定位技术”将在另一门专业课中介绍）,在不规则的地球自然表面上获得的距离、角度等观测数据必须先归算到规则的椭球面上，再按保角投影方法投影到高斯平面上才能进行平差计算。等级控制点的坐标都是高斯平面直角坐标系中的坐标。如果不对地面观测数据进行归算、改化就直接按坐标公式计算，则当范围较大时控制网无法拼接。

2、如公路测量等。归算和改化工作分两步进行。不难理解，椭球体实际上只是一个过渡体。在第一章中已经简介过参考椭球体的有关概念和参数。本章将比较系统、详细地介绍椭球体的参数、坐标系以及在椭球面上的测量计算问题。椭球面上的测量计算公式很多。因时间有限，不一定一一推导。课堂上讲过的主要公式，未推导部分请同学们课后尽量自学。,地球自然表面,参考椭球面,高斯平面,7-1 地球椭球的几何参数及其相互关系,一、地球椭球的基本几何参数地球椭球：在控制测量中用来代表地球的椭球，它是地球的数学表达模型。参考椭球：具有一定的几何参数，经过定位和定向，用以代表某一地区大地水准面的地球椭球叫做参考椭球。地面上一切观测元素都应

3、归算到参考椭球面上，并在该面上进行计算。参考椭球面是大地测量计算的基准面。,椭球体有关元素O为椭球中心；NS为旋转轴；EAE为赤道（NS）；a为长半轴；b为短半轴；NKAS为子午圈（经圈或子午椭圆）；QKQ 为平行圈（纬圈，与赤道平行）。决定椭球形状和大小的五个参数长半轴（长度元素）：a短半轴（长度元素）：b扁率（反映椭球体的扁平程度）：第一偏心率：第二偏心率：e和e是子午椭圆的焦点离开中心的距离与椭圆长短半径之比，它们也能反映椭球体的扁平程度。偏心率越大，椭球愈扁。,五个参数中，知道其中的两个就可决定椭球的形状和大小，但其中至少应有一个是长度元素（如a或b）。习惯上通常用a和。椭球面上的测量

4、计算公式很多。为简化书写，引入下列符号：式中B为大地纬度；W、V为辅助函数，其中W叫第一基本纬度函数，V叫第二基本纬度函数。自1738年布格（法国）推算出第一组椭球参数以来，200多年来各国大地测量工作者根据某一国或某一地区的资料，求出了数目繁多、数值各异的椭球参数，比较著名的就有30多组。,我国1954年北京坐标系应用的是克拉索夫斯基椭球参数，1980年西安坐标系应用的是1975年国际椭球参数，而GPS应用的是WGS-84椭球参数。涉及我国的这三组参数值见表71。,经国务院批准，根据中华人民共和国测绘法，我国自2008年7月1日起启用2000国家大地坐标系（简称“2000坐标系”，或CGS2

5、000，即China Geodetic System 2000，或CGCS2000，即China Geodetic Coordinate System 2000）。2000国家大地坐标系是一个地心坐标系，原点与地球质心重合。2000国家大地坐标系采用的地球椭球参数如下：,二、椭球参数间的关系（运用初等数学知识很容易推导，过程略）,7-2 椭球面上的常用坐标系及其相互关系,椭球面上的坐标系有：大地坐标系空间直角坐标系子午平面直角坐标系大地极坐标系地心纬度坐标系归化纬度坐标系前两种是大地测量中最常用的基本坐标系，而最后两种坐标系在实际应用中很少采用，仅用于大地测量公式推导和某些特殊测量计算中。,一

6、、椭球面常用坐标系1、大地坐标系过P点的子午面NPS与起始子午面NGS之间的夹角叫做P点的大地经度L，经度有东经、西经之分；P点的法线Pn与赤道面的夹角叫P点的大地纬度B，纬度有北纬、南纬之分。若P点不在椭球面上，还需增加第三参数大地高H。大地高与正高及正常高的关系为：在大地坐标系中，P点的位置用（L，B，H）表示。,2、空间直角坐标系以椭球中心O为原点，起始子午面与赤道面交线为X轴，在赤道面上与X轴正交的方向为Y轴，椭球体旋转轴为Z轴。X、Y、Z构成右手坐标系。在空间直角坐标系中，P点的位置用（X，Y，Z）表示。,3、子午面直角坐标系设P点的大地经度为L。在过P点的子午面上，以子午椭圆中心为

7、原点，建立平面直角坐标系（注意：该平面直角坐标系的x、y轴的方向与测量习惯不一致）。在该坐标系中，P点的位置用（L，x，y）表示。子午面直角坐标系仅用于大地测量公式推导。,4、大地极坐标系M为椭球面上任意一点。MN 为过M 点的子午线。S 为连结M 点和P点的大地线（最短曲线，将在75介绍）长度。A 为大地线在M 点的方位角。以M 为极点，MN 为极轴，S 为极径，A为极角，便构成大地极坐标系。在大地极坐标系中，P点的位置用（S,A）表示。,二、几种坐标系之间的关系1、子午面直角坐标与大地坐标的关系在这两个坐标系中，同点的L相同。过点P作子午椭圆的法线，该法线与x 轴、y 轴的交点分别为Q、n

8、，与x 轴的夹角为B（P点纬度）。可以证明（过程略），子午面直角坐标（x，y）与大地纬度B 有以下关系式：式中，a、b 为椭球长、短半径，W、V为前述辅助函数：,设Pn=N，由右图可以看出：xNcosB(7-18)比较(7-16)式，有：Na/W(7-19)于是 y N(1e2)sinB(7-20)又由图可知：yPQsinB(7-21)所以：PQN(1e2)(7-22)QnNPQNe2(7-23)由(7-22)、(7-23)可知P点法线Pn在赤道两侧的长度。,N,2、空间直角坐标与子午面直角坐标的关系从右边2图容易得出：,3、空间直角坐标与大地坐标的关系（1）由大地坐标（L，B，H）求空间直角

9、坐标（X，Y，Z）当P点在椭球面上，即H 0时，将(7-18)、(7-20)二式代入(7-24)式即得：当P点不在椭球面上，即H 0时，有,3、空间直角坐标与大地坐标的关系（2）由空间直角坐标（X，Y，Z）求大地坐标（L，B，H）关系式如下：（731）可直接由（725）得到。（732）可根据右图得到。OPx因等式右边也包含B，故需迭代计算，其初始值可设为0;N值也需逐次迭代。同样，（734）可根据右图得到。sinBZ/(H+PQ)教材7.2.3“地心站平坐标系”实际应用较少。,7-3 椭球面上的几种曲率半径,为了进行控制测量计算，必须了解椭球面上曲线的曲率半径。先介绍两个名词。法截面过椭球面上

10、的一点作垂直于椭球面的法线，包含这条法线的平面。法截线法截面与椭球面的交线。法截线又叫法截弧。由一条法线可作无数个法截面，相应有无数条法截线。在圆球面上，各点的法截线形状相同，曲率半径也相同。在椭球面上，不仅各点的法截线形状不同、曲率半径不同，而且同一点上不同方向的法截线的形状及曲率半径也各不相同。下面介绍椭球面上几种主要曲线的曲率半径。,一、子午圈曲率半径在子午椭圆上取微分弧长dS（dS=DK），相应地有坐标增量(子午面直角坐标系)dx、dy。设微分弧dS的曲率中心为n，则线段Dn及Kn即子午圈曲率半径，用M表示。由平面曲线的曲率半径定义公式知：进一步推导（过程略），可得：式中：,M,子午圈

11、曲率半径与纬度的大小有关。其变化规律如表72所示：M随纬度B的变化情况亦可从下图中看出。,二、卯酉圈曲率半径 过椭球面上P点可作无数个法截面。其中与子午面垂直的那个法截面同椭球面的交线称为卯酉圈。右图中的EPE即为过P点的卯酉圈，其曲率半径用N表示。过P点作以O为中心的平行圈PHK的切线PT，该切线位于垂直于子午面的平行圈平面内。卯酉圈EPE也垂直于子午面，故PT也是卯酉圈在P点处的切线，垂直于法线Pn。所以PT是平行圈PHK及卯酉圈EPE在P点处的公切线。卯酉圈EPE是子午椭圆上P点处的法截弧，而平行圈PHK可看作子午椭圆上P点处的“斜截弧”。,麦尼尔定理：如果过曲面上一点引两条截弧，一为法

12、截弧，一为斜截弧，且在该点处两截弧具有公切线，则斜截弧在该点的曲率半径等于法截弧的曲率半径乘以两截弧平面夹角的余弦。由右图可知，两截弧平面夹角为B。斜截弧（平行圈PHK）的曲率半径r等于P点的横坐标x（子午面直角坐标）。设法截弧（EPE）的曲率半径为N，则根据麦尼尔定理，有：在直角三角形PO n中，显然有：因此，N=Pn。说明卯酉圈的曲率中心位于椭球的旋转轴上，曲率半径等于椭球面与短轴之间的一段法线的长度。,将上节导出的子午面直角坐标与大地坐标的关系式（716）代入（762），容易得到：（764）式便是卯酉圈曲率半径计算公式。由（78）式中的第3和第7个关系式容易得到N的另一个计算式：,卯酉圈

13、曲率半径N也随纬度B而变,变化规律如表 73。,M和N是两个互相垂直的法截弧（子午圈和卯酉圈）的曲率半径，称为主曲率半径。计算主曲率半径的公式（760）和（764）式应熟记。,三、任意法截弧的曲率半径椭球面上P点处的子午法截弧与卯酉法截弧正交。子午法截弧的方位角为卯酉法截弧的方位角为前面已导出了M、N的计算公式。本段讨论方位角为的任意法截弧的曲率半径的计算问题。根据尤拉公式（又叫欧拉公式）知：由曲面上任意一点主曲率半径计算该点任意方位角A的法截弧的曲率半径RA的公式为：由（760）、（764）和（79）式可得：将（783）式右边分子分母同除以M，并顾及（784）式，有,0或180,90或270

14、,推导过程中利用了（74）中的关系式：该式为计算方位角为A的任意方向法截弧的曲率半径的理论公式（785）。将上式按级数展开，并经简化后得实用公式：式中R为P点的平均曲率半径（下段介绍）。R和可用计算器编程计算或查表得到。,四、平均曲率半径任意方向的法截弧曲率半径随方位而变，给测量计算带来不便。当精度要求不太高、范围不太大时，可将椭球面当作球面处理。找出测区的中心点。将中心点处的平均曲率半径（即所有法截弧的曲率半径的算术平均值的极限）作为测区近似球面的半径R。其计算公式为：（793）式说明，某点的平均曲率半径R等于该点M和N的几何平均值。,五、M、N、R的关系由前面介绍的公式可知，一般情况下，椭

15、球面上某点处的M、N、R值不相等：所以它们有如下关系：N R M 只有在极点上，它们才相等，且均等于极曲率半径c。,7-4 椭球面上的弧长计算,控制测量计算中需要用到子午线弧长及平行圈弧长。本节介绍这两种弧长的计算公式。,一、子午线弧长计算公式在子午线上取某微分弧PPdx设：P点的纬度为B，P点的纬度为 BdB，P点的子午线曲率半径为M。显然，有：从赤道开始到任意纬度B的子午线弧长为：式中M可用下式表达（由（770）式变形而得）：,不难发现，系数a0a8都是椭圆长半径a和第一偏心率e的函数。椭球元素确定时，它们都是常数；椭球元素不同时，a0a8之值也会不同。,经积分、整理后，得子午线弧长计算式

16、：式中的系数很有规律。对于克拉索夫斯基椭球，子午线弧长计算公式为：式中B表示以度为单位的纬度。对于1975年国际椭球，子午线弧长计算公式为：求子午线上两点间的弧长时，只需分别将两点的纬度值代入弧长公式，算出相应的弧长后相减即得。反过来，根据子午线弧长X(高斯平面直角坐标系中的纵坐标)也可以反算大地纬度B。反算得到的纬度叫“底点纬度Bf”（高斯平面上过已知点向纵坐标轴作垂线与纵坐标轴交点的大地纬度），一般采用迭代法和直接解法。（详见）,二、平行圈弧长公式椭球体的每一个平行圈都是圆，其半径就是圆上任意一点的子午面直角坐标x，即：对于同一平行圈上经差为l 的两点，其弧长为：,三、子午线弧长和平行圈弧

17、长变化的比较表76列出了在不同的纬度处，纬度每变化1、1、1时子午线弧长的变化情况以及经度每变化1、1、1 时平行圈弧长的变化情况。由表可以看出，单位纬差的子午线弧长随B的增大而缓慢地增大；而单位经差的平行圈弧长则随B的增大而急剧缩短。由表还知道：纬差1对应的子午线弧长约为110KM，1约为1.8KM，1约为30M；在15以下的低纬度地区，平行圈弧长与经差的关系也大致如此。（测量技术人员必须有此概念）,7-5 大地线,两点间的最短距离，在平面上是两点间的直线；在球面上是两点间的大圆弧。在椭球面上呢？是大地线。本节介绍大地线的有关概念。先引入“相对法截线”。,一、相对法截线设椭球面上的任意两点A

18、、B，其纬度分别为B1、B2，且B1 B2。过A、B两点分别作法线与椭球短轴交于na、nb点。由73知，na、nb分别为A、B两点的 卯酉圈曲率中心。因为OniNie2sinBi,所以当B1 B2 时，Ona Onb，亦即na、nb两点不重合。由此可见：（1）纬度不同的两点，法线交于椭球旋转轴上的不同点；（2）椭球面上点的纬度愈高，法线与旋转轴的交点就愈低；（3）当两点的纬度不同，又不在同一子午圈上时，这两点的法线将在空间交错而不相交。（4）既不在同一子午圈上也不在同一平行圈上的两点之间有两条法截线。,过A点作一个法截面通过B点，过B点作一个法截面通过A点。这两个法截面不重合，其与椭球面的交线

19、也不重合。,假设经纬仪的竖轴分别与A、B两点的法线重合（忽略垂线偏差），以其中的一点为测站照准另一点时，经纬仪的照准面就是一个法截面。在A点照准B点时，照准面A na B与椭球面的截线A aB叫做A点的正法截线，或B点的反法截线；同理，在B点照准A点时，照准面B nb A与椭球面的截线B bA叫做B点的正法截线，或A点的反法截线。A aB和B bA叫做A、B两点的相对法截线，二者一般不重合。只有当A、B两点位于同一子午圈上（经度相同）或同一平行圈上（纬度相同）时，正反法截线才会合二为一，彼此重合。需要注意的是，同一平行圈上的两点，正反法截线虽然重合，但它们并不在平行圈上。（北半球，法截线在平行

20、圈上方）,AB方向位于不同的象限时，正、反法截线之间的关系如右图所示。（某点的正法截线总是凹向该点处的子午线）图中实线表示仪器端，虚线表示目标端。一般情况下，正反法截线不重合。因此，在椭球面上A、B、C三点处所测得的角度（各点上正法截线之夹角）将不能构成闭合三角形，如右下图所示。为消除这种矛盾，需在相邻两点间另选一条单一截线代替相对法截线，从而得到由单一截线构成的闭合三角形。这种单一截线就是大地线。,二、大地线的定义及其性质先从空间曲线（不限于在同一平面内的曲线）的有关概念入手。图中，ADB为椭球面上的一条空间曲线。切线通过曲线上的D点和与它邻近的一点D1作一条直线D1D，当D1无限趋近D点时

21、，直线D1D的极限位置，就是该曲线在D点处的切线。密切平面通过曲线上的D点和与它邻近的两点D1和D2作一个平面D1D D2（过三点可作一平面），当D1和D2同时无限趋近D点时，平面D1D D2的极限位置就是该曲线在D点处的密切平面。,主法线通过曲线上的D点，垂直于曲线在该点的切线可以作无数条直线，这些直线都叫做曲线在D点的法线。其中位于密切平面内的一条法线叫做主法线。曲面的法线过曲面上的D点可以作一个切平面。通过该点且垂直于切平面的直线叫做曲面在该点的法线。椭球面上某点的法线是该点所有法截弧的公法线，但各法截弧的曲率中心的位置不同，即曲率半径不同。,大地线对于曲面上的一条曲线，如果其上每一点处

22、的密切平面都包含曲面在该点的法线，则这条曲线称为该曲面上的一条大地线。大地线又叫测地线。例如，椭球面上的子午线就是一条大地线，因为其上每一点处的密切平面都包含椭球面在该点的法线。但是，除赤道以外的平行圈却不是大地线，因为平行圈上每一点处的密切平面（赤道的平行面）不包含椭球面在该点的法线。大地线上各点的主法线与该点的曲面法线重合。因各点的曲面法线一般在空间交错，不相交、不平行，所以大地线是一条“空间曲面曲线”（不在同一平面）。,假设在椭球面上的A、B两点间拉紧一条弹性细绳，细绳无重量，与椭球面之间也没有摩擦力。则当细绳平衡静止时，其形状就是大地线的形状。由弹性力学可知，细绳平衡时，其上每一点的弹

23、力的合力将位于该点的密切平面内；而曲面对细绳上各点的反作用力则必然与曲面在该点的法线一致。因细绳处于平衡静止状态，故作用力与反作用力必在同一直线上，大小相等方向相反。意味着细绳上每一点处的密切平面包含了该点的曲面法线。故平衡状态下的细绳形状是大地线。因细绳是在两端有拉力的情况下平衡的，距离最短，故大地线是椭球面上两点间的最短曲线。A、B两点的大地线位于其相对法截线之间。如右上图。,在子午圈上，大地线与相对法截线重合，都位于子午线上。在除赤道以外的平行圈上，尽管大地线也与相对法截线重合，但它们都不在平行圈上。（北半球偏上）,不在同一子午圈或同一平行圈上的两点的正反法截线一般不重合。它们之间的夹角

24、用表示。在一等三角测量中，平均边长为25公里，约为0.004。大地线位于相对法截线之间，并靠近正法截线。大地线与正法截线间的夹角用表示：因和都很小，所以在一般的工程测量中，可以不考虑法截线与大地线之间的方位角差，即在方向观测值中无需施加此项改正。大地线与法截线的长度大约每公里相差百万分之一米，完全可以认为它们相等。在椭球面上进行测量计算时，以两点间的大地线为依据。,三、大地线的微分方程和克莱洛方程*（该部分内容主要用在椭球大地测量学中，用于大地主题解算。有兴趣的同学请自学）,7-6 将地面观测值归算至椭球面,我们知道，测量计算的基准面是参考椭球面。野外测量工作都是在地面上进行的，测站点和照准点

25、一般都超过参考椭球面一定高度。观测时的基准线不是各点相应的椭球面的法线，而是铅垂线。各点的垂线与法线间存在着垂线偏差。地表有较大起伏，故不能直接在地面上处理观测成果，须将地面观测的元素归算至椭球面上。归算时有两个基本要求（1）以椭球面的法线为基准；（2）将地面观测元素化为椭球面上大地线的相应元素。本节将介绍方向值和长度值的归算问题。,一、将地面观测的水平方向值归算至椭球面-三差改正将水平方向值归算至椭球面，需要进行三差改正，即垂线偏差改正、标高差改正和截面差改正。1、垂线偏差改正 地面上所有水平方向的观测都是以垂线为依据的，而在椭球面上则是以法线为依据。将以垂线为依据的地面观测的水平方向值归算

26、为以法线为依据的方向值时应加的改正叫垂线偏差改正，常用u表示。,M是地面观测目标（m）在球面上的投影。延长ZZ1，与过A点的水平圆相交于O点。如果AM位于垂直面AZZ1O内，则无论观测方向是以法线为准还是以垂线为准，照准面都是一个，没有垂线偏差改正。故可将 A O作为参考零方向。以垂线AZ1为准时，照准m点，得方向值OR1（R1是平面AZ1M与过A点的水平圆的交点。过A点的水平圆可看作水平度盘）；以法线AZ为准时，照准m点，得方向值OR（R是平面AZM与过A点的水平圆的交点）。,以测站A为中心作单位半径的辅助空心圆球。Z为对应于法线的大地天顶，Z1为对应于铅垂线的天文天顶。设u为垂线偏差（在辅

27、助球面上，uZZ1）。它在子午圈和卯酉圈上的分量分别为和。,水平方向观测值的垂线偏差改正计算公式为：式中，Am和分别为测站点至照准点的大地方位角和垂直角。垂线偏差的子午分量和卯酉分量可在测区“垂线偏差分量图”中查得。在地势平坦地区，和都很小。,很显然，垂线偏差对水平方向观测值的影响是OR与OR1之差，即：u RR1以上是垂线偏差的几何解释。,2、标高差改正前已述及，不在同一子午面或同一平行圈上的两点，法线在空间交错，不共面。如果照准点高出椭球面一定高度，则照准面就不能通过照准点的法线与椭球面的交点，由此引起的水平方向改正叫标高差改正，用 h表示。假设测站点A位于椭球面上，照准点B高出椭球面的高

28、度为H2；Ana、Bnb分别为A点、B点的法线；b为B点的法线与椭球面的交点。如右图所示。因Ana与 Bnb不在同一平面，故在A点照准B点时的法截线是Ab而不是Ab。将观测值Ab归算为正确的方向值Ab时需加的改正数就是标高差改正。标高差改正的计算公式为：H2、B2和M2分别为照准点的大地高、大地纬度和子午圈曲率半径；Am为测站点至照准点的大地方位角。,3、截面差改正前已述及，当两点的纬度不同时，测站点和照准点的相对法截线不重合，需用大地线代替相对法截线。将观测得到的法截线方向值化算为大地线方向值时应加的改正数叫截面差改正，用g表示。上节已介绍，当平均边长为25公里时，g约为0.0013，故在二

29、等以下的控制测量中一般不考虑施加截面差改正。截面差改正的计算公式见（7167）式。,4、施加三差改正的条件水平方向观测值经三差改正后，便是椭球面上相应大地线的方向值。但是，并非任何情况下都需对观测结果施加三差改正。规范规定，内业计算时水平方向值的取位与控制测量的等级有关 一等计算至0.001 二等计算至0.01 三等和四等计算至0.1如果三差改正数中的某一项小于相应的计算精度（取位精度），则可不考虑该项改正。对三差改正数的分析结果表明，一般情况下：一等三角测量应加三差改正；二等三角测量应加垂线偏差改正和标高差改正；三等和四等三角测量可不加三差改正，但当、超过10或H2000m时，则应分别考虑施

30、加垂线偏差改正和标高差改正。,二、天顶距观测值的垂线偏差改正（教材中无此部分内容）为了进行距离改化等椭球面上的计算，需要用到大地高高差。采用三角高程测量时，高差是根据以垂线为依据的天顶距观测值计算的，不是大地高高差。为获得大地高高差，须将以垂线为依据的天顶距观测值归算为以法线为依据。在右上图中,p1、p2为两个三角点；z1、z2分别为在p1、p2观测的天顶距；z1、z2为归化后以法线为依据的天顶距，即大地天顶距；u1、u2是两点的垂线偏差。两种天顶距之间的关系为：ui,j为测站i的垂线偏差ui在测线ij方向上的分量。由右图知：,所以，计算大地天顶距的公式为：上式与第五章介绍的公式在符号上是一致

31、的：此处是天顶距！垂线偏差子午分量和卯酉分量可根据天文经纬度（、）和大地经纬度（L、B）计算,也可从测区垂线偏差分量图上查取。=-B；=（-L）cos 式中A为测站点至照准点的大地方位角。根据大地天顶距求得的高差，为大地高高差。,三、将地面观测的长度归算至椭球面在地球自然表面上测得的长度需归算至参考椭球面上。控制测量中，测边方法有两种：铟钢基线尺量距和电磁波测距。上世纪80年代之前，铟钢基线尺量距是精密测距的主要手段。随着电磁波测距仪的普遍应用，现在已基本不采用基线尺量距这一方法。电磁波测距的归算问题在上册第四章中已经介绍。本段再简单回顾一下。,如右图所示，地面点Q1、Q2的大地高分别为H1、

32、H2，用电磁波测距仪测得两点间的斜距为D；Q1、Q2沿各自的法线投影到椭球面上的位置分别为Q1、Q2；S为椭球面上Q1、Q2之间的大地线长度。所谓长度归算就是根据斜距D求S。在工程测量中，可以认为大地线长度S等于Q1点处方位角为A12的法截弧曲率半径RA所对应的圆弧长。经推导，电磁波测距归算的严密公式为：式中：,一般情况下，（7177）式右边第二项非常小，D为5km时才0.13 mm，故可忽略不计。,（7177）式跟上册第四章（4141）式是一致的。区别在于一为弧长，一为弦长。【（7177）式比（4141）式多了一项。】电磁波测距归算的实用公式为：式中：精度要求很高时，上式可取至二次项,即,为

33、了解释电磁波测距归算的几何意义，（7179）式可以进一步化简为：上式右边第二项是用高差进行倾斜改正的主项 经过此项改正，斜距已化算为平距；第三项是由于两点的平均高程面高出参考椭球面应进行的归算改正，经过此项改正后，测线已化算为椭球面上的弦长d；第四项则是由弦长改化为弧长的改正项。该式没有（7177）式严密，也比（7179）式要差。实际工作中，应尽量采用（7177）式或（7179）式进行归算。,7-7 大地测量主题解算简介,一、大地主题解算的概念大地主题解算：根据某些已知大地元素推求另一些未知大地元素。大地元素：椭球面上点的大地经度L、大地纬度B，两点间的大地线长度S及其正、反大地方位角A12、

34、A21（二者之差不是刚好等于180），统称为大地元素。大地主题解算包括正算和反算。,已知P1点的大地坐标（L1、B1）、P1至P2点的大地线长S及大地方位角A12，计算P2点的大地坐标（L2、B2）和P2点处大地线的方位角A21的过程叫做大地主题正解（算）。已知P1和P2点的大地坐标（L1、B1）和（L2、B2），计算P1至P2点的大地线长S及正、反大地方位角A12和A21的过程叫做大地主题反解（算）。大地主题正、反算与平面直角坐标的正、反算类似，只是多了反方位角的计算。从解析意义来讲，大地主题解算就是研究大地极坐标（S，A）与大地坐标(L,B)间的相互变换关系。,二、大地主题解算的用途在天文

35、大地测量中推算一等点的经纬度及点间大地方位角；在空间技术和航空、航海、国防等领域进行科学计算。现代控制测量计算中已很少进行大地主题解算，相关计算工作均在高斯投影面上进行（见下章）。三、大地主题解算方法分类1、按解算距离分 短距离（400KM以内）、中距离（400-1000KM）和长距离（1000KM以上）三种。距离不同，解算方法也不同。2、按解算方法分直接解法（白塞尔Bessel法等）和间接解法（高斯平均引数公式等）。一百多年来，测量学者共提出了70余种解算公式和方法。教材介绍了其中主要的几种解法，有兴趣的同学可自学。,7-8 椭球面上三角形的解算,一、“椭球面上三角形解算”的实质“椭球面上三

36、角形解算”的实质就是在椭球面上解算三角形，求椭球面上三角形的边长。前面已介绍了将地面三角网的方向值、长度值归算至椭球面上。在椭球面上已得到由大地线组成的“三角网”。椭球面上三角网的各个角度为已知。除起算边外，其余各边的长度未知，需要解算。,二、用球面三角形代替椭球面三角形的条件大地线是空间曲线，其上各点处的曲率半径不相同。因此，在椭球面上解算由大地线组成的三角形非常复杂。数学研究结果表明：在半径小于140km的圆形区域（相当于边长约240km的等边三角形）内，椭球面三角形可当作球面三角形看待边长相等，角度差小于0.001。圆球的半径取三点曲率半径的平均值。所以，用球面三角形代替椭球面三角形的条

37、件是“三角形的边长小于240km”。,三、球面三角形与同边长平面三角形的关系由球面三角学知：球面三角形的三角之和 大于180，小于540。即 180A0+B0+C0540超过180的部分叫“球面角超”，用表示。勒让德（Legendre，法国数学家，1752-1833）定理 如果平面三角形和球面三角形对应边相等，则平面角等于对应球面角减去三分之一球面角超。,如图所示。设球面三角形A0B0C0的三边与平面三角形A1B1C1的三边对应相等，则据勒让德定理，它们的角度有如下关系：球面角超按下式计算：式中，R为圆球半径，F为平面三角形的面积。,在平面三角形中，根据正弦定理解算边长十分容易。,作业四（1）,我国主要采用哪几种参考椭球?其主要参数是什么?简要说明并图示地面某一点的大地高、正高、正常高、大地水准面差距和高程异常。试推证由一点的空间直角坐标（X，Y，Z）求大地坐标（L，B，H）的计算公式。什么叫法截线？什么叫大地线？大地线有何特性？将水平方向值归算至椭球面时需要进行三差改正。试述三差改正的定义，写出相应的计算公式，并说明各自的应用条件。垂线偏差对水平方向观测值和垂直角观测值各有什么影响？写出电磁波测距测得的斜距化算为大地线长度的计算公式，并说明式中各符号的意义。,