我正在探索 Brandon Rhodes 開發的出色軟件 Skyfield 的可能性。我制作了一個腳本來計算隨機對象之間的 Right Ascension 連詞。我使用以下腳本：from skyfield import almanacfrom skyfield.searchlib import find_maxima, find_minima, find_discretefrom skyfield.api import Star, loadfrom datetime import datetime, date,timedeltaimport pytzplanets = load('de430t.bsp')earth = planets['earth']x = [['Aldebaran',[4, 35, 55.2],[16, 30, 33]],['Regulus',[10, 8, 22.3],[11, 58, 2]],['Pollux',[7, 45, 18.9],[28, 1, 34]],['Antares',[16, 29, 24.4],[-26, 25, 55]],]ts = load.timescale(builtin=True)t = ts.now()tzn = 'Europe/Amsterdam'tz = pytz.timezone(tzn)now = datetime(2020, 1, 1, 0, 0, 0)t0 = ts.utc(tz.localize(now))t1 = ts.utc(tz.localize(now) + timedelta(days=+365))def difference(t):    e = earth.at(t)    ra11, dec11, distance = e.observe(object).radec()    ra12, dec12, distance2 = e.observe(planets[target1]).radec()    diff = ra11.hours - ra12.hours    return diff >= 0difference.rough_period = 1.0for count in range (len(x)):    object = Star(ra_hours=(x[count][1][0], x[count][1][1], x[count][1][2]),dec_degrees=(x[count][2][0], x[count][2][1], x[count][2][2]))    target1 = 'venus'    t, b = find_discrete(t0, t1, difference)    if len(t) > 0:        print (f"{x[count][0]} and {target1}")        for i, (a, b) in enumerate(zip(t, b)):            e = earth.at(a)            ra11, dec11, distance = e.observe(object).radec('date')            ra12, dec12, distance2 = e.observe(planets[target1]).radec('date')            print (f"Diff: {ra11.hours - ra12.hours}, ra_{x[count][0]}: {ra11.hours}, ra_{target1}: {ra12.hours}")            print(f"{a.utc_iso()},{dec11._degrees - dec12._degrees}")            print ("")我相信這是在計算兩個對象具有相同 RA 的時間實例。以“Diff”開頭的行是監視輸出有效性的行。Diff 代表計算出的 RA 差異。它應該接近于零。其他兩個值是兩個對象的赤經。它們應該非常相似。第二行是我想要的結果，它是計算的時間和對象之間的度數距離。正如您所看到的，出于某種原因，對于每組對象，我得到了時間實例的無效結果：2020-02-07T21:20:06Z，并且該實例的差值肯定不接近于零。如果我將對象 venus 更改為 moon 它會變得更糟，因為每秒結果都是無效的。我根據 Skychart / Cartes du Ciel 軟件和那些結帳檢查了其他結果。我不知道這里出了什么問題。有人能幫助我嗎？