GPS坐标系统和时间系统

1、第2章坐标系统和时间系统2.1概述2.1.1坐标系统的重要性坐标系统是一切测量工作的位置基准，是国民经济和国防建设的基础，贯穿于测 绘科学 与技术的各分支之中。GPS定位技术是通过安置于地球表面的 GPS接收机 接收GPS卫星信号 以测定地面点位置的技术。观测站固定在地球表面，其空间位置随同地球的自转而运动；GPS卫星围绕地球质心旋转，与地球自转无关，需要建立 卫星在其轨道上运动的坐标系。因此，GPS定位必须寻求两种坐标系之间的关系，以 实现坐标系之间的转换。2.1.2时间系统的重要性时间作为精确描述天体和人造卫星运行位置及其相互关系的重要基准，是GPS系统中的重要量测值，其精度将直接影响GP

2、S定位的精度。(1卫星位置：若欲使卫星位置误差小于1cm ,相应的时刻误差应小于2.6 10-6s ; (2星站距离:若欲使距离误差小于1cm，信号传播时间的测定误差应不超过 3 X0- 11s ; (3地球上点的位置：若欲使赤道上一点的位置误差小于1cm ,时间的测定误差应不超过2 X0-5s。(4 GPS授时测量对时间的测定精度要求更高。2.1.3坐标系统和时间系统的关系可以看出GPS定位的坐标系统和时间系统密不可分，且两者的融合构成四维大地测量学的基准。2.2坐标系统在GPS定位测量中,采用天球坐标系和地球坐标系两类坐标系：地球坐标系随同地球自转，可看作固定在地球上的坐标系,便于描述地面

3、观测站的空间位置；天球坐标系是一种惯性坐标系，与地 球自转无关，便于描述人 造地球卫星的位置。2.2.1协议天球坐标系一、天球及其相关概念(见图2-1，以研(1天球:以地球质心为中心，半径无穷大的理想球体。可作为天体的投影面究其位置、运动规律及其相互关系(2天轴和天极：前者指地球自转轴的延伸直线；天轴和天球表面的交点称为天极 P , P N和PS分别表示北天极和南天极。(3天球赤道面和天球赤道：前者指通过地球质心并与天轴垂直的平面，其与天球 表面的交线为天球赤道。(4黄道面、黄道和黄极：地球绕太阳公转时的轨道平面为黄道面；其与天球表明 相交的大圆为黄道，即我们看到的太阳运行轨迹；过天球中心并垂

4、直于黄道平面的 直线和天球表面的交点为黄极，分黄北极和黄南极。(5春分点:太阳在黄道上从天球南半球向北半球运行时，黄道与天球赤道的交点 Y (6黄赤交角:黄道平面和天球赤道面的夹角 &(7天球子午面与子午圈、时圈：包含天轴并通过天球面上任意一点的平面称为 天球子午面,天球子午面和天球表面相交的大圆称为天球子午圈，通过天轴的平面和天球表面相 交的半个大圆称时圈。二、天球坐标系原点位置、坐标轴指向和长度单位构成坐标系统 的三要素，因此，应从上述三个 方面描述坐标系统。常见的天球坐标系包括天球空间直角坐标系和天球球 面坐标 系两种形式。（1天球空间直角坐标系：以地心M为坐标原点，其Z轴指向北

5、天极，X轴指向春分点，丫轴垂直于XOZ轴并构成右手坐标系。该坐 标系中天体的位置为（x , y ,z。以地心M为坐标原点，赤经a为含天轴和春分点的天球子午面与过天体S的天球子午面之间的夹角，赤纬S为原点M至天体s的连线与天球赤道面之间的夹角；向径r为原点M至天体 s的距离，各轴正向如图2-2所示。岁差与章动日月对地球的引力产生力距，从而使地球自转轴的方向在惯性空间缓慢地移 动。可以将运动分解为一个长周期变化和一系列短周期变化的叠加。地球自转轴的长周期变化约25800N 2-3岁望和盘动示意图年绕黄极一周。使春分点产生每年约 50.26 '的长期变化,称之为日月岁差。一 系列短日、月周期

6、变化中幅值最大的约为 9.2 周期为18.6年，这些短周期变化统 称为章动。岁差使北天极绕黄北极以顺时针方向缓慢旋转，构成如图2-3所示的以 黄赤交角&为的小圆。这种有规律运动的北极称为平北天极，相应的天球赤道和春分点称为天球平赤 道和平春分点；若把观测时的北天极称为瞬时北天极（简称真北天极，相应瞬时天球 赤道和瞬时春分点（或称真天球赤道和真春分点，章动将使真北天极绕平北天极产 生旋转，轨迹见图2-3。四、协议天球坐标系（一协议天球坐标系在岁差和章动的影响下，瞬时天球坐标系的坐标轴指向是不断变化的。在这样的坐标系中不能直接使用牛顿第二定律，这对研究卫星的运动很不方便。因此需要 建立一个

7、三轴指向不变的天球坐标系，以便在这个坐标系内研究人造卫星的运动（计 算卫星的位置。而在这 个坐标系中所得到的卫星位置又可以方便地变换为瞬时天 球坐标系中的值，以便与地球坐标 系进行坐标变换。为此，选择某一个历元时刻，以 此瞬间的地球自转轴和春分点方向分别扣除此瞬间的章动值作为z轴和x轴指向，y轴按构成右手坐标系取向，坐标系原点与真天球坐标系相同。这样的坐标系称为 该历元时刻的平天球坐标系，也称协议天球坐标系。国际大地测量学协会（IAG和国际天文学联合会（IAU决定，从1984年1月1日后启用的协议天球坐标系，其坐标轴的指向是以2000年1月15日太阳质心力学时（TDB为标准历元（记为J2000

8、.0的赤道和春分点所定义的。（二天球坐标系之间的 坐标转换为了将协议天球坐标系的卫星坐标转换到观测历元t的瞬时天球坐标，可分两步进行:首先将协议天球坐标系中的坐标换算到瞬时平天球坐标系；然后将瞬时平天球坐标系的坐标转换到瞬时天球坐标系统。表2-1为三种天球坐标系统的定义与 缩写。人球坐标系原点Z轴X轴坐标系缩坷协议人球堰标系地卍标准历兀平天战标准历元平祥分点CIS平人球坐标系地心瞬时平天极瞬时平春分点瞬时天球坐标系地心瞬时真天极瞬时良春分点t（1将协议天球坐标系转换为瞬时平天球坐标系（岁差旋转（2-1公式中z、 a ?分别表示与岁差有关的三个旋转角，公式中的3 B矩阵为 岁差旋转矩阵。（2将瞬

9、时平天球坐标系转换为瞬时天球坐标系（章动旋转（2-2公式（2-2中&£? ?分别表示黄赤交角、交角章动和黄经章动。根据以上两个公式，可将协议天球坐标系转换为瞬时天球坐标系的公式简化为CISM 1000cos sin 0sin cos cos 0sin 010sin Ocos 1OOOcos sin 0sin cos ? ? ? ?cos sin 0sin cos 0001cos sin 0sin cos cos( sin(0sin( cos(00?_? ?+?+ _?+?+= ?x z y x t ££££££

10、63;££££CISzyz xzx t z y x R R z y x ?=? (2-32.2.2协议地球坐标系、地球坐标系为了表达地面观测站的位置，需采用固 联在地球上、随同地球自转的地球坐标 系。地球空间直角坐标系以地球质心为坐标原点，以地球自转轴作为Z轴的正向，与天球坐标系不同的是以地球赤道面与格林威治子午面交线的方向作为X轴的正 向。大地坐标 系定义为:地球椭圆中心与地球质心重合，椭球短轴与地球自转轴重 合,大地纬度B为过地面点的拖球法线与拖球赤道面的夹角，大地经度L为过地面 点的椭球子午面与格林 威治子午面之间的夹角，大地高H为地面点 沿椭球法

11、线至 椭球面的距离。二、协议地球坐标系(1极移:由于受到地球内部质量不均匀影响，地球自转轴相对于地球体产生运动 导致地极点在地球表面的位置随时间而变化，这种现象称为地极移动，简称极移。 为了定量 描述极移，可构造一平面直角坐标系，取平地极为原点，xp轴指向格林尼治 平子午圈,即指向经度为0的方向,yp轴指向经度为270 的方向。(2协议地球坐标系:1900年国际大地测量与地球物理联合会以 1900.00至 1905.05年地球自转轴瞬时位置的平均位置作为地球的固定极称为国际协议原点 CIO,以此作为协议 地极CTP。以协议地极为基准点的坐标系称为协议地球坐标 系;与瞬时极对应的地球坐标 系，则

12、称为瞬时地球坐标系。图极移分量和极移旋转(3地球瞬时坐标系与协议坐标系的转换：如图2-5所示,有瞬时地球坐标系到协议地球坐标系的转换可通过绕 Xt轴顺时针转动极移分量 Y P和绕Y t轴顺时针转动极移分量X P实现，转换模型见公式(2-4P P P PPPZYXY Y XZY?_ 均？? 101001cos sin 0sin cos 001cos 0sin 010sin 0cos 0（2-4三、天球坐标系与地球坐标系的转换根据协议天球坐标系（x , y , z CIS和协议地球坐标系（X , 丫 , Z CIS的定义，二者 坐标原点和纵轴指向均相同，x和X轴间夹角为春分点的格林威治恒星时（记为

13、 GAST ,则瞬时天球坐标系（x , y , z t转换为瞬时地球坐标系（X，丫 , Z t的公式可表 示为，tt z t z y x GAST GAST GAST GAST z y x GAST R Z 丫 X ? ?-=?=? 1000cos( sin(0 sin(cos( ( (2-5若记公式(2-4中地球瞬时坐标系与协议坐标系的转换矩 阵为M ,即(P X P Y Y R X R M -=,则瞬时天球坐标系(x , y , z t转换为协议地球 坐标系(X，丫 ,Z CIS的公式为tz CISzyxGAST MRZY X?=? ( (2-6则可根据协议天球坐标系转换为瞬时天球坐标系的

14、公式(2-3得到协议天球坐标系与协议地球坐标系之间的转换公式（2-7CISzyzxzx z CISz y x R R GAST MR ZY X?=? （ （2-7可以对协议天球坐标系（x , y , z CIS和协议地球坐标系（X，丫 , Z CIS之间的转换步骤进行归纳总结，见图2-6223国家坐标系与地方坐标系GPS卫星定位测量是用三维地心坐标为依据测定和表示点的空间位置，为充分利用已有测绘成果，需要将GPS测量成果纳入国家坐标系或地方独立坐标系；因此， GPS定位测量数据处理中，应考虑将GPS测量成果由世界地心坐标系转换至国家或 地方独立坐标系。一、地心坐标系和参心坐标系(1旋转椭球:与

15、大地水准面接近的规则的具有微小扁率的数学曲面，现代大地测 量中常以椭球的长半径a、地球重力场二阶带谐系数 J 2、地球引力常数与地球质 量的乘积GM和地球自转角速度宀4个参数描述其几何物理特性。(2地心坐标系：将椭球中心与地球质心重合，且与全球大地水准面最为密合的旋 转椭球。(3参心坐标系：为了研究局部球面的形状，且使地面测量数据归算至椭球的各项 改正数最小，各个国家和地区分别选择和某一局部区域的大地水准面最为密合的椭 球建立坐标系。这样选定和建立的椭球称为参考椭球，对应的坐标系称为参心坐标 系。二、国家坐标系(1 1954年北京坐标系：采用前苏联的克拉索夫斯基椭球体，由前苏联西伯利亚地 区的

16、一等锁，经我国的东北地区的呼玛、吉拉林、东宁三个基准网传算；基于1954年北京坐 标系的我国天文大地网未进行整体平差；高程异常是以前苏联1955年大地 水准面重新平差 的结果为起算值，按我国天文水准路线推算出来的，而高程又是以 1956年青岛验潮站的黄海平均海水面为基准。缺点:椭球参数与现代精确的椭球参数的差异较大，不包含表示地球物理特性的 参数,给理论和实际应用带来了许多的不便；椭球定向不十分明确，既不是指向CIO 极,也不是指向我国目前使用的JYD极;采用局部分区平差，参考椭球面与我国大地 水准面呈西高东低的系统性倾斜，东部高程异常最大达67米。椭球参数为(2 80西安国家坐标系：1978

17、年决定对我国天文大地网进行整体平差 ,重新选定 椭球,并进行椭球的定位、定向。椭球的短轴由地球质心指向 1968.0 JYD,起始子午 面平行于格 林尼治平均天文子午面，椭球面与大地水准面在 我国境内符合最好，高 程系统采用1956年黄海平均海水面为高程起算基准。椭球参数为其特点是:采用多点定位原理建立，理论严密,定义明确；椭球参数为现代精确的地球总椭球参数;椭球面与我国大地水准面吻合得较好 椭球短半轴指向明确；经过了整体平差，点位精度高。三、站心坐标系 如果测量工 作以测站为原点，则所构成的坐标系称为测站中心坐标系(简称站心坐标系，分为站 心地平直角坐标系和站心极坐标系。(1站心地平直角坐标

18、系:Z轴与过测站的椭球法线重合，X轴垂直于Z轴(与过 测站的大地子午线相切并指向椭球的短轴，丫轴与X轴和Z轴一起构成左手坐标 系，如图2-7所示。(2站心极坐标系：以测站的铅垂线为准，以测站点到某点j的空间距离D、天 顶距Z14325163781403.986005101.08263107.29211510a mG M m s J rad s-=? =? =? 63782451/298.3a m f =天和大地方位角A表示j点的位置。四、地方独立坐标系基于方便实用和限制变形目的，选取过测区中心的经线或某个起算点的经线作 为独立的中央子午线,以某个特定方便使用的点和方位为起算原点和方位，并选取当

19、 地平均高程面H M为投影面。其参考椭球半径a 1为?+=?+=0111 三 m H a a a a （28E（为该地区的平均高程异常。五、高斯平面直角坐标系和 UTM（墨卡托坐标系（略224 WGS-84坐标系、WGS-84坐标系WGS （世界大地坐标系统 属于协议地球坐标系，该系统以地球质心为原点,Z轴 指向国际时间局BIH1984.0定义的协议地极CTP，轴指向BIH1984.0定义的零子午 面与CTP相应的赤道的交点，丫轴构成右手坐标 系。WGS-84坐标系采用的地球 椭球称为WGS-84椭球,其常数为国际大地测量学与图2-S WGS-S4坐标地球物理学联合会（IUGG第17界大会的推

20、荐值：长半轴：a=6378137m±m地球引力常数:GM=3986005 108 ±.6 108 （m 3s -2正常化二阶带谐系数:C2.0=- 484.16685 10-6 ±30 10-6 地球自转角速度:w=7292115 10-11 0±5 10-11 rads-1二、椭球扁率：f84=1/298.257223563三、GPS定位测量中的坐标转换我们知道,GPS定位结果属于协议地球地心坐标系,即WGS-84,且通常以空间 直角坐标（X , Y , Z或椭球大地坐标（B , L , H的形式给出。而我们所需要的成果一 般都是 北京54（BJZ54

21、、西安80坐标系（GDZ80或地方独立坐标系，为此，必须要实现GPS坐标系与实用坐标系间的转换。GPS定位测量中的坐标转换包括以下两种 情形，-(X , Y , Z (1同一坐标系内空间直角坐标与大地坐标的换算,即(B丄，H图2-8 WGS-S4坐图2弋222( cos cos (cos sin (1 sin sin X N H B L Y N H B L Z N e H BaN H Bb =+=+? ? =-+? ?=?+? (2arcta n arcta n 1si n Y L X B Z H N e B?=?-a N =222222a b e f f=-a b f a(2不同的空间直角坐

22、标系之间的坐标换算常用七参数法，一般使用布尔莎模型和莫洛金斯基模型。1布尔莎-沃尔夫(Bursa-Wolf 模型在该模型中采用了 7个参数，分别是3个平移参数、3个旋转参数(3个欧拉角 和1个尺度差参数。该模型的基本转换可分解为平移变换、缩放变换和旋转变换三 个过程，I .从A X正向看向原点 A O ,以A O点为固定旋转点，将A O -A A A X 丫 Z 绕A X轴逆时针旋转A B角度，使经过旋转后的A 丫轴与B O -B B X 丫平面平行;n .从A 丫正向看向原点A O ,以A O点为固定旋转点，将A O -A A A X 丫 Z 绕A 丫轴逆时针旋转a平移转换b缩放转换c旋转转

23、换图2-9 布尔莎-沃尔人(Bursa-Wolf)模型A B丫角度，使经过旋转后的AX轴与B O -B B X 丫平面平行，显然,此时A Z轴也与B Z平行;川.从A Z正向看向原点 A O ,以A O点为固定旋转点，将A O -A A A X 丫 Z 绕A Z轴逆时针旋转Z 3角度，使经过旋转后的A X轴也与B X平行，显然，此时A O -A A A X 丫 Z 的三个坐标 轴已与B O -B B B X 丫 Z中相应的坐标轴平行；IV .将A O -A A A X 丫 Z中的长度单位缩放(1+m倍,使其长度单位与B O -B B B X 丫 Z的一致；V .将A O -A A A X 丫

24、Z 的原点分别沿 A X、A 丫和A Z轴移动A BXT -、A B丫 T -和 A BZ T -,使其与B O -B B B X 丫 Z的原点重合。由图2-9可知,任意点P i在两坐标系中的坐标之间有如下关系，(? ? +? ? ?=?AAAX丫 Z Z 丫 X B B B Z 丫 X R R R m T T T Z 丫 X£-S£ (1 (2考虑到两坐标轴定向的差别一般很小，因此欧拉角&X、£ Y £ Z通常都是微小量，有(? ? ? ? ? _=111( (XYX ZY ZXY Z R R R R &££

25、3;££££££2莫洛金斯基模型莫洛金斯基模型认为受尺度和旋转影响的只是任意点与参考点的坐标差，在该模型中也 是采用了 7个参数,分别是3个平移参数、3个旋转参数(也被称为3个欧 拉角和1个尺度参数,不过定义与布尔沙模型有所不同。 转换过程可描述为：将0 A -X A Y A Z A 的原点平移到某 点P，形成一个过渡坐标系P-X'Y'Z'将OP-X'Y'Z' 依次分别绕X、Y和Z轴旋转s X、coy和3 z三个角度后使其坐标轴与 0 B -X B 丫 B Z B中相应的坐标轴平行，旋转

26、方 式和次序与布尔沙-沃尔夫模型相似；再将P-X'Y'Z'中的长度单位缩放(1+m倍，使其 长度单位与0 B -X B 丫 B Z B的一致;最后将0 A -X A Y A Z A的原点分别沿X、Y和Z轴移动-T X、-T Y 和-T Z ,使其与O B -X B Y B Z B的原点重合。其转换模型为,两种模型的转换结果是等价的，但在实际应用过程中，还是有所差 异：布尔莎 模型在进行全球或 较大范围的基准转换时较为常用，但是,旋转参数与平移参数具 有较高的 相关性。对于小范围可以3参数（3个平移参数；或者是3个平移和1个尺 度参数（4参数;最好的情况 除了上述4个参数

27、外,可确定一个旋 转参数（5参数。2-10罠洛金斯基模型采用莫洛金斯基模型则可以 克服这一问题，因为其旋转中心可以 人为选定， 当网的规模不大时，可以选取网中任意一个点，当网的规模较 大时，则可选取网的重 心,然后以该点作为为固定旋转点进行旋转。 2.3时间系统时间作为精确描述天体和人造卫星运行位置及其相互关系的重要基准，是GPS系统中的重要量测值，其精度将直接影响GPS定位的精度。时间系统和其他测量 基准一样,要定义时间单位（尺度和原点（起始历元，只有把尺度与原点结合起来，才 能给出时刻的概念。理论上,任何一个周期运动，只要它的运动是连续的，其周期是恒 定的，并且是可观测和用 实验复现的，都

28、可以作为时间尺度(单位。实践中，由于所选 用的周期运动现象不同，便产生了不同的时间系统。与 GPS测量有关的时间系统包 括世界时、原子时和力学时。2.3.1世界时系统世界时系统是以地球自转为基准的一种时间系统。由于观察地球自转运动时，所选择的空间参考点不同，世界时系统又包括以下几种不同的形式。一、恒星时(ST1(1 11B PZY A PX B P Z XA P Y BPYX A P Z X X X X T Y Y m Y Y T Z Z Z Z Tcd(-?d?g?g? ?=+-+? ?-?以春分点为参考点，由春分点的周日视运动所定义的时间系统为恒星时系统。其时间尺 度为:春分点连续两次经过

29、本地子午圈的时间间隔为一恒星日 ，一恒星日 分为24个恒星时。恒星时以春分点通过本地上子午圈时刻为起算原点 ，所以恒星 时在数值上等于春分点相对于 本地子午圈的时角。同一瞬间对不同测站的恒星各 异,所以恒星时具有地方性，有时也称为 地方恒星时。恒星时是以地球自转为基础的。由于岁差、章动的影响，地球自转轴在空间的 指向是变化的,春分点在天球上的位置并不固定。 对于同一历元所相应的真天极和 平天极,有真春分点和平春分点之分。因此，相应的恒星时也有真恒星时和平恒星时 之分。由于地球围绕太阳的公转轨道为一椭圆，太阳的视运动速度是不均匀的。 若以 真太阳作为观察地球自转运动的参考点，将不符合建立时间系统

30、的基本要求。为此选取平太阳作为参 考点,其视运动速度为真太阳周年运动的平均速度。其时间尺度为:平太阳连续两次经过本 地子午圈的时间间隔为一平太阳日，一平太阳日分为24 平太阳时，平太阳时在数值上等于 平太阳通过本地上子午圈的时角。同样,平太阳时也具有地方性，故常称其为地方平太阳时 或地方时。平太阳日由平正午开始，即平正午为0时,平子夜为12时。1925年国际天文联 合会决定，改平太阳日由平子夜开始，即平子夜为0时，平正午为12时,简称平时或民 用时。以平子夜为零时起算的格林尼治平太阳时定义为世界时UT 0。世界时与平太阳时的尺度 相同，但起算点不同。这是以地球自转这一周期运动作为基础的时间 尺度。由于地球自转的 不稳定性，在UT 0中加入极移改正 ?即得到UT 1。UT1 加上地球自转速度季节性变化 T s后为UT 2,即入？+=01UT UT (210s T UT UT ?+=02 （2-11显然,世界时UT 1经过极移改正后仍含有地球自转速度 变化的影响，而UT 2虽经地球自转季节性变化的影响，但仍含有地球自转速度长期变化和不规则变化的影响 ，所以 世界时UT 2仍不是一个严格均匀的时间系统。世界时 UT 0、UT 1和UT 2之间关 系为s T UT UT ?+?+=入 02 <222.3.2原