как преобразовать градусы минуты секунды в десятичные градусы (питон / небо)

Question

Я хотел бы извлечь координату спутника (МКС) в десятичном виде, используя этот код:

from skyfield.api import EarthSatellite, Topos, load
import time

line1 = '1 25544U 98067A   14020.93268519  .00009878  00000-0  18200-3 0  5082'
line2 = '2 25544  51.6498 109.4756 0003572  55.9686 274.8005 15.49815350868473'

satellite = EarthSatellite(line1, line2, name='ISS (ZARYA)')

while True:
    ts = load.timescale()
    t = ts.now()
    geometry = satellite.at(t)

    subpoint = geometry.subpoint()
    print(subpoint.latitude)
    print('\n')
    print(subpoint.longitude)
    time.sleep(1)

Выходные данные - строка: -45deg 44' 13.5".

Какой самый простой способ преобразовать его в нечто вроде: -77.0089°?

Answer 1

К счастью, объекты latitude и longitude не являются простыми строками, а представляют собой объекты с необычными углами, которые просто печатаются в виде трехкомпонентных строк, чтобы их было легко читать на экране. Вы можете узнать о них больше, спросив у Python их документацию. В конце l oop попробуйте добавить:

help(subpoint.latitude)

Появится документация для класса Angle. Вы также можете найти его в Интернете здесь:

https://rhodesmill.org/skyfield/api-units.html#skyfield .units.Angle

Вы захотите использовать атрибут degrees, который выражает угол как число с плавающей запятой. Измените вызовы печати в вашей программе на:

print(subpoint.latitude.degrees)
print('\n')
print(subpoint.longitude.degrees)

Answer 2

Попробуйте это

from skyfield.api import EarthSatellite, Topos, load
import time
line1 = '1 25544U 98067A   14020.93268519  .00009878  00000-0  18200-3 0  5082'
line2 = '2 25544  51.6498 109.4756 0003572  55.9686 274.8005 15.49815350868473'

satellite = EarthSatellite(line1, line2, name='ISS (ZARYA)')

def convert(deg):
  d, m, s = str(deg).replace('deg', '').split(" ")
  ans = float(d) + (float(m.strip("'")) / 60) + (float(s.strip('"')) / 3600)
  return str(ans) + chr(176)

while True:
    ts = load.timescale()
    t = ts.now()
    geometry = satellite.at(t)

    subpoint = geometry.subpoint()
    lat = convert(subpoint.latitude)
    lng = convert(subpoint.longitude)
    print(lat, lng)
    time.sleep(1)

Вывод:

48.522305555555555° 133.80061111111112°
48.49586111111111° 133.89988888888888°
48.46933333333334° 133.99902777777777°
48.44269444444444° 134.09808333333334°
48.416° 134.19702777777778°