第二章地球坐标系和地球椭球ppt课件

1、第二章第二章 地球坐标系和地球椭球地球坐标系和地球椭球2.1 概 述大地丈量采用的坐标系：天球坐标系、地球坐标系地球坐标系：固定在地球上与地球一同自转和公转的 坐标系地球坐标系分类：参心坐标系、地心坐标系定义坐标系的要素：原点位置、尺度与坐标轴指向；还包括一些天文、物理、地球等参数，假设采用大地 坐标表述方式还需求椭球元素。2.2 地球椭球面的数学计算和有关计算2.2.1 地球椭球的几何、物理元素椭球方程：扁率：第一偏心率：第二偏心率：XYZO1222222bZaYaXabaaEabae222bEbbae2222.2.1 地球椭球的几何、物理元素续1几个关系式：几个关系式：2222222222

2、1 111 1 2eeeeeeeee21eba21 eab1954年北京坐标系，克拉索夫斯基椭球元素：年北京坐标系，克拉索夫斯基椭球元素：3 .2981 m 6378245a2.2.1 地球椭球的几何、物理元素续21980年大地坐标系采用第16届 IAGIUGG 椭球，其椭球元素为：257.2981 /10292115.7 10108263/103.986005GM m 63781405822314可求得扁率：sradJsma2.2.2 旋转椭球面的参数表示及数学性质1、经线和纬线的曲线方程在XOZ坐标面上的起始经线方程：OXYZM1M0MLLrARS0 12222YbZaXM0饶Z轴旋转，构

3、成纬圈平行圈，其半径：22YXr经度为L的经线方程：LXYbZaYaXtan 12222222.2.2 旋转椭球面的参数表示及数学性质续1OXYZM1M0MLLrARS纬圈方程：0222222Z 1ZbZaYaX2.2.2 旋转椭球面的参数表示及数学性质续22、椭球面法线与子午线主法线的同一性、经纬线的Frenet标架POQMPRTNA如图为过M点的子午面。子午线的主法线MP位于子午面内，且垂直于子午线切线T；R为过M点的平行圈切线，显然R垂直于M点的子午面，因此R垂直于MP。所以， MP垂直于椭球面在M点的切平面，因此它是椭球面的法线。Frenet标架：曲线上恣意一点处的三个相互正交的单位标

4、架：曲线上恣意一点处的三个相互正交的单位向量取切向、主法向和与该两个方向正交的第三个方向量取切向、主法向和与该两个方向正交的第三个方向构成的三维直角坐标系。向构成的三维直角坐标系。2.2.2 旋转椭球面的参数表示及数学性质续33、旋转椭球面及经纬线的参数方程1). 以大地经度L及归化纬度u为参数的方程uaXZOMMubZuaXsincos在XOZ子午面内，有在三维空间坐标系中：ubZLuaYLuaXsinsincoscoscos(2). 以大地经纬度L、B为参数的方程XZK0B90+ BOTM 0切线M0T的斜率的导数式：BBdXdZctg90tan0由椭圆方程求导得：12222bZaXZXe

5、ZaXbdXdZ2221代入第一式得：BeXZtan1212.2.2 旋转椭球面的参数表示及数学性质续42.2.2 旋转椭球面的参数表示及数学性质续5将将 代入椭圆方程，化简后得：代入椭圆方程，化简后得：1BeBeaZBeBaX22222sin1sin1 sin1cos引入辅助符号：WaNBeW sin122那么有： sin1 cos2BeNZBNX2.2.2 旋转椭球面的参数表示及数学性质续6以大地经纬度两个参数表示的椭球面上一点的三维坐标以大地经纬度两个参数表示的椭球面上一点的三维坐标为椭球面参数方程式：为椭球面参数方程式： sin1sincoscoscos2BeNZLBNYLBNX以大地

6、纬度为参数的经度为LC的子午线参数方程为:BeNZLBNYLBNXCCsin1sincoscoscos2在一点在一点BC ,LC 处的子午线切向量处的子午线切向量CcCCcCCcBMBZLBMBYLBMBXcossinsincossin子午线切线单位向量子午线切线单位向量CCCCCBLBLBcossinsincossin以大地经度为参数的大地纬度为以大地经度为参数的大地纬度为的的BC纬线的参数方程为纬线的参数方程为CBCBCBBeNZLBNYLBNXCCCsin1sincoscoscos2在一点在一点BC ,LC 处的平行圈切向量处的平行圈切向量0coscossincosLZLBNLYLBNL

7、XCCcCCc平行圈切线单位向量平行圈切线单位向量0coscossincosCCCCLBLB0coscossincosCCCCLBLBCCCCCCCCCCBLBLBBLBLBsinsincoscoscoscossinsincossin椭球面单位法向量为其矢量积：椭球面单位法向量为其矢量积：2.2.2 旋转椭球面的参数表示及数学性质续7(3). 以大地经度L及球心纬度为参数的方程XZOM 0 球心纬度，向径，那么对于XOZ平面上的椭圆有：sincosZX 在椭圆上，向径由球心纬度独一确定，将上式代入椭圆方程，得：222cos11eea2.2.2 旋转椭球面的参数表示及数学性质续8对于对于XOZ平

8、面上的椭圆有：平面上的椭圆有：222222cos11sin cos11coseeaZeeaX以经度、球心纬度两个参数表示的椭球面上一点的三维以经度、球心纬度两个参数表示的椭球面上一点的三维坐标为参数方程式为：坐标为参数方程式为：2222222221coscos 1cos1cossin1cos1sin1coseXaLeeYaLeeZae2.2.2 旋转椭球面的参数表示及数学性质续9 不难得出，不难得出，u, B, 的关系为：的关系为：BeeBuBeBu22222sin11sinsin sin1coscos因此有：BabBeutantan1tan2由球心纬度公式，得：BabBetantan1tan

9、2222.2.2 旋转椭球面的参数表示及数学性质续104、旋转椭球面的几何性质 a). 对称性 b). 有界性 c). 正那么性：曲面上每点都对应于独一 确定的非零法向量。 d). 不可展性BLBLBsinsincoscoscosn2.2.3 法截线曲率及曲率半径 1、空间曲线的曲率几曲率半径 假设以曲线的弧长s为参数，曲线上的点位用向量r(s)表示。那么曲线的曲率为： cccsssdssddssdsk22rT假设以t参数，那么曲线的曲率可表示为： 322ccctttdttddttddttdtkrrr2、椭球面法截线的曲率(1). 子午线曲率半径 不失普通性，以起始子午线为例推导。假设以归化纬

10、度u为子午线方程的参数，那么有：2.2.3 法截线曲率及曲率半径续1 ubuaduudubuaduudubuausin0coscos0sinsin0cos22rrr2.2.3 法截线曲率及曲率半径续2那么有：那么有： 0022abduudduudrr 232222322sincosuaubabduudduudduudukrrr同理，假设以大地纬度为参数，得： 232232211sin1eaWeaBeBk子午曲率半径M，就是曲率是倒数，即： 3211WeaBkM2.2.3 法截线曲率及曲率半径续3(2). 卯酉线曲率半径定义：与子午面切线正交的法截面与椭球面的交线为卯 酉线。 根据微分几何中的麦

11、尼尔定理，卯酉圈曲率kn与平行圈曲率kr的关系为：Bkkrncos平行圈半径为子午面XOZ 平面内的X坐标，即：BWaXrBcos那么有，上述两式得卯酉曲率半径N为：WaBrBkkNBrncoscos112.2.3 法截线曲率及曲率半径续 4(3). 恣意方向法截线的曲率半径 根据微分几何中的Euler公式，恣意方向法截线的曲率与子午、卯酉曲率半径的关系为：NAMAkA22sincos因此，恣意方向的曲率半径为：AMANMNkRAA22sincos1当A为0，/2， 3/2时，获得极值。2.2.3 法截线曲率及曲率半径续 5(4). 平均曲率半径 定义：一切方向法截线曲率半径的平均值。dAAN

12、MdtANMtdAAMANMNdARRA2202220cos1 ,tansincos22代入上式，得：0222112WeaMNtdtMNR2.2.3 法截线曲率及曲率半径续 6不难得到：不难得到：N R M引入辅助量：引入辅助量：222221cos11BeVeabac存在以下关系：BeBtVcRVcNVcMeVWeWVecaeac222232222cos tan 1 11 12.2.4 椭球面上第一根本方式及弧长面积计算1. 椭球面的第一根本方式dLLdBBdBeNLBNLBNrrrrrsin1sincoscoscos2椭球面上点的向量：椭球面上的微分弧长：22222GdLFdBdLEdBdd

13、ddSrrr其中：dLddLdGdLddBdFdBddBdErrrrrr 对于椭球面：BNGFME222cos 0 222222cos BdLNdBMdS2.2.4 椭球面上第一根本方式及弧长面积计算续12、子午线弧长子午线微分弧长：MdBdS 积分得：21212322221sin11BBBBdBBeeaMdBS11112211119229117227115225113223112212221cossincossin cossincossincossincossin cossincossincossincossin cossincossin1BBBBGBBBBFBBBBEBBBBDBBBBCB

14、BBBBarcBarcBAeaS用二项式展开，并逐项积分得：常数 A、B、C、D、E、F、G的计算公式见教材2.2.4 椭球面上第一根本方式及弧长面积计算续2对于小于对于小于400km的弧长，可采用以下简化式。的弧长，可采用以下简化式。3331221241BdBSdBdBdSSSSmm其中：1221 21BBBBBBm根据： 33dBdSdBddBddBSdMdBdS求出导数，代入上式并化简，得： 2cos812221BBeBMSmm对于小于40km的弧长，可进一步简化为： 21BMSm2.2.4 椭球面上第一根本方式及弧长面积计算续3知知B1和弧长和弧长S12求求B2称为反算，可采用叠代法计

15、算。称为反算，可采用叠代法计算。初值：初值：1221021arcBAeaSarcB叠代格式： kBkkkSBSSBB22122121212其中： 122281222612224222222221cos82cossin cos62cossin cos42cossin cos22cossin2cos12BBBFBBBEBBBDBBBCBBAeaSkkkkkkkkkkBk mBSSk001. 022121要求：2.2.4 椭球面上第一根本方式及弧长面积计算续43、平行圈的半径与弧长BeLLBaLLBNSBeBaBNrB2212122122sin1coscossin1coscos一样经差的平行圈弧长在

16、赤道最长，越接近两极越小。2.2.4 椭球面上第一根本方式及弧长面积计算续54、利用经纬格网计算椭球面的面积LL+dLBB+dBMdBNcosBdLd2222sin1cos1cosBeBdBdLeaBdBdLMNd21212222sin1cos1BBLLDBeBdBdLeadA2.2.4 椭球面上第一根本方式及弧长面积计算续6上式利用二项式展开并积分，得：上式利用二项式展开并积分，得：17276152541323212122sinsin74sinsin53 sinsin32sinsin1BBeBBeBBeBBLLeaA 取 L2-L1 = 2，B2 = /2，B1 = 0 算得半球面积，乘2可

17、以估算全球面积约为5.1亿平方公里习 题1、导出三种纬度、导出三种纬度 、u与与B的关系。的关系。2、导出子午曲率半径、导出子午曲率半径M与卯酉曲率半径与卯酉曲率半径N的计算公式。的计算公式。3、M、N、R的关系如何？在什么条件下三者一样？的关系如何？在什么条件下三者一样？4、某点到赤道的子午弧长、某点到赤道的子午弧长 ,求该点的求该点的纬度。纬度。 a=6378245, =1/298.35、知某点的纬度、知某点的纬度 ,求该点自赤道起的子求该点自赤道起的子午弧长。午弧长。 a =6378245, =1/298.32831.1628310B米193.3745682S2.2.5 大地线1、大地线

18、的定义与性质法截弧：由椭球面上A点的法线与B点所确定的法截面与椭球面相割得到的曲线称为A到B的法截弧。相对法截弧： A到B的法截弧与B到A的法截弧。 由相对法截弧构成的椭球面三角形不是闭合图形。2.2.5 大地线续1大地线的定义：大地线的主法线与曲面法线处处重合。大地线的性质：1、大地线上任何点的亲密平面就是该点 的法截面； 2、曲面上衔接任何两点的最短直线必为 大地线。 3、大地线的测地曲率等于0曲线的测地曲率：曲线的曲率在曲面切平面上的投影。大地线的曲率：大地线的挠率ANNAMAkg2222cos11sincosAANAAMNgcossincossin1122.2.5 大地线续22、大地坐

19、标系中大地线的微分方程(1). 大地线的二阶微分方程以u,v 为参数的普通曲面的大地线微分方程可表示为： 下标为相应的偏导数。 22 22 22 222222222322FEGEFEEFESFEGFFFEGEFEGFEEGRFEGFFFGEGFEGFGGEQFEGGFGGFGPSdudvRdudvQdudvPduvduvuuvuuvvuvuvvuv2.2.5 大地线续3对于椭球面，有：对于椭球面，有：代入前面公式，得：代入前面公式，得：那么旋转椭球面上大地线的微分方程为：那么旋转椭球面上大地线的微分方程为：BNGFME222cos 0 0 2 0 sincos222StVRQBBVP32222

20、2cossin2dBdLBBVdBdLVtdBLd2.2.5 大地线续4(2). 克莱劳定理克莱劳定理直角坐标系中的椭球面方程：01222222bZaYaXF椭球面法向量为：222222bZaYaXZFYFXFN以大地线弧长为参数的大地线主法线向量为：222222dSZddSYddSXdn两者指向一致，即：222222222222dSZdbZdSYdaYdSXdaX2.2.5 大地线续5由上式的前两个方程得：由上式的前两个方程得：CdSdYXdSdXYdSYdXdSXdY 02222积分得：将三维空间坐标与大地坐标的关系式及微分式代入：dSdLLrLdSdrdSdYdSdLLrLdSdrdSd

21、XLrLBNYLrLBNXBBBBBBcossin sincossinsincos coscoscos1代入 式，整理得：1CdSdLrB222.2.5 大地线续6将关系式：将关系式：即：大地线上各点的平行圈即：大地线上各点的平行圈半径与该点的大地线方位角半径与该点的大地线方位角正弦的乘积是常数。正弦的乘积是常数。AdSdLrBsin代入上式，即得克莱劳定理：CArBsinBrdLrBAdSdAAMdBdL2.2.5 大地线续7(3). 大地线的一阶微分关系式BrdLrBAdSdAAMdBdLBNAdSdLAdSBdLNMAdSdBAdSMdBcossin sincoscos cos由克莱劳定

22、理，微分得：BBBBBBdrrAdrArAdAAdArAdrtancossin:0cossin则2.2.5 大地线续8又如下图：又如下图：MdBdrBBMdBdrBsin代入上式，得：AdSNBdSMAABNBAMABNBdBAMdAsintancoscoscossinsin coscossinsin三个微分关系式可整理为：ANBdSdABNAdSdLMAdSdBsintan cossin cos 32.2.5 大地线续93、以弧长和大地方位角为参数的大地线方程 大地线始点坐标P0(B0,L0)，大地线上任何点的位置向量都可以展开成S,A的级数方式： 5554443332221201241 6

23、121sdsdsdsdsdsdsdsdsdsdsrrrrrrFrenet标架的坐标轴定义：x指向大地线的切向t， y指向大地线的主法向n，向内为正， z指向大地线的副法向b，构成左手系。000,LBPASLBP,xyz42.2.5 大地线续10显然有：显然有：dsdrt 根据曲线论中的Frenet公式：nbbtnntggggdSdkdSdkdSd 由以上两式可求出各阶导数：23223243242545 32gggggggggggggggggdkddddkkkkdSdSdSdSdSdkd kdkddkkkkdSdSdSdSdSdkdS rtrnntnbrtnbrt2.2.5 大地线续11 将上式

24、代入大地线展开式将上式代入大地线展开式 ，得，得Frenet标架下的三维坐标：标架下的三维坐标：443423223254432224161241612112018161SdSdkdSdkSkzSkkdSkdSdSdkSkySkSdSdkkSkSxgggggggggggggggg5ANNAMAkg2222cos11sincosAANAAMNgcossincossin112顾及公式：2.2.5 大地线续12和：和：ANBdSdAMAdSdBsintan cos 求导得：ANAtdSdAAdSdBBdSdAtNdSdAdSdkAdSdBdSdkBdSkdANAtdSdAAkdSdBBkdSdkggg

25、gggggg22222232222222cos31sincos3cos1cos3 2.2.5 大地线续13代入代入Frenet标架下的三维坐标公式标架下的三维坐标公式 ，得：，得：542332242223322222544323222sin2412sin121cos231241 cos21cos12112018cos3cos161AStNASNzStANAStNSANySNSNAtSANSx2.2.5 大地线续14 将坐标系饶将坐标系饶 y 逆时针旋转逆时针旋转A，得得x、y、z坐标系，那么有：坐标系，那么有：SA000, LBPLBP,xyzzx zyxAAAAzyxcos0sin010si

26、n0cos 以P0点为原点的地平坐标系(站心坐标系) x、y、z，与x、y、z坐标系的关系为：A000, LBPLBP,xyzzxyyzzyxx 2.2.5 大地线续15 最后得到地平坐标系站心系中的大地线方程，最后得到地平坐标系站心系中的大地线方程，称为称为Weingarten级数式。级数式。4222332222254432322254423232222cos231241 2coscos121120sin3cossin cos16sinsin120coscos8124 cos16coscosStANSNAtSANzSNASNAAtSANAASySNASANtSANAASx62.2.5 大地线

27、续16法截弧为平面曲线，其挠率为法截弧为平面曲线，其挠率为0，同理可推得地平坐标系中的计，同理可推得地平坐标系中的计算式为：算式为：42223322222544323222544323222cos231241 2coscos121120sin8cossin3 cos16sinsin120cos8cos3 cos16coscosStANSNAtSANzSNASNAAtSANAASySNASNAtSANAASx2.2.5 大地线续174、基于大地线的椭球面曲线坐标系(1). 大地线极坐标系大地圆：到极点具有一样大地线长 度的点所构成的轨迹。 由大地线长度和大地方位角可描画曲面点的位置 。As,rA

28、dAsdsmdAdS0P如下图：2222dAmdsdS对照第一根本方式，得：2 0 1mGFE由图中的微分直角三角形，得大地极坐标系中的微分关系式：sin1 sin1 cosdSdmmdSdmdSdAdSds2.2.5 大地线续18大地线的归化长度大地线的归化长度 m 的计算公式：的计算公式：44233222120cos3161NSAtNSNSSm222AzAyAxGm由 式求出偏导数代入得：62.2.7以长度量为坐标参数的新大地坐标系(1)(1)以长度量表示的椭球面上坐标系的由以长度量表示的椭球面上坐标系的由来来早在早在18101810年年SoldnerSoldner就提出了球面直角坐就提出

29、了球面直角坐标系统标系统 。以后。以后HelmertHelmert，GrossmannGrossmann，HeckHeck等德国丈量学者基于等德国丈量学者基于SoldnerSoldner球面球面直角坐标系推行提出了椭球面直角坐直角坐标系推行提出了椭球面直角坐标系。教材上所述的测地坐标系与其标系。教材上所述的测地坐标系与其有类似之处。有类似之处。(2)一种新型的大地坐标系一种新型的大地坐标系2019年提出了椭球面上一种新型的大地坐标年提出了椭球面上一种新型的大地坐标系：它仍以经纬线作为坐标曲线，且与大系：它仍以经纬线作为坐标曲线，且与大地坐标系之间能进展准确的坐标转换；它地坐标系之间能进展准确的

30、坐标转换；它所采用的坐标参数是以长度而不是以角度所采用的坐标参数是以长度而不是以角度为单位；可简化椭球面上的繁复计算。为单位；可简化椭球面上的繁复计算。施一民，朱紫阳，范业明施一民，朱紫阳，范业明.坐标参数为长度量坐标参数为长度量的一种新型的大地坐标系的一种新型的大地坐标系. 同济大学学报，同济大学学报，2019,3311: 1537-1540 新型大地坐标系的定义为构建新型的大地坐标系为构建新型的大地坐标系, ,可在区域中心附近可在区域中心附近选择一点作为其坐标原点选择一点作为其坐标原点, ,其在大地坐标系其在大地坐标系中的大地经纬度设为中的大地经纬度设为B0B0，L0L0。两族互。两族互为

31、正交的经纬线构成坐标系的坐标格网。为正交的经纬线构成坐标系的坐标格网。设过经纬度为设过经纬度为B B，L L的任一点的经线与的任一点的经线与起始纬线的交点为起始纬线的交点为PB0,PB0,点点P P至点至点PB0PB0的经线的经线上弧长作为纵坐标上弧长作为纵坐标sL,sL,点点P0P0至点至点PB0PB0的纬线的纬线上弧长取为横坐标上弧长取为横坐标sB,sB,如下图。如下图。新型大地坐标系的定义续新型大地坐标系的定义续过点P的纬线上的微分弧长ds与起始纬线上相应的微分弧长ds之比n称为按纬度变化方向的长度归化因子 n= ds/ds在大地坐标系中，与这两个平行圈上微分弧长相应的经差dl 相等，故

32、可准确求得 n= Nsin B/ N0sin B0式中N0，N分别为纬度B0，B处的卯酉线曲率半径。基于微分几何和椭球大地丈量的实际，n可用新大地坐标表示为(取至二次项) 202002/tan1NsNsBnLL4.新型大地坐标与大地坐标的转换关系BBLBMdBsLLBNs0000cos习 题1. 纬度一样的两个点的相对法截弧能否重合？此线能否纬度一样的两个点的相对法截弧能否重合？此线能否就是大地线？就是大地线？2. 推导大地线的三个微分式。推导大地线的三个微分式。3. 试述测地坐标系的定义？测地平行线能否等距？测地试述测地坐标系的定义？测地平行线能否等距？测地大地线能否等距？大地线能否等距？4

33、. 简述简述weingarten级数的推导步骤。级数的推导步骤。2.3 椭球面上大地坐标的计算2.3.1 程度方向、边长观测值归算到椭球面1、程度方向观测值归算到参考椭球面的矫正 包括三项矫正，称为三差矫正。(1). 垂线偏向矫正(2). 标高差矫正12121121cossinctgzAAu2221222cos2sin2HBAMeh)(2sincos1089. 021222KmHABh用椭球半径的近似值代入得：2.3.1 程度方向、边长观测值归算到椭球面(3). 法截弧方向归算到大地线方向的矫正121222121212212122sincos12cossin6ABeNSAANSg 该项矫正很小

34、，100公里约0.03“，只需一等控制网才估计此项矫正。2.3.1 程度方向、边长观测值归算到椭球面2、空间边长归算至参考椭球面的矫正 测线端点的大地高为：Dd1HS2H12AR21AR1P2PvhHihH222111椭球面上弦长 d 的计算公式211221212211AARHRHHHDd省略H/R的二次项，得： 211221 211212122AAAmAmRRRHHHRHHHDd2.3.1 程度方向、边长观测值归算到椭球面椭球面上的弧长为：椭球面上的弧长为：Dd1HS2H12AR21AR1P2P2232122212223331241 2448222sin2AAmAAAAAARHHDRHHHD

35、RddRdRdRRdRS2.3.1 程度方向、边长观测值归算到椭球面3. 工程控制网中的地面观测元素的归算 以平均高程面作投影面，范围小，可以用球替代椭球；球半径采用高斯平均曲率半径。计算公式为：22321222122241 RHHDRHHHDSm不难证明：椭球半径的误差对边长归算结果影响很小，R取6371km即可，但高差误差对边长归算比较敏感。2.3.2 椭球面上三角形解算1、球面角超AABBCC222222224224224RRFRRFRRF三块面积之和为：FRFFF222代入球面角超定义式，得：2RF2.3.2 椭球面上三角形解算按球面三角公式：按球面三角公式：222222222241s

36、in213241sin21RcbaabRRcbaabF当边长小于40公里时，第二项影响小于0.0004“，可略去sin212abR2.3.2 椭球面上三角形解算2、解算球面三角形的勒让德定理勒让德定理：对于较小的球面三角形，可用平面三角公式来解算，只需使三个平面角等于相应的球面角减去三分之一的球面角超，而边长坚持不变。ABCabccCbBaA3sin3sin3sin2.3.3 大地主题解算大地主题解算分类：正算：知(B1, L1)，A12，S12，计算(B2, L2)，A21反算：知(B1, L1)， (B2, L2)， 计算A12，S12 ，A21短间隔中间隔长间隔中短间隔解算方法：按级数展

37、开采用Gauss平均引数公式； 长间隔解算方法：贝塞尔公式KmS120KmSKm400120KmS4002.3.3 大地主题解算1、纬度差、经度差和方位角差展开为大地线长度的级数式626262303320220123033202201230332022012SdSAdSdSAdSdSdAAAaSdSLdSdSLdSdSdLLLlSdSBdSdSBdSdSdBBBb由大地线的微分公式，得其一阶导数为：ANBdSdABNAdSdLMAdSdBsintan cossin cos 2.3.3 大地主题解算二阶和三阶导数采用复合函数求导法计算：二阶和三阶导数采用复合函数求导法计算：dSdAdSBdAdS

38、dBdSBdBdSBddSdAdSdBAdSdBdSdBBdSBd22223322 同理可求出四阶以上的导数和L、A的高阶导数，代入展开式即可。2.3.3 大地主题解算2、高斯平均引数公式假设取大地线中点展开，得：4882488233322213332222SdSBdSdSBdSdSdBBBSdSBdSdSBdSdSdBBBMMMMMMMM两式相减，得：111,LBP222,LBPMN12AMA2S2S2433312SdSBdSdSdBbBBMM类似地，有：2433312SdSLdSdSdLlLLMM243331221SdSAdSdSdAaAAMM12.3.3 大地主题解算两式相加，得：两式相

39、加，得：8222SdSBdBBMMm类似地，有：8 ,8222222SdSAdAASdSLdLLMMmMMm其中：21212121 ,21 ,21AAALL LBBBmmm将 展开成级数，得：MdSdBMmmMmmmMAAdSdBABBdSdBBdSdBdSdB22.3.3 大地主题解算由大地线的微分公式：由大地线的微分公式：mmmANVMAdSdBcoscos 2求导，得：ANVdSdBAAtMdSdBBmmmmmmmsin cos322mMmMdSAddSAddSBddSBd22222222 可取：代入 式，得 的计算公式。并取2MdSdBmMdSBddSBd3333代入 式，求出各阶导数

40、后整理得：12.3.3 大地主题解算2222222222221241cos3 232sin241cosmmmmmmmmmmmmmttAtANSASNVBBb同理可得：222222221291cos sin241sincos1mmmmmmmmmmtAtANSASBNLLl4222222222125972cos 22sin241sinmmmmmmmmmmmmtAtANSASNtAAa以上3式具有4次方精度，可用于解算200公里下的大地主题。32.3.3 大地主题解算因计算因计算Bm , Lm要用到要用到B2 , L2，因此需求叠代计算。其初值为：，因此需求叠代计算。其初值为： 1211120121

41、10sintan21cos21ASBNAAASMBBmm叠代计算公式为： 2 2 12111kkmkkmaAAbBB直到 为止。 10.00 10.000 11 kmkmAkmkmBAABB最后计算纬度、经度和方位角：aAAlLLbBB12211212 2.3.3 大地主题解算3、高斯平均引数反算公式由正算公式，反解得：222222222222241cos3 232sin24coscosmmmmmmmmmmmmmmttAStASNASNbVAS22222222291cossin24sincossinmmmmmmmmmmmtAStASNASBlNAS右端第二项与第一项相比为小量，可以作近似：mm

42、mmmmBlNASNbVAScossin cos22.3.3 大地主题解算代入上式第二项，得：代入上式第二项，得：3222222222824332coscosbtNblttBNVbNASmmmmmmmmmmmm32222224sincos9124coscossinlBBNlbtBNlBNASmmmmmmmmmmm由此可求得平均方位角和大地线长度如下：mmmmmAAS SASASAsinsin cossintan12.3.3 大地主题解算由正算公式的第三式，计算由正算公式的第三式，计算a：4222222222125972cos 22sin241sinmmmmmmmmmmmmtAtANSASNtA

43、Aa最后得起终点的大地方位角为：2 22112aAAaAAmm2.3.3 新大地主题解算4、新大地坐标系与大地坐标系间的坐标转换(1). 由 B,L求解sL ,sB (2). 由sL ,sB求解B,L:B由下式迭代反解 而L由下式直接反解.BBLBMdBsLLBNs0000cosBBLMdBs0000cosLBNsLB000cosLBNsLB000cosLBNsLB000cosLBNsLB2.3.4 大地主题微分公式1、大地主题正解微分公式 终点的经纬度B2,L2和大地线方位角A21，与起点的经纬度B1,L1和大地线方位角A12，以及大地线长度S的微分关系。121112211222121222

44、1221212212112212122sinsincostansin0cossincoscoscossin1tansinsincos0cosdAdSdLdBAALLRBABLLBRASBRABLLRASRAAAdAdLdB2.3.4 大地主题微分公式2、大地主题反解微分公式 起点大地线方位角A12和大地线方位角A21，以及大地线长度S与起点和终点的经纬度B1,L1和B2,L2的微分关系。21211211121121212121222122212121212221222121212112coscoscoscossinsincoscoscoscossinsinsincossincoscoscosd

45、LdLdBdBASBNASBNASMASMASBNASBNASMASMABNABNAMAMdAdAdS习 题1、地面观测方向归算到椭球面上需求加哪几项矫正？、地面观测方向归算到椭球面上需求加哪几项矫正？2、地面观测间隔归算到椭球面上二步矫正的几何意义？、地面观测间隔归算到椭球面上二步矫正的几何意义？3、 P1与与P2与为控制点，知：与为控制点，知：计算归算到椭球面上的长度计算归算到椭球面上的长度4、知、知利用利用Gauss平均引数公式正反算。平均引数公式正反算。762.294175A ,814.249090953.2746115 ,3421.080040012120101mSLB,67.395

46、0 ,46.4130 , 6 .0338 , 6 .542366 .8356 , 3 .2238 ,456.2867821020210120111mHmHBAABmdPP2.4 空间大地直角坐标系及其转换模型2.4.1 空间直角坐标系与相应大地坐标系的关系1、X、Y、Z与B、L、H间的关系 空间坐标系的定义：Z自转轴，X位于赤道面，指格林尼治天文台，Y指东，构成右手系。大地坐标的定义：B为过一点的椭球面的法线与赤道面交角、L为过同一点的子午面与起始子午面二面角的平面角，H为点沿法线到椭球面的间隔。大地高与正高、正常高之间的关系：HNHHNXYZLBOPKPPQ2.4.1 空间直角坐标系与相应大

47、地坐标系的关系如下图：如下图：XYZLBOPKPPQBeNLBNLBNZYXPPOsin1sincoscoscos2rBHLBHLBHHPPsinsincoscoscosnrBHeNLBHNLBHNZYXPPPOOPsin1sincoscoscos2rrr12.4.1 空间直角坐标系与相应大地坐标系的关系2、由X、Y、Z计算B、L、H的迭代解法计算L：22arcsinarctanYXYXYL迭代计算B： 2221sinarctanYXBeNZBiii迭代初值为： 220arctanYXZB最后计算H：NBYXeNBZHsec1csc2222.4.1 空间直角坐标系与相应大地坐标系的关系3、X、

48、Y、Z与B、L、H间的微分关系由前面 式微分得；1dHdLdBdHdLdBBBHMLBLBHNLBHMLBLBHNLBHMdZdYdXAJ sin0cossincoscoscossinsincoscossincoscossinBBLBLLBLBLLBsin0cossincoscossinsincoscossincossinA1000cos000BHNHMJ其中：对角阵2.4.1 空间直角坐标系与相应大地坐标系的关系顾及顾及A是正交阵，是正交阵，J是对角阵，得：是对角阵，得：dZdYdXBLBLBBHNLBHNLHMBHMLBHMLBdZdYdXdZdYdXdHdLdBTsinsincoscos

49、cos0coscoscossincossinsincossin 11AJAJ2.4.1 空间直角坐标系与相应大地坐标系的关系4、B、L、H与椭球元素a, e2 之间的微分关系假设顾及椭球元素的变化，那么前面的微分公式变为：2dedadHdLdBdZdYdXBAJ其中：222222222coscossin12sincossinsincos2coscossincoscosWBWBNaBeNWLBBNaLBNWLBBNaLBNB2.4.1 空间直角坐标系与相应大地坐标系的关系由上式可得：由上式可得：211dedadZdYdXdHdLdBTTBAJAJ 假设空间坐标系的原点和坐标轴指向坚持不变，即椭球

50、的定位与定向不变，那么：0dZdYdX2.4.1 空间直角坐标系与相应大地坐标系的关系上式简化成大地坐标与椭球元素间的微分关系：上式简化成大地坐标与椭球元素间的微分关系：222222212sin002sin2sincossincos dedaBNWHMWBeBBNHMWBBededadHdLdBTBAJ2.4.2 空间直角坐标系之间的旋转变换 两个右手旋转坐标系之间的旋转角，取逆时针旋两个右手旋转坐标系之间的旋转角，取逆时针旋转为正，顺时针旋转为负，旋转矩阵为正交阵，可表转为正，顺时针旋转为负，旋转矩阵为正交阵，可表示为：示为：XXXXXXcossin0sincos0001RYYYYYYcos

51、0sin010sin0cosR1000cossin0sincosZZZZZZR2.4.2 空间直角坐标系之间的旋转变换方法一：XYZXYZOX Y Y 将将X、Y、Z转换到转换到X、Y 、 Z 坐 标 系 ：坐 标 系 ：先绕先绕Z将将X旋转到旋转到XOY平面与平面与XOY平面的交平面的交线线X ，再绕，再绕X 轴将轴将Z旋转到旋转到Z轴，最后再绕轴，最后再绕Z轴，将轴，将X 旋转到旋转到X轴方轴方向。由于三坐标轴的正向。由于三坐标轴的正交关系，经最后一次旋交关系，经最后一次旋转的转的Y必位于必位于Y轴上。轴上。 2.4.2 空间直角坐标系之间的旋转变换坐标变换公式为：坐标变换公式为：XYZX

52、YZOX Y Y ZYXZYXZZXXZZ12RRR12ZZXXZZ RRRR旋转矩阵：是正交矩阵。2.4.2 空间直角坐标系之间的旋转变换假设、分别表示X与X和Y与Y之间的夹角，那么有：XZZZZXZZZZcossinsincoscoscoscoscoscossinsincos21212121假设表示Z与Z之间的夹角，那么有：Xcoscos2.4.2 空间直角坐标系之间的旋转变换方法二：XYZXYZ OX Y ZYZXZXY将将X、Y、Z转换到转换到X、Y 、 Z 坐 标 系 ：坐 标 系 ：先绕先绕X将将Y旋转到旋转到YOZ平面与平面与YOZ平面的交线平面的交线Y ，再绕，再绕Y 轴将轴将

53、Z 旋转到旋转到Z轴，最后再绕轴，最后再绕Z轴，将轴，将X 旋转到旋转到X轴方轴方向。由于三坐标轴的正向。由于三坐标轴的正交关系，经最后一次旋交关系，经最后一次旋转的转的Y必位于必位于Y轴上。轴上。 2.4.2 空间直角坐标系之间的旋转变换坐标变换公式为：坐标变换公式为： ZYXZYXXXYYZZRRRXXYYZZ RRRR其中，旋转矩阵：是正交矩阵。XYZXYZ OX Y ZYZXZXY2.4.2 空间直角坐标系之间的旋转变换 假设、 分别表示X与X、Y与Y和Z与Z之间的夹角，那么有：YXZYXZXZYcoscoscoscossinsincoscoscoscoscoscos2.4.2 空间直

54、角坐标系之间的旋转变换当旋转角是小角度时，可略去其二次项，取：当旋转角是小角度时，可略去其二次项，取：1cos ,sin111XYXZYZR旋转矩阵简化为：ZYXZYXXYXZYZ111坐标转换模型简化为：2.4.2 空间直角坐标系之间的旋转变换当旋转角较大时当旋转角较大时,因旋转矩阵是正交阵，满足条件：因旋转矩阵是正交阵，满足条件：IRR T假设：333231232221131211rrrrrrrrrR那么根据正交条件，得：111000233232231223222221213212211332332223121331332123111231322122111rrrrrrrrrrrrrrrr

55、rrrrrrrrrrr2.4.2 空间直角坐标系之间的旋转变换 旋转矩阵中只需旋转矩阵中只需5个独立未知数。在进展坐标转换时，个独立未知数。在进展坐标转换时，可以直接以旋转矩阵中的可以直接以旋转矩阵中的9个元素为未知数，加上个元素为未知数，加上6个个约束条件直接解算。求得旋转矩阵元素后，进展坐标约束条件直接解算。求得旋转矩阵元素后，进展坐标转换，不用解算旋转角。转换，不用解算旋转角。 这样可防止大旋转角时，线性化过程的复杂方式。这样可防止大旋转角时，线性化过程的复杂方式。 习 题1、假设采用克拉索夫斯基椭球，知大地坐标：、假设采用克拉索夫斯基椭球，知大地坐标： 计算三维空间坐标，并反算检核。计

56、算三维空间坐标，并反算检核。2、在上题中，大地经纬度和大地高分别变化了、在上题中，大地经纬度和大地高分别变化了 用微分公式计算三维空间坐标的变化量。用微分公式计算三维空间坐标的变化量。3、在球近似下，给出球心经纬度和高程与三维空间坐标、在球近似下，给出球心经纬度和高程与三维空间坐标的微分关系式。的微分关系式。4、假设要求相对误差小于、假设要求相对误差小于10-7，那么当旋转角超越多少，那么当旋转角超越多少时，不能采用略去二次项的线性近似。时，不能采用略去二次项的线性近似。mHLB391.108 4015.5013121 ,1283.16823100 mHLB53 . 02 . 0 2.4 空间

57、大地直角坐标系及其转换模型2.4.3 站心地平坐标系及其运用1、站心地平直角坐标系与空间大地直角坐标系的转换关系定义：站心点的法线为z轴，向上为正在地平面上以子午线方向为x轴，y与x、z轴正交，指向以东为正。O000, LBPxzyXYZLBPKQ 将站心坐标轴 xyz 变换成与空间坐标系的指向一致，需求如下几步：(1). z 坐标轴反向；(2). 绕y轴90。+B；(3). 绕z轴旋转-L。2.4.3 站心地平坐标系及其运用 将站心系坐标轴变换到与三维空间直角坐标将站心系坐标轴变换到与三维空间直角坐标轴指向一致时的旋转矩阵为：轴指向一致时的旋转矩阵为：O000, LBPxzyXYZLBPKQ

58、000000000000000sin0cossincoscossinsincoscossincossin 10001000190BBLBLLBLBLLBBLyzRRR顾及，站心系原点在空间坐标系中的坐标为：002000000000sin1sincoscoscos000BHeNLBHNLBHNZYXPPP2.4.3 站心地平坐标系及其运用那么，站心系坐标到空间直角坐标系的变换公式为：那么，站心系坐标到空间直角坐标系的变换公式为：zyxBBLBLLBLBLLBBHeNLBHNLBHNzyxZYXZYXPPP000000000000002000000000sin0cossincoscossinsin

59、coscossincossin sin1sincoscoscos 000R2.4.3 站心地平坐标系及其运用由上式得，空间直角坐标系到站心系的变换公式为：由上式得，空间直角坐标系到站心系的变换公式为：000000 sinsincoscoscos0cossincossinsincossin 000000000000PPPPPPTZZYYXXBLBLBLLBLBLBZZYYXXzyxR2.4.3 站心地平坐标系及其运用2、站心极坐标系与站心地平直角坐标系的关系定义：以站心系原点到点的空间间隔、方位角和天顶距为坐标参数来确定三维点位，称为站心极坐标系。zxyoAZDPZDAZDAZDzyxcossi

60、nsincossin由上式，得：ZzZAyAxzAyAxxyDZAcossinsincossincosarctanarctan2.4.3 站心地平坐标系及其运用也可以用以下公式计算：也可以用以下公式计算：22222arctanarctanzyxzyxxyDZA 公式中的天顶距和方位角都归算到以法线为基准。丈量时以垂线为基准的，需求作垂线偏向矫正。矫正公式下面将讲到。2.4.3 站心地平坐标系及其运用3、空间直角坐标系与站心地平直角坐标系的旋转矢量之间的关系 假设x、y和z为空间坐标系的旋转矢量， x、 y和z为站心坐标系的旋转矢量。顾及旋转矢量是平移不变量，旋转关系与坐标矢量一样。zyxZYX