PyEphem 和 Skyfield HA/dec 结果之间的差异

Discrepancy between PyEphem and Skyfield HA/dec results

我正在重构一些使用 PyEphem 以使用 Skyfield 的旧代码,并且 body 的 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 — epoch-of-date"

的地心赤经

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

而使用Skyfield,你必须在一个位置上调用.apparent()方法来学习相应的表观位置;它不会自动发生:

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

看看该更改是否消除了坐标之间的大部分差异。