Die Distanz zwischen zwei Breitengrad-Längengrad-Punkten berechnen? (Haversine-Formel)

Question

Wie berechne ich die Entfernung zwischen zwei Punkten, die durch Breiten- und Längengrad angegeben sind?

Zur Klarstellung möchte ich die Entfernung in Kilometern; die Punkte verwenden das WGS84-System und ich möchte die relativen Genauigkeiten der verfügbaren Methoden verstehen.

Answer 1

Sie können es mit der Haversine-Formel berechnen, die lautet:

a = sin²(/2) + cos 1  cos 2  sin²(/2)
c = 2  atan2( a, (1a) )
d = R  c

Ein Beispiel zur Berechnung der Entfernung zwischen zwei Punkten finden Sie unten

Angenommen, ich muss die Entfernung zwischen Neu-Delhi und London berechnen, wie kann ich diese Formel nutzen:

Neu Delhi Koordinaten= 28.7041° N, 77.1025° E
London Koordinaten= 51.5074° N, 0.1278° W

var R = 6371e3; // Meter
var 1 = 28.7041.toRadians();
var 2 = 51.5074.toRadians();
var  = (51.5074-28.7041).toRadians();
var  = (0.1278-77.1025).toRadians();

var a = Math.sin(/2) * Math.sin(/2) +
        Math.cos(1) * Math.cos(2) *
        Math.sin(/2) * Math.sin(/2);
var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));

var d = R * c; // Meter
d = d/1000; // km

Answer 2

Es gibt ein gutes Beispiel hier, um die Entfernung mit PHP zu berechnen http://www.geodatasource.com/developers/php:

 function distance($lat1, $lon1, $lat2, $lon2, $unit) {

     $theta = $lon1 - $lon2;
     $dist = sin(deg2rad($lat1)) * sin(deg2rad($lat2)) +  cos(deg2rad($lat1)) * cos(deg2rad($lat2)) * cos(deg2rad($theta));
     $dist = acos($dist);
     $dist = rad2deg($dist);
     $miles = $dist * 60 * 1.1515;
     $unit = strtoupper($unit);

     if ($unit == "K") {
         return ($miles * 1.609344);
     } else if ($unit == "N") {
          return ($miles * 0.8684);
     } else {
          return $miles;
     }
 }

Answer 3

Hatte ein Problem mit math.deg in LUA... falls jemand eine Lösung kennt, bitte bereinigen Sie diesen Code!

In der Zwischenzeit hier eine Implementierung des Haversine in LUA (verwenden Sie dies mit Redis!)

function calcDist(lat1, lon1, lat2, lon2)
    lat1= lat1*0.0174532925
    lat2= lat2*0.0174532925
    lon1= lon1*0.0174532925
    lon2= lon2*0.0174532925

    dlon = lon2-lon1
    dlat = lat2-lat1

    a = math.pow(math.sin(dlat/2),2) + math.cos(lat1) * math.cos(lat2) * math.pow(math.sin(dlon/2),2)
    c = 2 * math.asin(math.sqrt(a))
    dist = 6371 * c      -- multipliziere mit 0.621371, um in Meilen umzurechnen
    return dist
end

Prost!

Answer 4

Hier ist die Erlang-Implementierung

lat_lng({Lat1, Lon1}=_Point1, {Lat2, Lon2}=_Point2) ->
  P = math:pi() / 180,
  R = 6371, % Radius der Erde in KM
  A = 0.5 - math:cos((Lat2 - Lat1) * P) / 2 +
    math:cos(Lat1 * P) * math:cos(Lat2 * P) * (1 - math:cos((Lon2 - Lon1) * P))/2,
  R * 2 * math:asin(math:sqrt(A)).

Answer 5

Wenn Sie Python verwenden; pip install geopy

from geopy.distance import geodesic

origin = (30.172705, 31.526725)  # (Breitengrad, Längengrad) nicht verwechseln
destination = (30.288281, 31.732326)

print(geodesic(origin, destination).meter)  # 23576.805481751613
print(geodesic(origin, destination).kilometers)  # 23.576805481751613
print(geodesic(origin, destination).miles)  # 14.64994773134371