目前,我正在按照 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 转换矩阵的函数,而且近期似乎没有这样做的计划。