由前面的大地坐标微分公式-同济大学.ppt

第二章地球坐标系和地球椭球,2,2.1概述,大地测量采用的坐标系：天球坐标系、地球坐标系地球坐标系：固定在地球上与地球一起自转和公转的坐标系地球坐标系分类：参心坐标系、地心坐标系定义坐标系的要素：原点位置、尺度与坐标轴指向；还包括一些天文、物理、地球等参数，若采用大地坐标表述形式还需要椭球元素。,3,2.2地球椭球面的数学计算和有关计算,2.2.1地球椭球的几何、物理元素椭球方程：扁率：第一偏心率：第二偏心率：,4,2.2.1地球椭球的几何、物理元素（续1）,几个关系式：,1954年北京坐标系，克拉索夫斯基椭球元素：,5,2.2.1地球椭球的几何、物理元素（续2）,1980年大地坐标系采用第16届IAGIUGG椭球，其椭球元素为：,6,2.2.2旋转椭球面的参数表示及数学性质,1、经线和纬线的曲线方程在XOZ坐标面上的起始经线方程：,M0饶Z轴旋转，形成纬圈（平行圈），其半径：,经度为L的经线方程：,7,2.2.2旋转椭球面的参数表示及数学性质（续1）,纬圈方程：,8,2.2.2旋转椭球面的参数表示及数学性质（续2）,2、椭球面法线与子午线主法线的同一性、经纬线的Frenet标架,如图为过M点的子午面。子午线的主法线MP位于子午面内，且垂直于子午线切线T；R为过M点的平行圈切线，显然R垂直于M点的子午面，因此R垂直于MP。所以，MP垂直于椭球面在M点的切平面，因此它是椭球面的法线。,Frenet标架：曲线上任意一点处的三个相互正交的单位向量取切向、主法向和与该两个方向正交的第三个方向构成的三维直角坐标系。,9,2.2.2旋转椭球面的参数表示及数学性质（续3）,3、旋转椭球面及经纬线的参数方程1).以大地经度L及归化纬度u为参数的方程,在XOZ子午面内，有,在三维空间坐标系中：,10,(2).以大地经纬度L、B为参数的方程,切线M0T的斜率的导数式：,由椭圆方程求导得：,代入第一式得：,1,2.2.2旋转椭球面的参数表示及数学性质（续4）,11,2.2.2旋转椭球面的参数表示及数学性质（续5）,将代入椭圆方程，化简后得：,1,引入辅助符号：,则有：,12,2.2.2旋转椭球面的参数表示及数学性质（续6）,以大地经纬度两个参数表示的椭球面上一点的三维坐标为椭球面参数方程式：,13,以大地纬度为参数的经度为LC的子午线参数方程为:,14,在一点BC,LC处的子午线切向量,15,子午线切线单位向量,16,以大地经度为参数的大地纬度为的BC纬线的参数方程为,17,在一点BC,LC处的平行圈切向量,18,平行圈切线单位向量,19,椭球面单位法向量为其矢量积：,20,2.2.2旋转椭球面的参数表示及数学性质（续7）,(3).以大地经度L及球心纬度为参数的方程,球心纬度，向径，则对于XOZ平面上的椭圆有：,在椭圆上，向径由球心纬度唯一确定，将上式代入椭圆方程，得：,21,2.2.2旋转椭球面的参数表示及数学性质（续8）,对于XOZ平面上的椭圆有：,以经度、球心纬度两个参数表示的椭球面上一点的三维坐标为参数方程式为：,22,2.2.2旋转椭球面的参数表示及数学性质（续9）,不难得出，u,B,的关系为：,因此有：,由球心纬度公式，得：,23,2.2.2旋转椭球面的参数表示及数学性质（续10）,4、旋转椭球面的几何性质a).对称性b).有界性c).正则性：曲面上每点都对应于唯一确定的非零法向量。d).不可展性,24,2.2.3法截线曲率及曲率半径,1、空间曲线的曲率几曲率半径若以曲线的弧长s为参数，曲线上的点位用向量r(s)表示。则曲线的曲率为：,若以t参数，则曲线的曲率可表示为：,25,2、椭球面法截线的曲率(1).子午线曲率半径不失一般性，以起始子午线为例推导。若以归化纬度u为子午线方程的参数，则有：,2.2.3法截线曲率及曲率半径（续1）,26,2.2.3法截线曲率及曲率半径（续2）,则有：,同理，若以大地纬度为参数，得：,子午曲率半径M，就是曲率是倒数，即：,27,2.2.3法截线曲率及曲率半径（续3）,(2).卯酉线曲率半径定义：与子午面切线正交的法截面与椭球面的交线为卯酉线。根据微分几何中的麦尼尔定理，卯酉圈曲率kn与平行圈曲率kr的关系为：,平行圈半径为子午面XOZ平面内的X坐标，即：,则有，上述两式得卯酉曲率半径N为：,28,2.2.3法截线曲率及曲率半径（续4）,(3).任意方向法截线的曲率半径根据微分几何中的Euler公式，任意方向法截线的曲率与子午、卯酉曲率半径的关系为：,因此，任意方向的曲率半径为：,当A为0，/2，3/2时，取得极值。,29,2.2.3法截线曲率及曲率半径（续5）,(4).平均曲率半径定义：所有方向法截线曲率半径的平均值。,代入上式，得：,30,2.2.3法截线曲率及曲率半径（续6）,不难得到：NRM引入辅助量：,存在下列关系：,31,2.2.4椭球面上第一基本形式及弧长面积计算,1.椭球面的第一基本形式,椭球面上点的向量：,椭球面上的微分弧长：,其中：,对于椭球面：,32,2.2.4椭球面上第一基本形式及弧长面积计算（续1）,2、子午线弧长子午线微分弧长：,积分得：,用二项式展开，并逐项积分得：,常数A、B、C、D、E、F、G的计算公式见教材,33,2.2.4椭球面上第一基本形式及弧长面积计算（续2）,对于小于400km的弧长，可采用以下简化式。,其中：,根据：,求出导数，代入上式并化简，得：,对于小于40km的弧长，可进一步简化为：,34,2.2.4椭球面上第一基本形式及弧长面积计算（续3）,已知B1和弧长S12求B2称为反算，可采用叠代法计算。初值：,叠代格式：,其中：,要求：,35,2.2.4椭球面上第一基本形式及弧长面积计算（续4）,3、平行圈的半径与弧长,相同经差的平行圈弧长在赤道最长，越靠近两极越小。,36,2.2.4椭球面上第一基本形式及弧长面积计算（续5）,4、利用经纬格网计算椭球面的面积,37,2.2.4椭球面上第一基本形式及弧长面积计算（续6）,上式利用二项式展开并积分，得：,取L2-L1=2，B2=/2，B1=0算得半球面积，乘2可以估算全球面积约为5.1亿平方公里,38,习题,1、导出三种纬度、u与B的关系。2、导出子午曲率半径M与卯酉曲率半径N的计算公式。3、M、N、R的关系如何？在什么条件下三者相同？4、某点到赤道的子午弧长,求该点的纬度。a=6378245,=1/298.35、已知某点的纬度,求该点自赤道起的子午弧长。a=6378245,=1/298.3,39,2.2.5大地线,1、大地线的定义与性质法截弧：由椭球面上A点的法线与B点所确定的法截面与椭球面相割得到的曲线称为A到B的法截弧。相对法截弧：A到B的法截弧与B到A的法截弧。由相对法截弧构成的椭球面三角形不是闭合图形。,40,2.2.5大地线（续1）,大地线的定义：大地线的主法线与曲面法线处处重合。大地线的性质：1、大地线上任何点的密切平面就是该点的法截面；2、曲面上连接任何两点的最短直线必为大地线。3、大地线的测地曲率等于0曲线的测地曲率：曲线的曲率在曲面切平面上的投影。大地线的曲率：大地线的挠率,41,2.2.5大地线（续2）,2、大地坐标系中大地线的微分方程(1).大地线的二阶微分方程以u,v为参数的一般曲面的大地线微分方程可表示为：下标为相应的偏导数。,42,2.2.5大地线（续3）,对于椭球面，有：代入前面公式，得：则旋转椭球面上大地线的微分方程为：,43,2.2.5大地线（续4）,(2).克莱劳定理,直角坐标系中的椭球面方程：,椭球面法向量为：,以大地线弧长为参数的大地线主法线向量为：,两者指向一致，即：,44,2.2.5大地线（续5）,由上式的前两个方程得：,将三维空间坐标与大地坐标的关系式及微分式代入：,1,代入式，整理得：,1,2,45,2.2.5大地线（续6）,将关系式：即：大地线上各点的平行圈半径与该点的大地线方位角正弦的乘积是常数。,代入上式，即得克莱劳定理：,46,2.2.5大地线（续7）,(3).大地线的一阶微分关系式,由克莱劳定理，微分得：,47,2.2.5大地线（续8）,又如图所示：,代入上式，得：,三个微分关系式可整理为：,3,48,2.2.5大地线（续9）,3、以弧长和大地方位角为参数的大地线方程,大地线始点坐标P0(B0,L0)，大地线上任何点的位置向量都可以展开成S,A的级数形式：,Frenet标架的坐标轴定义：x指向大地线的切向t，y指向大地线的主法向n，向内为正，z指向大地线的副法向b，构成左手系。,4,49,2.2.5大地线（续10）,显然有：,根据曲线论中的Frenet公式：,由以上两式可求出各阶导数：,50,2.2.5大地线（续11）,将上式代入大地线展开式，得Frenet标架下的三维坐标：,4,5,顾及公式：,51,2.2.5大地线（续12）,和：,求导得：,52,2.2.5大地线（续13）,代入Frenet标架下的三维坐标公式，得：,5,53,2.2.5大地线（续14）,将坐标系饶y逆时针旋转A，得x”、y”、z”坐标系，则有：,以P0点为原点的地平坐标系(站心坐标系)x、y、z，与x”、y”、z”坐标系的关系为：,54,2.2.5大地线（续15）,最后得到地平坐标系（站心系）中的大地线方程，称为Weingarten级数式。,6,55,2.2.5大地线（续16）,法截弧为平面曲线，其挠率为0，同理可推得地平坐标系中的计算式为：,56,2.2.5大地线（续17）,4、基于大地线的椭球面曲线坐标系(1).大地线极坐标系大地圆：到极点具有相同大地线长度的点所构成的轨迹。由大地线长度和大地方位角可描述曲面点的位置。,如图所示：,对照第一基本形式，得：,由图中的微分直角三角形，得大地极坐标系中的微分关系式：,57,2.2.5大地线（续18）,大地线的归化长度m的计算公式：,由式求出偏导数代入得：,6,58,2.2.7以长度量为坐标参数的新大地坐标系,(1)以长度量表示的椭球面上坐标系的由来早在1810年Soldner就提出了球面直角坐标系统。此后Helmert，Grossmann，Heck等德国测量学者基于Soldner球面直角坐标系推广提出了椭球面直角坐标系。教材上所述的测地坐标系与其有类似之处。,59,(2)一种新型的大地坐标系,2005年提出了椭球面上一种新型的大地坐标系：它仍以经纬线作为坐标曲线，且与大地坐标系之间能进行精确的坐标转换；它所采用的坐标参数是以长度而不是以角度为单位；可简化椭球面上的繁复计算。施一民，朱紫阳，范业明.坐标参数为长度量的一种新型的大地坐标系.同济大学学报，2005,33（11）:1537-1540,60,新型大地坐标系的定义,为构建新型的大地坐标系,可在区域中心附近选择一点作为其坐标原点,其在大地坐标系中的大地经纬度设为（B0，L0）。两族互为正交的经纬线构成坐标系的坐标格网。设过经纬度为（B，L）的任一点的经线与起始纬线的交点为PB0,点P至点PB0的经线上弧长作为纵坐标sL,点P0至点PB0的纬线上弧长取为横坐标sB,如图所示。,61,新型大地坐标系的定义（续）,62,新型大地坐标系的定义（续）,过点P的纬线上的微分弧长ds与起始纬线上相应的微分弧长ds之比n称为按纬度变化方向的长度归化因子n=ds/ds在大地坐标系中，与这两个平行圈上微分弧长相应的经差dl相等，故可精确求得n=NsinB/N0sinB0式中N0，N分别为纬度B0，B处的卯酉线曲率半径。基于微分几何和椭球大地测量的理论，n可用新大地坐标表示为(取至二次项),63,4.新型大地坐标与大地坐标的转换关系,64,习题,1.纬度相同的两个点的相对法截弧是否重合？此线是否就是大地线？2.推导大地线的三个微分式。3.试述测地坐标系的定义？测地平行线是否等距？测地大地线是否等距？4.简述weingarten级数的推导步骤。,65,2.3椭球面上大地坐标的计算,2.3.1水平方向、边长观测值归算到椭球面1、水平方向观测值归算到参考椭球面的改正包括三项改正，称为三差改正。(1).垂线偏差改正(2).标高差改正,用椭球半径的近似值代入得：,66,2.3.1水平方向、边长观测值归算到椭球面,(3).法截弧方向归算到大地线方向的改正,该项改正很小，100公里约0.03“，只有一等控制网才估计此项改正。,67,2.3.1水平方向、边长观测值归算到椭球面,2、空间边长归算至参考椭球面的改正测线端点的大地高为：,椭球面上弦长d的计算公式,省略H/R的二次项，得：,68,2.3.1水平方向、边长观测值归算到椭球面,椭球面上的弧长为：,69,2.3.1水平方向、边长观测值归算到椭球面,3.工程控制网中的地面观测元素的归算,以平均高程面作投影面，范围小，可以用球代替椭球；球半径采用高斯平均曲率半径。计算公式为：,不难证明：椭球半径的误差对边长归算结果影响很小，R取6371km即可，但高差误差对边长归算比较敏感。,70,2.3.2椭球面上三角形解算,1、球面角超,三块面积之和为：,代入球面角超定义式，得：,71,2.3.2椭球面上三角形解算,按球面三角公式：,当边长小于40公里时，第二项影响小于0.0004“，可略去,72,2.3.2椭球面上三角形解算,2、解算球面三角形的勒让德定理勒让德定理：对于较小的球面三角形，可用平面三角公式来解算，只需使三个平面角等于相应的球面角减去三分之一的球面角超，而边长保持不变。,73,2.3.3大地主题解算,大地主题解算分类：正算：已知(B1,L1)，A12，S12，计算(B2,L2)，A21反算：已知(B1,L1)，(B2,L2)，计算A12，S12，A21短距离中距离长距离中短距离解算方法：按级数展开采用Gauss平均引数公式；长距离解算方法：贝塞尔公式,74,2.3.3大地主题解算,1、纬度差、经度差和方位角差展开为大地线长度的级数式,由大地线的微分公式，得其一阶导数为：,75,2.3.3大地主题解算,二阶和三阶导数采用复合函数求导法计算：,同理可求出四阶以上的导数和L、A的高阶导数，代入展开式即可。,76,2.3.3大地主题解算,2、高斯平均引数公式若取大地线中点展开，得：,两式相减，得：,类似地，有：,1,77,2.3.3大地主题解算,两式相加，得：,类似地，有：,其中：,将展开成级数，得：,2,78,2.3.3大地主题解算,由大地线的微分公式：,求导，得：,可取：,代入式，求出各阶导数后整理得：,1,79,2.3.3大地主题解算,同理可得：,以上3式具有4次方精度，可用于解算200公里下的大地主题。,3,80,2.3.3大地主题解算,因计算Bm,Lm要用到B2,L2，因此需要叠代计算。其初值为：,叠代计算公式为：,直到为止。,最后计算纬度、经度和方位角：,81,2.3.3大地主题解算,3、高斯平均引数反算公式由正算公式，反解得：,右端第二项与第一项相比为小量，可以作近似：,82,2.3.3大地主题解算,代入上式第二项，得：,由此可求得平均方位角和大地线长度如下：,83,2.3.3大地主题解算,由正算公式的第三式，计算a：,最后得起终点的大地方位角为：,84,2.3.3新大地主题解算,4、新大地坐标系与大地坐标系间的坐标转换(1).由（B,L）求解（sL,sB）(2).由（sL,sB）求解（B,L）:B由下式迭代反解,而L由下式直接反解.,85,2.3.4大地主题微分公式,1、大地主题正解微分公式终点的经纬度（B2,L2）和大地线方位角A21，与起点的经纬度（B1,L1）和大地线方位角A12，以及大地线长度S的微分关系。,86,2.3.4大地主题微分公式,2、大地主题反解微分公式起点大地线方位角A12和大地线方位角A21，以及大地线长度S与起点和终点的经纬度（B1,L1）和（B2,L2）的微分关系。,87,习题,1、地面观测方向归算到椭球面上需要加哪几项改正？2、地面观测距离归算到椭球面上二步改正的几何意义？3、P1与P2与为控制点，已知：计算归算到椭球面上的长度4、已知利用Gauss平均引数公式正反算。,88,2.4空间大地直角坐标系及其转换模型,2.4.1空间直角坐标系与相应大地坐标系的关系1、X、Y、Z与B、L、H间的关系空间坐标系的定义：Z自转轴，X位于赤道面，指格林尼治天文台，Y指东，构成右手系。大地坐标的定义：B为过一点的椭球面的法线与赤道面交角、L为过同一点的子午面与起始子午面二面角的平面角，H为点沿法线到椭球面的距离。大地高与正高、正常高之间的关系：,89,2.4.1空间直角坐标系与相应大地坐标系的关系,如图所示：,90,2.4.1空间直角坐标系与相应大地坐标系的关系,2、由X、Y、Z计算B、L、H的迭代解法,计算L：,迭代计算B：,迭代初值为：,最后计算H：,91,2.4.1空间直角坐标系与相应大地坐标系的关系,3、X、Y、Z与B、L、H间的微分关系,由前面式微分得；,其中：对角阵,92,2.4.1空间直角坐标系与相应大地坐标系的关系,顾及A是正交阵，J是对角阵，得：,93,2.4.1空间直角坐标系与相应大地坐标系的关系,4、B、L、H与椭球元素a,e2之间的微分关系若顾及椭球元素的变化，则前面的微分公式变为：,其中：,94,2.4.1空间直角坐标系与相应大地坐标系的关系,由上式可得：,若空间坐标系的原点和坐标轴指向保持不变，即椭球的定位与定向不变，则：,95,2.4.1空间直角坐标系与相应大地坐标系的关系,上式简化成大地坐标与椭球元素间的微分关系：,96,2.4.2空间直角坐标系之间的旋转变换,两个右手旋转坐标系之间的旋转角，取逆时针旋转为正，顺时针旋转为负，旋转矩阵为正交阵，可表示为：,97,2.4.2空间直角坐标系之间的旋转变换,方法一：,将X、Y、Z转换到X、Y、Z坐标系：先绕Z将X旋转到XOY平面与XOY平面的交线X”，再绕X”轴将Z旋转到Z轴，最后再绕Z轴，将X”旋转到X轴方向。由于三坐标轴的正交关系，经最后一次旋转的Y”必位于Y轴上。,98,2.4.2空间直角坐标系之间的旋转变换,坐标变换公式为：,旋转矩阵：,是正交矩阵。,99,2.4.2空间直角坐标系之间的旋转变换,若、分别表示X与X和Y与Y之间的夹角，则有：,若表示Z与Z之间的夹角，则有：,100,2.4.2空间直角坐标系之间的旋转变换,方法二：,将X、Y、Z转换到X、Y、Z坐标系：先绕X将Y旋转到YOZ平面与YOZ平面的交线Y”，再绕Y”轴将Z”旋转到Z轴，最后再绕Z轴，将X”旋转到X轴方向。由于三坐标轴的正交关系，经最后一次旋转的Y”必位于Y轴上。,101,2.4.2空间直角坐标系之间的旋转变换,坐标变换公式为：,其中，旋转矩阵：,是正交矩阵。,102,2.4.2空间直角坐标系之间的旋转变换,若、分别表示X与X、Y与Y和Z与Z之间的夹角，则有：,103,2.4.2空间直角坐标系之间的旋转变换,当旋转角是小角度时，可略去其二次项，取：,旋转矩阵简化为：,坐标转换模型简化为：,104,2.4.2空间直角坐标系之间的旋转变换,当旋转角较大时,因旋转矩阵是正交阵，满足条件：,若：,则根据正交条件，得：,105,2.4.2空间直角坐标系之间的旋转变换,旋转矩阵中只有5个独立未知数。在进行坐标转换时，可以直接以旋转矩阵中的9个元素为未知数，加上6个约束条件直接解算。求得旋转矩阵元素后，进行坐标转换，不必解算旋转角。这样可避免大旋转角时，线性化过程的复杂形式。,106,习题,1、若采用克拉索夫斯基椭球，已知大地坐标：计算三维空间坐标，并反算检核。2、在上题中，大地经纬度和大地高分别变化了用微分公式计算三维空间坐标的变化量。3、在球近似下，给出球心经纬度和高程与三维空间坐标的微分关系式。4、若要求相对误差小于10-7，则当旋转角超过多少时，不能采用略去二次项的线性近似。,107,2.4空间大地直角坐标系及其转换模型,2.4.3站心地平坐标系及其应用1、站心地平直角坐标系与空间大地直角坐标系的转换关系定义：站心点的法线为z轴，向上为正在地平面上以子午线方向为x轴，y与x、z轴正交，指向以东为正。,将站心坐标轴xyz变换成与空间坐标系的指向一致，需要如下几步：(1).z坐标轴反向；(2).绕y轴90。+B；(3).绕z轴旋转-L。,108,2.4.3站心地平坐标系及其应用,将站心系坐标轴变换到与三维空间直角坐标轴指向一致时的旋转矩阵为：,顾及，站心系原点在空间坐标系中的坐标为：,109,2.4.3站心地平坐标系及其应用,则，站心系坐标到空间直角坐标系的变换公式为：,110,2.4.3站心地平坐标系及其应用,由上式得，空间直角坐标系到站心系的变换公式为：,111,2.4.3站心地平坐标系及其应用,2、站心极坐标系与站心地平直角坐标系的关系定义：以站心系原点到点的空间距离、方位角和天顶距为坐标参数来确定三维点位，称为站心极坐标系。,由上式，得：,112,2.4.3站心地平坐标系及其应用,也可以用以下公式计算：,公式中的天顶距和方位角都归算到以法线为基准。测量时以垂线为基准的，需要作垂线偏差改正。改正公式下面将讲到。,113,2.4.3站心地平坐标系及其应用,3、空间直角坐标系与站心地平直角坐标系的旋转矢量之间的关系若x、y和z为空间坐标系的旋转矢量，x、y和z为站心坐标系的旋转矢量。顾及旋转矢量是平移不变量，旋转关系与坐标矢量相同。,114,2.4.3站心地平坐标系及其应用,4、站心地平直角坐标系的应用(1).计算基线向量的大地方位角,其中，B0,L0为基线始端的纬度和经度。,(2).绕站心系坐标轴的旋转向量有特殊意义,z为方位旋转角。因涉及大地方位角起始方向的变动,115,2.4.3站心地平坐标系及其应用,(4).计算卫星的高度角和方位角卫星Q的方位角和高度角可用其站心坐标xQ、yQ计算。,116,2.4.4两个空间大地直角坐标系间的转换模型,1、Bursa-Wolf模型转换参数包括三个平移参数、三个旋转参数与一个尺度参数。,R为前面所述的旋转矩阵。当旋转角为小角度时，上式可简化为：,117,2.4.4两个空间大地直角坐标系间的转换模型,略去尺度参数和旋转参数的乘积项，上式可进一步简化为：,上式第二式常用于转换参数未知时，利用同名点在两个坐标系中的坐标计算转换参数。,118,2.4.4两个空间大地直角坐标系间的转换模型,上式再应用于Pj，并与上式相减，得Pi与Pj两点坐标差的坐标变换模型如下：,119,2.4.4两个空间大地直角坐标系间的转换模型,2、Molodensky模型,如果旋转与尺度是相对于参考点PK，即以参考点PK作变换中心。则有Molodensky模型。,旋转角为小角度时，上式可简化为：,120,2.4.4两个空间大地直角坐标系间的转换模型,上式同样可以简化为求解转换参数的形式如下：,其中，,相应于Molodensky模型的坐标差的转换模型与Bursa-Wolf模型相同。,121,2.4.4两个空间大地直角坐标系间的转换模型,3、范士转换模型,若旋转角是围绕参考点的站心地平坐标系的坐标轴，即为范士转换模型。将三维空间坐标系的旋转角与站心系旋转角的关系代入Molodensky模型，即得范士转换模型如下：,122,2.4.4两个空间大地直角坐标系间的转换模型,4、卫星网与地面网之间的转换卫星网精度高，地面网平面坐标与高程点不重合。,123,2.4.5大地坐标的微分公式,根据大地坐标与三维空间直角坐标间的微分公式：,大地直角坐标的变动是由于原点平移、坐标轴旋转和尺度变化引起。即：,代入上式，得大地坐标微分公式。,124,2.4.5大地坐标的微分公式,大地坐标微分公式的矩阵形式可表示为：,125,习题,1、给出站心坐标系的定义。2、经过哪几步旋转和平移变换，可将站心系坐标变换到三维空间直角坐标系中。3、导出两点的大地方位角、距离和天顶距与站心坐标的关系。4、三维空间坐标变换有哪几种模型？各种模型间的差异在哪里？5、范士变换模型的旋转参数有什么意义？,126,2.5参心坐标系和参考椭球,2.5.1垂线偏差与Laplace方程1、天文经度、天文纬度和天文方位角天文经度：包含测站垂线的子午面与起始子午面的夹角；天文纬度：测站垂线的与赤道面的夹角；天文方位角：包含测站垂线的子午面与测站垂线和照准面所张成的垂直面的夹角；天文天顶距：测站垂线与观测方向的夹角,127,2.5.1垂线偏差与Laplace方程,因地极移动(x,y)，观测的天文经纬度、方位角需要归算到地极原点，称为极移改正，其公式如下：,观测值在地面取得，归算到椭球面上时，天文纬度和方位角需要作如下改正：,128,2.5.1垂线偏差与Laplace方程,2、垂线偏差和大地水准面差距,数值积分，得：,129,2.5.1垂线偏差与Laplace方程,3、垂线偏差公式和Laplace方位角,如图所示：xyz为大地站心坐标系，x1y1z1为天文站心坐标系。两者的关系为：,1,130,2.5.1垂线偏差与Laplace方程,天文和大地坐标系分别与原点在站心，坐标轴与三维空间直角坐标系指向相同的坐标系的关系如下：,131,2.5.1垂线偏差与Laplace方程,由上面第一式代入第二式，略去高次项，整理得：,132,2.5.1垂线偏差与Laplace方程,并得出Laplace方程：,顾及天文站心系（x1，y1，z1）与大地站心系（x，y，z）的关系：,和天顶距、方位角和站心坐标的关系：,133,2.5.1垂线偏差与Laplace方程,将第二式代入第一式，得：,