1、第四章第四章 地球椭球数学投影变换的基本理论地球椭球数学投影变换的基本理论 4.1地球椭球基本参数及其互相关系地球椭球基本参数及其互相关系 4.2 椭球面上常用坐标系及其关系椭球面上常用坐标系及其关系 4.3 椭球面上的几种曲率半径椭球面上的几种曲率半径 4.4 椭球面上的弧长计算椭球面上的弧长计算 4.5 大地线大地线 4.6 将地面观测值归算至椭球面将地面观测值归算至椭球面 4.7 大地测量主题解算概述大地测量主题解算概述 4.8 地图数学投影变换的基本概念地图数学投影变换的基本概念 4.9 高斯平面直角坐标系高斯平面直角坐标系 4.10 横轴墨卡托投影和高斯投影簇的概念横轴墨卡托投影和高
2、斯投影簇的概念 4.11 兰勃脱投影概述兰勃脱投影概述1、球面三角形、球面三角形(spherical triangle)2、球面角超、球面角超(spherical excess)180CBAn 定义定义2RS 球面三角学的基本知识球面三角学的基本知识 Basic of Spherical TrigonometryOabcBCAn 计算公式计算公式3、球面三角形公式（单位球）球面三角形公式（单位球）n 正弦公式（四元素）正弦公式（四元素）CcBbAasinsinsinsinsinsinBDOBBFBEBFBDOBBECcsinsinBDOBBFBEBEBDOBBFAasinsin单位球单位球)(
3、sinsinsinsinsinsinRCRcBRbARa半径为n 边的余弦公式（四元素）边的余弦公式（四元素）Acbcbacossinsincoscoscosn 角的余弦公式（四元素）角的余弦公式（四元素）aCBCBAcossinsincoscoscosn 正余弦公式（五元素）正余弦公式（五元素）cCBCBcAaBCBCbAAcbcbCaAbcbcBacossincoscossincossincossincoscossincossincossincoscossincossincossincoscossincossinOabcBCA3、球面三角形公式（单位球）球面三角形公式（单位球）n 余切公式
4、（四元素）余切公式（四元素）acCBBAabbCCAACCbbaABBccacotsincoscossincotcotsincoscossincotcotsincoscossincotcotsincoscossincotn 正切公式（四元素）正切公式（四元素）)(21tan)(21tan)(21tan)(21tanbabaBABA3、球面三角形公式（单位球）、球面三角形公式（单位球）OabcBCA球面直角三角形公式的纳白尔规则（单位球）球面直角三角形公式的纳白尔规则（单位球）纳白尔规则纳白尔规则 环形上任一元素的正弦等于：环形上任一元素的正弦等于：1）相邻两元素正切的积；）相邻两元素正切的积；
5、2）相对两元素余弦的积。）相对两元素余弦的积。4.1地球椭球基本参数及其互相关系地球椭球基本参数及其互相关系 地球椭球是选择的旋转椭球,旋转椭球的形状和大小常用子午椭圆的五个基本几何参数(或称元素):长半轴 短半轴 椭圆的扁率 椭圆的第一偏心率 椭圆的第二偏心率 通常用a,aba abae22abae22为简化书写，还常引入以下符号2222,tan,cosactBeBbBeVBeW2222cos1sin1221,11,11,11,12222222222eeVWeWVeeeeeeecaeaceabeba222222222221()1()1sin(1)1(1)bWeVVaaVeWWbWeBe VV
6、eW 4.2 椭球面上常用坐标系及其关系椭球面上常用坐标系及其关系4.2.1 各种坐标系的建立各种坐标系的建立1、大地坐标系、大地坐标系大地经度大地经度B 大地纬度大地纬度L 大地高大地高H2、空间直角坐标系空间直角坐标系定义：1、坐标原点位于总地球椭球(或参考椭球)质心；2、Z轴与地球平均自转轴相重合，亦即指向某一时刻的平均北极点；3、X轴指向平均自转轴与平均格林尼治天文台所决定的子午面与赤道面的交点G；4、Y轴与此平面垂直，且指向东为正。地心空间直角系与参心空间直角坐标系之分。3、子午面直角坐标系子午面直角坐标系 设P点的大地经度为L，在过P点的子午面上，以子午圈椭圆中心为原点，建立x,y
7、平面直角坐标系。在该坐标系中，P点的位置用L,x,y表示。(X,Y,Z)(L,x,y)4、地心纬度坐标系及归化纬度坐标系、地心纬度坐标系及归化纬度坐标系 设椭球面上P点的大地经度L，在此子午面上以椭圆中心O为原点建立地心纬度坐标系地心纬度坐标系;以椭球长半径a为半径作辅助圆，延长与辅助圆相交点，则OP与x轴夹角称为P点的归化纬度归化纬度u。(L,)(L,u)5 5、大地极坐标系、大地极坐标系 M是椭球面上一点，MN是过M的子午线，S为连接MP的大地线长，A为大地线在M点的方位角。以M为极点；MN为极轴；P点极坐标为（S,A)4.2.2 坐标系之间的相互关系坐标系之间的相互关系子午平面坐标系同大
8、地坐标系的关系ctgBBdxdy)90tan(022221(1)xy abyxabdxdy22222c(1)(2)bxxtgBeayyWBaBeBaxcossin1cos22Bexytan)1(2VBbBeWaBeBeaysinsin)1(sin1sin)1(2222WBaBeBaxcossin1cos22令:pn=NcosxNBWaNBeNysin)1(2BPQysin)1(2eNPQ2NeQn 空间直角坐标同子午面直角坐标系的关系空间直角坐标同子午面直角坐标系的关系c o s,s in,XxLYxLZy空间直角坐标系同大地坐标系空间直角坐标系同大地坐标系 在椭球面上在椭球面上的点：的点：不
9、在椭球面上不在椭球面上的点：的点：2coscoscossincossin(1)sinXxLNBLYxLNBLZyNeBnH0BHeNLBHNLBHNZYXsin)1(sincos)(coscos)(2n由空间直角坐标计算相应大地坐标 或2222arccosarcsinarctanYXXLYXYLXYL222sintanYXBNeZBNBYXHcos222(1)sinzHNeB B、u、之间的关系之间的关系 B和u之间的关系 2cos,sinsincos ,(1)sinxau ybuaabBxByeBWWVBWeusin1sin2BWucos1cosuVBsinsinuWBcoscosU、之间的
10、关系、之间的关系 大地纬度、地心纬度、归化纬度之间的差异很小，经过计算，当B=45时xytanuexytan12ue tan1tan2Betan)1(tan28.11)(9.5)(9.5)(maxmaxmaxBuuB uB大地测量数据处理过程大地测量数据处理过程地面观测元素地面观测元素天文方位角天文方位角距离距离方向方向归算归算椭球面上的元素椭球面上的元素大地方位角大地方位角距离距离方向方向高斯平面元素高斯平面元素归算归算距离距离方向方向坐标方位角坐标方位角大地坐标（大地坐标（L，B）平面坐标（平面坐标（x，y）平差平差平差平差4.3 椭球面上的几种曲率半径椭球面上的几种曲率半径 过椭球面上任
11、意一点可作一条垂直于椭球面的法线，包含这条法线的平面叫作 法截面法截面，法截面与椭球面的交线叫法截线法截线。子午圈曲率半径dBdSMBdxdSsinBdBdxMsin1WBaxcos2cossinWdBdWBBWadBdxWBBedBBeddBdWcossinsin1222)1(sin23eWBadBdx23(1)aeMW3VcM卯酉圈曲率半径(N)卯酉圈卯酉圈:过椭球面上一点的法线，可作无限个法截面，其中一个与该点子午面相垂直的法截面同椭球面相截形成的闭合的圈称为卯酉圈。麦尼尔定理麦尼尔定理:假设通过曲面上一点引两条截弧，一为法截弧，一为斜截弧，且在该点上这两条截弧具有公共切线，这时斜截弧在
12、该点处的曲率半径等于法截弧的曲率半径乘以两截弧平面夹角的余弦。BNrcosWBarxcosWaN VcN BrBPONPncoscos卯酉圈曲率半径的特点卯酉圈曲率半径的特点:卯酉圈曲率半径恰好等于法线介于椭球面和短轴之间的长度，亦即卯酉圈的曲率中心位在椭球的旋转轴上。主曲率半径的计算主曲率半径的计算 以上讨论的子午圈曲率半径M及卯酉圈曲率半径N，是两个互相垂直的法截弧的曲率半径，这在微分几何中统称为主曲率半径。23222)sin1)(1(BeeaM2122)sin1(BeaNBmBmBmBmmM886644220sinsinsinsinBnBnBnBnnN886644220sinsinsin
13、sin6284262240222089674523)1(memmemmemmemeam628426224022087654321nennennennenan2322)cos1(BecM2122)cos1(BecNBmBmBmBmmM886644220coscoscoscosBnBnBnBnnN886644220coscoscoscos1011)(89674523)1(/821062842622402220memmemmemmemmemeacm109)(876543211/821062842622402220nennennennenneneacn任意法截弧的曲率半径任意法截弧的曲率半径NAMAR
14、A22sincos1AMANMNRA22sincos21VMNABeNANRA2222coscos1cos1)coscos1(4422AANRA任意法截弧的曲率半径的变化规律:不仅与点的纬度B有关，而且还与过该点的法截弧的方位角A有关。当时，变为计算子午圈曲率半径的，即；当90时，为卯酉圈曲率半径，即。主曲率半径M及N分别是的极小值和极大值。当A由090时，之值由，当A由90180时，值由N，可见值的变化是以90为周期且与子午圈和卯酉圈对称的。平均曲率半径平均曲率半径 椭球面上任意一点的平均曲率半径 R 等于该点子午圈曲率半径M和卯酉圈曲率半径N的几何平均值。MNR 22221eWaVNVcW
15、bRM，N，R的关系 MRNcMRN909090对于克拉索夫斯基椭球4.4 椭球面上的弧长计算椭球面上的弧长计算子午线弧长计算公式子午线弧长计算公式 MdBdx BMdBX0BmBmBmBmmM886644220sinsinsinsinBBBBBBBBBBBBBB8cos12816cos1614cos3272cos16712835sin6cos3214cos1632cos3215165sin4cos812cos2183sin2cos2121sin8642BaBaBaBaBaX8sin86sin64sin42sin286420BaBaBaBaaM8cos6cos4cos2cos864201281
16、6323271638167321522128351653288866864486422864200mammammmammmmammbmmma如果以B90代入，则得子午椭圆在一个象限内的弧长约为10 002 137m。旋转椭球的子午圈的整个弧长约为40 008 549.995m。即一象限子午线弧长约为10 000km，地球周长约为40 000km。为求子午线上两个纬度B及间的弧长，只需按(11.42)式分别算出相应的X及X，而后取差：，该即为所求的弧长。当弧长甚短(例如X40km，计算精度到0.001m)，可视子午弧为圆弧，而圆的半径为该圆弧上平均纬度点的子午圈的曲率半径M由子午弧长求大地纬度
17、迭代解法:平行圈弧长公式 01/aXBf01/)(aBFXBifififififififBaBaBaBaBF8sin86sin64sin42sin2)(8642cos1lblBNS子午线弧长和平行圈弧长变化的比较子午线弧长和平行圈弧长变化的比较四、梯形图幅面积四、梯形图幅面积2、梯形图幅面积、梯形图幅面积2121cosBBLLBdBdLMNPBdBdLMNdPcos椭球面积元图幅ABCD的椭球面积BeW22sin132)1(,WeaMWaN212121222222222)sin1(cos cos)sin1()1(BBBBLLdBBeBlbBdBdLBeeaP212121222222222)si
18、n1(coscos)sin1()1(BBBBLLdBBeBlbBdBdLBeeaPmmmmBBBDBBBCBBBBBBBAblP7cos)(27sin5cos)(25sin3cos)(23sincos)(21sin901212121229cos)(29sin12mBBBE得略去具体推导，最后可为椭球短半径，为图幅的大地经度之差其中，bl)(2121BBBm8868648642864223045 25651121 645161803 1923516516361 1283516583211eeeeeeeeeeeeeeEDCBA4.5 大地线大地线 两点间的最短距离，在平面上是两点间的直线，两点间的
19、最短距离，在平面上是两点间的直线，在球面上是两点间的大圆弧，那么在椭球面上又是怎样在球面上是两点间的大圆弧，那么在椭球面上又是怎样的一条线呢的一条线呢?它应是大地线。它应是大地线。相对法截线相对法截线2211sinsinBnQOnBnQOnbbaa222121sinsinBeNOnBeNOnba相对法截线相对法截线相对法截线的特点相对法截线的特点:当当A，B两点位于同一子午圈或同一平行圈上两点位于同一子午圈或同一平行圈上时，正反法截线则合二为一。时，正反法截线则合二为一。在通常情况下，正反法截线是不重合的。因此在通常情况下，正反法截线是不重合的。因此在椭球面上在椭球面上A，B，C三个点处所测得
20、的角度三个点处所测得的角度(各各点上正法截线之夹角点上正法截线之夹角)将不能构成闭合三角形。将不能构成闭合三角形。为了克服这个矛盾，在两点间另选一条单一的为了克服这个矛盾，在两点间另选一条单一的大地线代替相对法截线，从而得到由大地线构大地线代替相对法截线，从而得到由大地线构成的单一的三角形。成的单一的三角形。大地线的定义和性质大地线的定义和性质椭球面上两点间的最短椭球面上两点间的最短程曲线叫做程曲线叫做大地线大地线。大地线的性质大地线的性质:大地线是两点间惟一最短线，而且位于相对法截线之间，大地线是两点间惟一最短线，而且位于相对法截线之间，并靠近正法截线，它与正法截线间的夹角并靠近正法截线，它
21、与正法截线间的夹角 在椭球面上进行测量计算时，应当以两点间的大地线为在椭球面上进行测量计算时，应当以两点间的大地线为依据。在地面上测得的方向、距离等，应当归算成相应依据。在地面上测得的方向、距离等，应当归算成相应大地线的方向、距离。大地线的方向、距离。长度差异可忽略长度差异可忽略,方向差异需改化。方向差异需改化。31补充内容补充内容：球面直角三角形的球面三角公式球面直角三角形的球面三角公式 任一元素的余弦等于不相邻两元素的正弦之积任一元素的余弦等于不相邻两元素的正弦之积 coscos)90sin()90sin(cossinsinsinsinsinsinsincossin)90sin(cosco
22、ssin)90sin(sincos0000baabccBbcAaAbAbBaBaBA任一元素的余弦等于相邻两元素的余切之积任一元素的余弦等于相邻两元素的余切之积ctgAabBbaBAccgacaBcbcbAtgsinctgtgsinctgctgcos ctgtctg)90(ctgcosctgtgctg)90(ctgcos00大地线的微分方程和克莱劳方程大地线的微分方程和克莱劳方程AdSMdBcosdSMAdBcosAdSBdLNsincosdSBNAdLcossincos(90)sinsin(90(90)dAdLBdB)sin(sinsindBBdLdABdLdAsinBdSNAdAtansi
23、n大地线的克莱劳方程大地线的克莱劳方程 在旋转椭球面上，大地线各点的平行圈半径与大地线在该在旋转椭球面上，大地线各点的平行圈半径与大地线在该点的大地方位角的正弦的乘积等于常数。式中常数点的大地方位角的正弦的乘积等于常数。式中常数C C也叫也叫大地线常数大地线常数 BdSNAdAtansindSMAdBcosBNBdBMAAdAcossincossincossinrNB MBdBdrCrAlnlnsinlnCAr sin当大地线穿越赤道时当大地线穿越赤道时当大地线达极小平行圈时当大地线达极小平行圈时由克莱劳方程可以写出由克莱劳方程可以写出 0sin AaC 0090sinrrC2112sinsi
24、nAArrCAuasincosCABNsincos4.6 将地面观测值归算至椭球面将地面观测值归算至椭球面 观测的基准线不是各点相应的椭球面的法线，而是各观测的基准线不是各点相应的椭球面的法线，而是各点的垂线，各点的垂线与法线存在着垂线偏差。点的垂线，各点的垂线与法线存在着垂线偏差。归算的两条基本要求：归算的两条基本要求：以椭球面的法线为基准；以椭球面的法线为基准；将地面观测元素化为椭球面上大地线的相应元素将地面观测元素化为椭球面上大地线的相应元素。将地面观测的水平方向归算至椭球面将地面观测的水平方向归算至椭球面 将水平方向归算至椭球面上，包括垂线偏差改正、标将水平方向归算至椭球面上，包括垂线
25、偏差改正、标高差改正及截面差改正，习惯上称此三项改正为高差改正及截面差改正，习惯上称此三项改正为三差改正三差改正。垂线偏差改正垂线偏差改正 以测站A为中心作出单位半径的辅助球,u是垂线偏差，它在子午圈和卯酉圈上的分量分别以,表示，M是地面观测目标m在球面上的投影。垂线偏差对水平方向的影响是(R-R1)11tan)cossin(cot)cossin(mmmmuAAZAA标高差改正222212cossin22heHBAMaHH常2截面差改正截面差改正2222111()cossin 212geSBAN 将地面观测的长度归算至椭球面将地面观测的长度归算至椭球面 基线尺量距的归算基线尺量距的归算 将基线
26、尺量取的长度加上测段倾斜改正后，将基线尺量取的长度加上测段倾斜改正后，可以认为它是基线平均高程面上的长度，以可以认为它是基线平均高程面上的长度，以表示，现要把它归算至参考椭球面上的大地表示，现要把它归算至参考椭球面上的大地线长度线长度S。1.1.垂线偏差对长度归算的影响垂线偏差对长度归算的影响 )(22122121HHuuhuuSu2.高程对长度归算的影响高程对长度归算的影响 RHRHRSSmm10101RHSSm2201RHRHSSmm2200RHSRHSSmmH)(21122110HHuuRHSSm电磁波测距的归算电磁波测距的归算12212coscos1ABeNRA)(2)()(cos21
27、22221HRHRDHRHRAAAAAARSRS2sin21coscos2)(4)(2sin1221222HRHRHHDRSAAA)1)(1()(1arcsin221212AAAARHRHDHHRDRS232121224)1)(1()(1AAARDRHRHDHHDS232242AAmRDRHDDhDS)1)(1()(121212AARHRHDHHDd2322241AAmRDRHhDS作业：作业：1为什么要引入正常地球，讨论正常位与正常重力？为什么要引入正常地球，讨论正常位与正常重力？目前一般实用的正常地球是什么形状？目前一般实用的正常地球是什么形状？2什么是水准测量的理论闭合差？水准测量的观测
28、什么是水准测量的理论闭合差？水准测量的观测高差一般要加入那些改正？高差一般要加入那些改正？3正常重力公式是用来计算什么地方的正常重力正常重力公式是用来计算什么地方的正常重力?4.为什么水准测量会产生多值性？为什么水准测量会产生多值性？4.7 大地主题解算概述 大地元素：大地经度L，大地纬度B，大地线长度S，正反大地方位角A12、A21大地主题解算：已知某些大地元素推求另一些大地元素。大地方位角：大地方位角是大地坐标系中表示方向的角量，是参考椭球面上过某点的子午圈与过该点某一方向的大地线间的夹角，大地方位角由子午圈北方向起按顺时针方向计算,通常用A表示。它不能直接测得,而是由天文方位角按拉普拉斯
29、方程换算而得。4.7 大地测量主题解算概述大地测量主题解算概述4.7.1 大地主题解算的一般说明大地主题解算的一般说明 主题解算分为主题解算分为:短距离短距离(400km)(400km)中距离中距离(1000km)(1000km)长距离长距离(1000km(1000km以上以上)1112122221(,),(,),P B L S A P B L A12正算:已知 求1122121221(,),(,),P B L P B L S A A12反算:已知，求1.以大地线在大地坐标系中的微分方程为基础，直以大地线在大地坐标系中的微分方程为基础，直接在地球椭球面上进行积分运算。接在地球椭球面上进行积分运
30、算。主要特点：解算精度与距离有关，距离越长，收敛主要特点：解算精度与距离有关，距离越长，收敛越慢，因此只适用于较短的距离越慢，因此只适用于较短的距离 典型解法：典型解法：高斯平均引数法高斯平均引数法c o ss inc o sta ns ind BAd SMd LAd SNBd ABAd SN212121212121co ssinco stansinPPPPPPABBd SMALLd SNBBAAA d SN2.以白塞尔大地投影为基础以白塞尔大地投影为基础1)1)按椭球面上的已知值计算球面相应值，即实现椭球面按椭球面上的已知值计算球面相应值，即实现椭球面 向球面的过渡；向球面的过渡；2)2)在
31、球面上解算大地问题；在球面上解算大地问题；3)3)按球面上得到的数值计算椭球面上的相应数值，即实现从圆按球面上得到的数值计算椭球面上的相应数值，即实现从圆球向椭球的过渡。球向椭球的过渡。典型解法：典型解法：白塞尔大地主题解算白塞尔大地主题解算 特点：计算公式可展为特点：计算公式可展为 的幂级数，的幂级数，解算精度与距离解算精度与距离长短无关，它既适用于短距离解算，也适用于长距离长短无关，它既适用于短距离解算，也适用于长距离解算。可适应解算。可适应20 000km20 000km或更长的距离，这对于国际联或更长的距离，这对于国际联测，精密导航，远程导弹发射等都具有重要意义。测，精密导航，远程导弹
32、发射等都具有重要意义。,SALB2e2221BB S LL S AA S(),(),()1112000BB LL AA(),(),()可在已知点可在已知点P1（S=0)，按照麦克，按照麦克劳林级数将劳林级数将P1和和P2点的纬度差、经点的纬度差、经度差和方位角之差展开为大地线长度差和方位角之差展开为大地线长度度S的幂级数。的幂级数。4.7.2 勒让德级数式勒让德级数式 为了计算为了计算 的级数展开式，关键问题是推求各阶导的级数展开式，关键问题是推求各阶导数。数。2221BB S LL S AA S(),(),()1112000BB LL AA(),(),()22332111112323nnnd
33、 BSdBd BSd BSBBBSdSndSdSdS()()()()!22332111112323nnnd LSdLd LSd LSLLLS dSndSdSdS()()()()!22332111112318023nnnd A SdAd A Sd A SAAASdSndSdSdS()()()()!B L A,3coscossinsecsincostansintansindBAVAdSMcdLAVBBdSNBcdABVABAdSNc 242222234 190()()(cossin)()d BdB dBdB dAVtAA dSB dS dSA dS dSc 22222()()secsincos(4
34、192)d LdL dBdL dAVtBAA dSB dS dSA dS dSc 2222221 24 194()()sincos()()d AdA dBdA dAVAAt dSB dS dSA dS dSc 一阶导数：一阶导数：3coscossinsecsincostansintansindBAVAdSMcdLAVBBdSNBcdABVABAdSNc 二阶导数：二阶导数：242222234 190()()(cossin)()d BdB dBdB dAVtAA dSB dS dSA dS dSc 22222()()sec sincos(4 192)d LdL dBdL dAVtBAA dSB
35、dS dSA dS dSc 2222221 24 194()()sin cos()()d AdA dBdA dAVAAt dSB dS dSA dS dSc 三阶导数三阶导数352222 22222331 39312 2d BVAAttAt5 t dSccos sin()cos()32322d LVtBAA dScsec sin cos 332222233213d LVBAA+ttA dScsecsincos()sin 1cosuSA 1sinvSA 4.7.3 高斯平均引数正算公式高斯平均引数正算公式 高斯平均引数正算公式推导的基本思想：高斯平均引数正算公式推导的基本思想：首先把勒让德级数在
36、首先把勒让德级数在 P P点展开改在大地线长度中点点展开改在大地线长度中点M M展开，以使级数公式项数减少，收敛快，精度高；其次，展开，以使级数公式项数减少，收敛快，精度高；其次，考虑到求定中点考虑到求定中点 M M 的复杂性，将的复杂性，将 M M 点用大地线两端点平均点用大地线两端点平均纬度及平均方位角相对应的纬度及平均方位角相对应的 m m 点来代替，并借助迭代计算点来代替，并借助迭代计算便可顺利地实现大地主题正解。便可顺利地实现大地主题正解。21,22SSMP MP(1)建立级数展开式建立级数展开式:21,22SSMP MP 223322311()()()(4200)22468MMdB
37、Sd B Sd B SBB dSdSdS223312311()()()(4201)22468MMMMdBSd BSd BSBB dSdSdS 33213()()(4202)24MMdBd BBBBSS dSdS同理可得同理可得:3321324MMdLd LLLLSSdSdS()()332112324MMdAd AAAASSdSdS()()MMmm BABA,(2)2121121118022mmBBB AAA(),()mMmMBB AA,MMMmMmmMmdBf BAF BBB AAAdS()(,)(,)+MmmMmMmmmdBfff BABBAAdSBA()(,)()()()()+MmmMmM
38、mmmdBdBdBdSdSf BABBAAdSBA()()()(,)()()()()+22222288MmMmSd BBBdSSd B dS()()22222288MmMmSd AAAdSSd A dS()()(3)由大地线微分方程依次求偏导数由大地线微分方程依次求偏导数:32mmmmmmmmAVVdBAAdSMcNcos()coscos323mmmmmmmVdBAdSctABBN()(cos)()cos 32mmmmmmVdBAVdScAAAN()(cos)()sin 222222222388mMmmmmmmmmS VSd BBBtAtAdSN()(sincos)22222221288Mmm
39、mmmmmSd ASAAAAt dSN()sincos()222222223223333812mmMmm mmmmmmm22mmmmmmVVdBSSAA tAAS+dSNNV AAt+S+5 8N()coscos(sincos)sincos()次23322222332222313924243155mMmmmmmmm22mmmmmmVSd BAAtt+dSN At+tS+()cossin()cos()次22222212222221232243195mmmmmmmmmmm mVSBBBSAAtNN At()cossin()cos()次同理可得：同理可得：22222222124195mmmmmmmm
40、m mSLSBAAtNN Atsecsinsincos()次22222242221279245225m mmmmmmmmmmmSASA tAtNN Atsincos()sin()次21212112,180BBB LLL AAA21211111()222mBBBBBBBB1212m AAA注意：从公式可知，欲求，及，必先有及。但由于2和21未知，故精确值尚不知，为此须用逐次趋近的迭代方法进行公式的计算。除此之外，此方法适合与200公里以下的大地问题解算，其计算经纬计算精度可达到0.0001”,方位角计算精度可达到0.001”。4.7.4 高斯平均引数反算公式高斯平均引数反算公式 高斯平均引数反算
41、公式可以依正算公式导出:上述两式的主式为:222222222sinsincossin24cos(19)mmmmmmmmmm mSALSANBS tAN SAt2222222222222coscossin(232)243cos(14)mmmmmmmmmmmmm mNSABSASAtVN SAtt 2sincos,cosmmmmmmNLBSANB SAV2301210323101230mmSArLrBLrLSAsBsBLsBsincos 230 12 10 3AtLtBLtL 3222201210333192424mmmmmmmm mmNNBNBrB rt rtcoscoscos,(),22222
42、22101230233233248mmmmmmmmmmmNNBNs stt stVcos,(),()2432012103221132212412mmm mmmm mmttB tB t tB tcos,cos(),cos()sintancosmmmSAASAsinsinmmSASA122111,18022mmAAA AAA已知：求得：147 46 52.6470B 135 49 36.3300L 1244 12 13.6640A 44 797.2826S m248 04 09.6384B 23614 45.0004L 2122430.550A 53计算范例计算范例4.7.5 白塞尔大地主题解算方
43、法白塞尔大地主题解算方法白塞尔法解算大地主题的基本思想白塞尔法解算大地主题的基本思想:以辅助球面为基础以辅助球面为基础,将椭球面三角形转换为辅助球面的将椭球面三角形转换为辅助球面的相应三角形相应三角形,由三角形对应元素关系由三角形对应元素关系,将椭球面上的大地元素按将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上，然后在球面上进行大地照白塞尔投影条件投影到辅助球面上，然后在球面上进行大地主题解算，最后再将球面上的计算结果换算到椭球面上。主题解算，最后再将球面上的计算结果换算到椭球面上。这种方法的关键问题这种方法的关键问题1 1、找出椭球面上的大地元素与球面上相应元素之间的关系式、找出椭球
44、面上的大地元素与球面上相应元素之间的关系式2 2、解决在球面上进行大地主题解算的方法。、解决在球面上进行大地主题解算的方法。12121212,BBAAL S 在球面上进行大地主题解算在球面上进行大地主题解算 球面上大地主题正算球面上大地主题正算:已知已知 求解求解 球面上大地主题反算球面上大地主题反算:已知已知 求解求解11 ,22 ,12 ,12 ,球面三角元素间的球面三角元素间的相互关系相互关系12211121221212sinsinsincos()sinsinsincos()sincoscossinsincoscos()sincossincoscossinc a b c 12122111
45、221112211os()cossinsincoscoscos()coscoscoscossinsincos()coscossinsincoscoscos()cossincossin d e f g 2111()sinsincoscossincos()h i球面上大地主题正解112,2已知 求,2111sinsincoscossincos()i1111sinsintan()()coscossinsincos af112111cossintan()()coscoscossinsin hg球面上大地主题反解方法1212,已 知，求,211212sin costan()()cos sinsin co
46、scosp acq121212sin costan()()coscoscossincosu bduuuu111212sincostansinsincoscoscospq2 椭球面和球面上坐标关系式dddddAdcossincostansin AdBdSMAdLdSNBBdAAdSNcossincostansin 在椭球面上与单位球面上的大地线微分方程为在椭球面上与单位球面上的大地线微分方程为:AdBdSMAdLdSNBBdAAdSNcossincostansin dddddAdcossincostansin 423842394240dBBdS dMddLAdS dNBddABA dS dNdc
47、os()coscossin()cossintansin()tansin 白塞尔提出如下三个投影条件：1.椭球面大地线投影到球面上为大圆弧2.大地线和大圆弧上相应点的方位角相等；3.球面上任意一点纬度等于椭球面上相应点的归化纬度。121212,Aa Aa2111111tan1tantan(1)tanueBuBu2222122tan1tantan(1)tanueBuBu111222121212(),(),uB uB A A,?S L ,?S L 423842394240dBBdS dMddLAdS dNBddABA dS dNdcos()coscossin()cossintansin()tansi
48、n 2211dSNuANucaNeedBABVVtansintantansintan 2222222222222111111eVeBWue VueueVcoscoscoscos 白塞尔提出的第三个投影条件：3.球面上任意一点纬度等于椭球面上相应点的归化纬度。2111111tan1tantan(1)tanueBuBu2222122tan1tantan(1)tanueBuBuBNddStantanBuNtantan21 eWaBddLsinsinV1BeV222cos1uWee2222cos11uVe222cos1ueV222cos11Va）公式（）公式（253425242122211PPLLLe
49、udcos 221cosdLeud2122221cos1cosppdSaeuSaeudd 以上为白塞尔微分方程以上为白塞尔微分方程.3 白塞尔微分方程的积分白塞尔微分方程的积分21221cosppSaeud2201Sbkd(sin)21221cosppSaeu d1019090uAcos()sin()sin 2221011uAcoscossin 22112222222001111ppppSaeAdaeeAd(cossin)cossin 2220keAcos 246221 2246112816kkkk(sin)sinsinsin 246112223 11248285153124616321632
50、xxxxxxxxxsincossincoscossincoscoscos 积分得到下式：积分得到下式：1111122222sin 2sin 2cos 2sin 2sin 2cos 2SABCSABC11222222SABCBCsin(cos)sin(cos)2462464635146 42 5 61 583 21 0 2 431 2 85 1 6kkkAbkkBbkkCbk()()()适合于反算:适合于正算:迭代法:直接法:1122sin 2(cos 2)sin 2(cos 2)SABCBC111112222SBCBCAsin(cos)sin()(cos()011122SBCAsin(cos)
