astronomy - 为什么 Skyfield 提供的计算与实际数据不符？

Question

我想将 Skyfield 计算的半长轴与 TLE 的半长轴进行比较。

from skyfield.elementslib import osculating_elements_of as OE
from skyfield.api import EarthSatellite, load, utc

ts = load.timescale()

sat = EarthSatellite(TLE1, TLE2)
date = sat.epoch.utc
test = sat.at(ts.utc(*date))
print(OE(test).semi_major_axis.km)
print(sat_df['semimajor_axis_km'].values[0])

输出：

7087.916048058872
7080.642

为什么他们不一样或至少更接近？