擅长:python、mysql、java
<p><strong>编辑:</strong>是!Skyfield现在直接支持transits,如其文档的新部分所述:</p>
<p><a href="https://rhodesmill.org/skyfield/almanac.html#meridian-transits" rel="nofollow noreferrer">https://rhodesmill.org/skyfield/almanac.html#meridian-transits</a></p>
<p>示例脚本:</p>
<pre><code>from skyfield import api
from skyfield import almanac
bluffton = api.Topos('40.8939 N', '83.8917 W')
ts = api.load.timescale()
t0 = ts.utc(2020, 11, 6)
t1 = ts.utc(2020, 11, 7)
eph = api.load('de421.bsp')
f = almanac.meridian_transits(eph, eph['Mars'], bluffton)
t, y = almanac.find_discrete(t0, t1, f)
print(t.utc_strftime('%Y-%m-%d %H:%M'))
print(y)
print([almanac.MERIDIAN_TRANSITS[yi] for yi in y])
</code></pre>
<p>其输出:</p>
<pre><code>['2020-11-06 03:32', '2020-11-06 15:30']
[1 0]
['Meridian transit', 'Antimeridian transit']
</code></pre>