Как вычислить румб по двум координатам?

Question

[Данил Екимов]

Всем привет, стоит задача вычислить румб по двум координатам

def haversine(lat1, lon1, lat2, lon2):
    # Конвертируем в список
    lon1, lat1, lon2, lat2 = map(radians, (lon1, lat1, lon2, lat2))

    dlon = lon2 - lon1
    dlat = lat2 - lat1
    a = sin(dlat / 2) ** 2 + cos(lat1) * cos(lat2) * sin(dlon / 2) ** 2
    c = 2 * asin(sqrt(a))
    #радиус земли(под вопросом, везде по разному указан)
    km =6371 * c
    km = km*1000
    km = toFixed(km, '2')
    r = (atan((dlat)/(dlon)))
    print ('Расстояние в метрах:', km)
    print ('Румбы(ещё не работает):', r)

С расстоянием всё отлично работает, выводится нормально. А вот что касается румбов, получается не очень...

Вот к примеру 2 точки:
54.5244194444444 61.2879
54.5246944444444 61.2872138888889

Расстояние в метрах: 104.62
Румбы(ещё не работает): 0.27477411336095403
Хотя румб на этих координатах должен быть: ~~56...Выручайте)

Comments

[hint000]

Вы уверены, что правильно пользуетесь терминологией? (румб это...)
Возможно, вам нужен азимут или пеленг?

Ну и atan у вас выдаст результат в радианах, разумеется. Надо перевести в градусы (разделить на Пи и умножить на 180).
(56 всё равно не получится)

[Данил Екимов]

hint000, Хотя бы азимут получить)

[Данил Екимов]

hint000, Азимут тут 304 должен быть, сейчас попробую его высчитать...Я в данной сфере не так давно, ещё разбираюсь, спасибо за ответ. Отпишусь, если всё получится

[hint000]

Данил Екимов, возможно, ваша фомула даст почти правильный результат, когда обе точки находятся близко к экватору.
Там длина дуги меридиана почти равна длине дуги параллели.
В вашей формуле r = (atan((dlat)/(dlon))) вместо dlat и dlon в градусах должны быть расстояния в метрах в соответствующих направлениях (по меридиану и по параллели).

[hint000]

это, конечно, грубое приближение, как будто Земля - идеальный шар, а не эллипсоид или геоид. Вот еще хорошая табличка.

[Данил Екимов]

hint000, Хорошо понял, кажется уже разобрался, функцию из паблика, мне выдаётся расстояние и азимут, а зная азимут не сложно румб получить. Большое спасибо за Вашу помощь и ответы, сейчас потихоньку разбираюсь во всех этих еденицах)

Answer 1

[lz961]

atan дает ответы в радианах, а вам, наверное, градусы нужны. Хотя все равно не спасет, поскольку для небольших расстояний нужно использовать поправку на длину параллели:
dx = dlon/cos(math.radians(lat))
dy = dlat
azimuth = 90-math.degrees(math.atan2(dy, dx))
для больших расстояний нужны формулы
gis-lab.info/qa/great-circles.html
gis-lab.info/qa/angles-sphere.html