J2000与WGS84坐标及转换

摘要

       在遥感卫星、航天器轨道计算、导航等领域，WGS84 和 J2000（J2000.0 惯性坐标系） 是两个最常用的坐标系。它们分别属于 地固坐标系（Earth-Fixed） 和 惯性坐标系（Inertial），适用于不同的应用场景。这份笔记详细介绍了两个坐标系的应用场景和相互转换方法。

一、 J2000与WGS84坐标系可视化
       

        WGS84（地固坐标系）随地球转动，J2000（J2000.0 惯性坐标系）在空间内固定，如下图所示。

图1  WGS84与J2000坐标系

        WGS84 是 地球固定坐标系（ECEF），随地球自转，适用于 地面或近地空间目标 的定位和导航。应用与卫星对地观测、全球定位系统（GPS/GNSS）和导弹/无人机导航等等。

        J2000 是 惯性坐标系，不随地球自转，适用于 航天器轨道计算、深空探测和天体测量。应用于卫星轨道计算、卫星姿态控制和天体测量与航天动力学等等。

二、 J2000与WGS84坐标系转换
     

  坐标转换过程如下图所示。

图2  J2000坐标系转WGS84坐标系过程

       从上图可知，要从惯性坐标系J2000转成地固坐标系WGS84需要四步。详细介绍如下所示。

2.1 岁差矩阵
   

   岁差（Precession）是地球自转轴方向在空间中的长期缓慢变化，类似于旋转陀螺的轴在重力作用下的摆动。这种运动主要由太阳和月球的引力对地球赤道隆起部分的扭矩作用引起。可视化如下所示(夸大了一点变率，看着明显)。

图3  地球岁差可视化

岁差矩阵的构造
设初始历元为J2000.0（儒略日2451545.0），目标历元为时刻。岁差矩阵通过三个欧拉旋转实现：
1．绕 z 轴旋转：将赤道面转到中间位置。
2．绕 y 轴旋转：调整赤道面倾角。
3．绕 z 轴旋转：完成最终对齐。
矩阵形式：

，

，

其中：
-  和为绕 z 轴和 y 轴的旋转矩阵。
- 角度 ,,是时间的多项式函数（以世纪为单位）：
                     

其中，为从J2000．0起算的儒略世纪数，。

2.2 章动矩阵
     

 章动（Nutation）是指地球自转轴在岁差运动基础上的周期性微小摆动。它主要由月球轨道平面相对于黄道面的变化引起，导致地球自转轴在空间中的方向出现周期性波动。可视化如下图所示(夸大了一点幅度，看着明显)。

图3  地球章动可视化

章动的分量

章动可分为两个主要分量（在赤道坐标系中）：
1．黄经章动 ：春分点在黄道上的周期性移动（影响赤经）。
2．交角章动（）：黄赤交角（）的周期性变化（影响赤纬）。
数学表达式：
                                      
其中：
- 为振幅系数（由天文观测或理论模型确定）。
-  为天文参数组合（如月球平近点角，太阳平近点角，月球交点经度等）。

参数具体定义可以参考IERS Conventions 2010。

章动矩阵的构造

（1）旋转步骤
1．绕 x 轴旋转 ：使黄道与瞬时真赤道对齐。
2．绕 z 轴旋转 ：修正春分点的黄经章动。
3．绕 x 轴旋转 ：回到初始参考框架。
（2）矩阵形式
章动矩阵可表示为：
                                         

，，

其中：
-  和  为绕 x 轴和 z 轴的旋转矩阵。
-  为平黄赤交角，。
-  和  为黄经章动和交角章动。

2.3 时角矩阵
   

   时角（Hour Angle, HA）是描述天体在赤道坐标系中相对于本地子午圈（观测者所在地的经线）的位置的角度，用于表示天体在当前时刻的方位变化。

时角矩阵表示为：

：绕Z轴旋转 格林尼治视恒星时（GAST） 的矩阵（对应时角部分）。

式中，GMST为格林尼治平恒星时，计算式为：

其中，GMST 0 为当天午夜的格林尼治平恒星时，r 为变化率，UT1 是极移修正世界时，为在儒略 （ 36525 天）中从 J 2000.0 到测量日的UT1零时经过的时间。

2.4 极移矩阵

       极移（Polar Motion）是指地球自转轴相对于地球表面的周期性微小摆动现象，表现为地理北极在地球表面围绕一个平均位置的运动。

时角矩阵表示为：

 和可以在IGS提供的ERP文件中获取。

三、总结
   

   J2000与WGS84坐标系转换过程如下所示。

其中为WGS84坐标系位置。为J2000坐标系位置，其余矩阵在第二章中详细介绍。