在 Python 中计算 ECI 到 ECEF 的转换矩阵

如何解决在 Python 中计算 ECI 到 ECEF 的转换矩阵

目前，我正在按照 MODIS Level 1A Earth Location: Algorithm Theoretical Basis Document 尝试将卫星观看矢量从 ECI (J2000) 转换为 ECR（例如 WGS84）参考系。

要做这个转换，首先要计算转换矩阵，定义为：

T = ABCD

其中A是极地运动，B是恒星时，C是天文章动，D > 是岁差。基本上，我正在寻找一个将 B 作为输入并输出 T 的 python 函数：

f(B) = T

不过，我好像在网上找不到这样的功能。因此，我很想自己编写它，但更愿意使用其他人尝试和测试过的东西。

编辑：

基本上，我想在 Python 中复制 here 中描述的 rECEFtoECI 函数。

编辑 2：

SOFA C/FORTRAN 库具有完成上述任务的函数（例如 iauC2t06a）。有 2 个项目试图为函数编写 Python 包装器：pysofa 和 pysofa2。但是，pysofa 不适用于 Python3，pysofa2 还没有封装计算 ECI 到 ECEF 转换矩阵的函数，而且近期似乎没有这样做的计划。