python 使用skyfield进行轨道计算
pip install skyfield
依赖numpy，scipy， 建议安装anaconda
计算火星在天空中的位置：
from skyfield.api import load
# Create a timescale and ask the current time.
ts = load.timescale()
t = ts.now()
# Load the JPL ephemeris DE421 (covers 1900-2050).
planets = load('de421.bsp')
earth, mars = planets['earth'], planets['mars']
# What's the position of Mars, viewed from Earth?
astrometric = earth.at(t).observe(mars)
ra, dec, distance = astrometric.radec()
print(ra)
print(dec)
print(distance)
输出为：
10h 47m 56.24s +09deg 03' 23.1" 2.33251 au
相比较PyEphem 而言，skyfiled 全部由python编写，更容易发挥python的优势，在后续的版本升级和迭代支持优于PyEphem，毕竟PyEphem的官网也写了，如果是新研发的项目，建议使用skyfield。
参见Skyfield — documentationhttps://rhodesmill.org/skyfield/
