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

第二章坐标系统和时间系统第一节天球坐标系与地球坐标系统一、坐标系概述坐标系统是描述点空间位置的参照系统，空间任何一点的位置均可用最少3个坐标参数唯一表示，三个参数随着坐标系或坐标系的组合形式不同而不同。常用的点的坐标表示形式：空间直角坐标（3维形式）--（XYZ）大地坐标（曲线坐标、3维形式）-（BLH）极坐标（3维形式）--（rαδ）高斯平面投影坐标+正常高（2维+1维组合）--（xy）、HN，点的高斯坐标由大地坐标（B、L）投影后形成，高程为正常高把一个坐标系统下同时拥有3个坐标参数的坐标系称3维坐标系；把仅有两个坐标参数（一般为B,L）的坐标系称二维坐标系二、GPS定位中常用到的坐标系天球坐标系CIS空间直角坐标系球面坐标系地球坐标系地球地心地固坐标系CTS坐标系、瞬时地球坐标系3维（BLH）（XYZ）地球参心地固坐标系2维（BL）-（x,y）ITRS-WGS84-GPSE2000-GALILEO？PZ90-GLONASSCGCS2000-BDS1954北京坐标系Beijing541980西安坐标系GDZ80地方独立坐标系国际通用，也是各国坐标系的母体坐标地球站心坐标系3维（XYZ）站心赤道直角坐标系站心地平直角坐标系站心极坐标系1.天球坐标系以圆球为坐标基准，用来描述天体位置的坐标系；主要有空间直角坐标和球面坐标两种等效坐标形式以S点的坐标为例两种坐标形式为：两种天球等效坐标形式的转换关系为：由于天轴（地轴的延伸）指向随时间变化，其变化表现为岁差和章动等天文现象，天极随时间变化形成协议极、瞬时平天极、瞬时真天极三种；天球坐标系也随时间分为CIS、CISM、CIST三种时空形式三者之间的关系由岁差矩阵、章动矩阵决定—《大地测量学基础》天球坐标系协议天球坐标（CIS）瞬时平天球坐标（CISM)瞬时真天球坐标（CIST)P---岁差改正矩阵1）CISCISM岁差参数2）瞬时平天球坐标系CISM瞬时真天球坐标系CIST3）综合上两步GISGIST章动参数2.地球地固坐标系：点在地球坐标系中的主要坐标形式有：空间直角坐标与大地坐标两种形式以T点的坐标为例，两种坐标形式为：两种地球等效地球坐标形式的转换关系为：由于地轴指向随时间变化，其变化表现为地轴北极的极移，使地球极随时间变化形成瞬时极、平极（协议极），相应的赤道又分为平赤道面、协议赤道面；极移又称S.C.Chandler运动；地球坐标系又随时间分为CTS、CTST两种时空形式1)瞬时真地球坐标（CTST）---协议地球坐标(CTS)转换旋转矩阵—极移矩阵3.瞬时真天球坐标---瞬时真地球坐标（天地转换）旋转矩阵，其中GAST为格林尼治视恒星时角4.天地坐标转换过程CIS-CISM-CIST-CTST-CTS上述转换过程可用于GNSS卫星（飞行中的已知点）坐标的转换如GPS卫星转换：由卫星星历提供Kepler根数计算卫星的轨道坐标--协议天球坐标系中的坐标---瞬时平天球坐标系---瞬时真天球坐标---瞬时地球坐标---协议地球坐标（美国的WGS84）第二节GNSS系统的大地坐标系及我国大地坐标系一、GPS坐标系-

WGS-84WGS-84定义：原点位于地球质心，Z轴指向BIH1984.0定义的协议地球极(CTP-ConventionalTerrestrialSystem)方向，X轴指向BIHl984.0的零子午面和CIP-Conventionalinternationalpole赤道的交点，Y轴与Z-O-X平面构成右手坐标系。WGS-84椭球参数:长半轴a＝6378137土±2m；地心引力常数(含大气层)GM＝(3986005l±0.6)l08（m3s-2）正常化二阶带谐系数484．16685*10-6±1.3×l0-9地球自转角速度

ω＝7292115×10-11土0.1500×10-11(rad·s-1)利用以上4个基本常数，可以计算出其它的椭球常数，如第一、第二偏心率e2、e'2和扁率α二、GLONASS坐标系-PZ-90定义原点位于地球质心，Z轴指向IERS推荐的协议地球极(1900～1905平均极)方向，X轴指向地球赤道与BIH定义的零子午面交点，Y轴与Z-O-X平面构成右手坐标系。PZ-90坐标基准参数地球平均半径Re=6378136m地球自转角速度ω

=0.00007292115rad/s地球引力常数GMe=398600.4km3/s-2地球重力场系数-0.00108263

三、BDS坐标系-CGCS2000定义原点位于包括和大气的整个地球质心，Z轴指向初始为(IRP)

BHI1984.0方向，经推算指向至参考历元为2000.0地球参考极，X轴指向地球赤道与IERS定义的零子午面(IRM)垂直于Z轴的方向，Y轴与Z-O-X平面构成右手坐标系；2008年7月1日启用。两者在空间三维坐标状态下之差为cm级2-1WGS-84/GPS与PZ-90/

GLONASS坐标系转换模型：CGCS2000坐标基准参数地球平均半径Re=6378137m地球自转角速度ω

=7.292115×10-5

rad/s-1地球引力常数GMe=3.986004418×104km3/s-2地球动力形状因子J2=-0.001082629832258正常椭球与参考椭球一致CGCS2000由三个层次的站网与坐标实现与维持第1层次：A级GNSS卫星定位连续运行参考网，用于全球地球动力学、地壳形变测量、卫星精密定轨测量，用GPS、GLONASS、VLBI、SLR等技术手段测量。第2层次：B级GNSS大地网，国家二等大地控制网，用于建立地方或城市坐标基准框架，区域地球动力学和地壳形变测量、精密工程测量等2000椭球几何参数a＝6378137mf=1/298.257222101

第3层次：GNSSC、D、D、E级网和天文大地网，用于建立三等及以下控制网测量和工程测量四、WGS-84/PZ-90/

CGCS-2000坐标系特点都属于CTS坐标，三维坐标；都至少有（XYZ）（BLH）两种坐标形式和投影坐标（x,y）形式都源于ITRS(或其实现ITRF)都是本系统单点定位或称绝对定位的坐标，以及相对定位中解算的基线向量的本坐标系同属地球地心、地固坐标系三种坐标系在三维层面的坐标之差为cm级，三维层面可有全球统一的转换参数，如PZ-90/WGS-84之转换关系式（2-1）式五、我国国家大地坐标系（二维、参心坐标系）1.1954北京坐标系坐标基准：克拉索夫斯基椭球，仅几何参数，无明确物理参数，正常重力计算采用赫尔黙特公式；几何参数大地原点：前苏联普尔科沃镇2.1980西安坐标系，GDZ80坐标基准：IUGG1975国际椭球大地原点：陕西省泾阳县永乐镇1980西安坐标系（参心坐标系）还进行了如下定义：椭球短轴Z轴平行于由地球地心指向1968.0地极原点(JYD)的方向；大地起始子午面平行于格林尼治平均天文台子午面，X轴在大地起始子午面与Z轴垂直指向经度零方向；y轴与Z，X轴成右手坐标系。从而确保了GDZ80坐标系椭球定向、定位的明确性，克服了1954北京坐标缺点。IUGG1975国际椭球赋予了明确物理意义，给出了相应椭球参数：可导出如下几何参数3.国家大地原点大地原点：国家水平大地控制网点大地坐标的起算点。大地起算数据：体现在大地原点上，共有四个个，包括大地原点的大地坐标值L0、B0、H0，以及大地原点至某一方向的大地方位角A0，这些数据构成了经典大地测量的基准。陕西泾阳县永乐镇大地原点外部内景第三节坐标系统之间的转换一、转换的主要形式（复习《大地测量基础》内容）1.同一坐标系统内不同坐标形式的变换1）大地坐标与空间直角坐标间等价变换2）球面（大地）坐标与平面坐标间投影变换3）空间直角坐标与极坐标间等价变换有关模型在《大地测量基础》和本章1，2节已述及2.不同坐标系统间的转换统一规定：下标A—目标坐标系；下标B—源坐标系1）不同空间直角坐标间坐标转换平移矩阵尺度系数旋转矩阵两空间直角坐标系坐标转换模型中有7个转换参数。其中平移参数3个、旋转参数3个、尺度参数1个。除上述布尔萨模型外，还有莫洛金斯基、武测等模型2）不同球面（大地）坐标间坐标转换3）不同平面坐标间坐标转换4）GNSS测量中，不同坐标间转换主要是下列两种另有三维七参、二维七参、二维四参等多种转换模型；参见《大地测量学基础》和其它专业书籍二、不同坐标间转换基本步骤（1）模型选择大范围坐标转换选择二维七参数转换模型；小区域坐标转换可选择三维四参数模型或平面四参数模型。对于相对独立的平面坐标系统与GNSS坐标系的联系可采用平面四参数模型或多项式回归模型。（2）重合点选取布GNSS网时应初步选用具有目标坐标系坐标的点做重合点，这些点的GNSS坐标（源坐标系）一经测定就有两套坐标可用。最终选定重合点时还要根据所确定的转换参数，计算重合点坐标残差，根据其残差值的大小来最后确定，若残差大于3倍中误差则剔除，重新计算坐标转换参数，直到满足精度要求为止；用于求二维平面四参数转换的最终重合点数量不少于2个，七参数转换的不少于3个。（3）模型参数计算以转换模型为基础方程，将重合点的两套坐标数据代入模型求转换参数，有多余重合点时应应用最小二乘原理求转换参数的最或然解。

（4）精度评估与检核选择部分重合点作为外部检核点，不参与转换参数计算，用转换参数计算这些点的转换坐标与已知坐标进行比较进行外部检核。应选定至少6个均匀分布的重合点对坐标转换精度进行检核三、附标准转换模型（供实用）（1）二维七参数转换模型

同一点位在两个坐标系下的纬度差经度差，单位为弧度椭球长半轴差（单位米）、扁率差（无量纲），平移参数，单位为米，旋转参数，单位为弧度，尺度参数（无量纲）（2）平面四参数转换模型

其中，x0，y0为平移参数，α为旋转参数，m为尺度参数。x2，y2为2000国家大地坐标系下的平面直角坐标，x1，y1为原坐标系下平面直角坐标。坐标单位为米。（3）Bursa七参数坐标转换模型（4）三维七参数坐标转换模型

（5）坐标转换精度评定和评估方法

依据计算坐标转换模型参数的重合点（没有参与转换参数计算的重合点）的残差中误差评估坐标转换精度。对于n个点，坐标转换精度估计公式如下：

①V（残差）=重合点转换坐标-重合点已知坐标②空间直角坐标X残差中误差③空间直角坐标Y残差中误差④空间直角坐标Z残差中误差点位中误差⑤平面坐标x残差中误差⑥平面坐标y残差中误差⑦大地高H残差中误差

平面点位中误差为（6）决定不同坐标间转换精度的因素重合点在网中的分布状况转换模型精度与适用性重合点的坐标分量及点位精度三、GNSS测量中涉及的其它坐标1.站心坐标系概述有站心赤道直角坐标系、站心地平直角坐标系和站心地平极坐标系三种坐标系形式。GNSS坐标计算中，常用的为站心赤道直角坐标系与站心地平直角坐标系：站心地平、站心赤道坐标和地球空间直角坐标关系图2.站心赤道直角坐标系定义：以测站P1为原点，且与地球空间直角坐标系，O-XYZ相应坐标轴平行的坐标系站心赤道坐标系与地球空间直角坐标间关系为简单的平移关系3.站心地平直角坐标系P1-xyz定义：以P1为原点，过测站P1的法线为z轴，指向天顶为正，以子午线方向为x轴，向北为正，以与x、z垂直方向为y轴，东为正。地平直角坐标系与站心赤道直角坐标系关系可通过旋转变换实现，关系式如下：

4.GNSS空间直角坐标与地平直角坐标的关系5.站心地平极坐标系P1—rAh（用于描述测站与卫星的相对位置）定义：以测站P1为原点，坐标参数为：

r--测站P1至卫星s的距离A--卫星的方位角，zox平面与zos平面的夹角，左旋为正h--卫星的高度角os与xoy平面的夹角。1）站心地平极坐标系与站心地平直角坐标系关系：A=arctan(y/x)h=arctan(z/(x2+y2)1/2)

四、GNSS卫星（飞行中的已知点）坐标计算流程GNSS定位坐标计算流程地面监测站监测卫星运动状态主控站处理注入站卫星用户获得导航电文导航电文计算卫星在协议天球坐标系中的坐标计算卫星在瞬时平天球坐标系中的坐标计算卫星在瞬时地球坐标系中的坐标通过星地距离方程解算基线将接收机坐标转换为北京54或西安80坐标计算卫星在协议地球坐标系CTS中的坐标计算卫星在轨道平面坐标系中的坐标计算卫星在瞬时真天球坐标系中的坐标卫星参考历元的Kepler参数+摄动改正参数=实际观测历元的Kepler参数计算卫星在协议天球坐标系中的坐标第四节坐标系统间的转换一、GNSS定位中时间的意义GNSS定位过程是在严格的时间定义下进行的，无论是坐标系的定义、GNSS卫星信号接收、星地距离测量，还是坐标转换都需在严格的时间定义下完成二、时间系统各类与构成时间系统类型：基于各种用途，时间系统有多种：ST恒星时、MT平太阳时、UT世界时、ATI原子时、UTC协调世界时。时间系统构成：由时间原点(历元)，时间尺度(时间单位)构成。1.恒星时ST(SiderealTime)原点：以春分点为参考点，由春分点的周日视运动所定义的时间系统为恒星时系统。ST尺度：春分点连续两次经过本地子午圈的时间间隔为一恒星日，一恒星日分为24个恒星时2.平太阳时MT(MeanSOlarTime)平太阳：以真太阳周年运动的平均速度在天球赤道上作周年视运动的太阳，其周期与真太阳一致。时间原点：以平太阳为参考点，由平太阳的周日视运动所定义的时间系统为平太阳时系统。时间尺度为：平太阳连续两次经过本地子午圈的时间间隔为一平太阳日，一平太阳日分为24平太阳时。

3.世界时UT(UniversalTime)时间原点：以平子夜为零时起算的格林尼治平太阳时定义为世界时UT。时间尺度：与平太阳时的尺度相同4.原子时ATI（InternationalAtomicTIME）时间原点：原子时的起点，按国际协定取为1958年1月1日0时0秒(UT2)(事后发现在这一瞬间AT1与UT2相差0.0039秒)。时间尺度：原子时的秒长被定义为铯原子C813基态的两个超精细能级间跃迁辐射振荡9192631170周所持续的时间。5.协调世界时UTC(CoordinateuniversalTime)时间原点：同UT相同时间尺度：与ATI相同考虑到原子时比世界时每年快约1秒，两者之差逐年积累，便采用跳秒(闰秒)的方法使协调时与世界时的时刻相接近，其差不超过1秒。跳秒一般在每年的6