Question

我正在使用skyfield计算行星与地球的相对距离作为时间的函数（如Skyfield主页中所述）。它工作得很好，现在我正在尝试实现Earth =＆gt;彗星距离（例如67P / Tchouri）。

我在NASA JPL上找到了一种为彗星（here）创建Spice SPK文件的方法，但它产生的xsp文件似乎无法用来自skyfield的load命令读取。

我考虑的另一种可能性是使用pyephem建议的轨道信息（参见here），但我不知道如何在Skyfield中阅读它们。我还看到彗星是Skyfield编码冲刺的路线图，所以这可能是我的答案，但如果你知道一种方法可以使它与当前版本一起工作，那将非常有帮助。

谢谢！

Answer 1

Skyfield确实获得了彗星的支持！您可以在此处找到详细信息：

https://rhodesmill.org/skyfield/kepler-orbits.html#comets

调整文档中的代码，以下是从小行星中心数据库到彗星的距离：

.findOne({
  $or: [ 
    { "slug.fr-FR": "test-1" },
    { "slug.nl-NL": "test-1" },
    { "slug.en-US": "test-1" }
  ]
})

我看到结果：

from skyfield.api import load
from skyfield.constants import GM_SUN_Pitjeva_2005_km3_s2 as GM_SUN
from skyfield.data import mpc

ts = load.timescale()

eph = load('de421.bsp')
sun, earth = eph['sun'], eph['earth']

with load.open(mpc.COMET_URL) as f:
    comets = mpc.load_comets_dataframe(f)
comets = comets.set_index('designation', drop=False)
row = comets.loc['1P/Halley']
comet = sun + mpc.comet_orbit(row, ts, GM_SUN)

t = ts.utc(2020, 10, 17)
ra, dec, distance = earth.at(t).observe(comet).radec()

print('Distance in AU:', distance.au)

使用Skyfields

1 个答案: