Как быстро оценить расстояние между двумя точками (широта, долгота)?

Question

Как быстро оценить расстояние между двумя точками (широта, долгота)?

Я хочу иметь возможность получить оценку расстояния между двумя точками (широта, долгота). Я хочу пропустить, так как это будет для поиска * graph, и я хочу, чтобы это было быстро. Точки будут находиться не более чем в 800 км друг от друга.

588 5

python math mapping latitude-longitude

5 ответов:

Comments

Ничего не найдено.

Aaron D · Accepted Answer · 2017-10-19 02:40:06

Ответы на Формулу Хаверсина в Python (Пеленг и расстояние между двумя точками GPS)предоставляют реализации Python, которые отвечают на ваш вопрос.

Используя реализацию ниже, я выполнил 100 000 итераций менее чем за 1 секунду на старом ноутбуке. Я думаю, для ваших целей этого должно быть достаточно. Тем не менее, вы должны профилировать что-либо, прежде чем оптимизировать производительность. {[0]}

Недооценивать haversine(lat1, long1, lat2, long2) * 0.90 или любой другой фактор, который вы хотите. Я не посмотрите, как полезно вводить ошибку в вашу недооценку.

TreyA · Accepted Answer · 2013-04-01 14:06:20

Поскольку расстояние относительно мало, можно использовать приближение равновеликого расстояния. Это приближение выполняется быстрее, чем при использовании формулы Хаверсина. Итак, чтобы получить расстояние от вашей контрольной точки (lat1/lon1) до точки, которую вы тестируете (lat2/lon2), используйте приведенную ниже формулу. Важное примечание: вам нужно преобразовать все точки lat/lon в радианы:
R = 6371  // radius of the earth in km
x = (lon2 - lon1) * cos( 0.5*(lat2+lat1) )
y = lat2 - lat1
d = R * sqrt( x*x + y*y )
Поскольку 'R' находится в км, расстояние ' d ' будет в км.

Ссылка: http://www.movable-type.co.uk/scripts/latlong.html

Raymond Hettinger · Accepted Answer · 2013-04-01 05:56:34

Одна из идей скорости состоит в том, чтобы преобразовать координаты long/lat в трехмерные (x,y,z) координаты. После предварительной обработки точек используйте евклидово расстояние между точками в качестве быстро вычисляемого недоразвития фактического расстояния.

Aaron D · Accepted Answer · 2013-04-01 06:55:39

Для максимальной скорости можно создать что-то вроде радужной таблицы для координатных расстояний. Похоже, вы уже знаете область, с которой работаете, поэтому кажется, что предварительное вычисление их может быть осуществимо. Затем вы можете загрузить ближайшую комбинацию и просто использовать ее.
Например, в континентальных Соединенных Штатах долгота составляет 55 градусов, а широта-20, что составляет 1100 целых точек. Расстояние между всеми возможными комбинации-это Задача рукопожатия , на которую отвечают (n-1)(n)/2 или около 600k комбинаций. Это кажется вполне осуществимым для хранения и извлечения. Если вы предоставите больше информации о своих требованиях, я могу быть более конкретным.

uher · Accepted Answer · 2016-08-02 10:26:23

Пожалуйста, используйте следующий код.

def distance(lat1, lng1, lat2, lng2):
    #return distance as meter if you want km distance, remove "* 1000"
    radius = 6371 * 1000 

    dLat = (lat2-lat1) * math.pi / 180
    dLng = (lng2-lng1) * math.pi / 180

    lat1 = lat1 * math.pi / 180
    lat2 = lat2 * math.pi / 180

    val = sin(dLat/2) * sin(dLat/2) + sin(dLng/2) * sin(dLng/2) * cos(lat1) * cos(lat2)    
    ang = 2 * atan2(sqrt(val), sqrt(1-val))
    return radius * ang