使用skyfield.api

Question

我尝试使用 tkinter ui 制作一个“简单”的程序，它可以在特定时间返回给定 TLE 的正确上升和下降：

from skyfield.api import EarthSatellite, Topos, load
import numpy
import tkinter as tk
from tkinter import ttk
import threading
import datetime
import pyfiglet
import urllib.request
import os


def localisation(): #I take this function from a website, but I adapted it to my usage and UI
  ts = load.timescale()
  yearstring = yearentry.get()+"/"+monthentry.get()+"/"+dayentry.get()
  daystring = hoursentry.get()+"h:"+minutentry.get()+"m:"+secondentry.get()+"s"
  timestring= yearstring+" "+daystring
  print(timestring)
  t = ts.utc(int(yearentry.get()), int(monthentry.get()), float(dayentry.get()), float(hoursentry.get()), float(minutentry.get()), float(secondentry.get())) # get time from Tkinter UI

  line1,line2=parsedData[satelliteselector.get()]
  loc = Topos(latentry.get(), logentry.get(), elevation_m=altentry.get()) 
  satellite = EarthSatellite(line1, line2)

  # Geocentric
  geometry = satellite.at(t)

  # Geographic point beneath satellite
  subpoint = geometry.subpoint()
  latitude = subpoint.latitude
  longitude = subpoint.longitude
  elevation = subpoint.elevation

  # Topocentric
  difference = satellite - loc
  geometry = difference.at(t)
  topoc= loc.at(t)
  #
  topocentric = difference.at(t)
  geocentric = satellite.at(t)
  # ------ Start outputs -----------
  print ('\n Heures locale :', timestring)
  # print (' Temps Solaire : ',t) pb
  print ('',loc)
  print ('\n Subpoint Longitude= ', longitude )
  print (' Subpoint Latitude = ', latitude )
  print (' Subpoint Elevation=  {0:.3f}'.format(elevation.km),'km')
  # ------ Step 1: compute sat horizontal coords ------
  alt, az, distance = topocentric.altaz()
  if alt.degrees > 0:
    print('\n',satellite, "\n est au dessu de l'horizon")
  print ('\n Altitude= ', alt )
  print (' Azimute = ', az )
  print (' Distance=  {0:.3f}'.format(distance.km), 'km')

在基础 function 中，在 loc = Topos(latentry.get()、logentry.get()、elevation_m=altentry.get()) 中，只有纬度和经度，但我想考虑海拔

但是当我尝试计算 RA 和 D 时，我有这个错误

Exception in Tkinter callback
Traceback (most recent call last):
  File "C:\Users\Utilisateur\AppData\Local\Programs\Python\Python39\lib\tkinter\__init__.py", line 1884, in __call__
    return self.func(*args)
  File "g:\Calamar Industries\Astrée\python\Astrée\Astrée.py", line 32, in localisation
    loc = Topos(latentry.get(), logentry.get(), elevation_m=altentry.get())    # location coords
  File "C:\Users\Utilisateur\AppData\Local\Programs\Python\Python39\lib\site-packages\skyfield\toposlib.py", line 257, in __init__
    iers2010.latlon(latitude_degrees, longitude_degrees, elevation_m,
  File "C:\Users\Utilisateur\AppData\Local\Programs\Python\Python39\lib\site-packages\skyfield\toposlib.py", line 142, in latlon
    elevation = Distance(m=_to_array(elevation_m))
  File "C:\Users\Utilisateur\AppData\Local\Programs\Python\Python39\lib\site-packages\skyfield\units.py", line 54, in __init__
    self.au = m / AU_M
TypeError: ufunc 'true_divide' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe''

我知道代码没有经过优化，但它应该可以工作。 你可以帮帮我吗？

奖金：

 if alt.degrees > 0:
    print('\n',satellite, "\n est au dessu de l'horizon")

这部分基于 function 也不起作用： AttributeError: 'Angle' object has no attribute 'degree' （当我尝试打印它时）