太原理工大学GNSS卫星定位原理及应用-第二章课件

1、第二章第二章 坐标系统和时间系统坐标系统和时间系统第一节第一节 天球坐标系与地球坐标系统天球坐标系与地球坐标系统一、坐标系概述一、坐标系概述坐标系统是描述坐标系统是描述点点空间位置的参照系统，空间任何一点的位空间位置的参照系统，空间任何一点的位置均可用最少置均可用最少3个坐标参数唯一表示，三个参数随着坐标系或坐个坐标参数唯一表示，三个参数随着坐标系或坐标系的组合形式不同而不同。标系的组合形式不同而不同。常用的点的坐标表示形式：常用的点的坐标表示形式：空间直角坐标（空间直角坐标（3维形式）维形式）-（X Y Z）大地坐标（曲线坐标、大地坐标（曲线坐标、3维形式）维形式）-（B L H）极坐标（极

2、坐标（3维形式）维形式）-（r ）高斯平面投影坐标高斯平面投影坐标+正常高（正常高（2维维+1维组合）维组合）-（x y）、）、HN，点的高斯坐标由大地坐标（点的高斯坐标由大地坐标（B、L）投影后形成，高程为正常高）投影后形成，高程为正常高把一个坐标系统下同时拥有把一个坐标系统下同时拥有3个坐标参数的坐标系称个坐标参数的坐标系称3维坐标维坐标系；把仅有两个坐标参数（一般为系；把仅有两个坐标参数（一般为B,L）的坐标系称二维坐标系）的坐标系称二维坐标系二、二、GPS定位中常用到的坐标系定位中常用到的坐标系天球坐标系天球坐标系CIS空间直角坐标系空间直角坐标系球面坐标系球面坐标系地球坐标系地球坐标

3、系地球地心地固坐标系地球地心地固坐标系CTS坐标系、瞬时地球坐标系坐标系、瞬时地球坐标系3维（维（BLH）（）（XYZ）地球参心地固坐标系地球参心地固坐标系2维（维（BL）-（x,y）ITRS-WGS84-GPSE2000-GALILEO？PZ90-GLONASSCGCS2000-BDS1954北京坐标系北京坐标系Beijing 541980西安坐标系西安坐标系GDZ80地方独立坐标系地方独立坐标系国际通用，也是各国国际通用，也是各国坐标系的母体坐标坐标系的母体坐标地球站心坐标系地球站心坐标系3维（维（XYZ）)(h站心赤道直角坐标系站心赤道直角坐标系站心地平直角坐标系站心地平直角坐标系站心极

4、坐标系站心极坐标系1.天球坐标系天球坐标系以圆球为坐标基准，用来描述天体位置的坐标系；主以圆球为坐标基准，用来描述天体位置的坐标系；主要有空间直角坐标和球面坐标两种等效坐标形式要有空间直角坐标和球面坐标两种等效坐标形式以以S点的坐标为例点的坐标为例两种坐标形式为：两种坐标形式为：),(),(rSZYXS22222tansintancossincoscosYXZarcrZXYarcrYZYXrrX两种天球等效坐标形式的转换关系为：两种天球等效坐标形式的转换关系为：由于天轴（地轴的延伸）指向随时间变化，其变化由于天轴（地轴的延伸）指向随时间变化，其变化表现为岁差和章动等天文现象，天极随时间变化形成

5、表现为岁差和章动等天文现象，天极随时间变化形成协议极、瞬时平天极、瞬时真天极三种；天球坐标系协议极、瞬时平天极、瞬时真天极三种；天球坐标系也随时间分为也随时间分为CIS、CISM、CIST三种时空形式三种时空形式三者之间的关系由岁差矩阵、章动矩阵决定三者之间的关系由岁差矩阵、章动矩阵决定大地测量学基础大地测量学基础天球坐标系天球坐标系协议天球坐标（协议天球坐标（CIS）瞬时平天球坐标（瞬时平天球坐标（CISM)瞬时真天球坐标（瞬时真天球坐标（CIST)P-岁差改正矩阵岁差改正矩阵1） CIS CISMCISMtzyxpzyxAAAAYAXAZZRRZRP,)()()90(岁差参数2）瞬时平天球

6、坐标系CISM 瞬时真天球坐标系CISTCIStzyxNPzyx3）综合上两步GIS GISTMttzyxNzyx?-,)()()(YZXRRRN章动参数2.地球地固坐标系：地球地固坐标系：点在地球坐标系中的主要坐标形式有：空间直角坐点在地球坐标系中的主要坐标形式有：空间直角坐标与大地坐标两种形式标与大地坐标两种形式以以T点的坐标为例，两种坐标点的坐标为例，两种坐标形式为：形式为：),(),(HLBTZYXT两种地球等效地球坐标形式的转换关系为：两种地球等效地球坐标形式的转换关系为：BHeNZLBHNYLBHNXsin)1(sincos)(coscos)(2)1(sin/)1(/ )(arct

7、an/arctan2222eNBZHHeNYXHNZBXYL由于地轴指向随时间变化，其变化表现为地轴北极的极移，由于地轴指向随时间变化，其变化表现为地轴北极的极移，使地球极随时间变化形成瞬时极、平极（协议极），相应的赤使地球极随时间变化形成瞬时极、平极（协议极），相应的赤道又分为平赤道面、协议赤道面；极移又称道又分为平赤道面、协议赤道面；极移又称S.C.Chandler运动；运动；地球坐标系又随时间分为地球坐标系又随时间分为CTS、CTST两种时空形式两种时空形式1)瞬时真地球坐标（瞬时真地球坐标（CTST）-协议地球坐标协议地球坐标(CTS)转换转换tCTSzyxMzyx11001ppppy

8、xyxM旋转矩阵旋转矩阵极移矩阵极移矩阵3.瞬时真天球坐标瞬时真天球坐标-瞬时真地球坐标（瞬时真地球坐标（天地转换天地转换）TZtZYXGASTRzyx)(1000cossin0sincosGASTGASTGASTGASTRZ旋转矩阵，其中旋转矩阵，其中GAST为格林尼治视恒星时角为格林尼治视恒星时角4.天地坐标转换过程天地坐标转换过程CIS-CISM-CIST-CTST-CTS上述转换过程可用于上述转换过程可用于GNSS卫星卫星（飞行中的已知点）（飞行中的已知点）坐坐标的转换标的转换如如GPS卫星转换：由卫星星历提供卫星转换：由卫星星历提供Kepler根数计算卫根数计算卫星的轨道坐标星的轨道

9、坐标-协议天球坐标系中的坐标协议天球坐标系中的坐标-瞬时平天瞬时平天球坐标系球坐标系-瞬时真天球坐标瞬时真天球坐标-瞬时地球坐标瞬时地球坐标-协议地协议地球坐标（美国的球坐标（美国的WGS84）WGS-84定义：原点位于地球质心，Z轴指向BIH1984.0定义的协议地球极(CTP-Conventional Terrestrial System)方向，X轴指向BIHl984.0的零子午面和CIP-Conventional international pole赤道的交点，Y轴与Z-O-X平面构成右手坐标系。WGS-84椭球参数:长半轴a6378137土2m； 地心引力常数(含大气层)GM(3986

10、005l 0.6)l08（m3s-2）正常化二阶带谐系数 48416685*10-61.3 l0-9地球自转角速度 729211510-11土0.150010-11(rads-1)利用以上4个基本常数，可以计算出其它的椭球常数，如第一、第二偏心率e2、e2和扁率定义定义原点位于地球质心，Z轴指向IERS推荐的协议地球极(19001905平均极)方向，X轴指向地球赤道与BIH定义的零子午面交点，Y轴与Z-O-X平面构成右手坐标系。PZ-90坐标基准参数地球平均半径Re=6378136m地球自转角速度 =0.000 072 921 15rad/s地球引力常数GMe=398 600.4km3/s-2

11、地球重力场系数 -0.001 082 63 02C定义定义原点位于包括和大气的整个地球质心，Z轴指向初始为(IRP) BHI1984.0方向，经推算指向至参考历元为2000.0 地球参考极，X轴指向地球赤道与IERS定义的零子午面(IRM)垂直于Z轴的方向，Y轴与Z-O-X平面构成右手坐标系；2008年7月1日启用。0 . 05 . 20 . 00 . 10 . 00 . 00 . 00 . 1109 . 10 . 0109 . 10 . 1906684PZWGSzyxzyx两者在空间三维坐标状态下之差为两者在空间三维坐标状态下之差为cm级级2-1WGS-84/GPS与与PZ-90/ GLON

12、ASS坐标系转换模型坐标系转换模型：CGCS2000坐标基准参数地球平均半径Re=6378137m地球自转角速度 =7.2 92115 10-5 rad/s-1地球引力常数GMe=3.98 6004418 104km3/s-2地球动力形状因子J2=-0.001082629832258正常椭球与参考椭球一致CGCS2000由三个层次的站网与坐标实现与维持第1层次：A级GNSS卫星定位连续运行参考网，用于全球地球动力学、地壳形变测量、卫星精密定轨测量，用GPS、GLONASS、VLBI、SLR等技术手段测量。第2层次：B级GNSS大地网，国家二等大地控制网，用于建立地方或城市坐标基准框架，区域地球

13、动力学和地壳形变测量、精密工程测量等2000椭球几何参数椭球几何参数a6378137mf=1/298.257222101 第3层次：GNSS C、D、D、E级网和天文大地网，用于建立三等及以下控制网测量和工程测量都属于CTS坐标，三维坐标；都至少有（XYZ）（BLH）两种坐标形式和投影坐标（x,y）形式都源于ITRS(或其实现ITRF)都是本系统单点定位或称绝对定位的坐标，以及相对定位中解算的基线向量的本坐标系同属地球地心、地固坐标系三种坐标系在三维层面的坐标之差为cm级，三维层面可有全球统一的转换参数，如PZ-90/WGS-84之转换关系式（2-1）式1.1954北京坐标系坐标基准：克拉索夫

14、斯基椭球，仅几何参数，无明确物理参数，正常重力计算采用赫尔黙特公式；几何参数大地原点：前苏联普尔科沃镇2.1980西安坐标系，GDZ80坐标基准：IUGG1975国际椭球大地原点：陕西省泾阳县永乐镇3 .298/16378245a1980西安坐标系（参心坐标系）还进行了如下定义：椭球短轴Z轴平行于由地球地心指向1968.0地极原点(JYD)的方向；大地起始子午面平行于格林尼治平均天文台子午面，X轴在大地起始子午面与Z轴垂直指向经度零方向；y轴与Z，X轴成右手坐标系。从而确保了GDZ80坐标系椭球定向、定位的明确性，克服了1954北京坐标缺点。IUGG1975国际椭球赋予了明确物理意义，给出了相

15、应椭球参数：1118223810729211510) 1108263(10) 33986005(56378140sradJsmGmaM257.298/16378140a可导出如下几何参数3.国家大地原点国家大地原点大地原点：国家水平大地控制网点大地坐标的起算点。大地起算数据：体现在大地原点上，共有四个个，包括大地原点的大地坐标值L0、B0、H0，以及大地原点至某一方向的大地方位角A0，这些数据构成了经典大地测量的基准。陕西泾陕西泾阳县永阳县永乐镇大乐镇大地原点地原点外部内景第三节第三节 坐标系统之间的转换坐标系统之间的转换一、转换的主要形式（复习一、转换的主要形式（复习大地测量基础大地测量基础

16、内容）内容）1.同一坐标系统内不同坐标形式的变换同一坐标系统内不同坐标形式的变换1）大地坐标与空间直角坐标间等价变换）大地坐标与空间直角坐标间等价变换)()(XYZBLH 2）球面（大地）坐标与平面坐标间投影变换）球面（大地）坐标与平面坐标间投影变换),()(yxBL 3）空间直角坐标与极坐标间等价变换）空间直角坐标与极坐标间等价变换)()(hrXYZ有关模型在有关模型在大地测量基础大地测量基础和本章和本章1，2节已述及节已述及NHH 2.不同坐标系统间的转换不同坐标系统间的转换统一规定：下标统一规定：下标A目标坐标系；下标目标坐标系；下标B源坐标系源坐标系1）不同空间直角坐标间坐标转换）不同

17、空间直角坐标间坐标转换BAXYZXYZ)()(AXYXZYZBZYXmZYXZYX111)1 (000平移矩阵尺度系数旋转矩阵两空间直角坐标系坐标转换模型中有7个转换参数。其中平移参数3个、旋转参数3个、尺度参数1个。除上述布尔萨模型外，还有莫洛金斯基、武测等模型除上述布尔萨模型外，还有莫洛金斯基、武测等模型2）不同球面（大地）坐标间坐标转换）不同球面（大地）坐标间坐标转换BABLHBLH)()(000sincos0coscossincossinsincoscoscoscossinsinXYLLNHBZNHBdLBLBLBdBMHMHMHdHBLBLB222220cossin1sincos0s

18、incossincossinsincoscos1sin0XYZtgBLtgBLNLLeBBMNeBBLNeBBLNeBm 22222222002sinsincossincos11sin1sinsin1MeBNeBBBBMH aMHNMeBeBBadad3）不同平面坐标间坐标转换）不同平面坐标间坐标转换BAyxyx),(),(BooAyxmyxyxcossinsincos)1 (4）GNSS测量中，不同坐标间转换主要是下列两种测量中，不同坐标间转换主要是下列两种BAyxyx),(),()2(BAXYZXYZ)()() 1 (另有三维七参、二维七参、二维四参等多种转换模型；参见大地测量学基础和其它

19、专业书籍二、不同坐标间转换基本步骤二、不同坐标间转换基本步骤 （1）模型选择 大范围坐标转换选择二维七参数转换模型；小区域坐标转换可选择三维四参数模型或平面四参数模型。对于相对独立的平面坐标系统与GNSS坐标系的联系可采用平面四参数模型或多项式回归模型。（2）重合点选取布GNSS网时应初步选用具有目标坐标系坐标的点做重合点，这些点的GNSS坐标（源坐标系）一经测定就有两套坐标可用。最终选定重合点时还要根据所确定的转换参数，计算重合点坐标残差，根据其残差值的大小来最后确定，若残差大于3倍中误差则剔除，重新计算坐标转换参数，直到满足精度要求为止；用于求二维平面四参数转换的最终重合点数量不少于2个，

20、七参数转换的不少于3个。（3）模型参数计算以转换模型为基础方程，将重合点的两套坐标数据代入模型求转换参数，有多余重合点时应应用最小二乘原理求转换参数的最或然解。 （4）精度评估与检核选择部分重合点作为外部检核点，不参与转换参数计算，用转换参数计算这些点的转换坐标与已知坐标进行比较进行外部检核。应选定至少6个均匀分布的重合点对坐标转换精度进行检核 三、附标准转换模型（供实用）三、附标准转换模型（供实用）（1 1）二维七参数转换模型）二维七参数转换模型 2222sincos0coscossincossinsincos0cossin1sincos0sincos00(2sin)sincos1xyzLL

21、XLNBNBYBBLBLBZMMMtgBLtgBLmNLLeBBMNeBeBBMasincosaBBff 同一点位在两个坐标系下的纬度差经度差，单位为弧度 椭球长半轴差（单位米）、扁率差（无量纲）， 平移参数，单位为米， 旋转参数，单位为弧度， 尺度参数（无量纲）,BL, af,XYZ,xyzm（2 2）平面四参数转换模型）平面四参数转换模型 其中，x0，y0为平移参数，为旋转参数，m为尺度参数。x2，y2为2000国家大地坐标系下的平面直角坐标，x1，y1为原坐标系下平面直角坐标。坐标单位为米。110022cossinsincos)1 (yxmyxyx（3 3）BursaBursa七参数坐标

22、转换模型七参数坐标转换模型SSSSSSZYXSSSSSSTTTZYXZYXmXYXZYZZYXZYX0002222sincos0()cos()cossincossinsincos()()()coscossinsinsin(1)(1)cossin1()sinLLXLNHBNHBYBBLBLBMHMHMHZHBLBLBNeHNeHtgBLtgBLNHNHNHNeB222222222222()sin0sincos0sincossinsincoscos0sincos()sin00(2sin)sincossincos1(1sin)(11xyzBNHNeBLLMHMHNeBBLNeBBLNeBBmMNHN

23、eeBNBBeBBfMaNMeBaa222sin)sinafeBB（4 4）三维七参数坐标转换模型）三维七参数坐标转换模型 （5 5）坐标转换精度评定和评估方法）坐标转换精度评定和评估方法 依据计算坐标转换模型参数的重合点（没有参与转换参数计算的重合点）的残差中误差评估坐标转换精度。对于n个点，坐标转换精度估计公式如下： V（残差）=重合点转换坐标-重合点已知坐标空间直角坐标X残差中误差空间直角坐标Y残差中误差1XXvvMn 1YYvvMn 空间直角坐标Z 残差中误差点位中误差平面坐标x残差中误差1ZZvvMn 222ZYXpMMMM1xxvvMn 平面坐标y残差中误差大地高H残差中误差 平面

24、点位中误差为1yyvvMn 1HHvvMn 22yxpMMM（6）决定不同坐标间转换精度的因素）决定不同坐标间转换精度的因素重合点在网中的分布状况转换模型精度与适用性重合点的坐标分量及点位精度三、三、GNSS测量中涉及的其它坐标测量中涉及的其它坐标1.站心坐标系概述站心坐标系概述 有站心赤道直角坐标系、站心地平直角坐标系和站心地平极坐标系三种坐标系形式。 GNSS坐标计算中，常用的为站心赤道直角坐标系与站心地平直角坐标系：站心地平、站心赤道坐标和地球空间直角坐标关系图2.站心赤道直角坐标系站心赤道直角坐标系定义：以测站P1为原点，且与地球空间直角坐标系，O-XYZ相应坐标轴平行的坐标系站心赤道

25、坐标系与地球空间直角坐标间关系为简单的平移关系ZYXP 1BHeNLBHNLBHNZYXZYXsin)1(sincos)(coscos)(23.站心地平直角坐标系P1-xyz定义：以P1为原点，过测站P1的法线为z轴，指向天顶为正，以子午线方向为x轴，向北为正，以与x、z垂直方向为y轴，东为正。地平直角坐标系与站心赤道直角坐标系关系可通过旋转变换实现，关系式如下： 地平站心赤道zyxPBRLRZYXxYz)90()180(地平zyxBBLBLLBLBLLBsin0cossincoscossinsincoscossincossin4.GNSS空间直角坐标与地平直角坐标的关系BHeNLBHNLBH

26、NzyxBBLBLLBLBLLBZYXsin)1 (sincos)(coscos)(sin0cossincoscossinsincoscossincossin2?5.站心地平极坐标系P1rAh（用于描述测站与卫星的相对位置） 定义：以测站P1为原点，坐标参数为： r-测站P1至卫星s的距离 A-卫星的方位角，zox平面与zos平面的夹角，左旋为正h-卫星的高度角os与xoy平面的夹角。1 1）站心地平极坐标系与站心地平直角坐标系关系）站心地平极坐标系与站心地平直角坐标系关系：coshcos Arx coshsin Ary sinhrz 222zyxrA=arc tan(y/x) h=arcta

27、n(z/(x2+y2)1/2) 四、四、GNSS卫星（飞行中的已知点）坐标计算流程卫星（飞行中的已知点）坐标计算流程GNSS定位坐标计算流程地面监测站监测地面监测站监测卫星运动状态卫星运动状态主控站处理主控站处理注入站注入站卫星卫星用户获得用户获得导航电文导航电文导航导航电文电文计算卫星在协议天球坐标系中的坐标计算卫星在瞬时平天球坐标系中的坐标计算卫星在瞬时平天球坐标系中的坐标计算卫星在瞬时地球坐标系中的坐标计算卫星在瞬时地球坐标系中的坐标通过星地距离方程解算基线通过星地距离方程解算基线将接收机坐标转换为北京将接收机坐标转换为北京54或西安或西安80坐标坐标计算卫星在协议地球坐标系计算卫星在协

28、议地球坐标系CTS中的坐标中的坐标计算卫星在轨道平面坐标系中的坐标计算卫星在轨道平面坐标系中的坐标计算卫星在瞬时真天球坐标系中的坐标计算卫星在瞬时真天球坐标系中的坐标卫星参考历元的卫星参考历元的Kepler参参数数+摄动改正参数摄动改正参数=实际实际观测历元的观测历元的Kepler参数参数计算卫星在协议天球坐标系中的坐标计算卫星在协议天球坐标系中的坐标第四节 坐标系统间的转换一、GNSS定位中时间的意义GNSS定位过程是在严格的时间定义下进行的，无论是坐标系的定义、GNSS卫星信号接收、星地距离测量，还是坐标转换都需在严格的时间定义下完成二、时间系统各类与构成时间系统类型：基于各种用途，时间系

29、统有多种： ST恒星时、MT平太阳时、UT世界时、ATI原子时、UTC协调世界时。时间系统构成： 由时间原点(历元)，时间尺度(时间单位)构成。1.恒星时ST(Sidereal Time)原点：以春分点为参考点，由春分点的周日视运动所定义的时间系统为恒星时系统。 ST尺度：春分点连续两次经过本地子午圈的时间间隔为一恒星日，一恒星日分为24个恒星时2.平太阳时MT(Mean SOlar Time) 平太阳：以真太阳周年运动的平均速度在天球赤道上作周年视 运动的太阳，其周期与真太阳一致 。时间原点：以平太阳为参考点，由平太阳的周日视运动所定义的时间系统为平太阳时系统。时间尺度为：平太阳连续两次经过本地子午圈的时间间隔为一平太阳日，一平太阳日分为24平太阳时。 3.世界时UT(Universal Time)时间原点：以平子夜为零时起