为什么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

为什么它们不相同或至少更接近？