ITRS,GCRS和J2000坐标系的相互转化

ITRS/GCRS/J2000坐标系的相互转换李云飞(blitheli@)南京大学天文系/上海航天技术研究院2008-12-6 本文主要阐述了目前国际天文界最新规定的岁差章动模型(IAU 2000A/B)，并根据此模型给出ITRS,GCRS和J2000平赤道地心系相互转换的详细步骤。

本文主要依据《IERS Conventions 2003》(IERS Technical Note No. 32)而来，由于文章为英文，且详细叙述了ITRS,GCRS参考系以及岁差章动等模型，会使得读者初次阅读（或天文背景知识不够）产生阅读障碍。

笔者也是在多次反复阅读的基础上并参考其它书籍和文献才稍微弄清楚，因此本文是笔者关于此文献的一个总结，希望能够对读者有所帮助。

本文不打算对坐标系转换的理论（岁差章动，极移等）进行详细讲解，仅仅给出坐标转换的基本步骤和必要常识，以便读者根据此文能够迅速掌握ITRS,GCRS和J2000平赤道地心系相互转换的基本原理，并能够依据IERS提供的Fortran源程序进行实际编程应用。

0 前言0.1 名词缩写和解释BCRS barycentric celestial reference systemCEO celestial intermediate originCIP celestial intermediate poleCIRS celestial intermediate reference systemparametersorientationearthEOPGCRS geocentric celestial reference systemGMST Greenwich mean sidereal timeGAST Greenwich apparent sidereal timeIAU international astronomical unionICRS international celestial reference systemIERS international Earth rotation and reference systems serviceITRS international terrestrial reference systemJ2000 2000 January 1.5SOFA standards of fundamental astronomyTEO terrestrial intermediate originTIRS terrestrial intermediate reference systemtimeterrestrialTTUT1timeuniversalUTUTC coordinated universal time0.001MAS milliarcsecond () ''0.2 ITRS就是我们常说的地固坐标系，其原点在地球质心（包含大气海洋等质量），坐标系xy平面为地球赤道面，z轴指向北极CIO处，x轴指向格林威治子午线与赤道面交点处。

此坐标系固定在地球上，地面站测控，以及地球引力场系数等都在此坐标系下定义。

0.3 J2000.0 Mean equator, mean equinox dynamical system常被称为J2000平赤道地心坐标系。

其原点也是在地球质心，xy 平面为J2000时刻的地球平赤道面，x 轴指向J2000时刻的平春分点（J2000时刻平赤道面与平黄道面的一个交点）。

此坐标系常被作为地球卫星的惯性坐标系，卫星运动积分等都在此坐标系计算。

0.4 GCRSJ2000地心天球坐标系，其定义与J2000平赤道地心坐标系仅有一个常值偏差矩阵B 。

目前IAU推荐用此坐标系逐渐取代J2000平赤道地心坐标系。

0.5 一些时间量1(132.184UT UTC UT UTC DAT TAI UTCTAI UTC DATTT TAI s)=+−=−=+=+ (1)给定UTC 时刻，可求得其距离J2000的世纪数，即 t = (TT - 2000 January 1d 12h TT) in days/36525.上式中2000 January 1d 12h TT 对应的儒略日(Julian day)为2451545.0天。

0.6 坐标转换当一坐标系绕其3个轴旋转θ角时，则坐标旋转矩阵可表述为：100()0cos sin 0sin cos cos 0sin ()010sin 0cos cos sin 0()sin cos 000x y z R R R 1θθθθθθθθθθθθθθθ⎛⎞⎜⎟=⎜⎟⎜⎟−⎝⎠−⎛⎞⎜⎟=⎜⎟⎜⎟⎝⎠⎛⎞⎜⎟=−⎜⎟⎜⎟⎝⎠其旋转方向符合右手螺旋法则，即逆时针旋转为正方向。

另外坐标旋转矩阵具备如下性质：1()()()T R R R θθθ−==−1 背景知识我们知道，地球的自转轴在惯性空间中不是固定的，而是不断摆动的。

此摆动造成地轴绕北黄级顺时针运动，夹角约为23.5度。

于此同时，地轴还在做微小的抖动，见下图。

前者的运动称为岁差(Precession)，后者运动成为章动(Nutation)。

岁差章动的原因主要有两个方面。

其一是太阳系行星对地球绕日轨道所产生的摄动影响；其二是太阳和月球对地球赤道隆起部分的摄动影响。

关于岁差章动的计算，此前一直采用IAU 1976岁差模型和1980章动模型。

随着时间的推移，此模型的精度逐渐跟不上需要。

因此，IAU 又规定，自2003年1月1日起，采用新的岁差章动模型，即IAU 2000A 模型（精度到达0.2mas ），或者IAU 2000B 模型（精度达到1mas ）。

岁差章动模型描述了地球自转轴的运动，另外，由于地球表面的海洋，大气运动以及地核内部液体的运动造成地球自转轴并不是相对地球不动的；相对地球北极CIO 点来说有个小范围的运动，此种现象成为极移。

此外，地球的自转也不是均匀的，也很复杂。

综上所述，从地固坐标系(ITRS)到地心惯性系(GCRS 或J2000平赤道地心系)的坐标转换矩阵由极移，自转和岁差章动组成。

图表 1 地轴的岁差章动现象2 ITRS 到GCRS 的转换矩阵本文采用IAU 2000A/B 岁差章动模型，在某历元UTC 时刻，ITRS 到GCRS 的转换矩阵可写成：(2)()()()()T GCRS ITRS ITRS r Q t R t W t r HG t r =⋅⋅⋅=⋅K K K 其中，和分别对应同一向量在ITRS 和GCRS 坐标系中的坐标。

ITRS r K GCRS r K上式中，，()W t ()R t 和分别对应极移，自转和岁差章动转换矩阵。

()Q t 在计算()R t 和的时候，会有两种计算方法，我们分别称为CEO-based 转换方法和equinox-based 转换方法。

其中前者为IAU 提出的新的计算方法。

()Q t 下面分别给出上述三个转换矩阵的求解过程。

2.1 极移矩阵()W t()(')()()z y p x W t R s R x R y p =−⋅⋅ (3)上式中，为： 's'-0.047 s mas t =⋅极移量(,)p p x y 的求解为：(,)(,)(,)(,)p p IERS tidal nutation x y x y x y x y =+ΔΔ+ΔΔ极移量主要是由IERS 根据天文观测给出的，每周都有新的观测数据，此外，由于地球潮汐和章动的影响，会对极移有微小的修正(,)IERS x y (,)tidal x y ΔΔ和(,)nutation x y ΔΔ。

上式中，(,由IERS 给出的观测数据计算求得，(,)IERS x y )(,)tidal nutation x y x y ΔΔ+ΔΔ可由公式计算得到，IERS 提供此fortran 源程序。

2.2 地球自转矩阵()R t()()z R t R θ=− (4)地球自转角θ的求解根据转换方法的不同有不同的求解方式（CEO-based 或者equinox-based ），具体求解IERS 给出了fortran 源程序。

2.3 岁差章动矩阵()Q t 前面提到，计算此矩阵有两种方法：1.) CEO-based 方法：这种方法是IAU 最新提出的并极力倡导的，其计算公式为： 2222221()1()1()121)z aX aXY X Q t aXY aY Y R s X Y a X Y a X Y ⎛⎞−−⎜⎟=−−⋅⎜⎜−−−+⎝⎠=++其中：⎟⎟ (5) 上式中：2000(,)(,)(,)IAU IERS X Y X Y dX dY =+ (6)2000(,)IAU X Y 和可根据IAU 2000A/B 岁差章动模型求解出，IERS 同样给出求解的fortran 源程序，另外，由于IAU2000A/B 岁差章动模型没有包含地轴的高频率运动，所以要加上IERS 通过观测数据给出的高频率修正项。

s (,)IERS dX dY 2.) Equinox-based 方法：()()()Q t B P t N t =⋅⋅ (7)其中，常值偏差矩阵B ，岁差矩阵和章动矩阵如下：()P t ()N t 0000()()()()()()()()()()()()z y x x z A x A z A x z x B R R R P t R R R R N t R R R δαξηεψωχεψεε=−⋅−⋅⎧⎪=−⋅⋅⋅−⎨⎪=−⋅Δ⋅+Δ⎩(8)章动量(,)ψεΔΔ为：2000(,)(,)(,)IAU IERS ψεψεδψδεΔΔ=ΔΔ+ΔΔ (9)上式中，2000(,)IAU ψεΔΔ(,dX )IERS 由IERS 2000A 章动模型给出。

前面提到过，IAU 2000A/B 模型提供的岁差章动不包含高频率项，而是由IERS 的观测数据提供（上式右端最后一项），但是在IERS 给出的观测数据中仅仅给出，我们可以通过IERS 提供的fortran 源程序将转换为)IERS dY (,)IERS dX dY (,δψδΔΔε。

其余参数皆为岁差参数，可以通过公式求出，此处从略。

值得一提的是常值偏差矩阵B 中的参数也是给定的，在CEO-based 方法求解中，此偏差是包含在中的。

2000(,)IAU X Y 3 利用IERS 提供的Fortran 源程序进行转换上节中详细讲述了ITRS 到GCRS 转换矩阵的求解过程，在实际应用中，如果是自己编写源程序的话是件非常琐碎的事情，因为IERS 2000A/B 章动模型的参数多达1000多项。

幸而这些基本子程序IERS 都提供了，我们所做的就是如何正确的运用这些源程序，并将它们组合起来。