使用Skyfield计算MODIS AQUA卫星位置与NASA LAADS数据不符的问题排查

我尝试用Skyfield工具包获取任意时间点MODIS AQUA卫星的位置，但计算结果和NASA LAADS页面的数据差异极大。比如2015年5月9日12:25时，LAADS显示卫星在瑞士/德国边境上空，但我的代码计算出来它在南太平洋上空，想知道是Skyfield使用有误，还是有其他原因导致分歧？

我的代码如下：

from skyfield.api import load, wgs84

# 获取TLE数据
satUrl = 'http://celestrak.com/NORAD/elements/active.txt'
sats = load.tle_file(satUrl)
sats = {sat.name: sat for sat in sats}
sat = sats['AQUA']

t = load.timescale().utc(2015, 5, 9, 12, 25)
sat_loc = wgs84.latlon_of(sat.at(t))
[x.degrees for x in sat_loc]  
# 输出结果：[-26.133974097154073, -172.11927450179587]

核心原因分析及解决办法 #

• TLE数据时效性不匹配：你用的是Celestrak的active.txt，里面是当前最新的TLE数据。卫星轨道会受大气阻力、引力摄动等因素长期漂移，用2024年的TLE去计算2015年的位置，结果必然完全错误。必须获取2015年5月9日前后的历史TLE数据才能得到准确结果。
• 时间基准可能混淆：确认LAADS页面的时间是UTC还是当地时间。2015年5月欧洲实行夏令时（CEST），12:25 CEST等于UTC 10:25。如果LAADS显示的是当地时间，你用UTC 12:25计算自然会差2小时，导致位置偏差。
• LAADS数据定义差异：LAADS显示的“卫星位置”可能不是卫星自身的实时经纬度，而是该时间点MODIS观测数据对应的地面覆盖区域中心。MODIS卫星过境时会扫描一条宽约2330公里的带状区域，LAADS可能展示的是数据条带的中心位置，而非卫星当时的位置。

内容的提问来源于stack exchange，提问作者AwkwardAttempts