我正在重构一些使用PyEphem来使用Skyfield的旧代码，并且我在主体的GHA/dec的结果中得到了轻微的差异。
def sf(year):
from skyfield.api import N, W, wgs84
from skyfield.api import load
from skyfield.units import Angle
ts = load.timescale()
eph = load('de421.bsp')
greenwich = wgs84.latlon(58.47722 * N, 0.0 * W)
t = ts.utc(year)
e = eph['earth'] + greenwich
v = eph['venus']
ha = e.at(t).observe(v).hadec()
gha = Angle(degrees=360.0 + ha[0]._degrees)
print(f'{gha}, {ha[1]}')
def pe(year):
import ephem
greenwich = ephem.Observer()
greenwich.lon = 0.0
greenwich.lat = ephem.degrees('51:28:38')
greenwich.pressure = 0.0
greenwich.horizon = '-0:34'
t = ephem.date(str(year))
v = ephem.Venus()
greenwich.date = t
greenwich.epoch = t
v.compute(greenwich)
gha = ephem.twopi - v.g_ra + greenwich.sidereal_time()
print(f'{ephem.degrees(gha)}, {ephem.degrees(v.g_dec)}')
if __name__ == '__main__':
sf(2016)
pe(2016)
产生结果：
219deg 41' 57.5", -18deg 37' 03.9"
219:42:15.7, -18:36:56.0
旧的PyEphem代码与我面前的航海年鉴一致。
这绝对是一个PEBKAC，但我在挠头，想知道我错过了什么时间或空间变换。
看起来你在向Skyfield询问天体测量位置，但向PyEphem询问明显位置。根据PyEphem文件：
"g_ ra和ra——日期时代的明显地心赤经"；
https://rhodesmill.org/pyephem/quick.html#body-计算日期
而对于Skyfield，你必须在一个位置上调用.apparent()方法来学习相应的表观位置；它不会自动发生：
https://rhodesmill.org/skyfield/positions.html#barycentric-天文视在
看看这种变化是否消除了坐标之间的大部分差异。