PyEphem与Skyfield HA/dec结果的差异

Question

我正在重构一些使用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，但我正在抓我的头来找出我错过的时间或空间转换。

Answer 1

看起来你是在向Skyfield询问占星术的位置，但PyEphem需要的是明显的位置。根据PyEphem文档：

"g_ra和ra -表观地心右提升为日期“

https://rhodesmill.org/pyephem/quick.html#body-compute-date

然而，对于Skyfield，您必须在一个位置上调用.apparent()方法来学习相应的表观位置；它不会自动发生：

https://rhodesmill.org/skyfield/positions.html#barycentric-astrometric-apparent

看看这种变化是否消除了坐标之间的大部分差异。