두 위도와 경도 지리 좌표 사이의 거리 계산


139

두 개의 GeoCoordinates 사이의 거리를 계산하고 있습니다. 3-4 개의 다른 앱에 대해 앱을 테스트하고 있습니다. 거리를 계산할 때 계산을 위해 평균 3.3 마일을 얻는 반면 다른 앱은 3.5 마일을 얻는 경향이 있습니다. 내가 수행하려는 계산에 큰 차이가 있습니다. 거리를 계산하기 위해 좋은 클래스 라이브러리가 있습니까? C #에서 다음과 같이 계산하고 있습니다.

public static double Calculate(double sLatitude,double sLongitude, double eLatitude, 
                               double eLongitude)
{
    var radiansOverDegrees = (Math.PI / 180.0);

    var sLatitudeRadians = sLatitude * radiansOverDegrees;
    var sLongitudeRadians = sLongitude * radiansOverDegrees;
    var eLatitudeRadians = eLatitude * radiansOverDegrees;
    var eLongitudeRadians = eLongitude * radiansOverDegrees;

    var dLongitude = eLongitudeRadians - sLongitudeRadians;
    var dLatitude = eLatitudeRadians - sLatitudeRadians;

    var result1 = Math.Pow(Math.Sin(dLatitude / 2.0), 2.0) + 
                  Math.Cos(sLatitudeRadians) * Math.Cos(eLatitudeRadians) * 
                  Math.Pow(Math.Sin(dLongitude / 2.0), 2.0);

    // Using 3956 as the number of miles around the earth
    var result2 = 3956.0 * 2.0 * 
                  Math.Atan2(Math.Sqrt(result1), Math.Sqrt(1.0 - result1));

    return result2;
}

내가 뭘 잘못하고 있니? 먼저 km 단위로 계산 한 다음 마일로 변환해야합니까?


1
지구 평균 반경 = 6,371km = 3958.76 마일
미치 밀


이 gis.stackexchange.com에 안
다니엘 파월

그것은 가능하지만 내 질문은 조금 다른 Windows Phone에서 이것을 계산하는 것과 관련이 있습니다. 수식은 동일하지만 DistanceTo 메서드와 같은 최신 메서드 호출을 반드시 사용할 필요는 없습니다.
Jason N. Gaylord

1
계산을 계속 반복 할 필요가 없도록 pi / 180을 저장하도록 제안하십시오.
Chris Caviness

답변:


313

정보 좌표 클래스 (.NET 프레임 워크 4 이상) 이미이 GetDistanceTo방법을.

var sCoord = new GeoCoordinate(sLatitude, sLongitude);
var eCoord = new GeoCoordinate(eLatitude, eLongitude);

return sCoord.GetDistanceTo(eCoord);

거리는 미터입니다.

System.Device를 참조해야합니다.


Nigel, DistanceTo 방법이 전화에서 작동합니까? WP7에 2.0 버전의 GeoCoordinate를 사용한다고 생각했습니다.
Jason N. Gaylord

1
나는 이것을 확인했고 장치의 GeoCordinate에는 당신이 참조 한 GetDistanceTo 메소드가 있습니다 (그러나 위에서 언급 한 것은 아님). 별거 아냐 내장 계산이 더 나은지 확인하기 위해 이것을 테스트 할 것입니다. 감사합니다 나이젤!
Jason N. Gaylord

1
나는 잘못된 질문을 할 수 있지만 결과는 어떤 단위입니까? 마일 또는 킬로미터입니다. 어디서나 찾을 수 없습니다.
Saeed Neamati

3
msdn.microsoft.com/en-us/library/ 에 따르면 @SaeedNeamati-이 역시 찾고있었습니다 .
Andy Butland

예, GeoCoordinate.GetDistanceTo ()는 값을 미터 단위로 반환합니다. 저에게 미국의 경우 1610보다 작 으면 피트 (미터 * 3.28084)로 변환하고 그렇지 않으면 마일 (미터 * 0.000621371)로 변환합니다. 내 목적에 정확도가 충분합니다.
user3235770

110

GetDistance가 가장 좋은 솔루션 이지만 대부분의 경우이 방법 (예 : Universal App)을 사용할 수 없습니다

  • 코디네이션 사이 의 거리계산하기 위한 알고리즘의 의사 코드 :

    public static double DistanceTo(double lat1, double lon1, double lat2, double lon2, char unit = 'K')
    {
        double rlat1 = Math.PI*lat1/180;
        double rlat2 = Math.PI*lat2/180;
        double theta = lon1 - lon2;
        double rtheta = Math.PI*theta/180;
        double dist =
            Math.Sin(rlat1)*Math.Sin(rlat2) + Math.Cos(rlat1)*
            Math.Cos(rlat2)*Math.Cos(rtheta);
        dist = Math.Acos(dist);
        dist = dist*180/Math.PI;
        dist = dist*60*1.1515;
    
        switch (unit)
        {
            case 'K': //Kilometers -> default
                return dist*1.609344;
            case 'N': //Nautical Miles 
                return dist*0.8684;
            case 'M': //Miles
                return dist;
        }
    
        return dist;
    }
  • 확장 메서드를 사용하는 Real World C # 구현

    용법:

    var distance = new Coordinates(48.672309, 15.695585)
                    .DistanceTo(
                        new Coordinates(48.237867, 16.389477),
                        UnitOfLength.Kilometers
                    );

    이행:

    public class Coordinates
    {
        public double Latitude { get; private set; }
        public double Longitude { get; private set; }
    
        public Coordinates(double latitude, double longitude)
        {
            Latitude = latitude;
            Longitude = longitude;
        }
    }
    public static class CoordinatesDistanceExtensions
    {
        public static double DistanceTo(this Coordinates baseCoordinates, Coordinates targetCoordinates)
        {
            return DistanceTo(baseCoordinates, targetCoordinates, UnitOfLength.Kilometers);
        }
    
        public static double DistanceTo(this Coordinates baseCoordinates, Coordinates targetCoordinates, UnitOfLength unitOfLength)
        {
            var baseRad = Math.PI * baseCoordinates.Latitude / 180;
            var targetRad = Math.PI * targetCoordinates.Latitude/ 180;
            var theta = baseCoordinates.Longitude - targetCoordinates.Longitude;
            var thetaRad = Math.PI * theta / 180;
    
            double dist =
                Math.Sin(baseRad) * Math.Sin(targetRad) + Math.Cos(baseRad) *
                Math.Cos(targetRad) * Math.Cos(thetaRad);
            dist = Math.Acos(dist);
    
            dist = dist * 180 / Math.PI;
            dist = dist * 60 * 1.1515;
    
            return unitOfLength.ConvertFromMiles(dist);
        }
    }
    
    public class UnitOfLength
    {
        public static UnitOfLength Kilometers = new UnitOfLength(1.609344);
        public static UnitOfLength NauticalMiles = new UnitOfLength(0.8684);
        public static UnitOfLength Miles = new UnitOfLength(1);
    
        private readonly double _fromMilesFactor;
    
        private UnitOfLength(double fromMilesFactor)
        {
            _fromMilesFactor = fromMilesFactor;
        }
    
        public double ConvertFromMiles(double input)
        {
            return input*_fromMilesFactor;
        }
    } 

1
이 미적분학에 사용 된 공식을 제공 할 수 있습니까? 변환하지 않고 마일 대신 결과 거리를 Km으로 직접 변경하려면 무엇을 변경해야합니까?
AlbertoFdzM

좋은 솔루션에 감사드립니다. 이제 데스크톱 응용 프로그램에서 사용할 수 있습니다.
Jamshaid Kamran

GeoCoordinate를 사용할 수없는 UWP 앱에서 훌륭하게 작동했습니다.
Zach Green

1
계산은 95 %입니다. 아래 함수는 100 % 정확합니다. stackoverflow.com/a/51839058/3736063
Malek Tubaisaht

31

그리고 여전히 만족스럽지 않은 사람들을 위해 .NET-Frameworks GeoCoordinate클래스 의 원본 코드 는 독립형 메서드로 리팩토링되었습니다.

public double GetDistance(double longitude, double latitude, double otherLongitude, double otherLatitude)
{
    var d1 = latitude * (Math.PI / 180.0);
    var num1 = longitude * (Math.PI / 180.0);
    var d2 = otherLatitude * (Math.PI / 180.0);
    var num2 = otherLongitude * (Math.PI / 180.0) - num1;
    var d3 = Math.Pow(Math.Sin((d2 - d1) / 2.0), 2.0) + Math.Cos(d1) * Math.Cos(d2) * Math.Pow(Math.Sin(num2 / 2.0), 2.0);

    return 6376500.0 * (2.0 * Math.Atan2(Math.Sqrt(d3), Math.Sqrt(1.0 - d3)));
}

8
아름다운 대답, 결과 거리가 미터 단위임을 지적하고 싶습니다. 공식 문서에
LeviathanCode

감사! GeoCoordinate 클래스에서 사용되는 실제 지구 반경을 찾고있었습니다.
KRoy

사소한 최적화 또는 더 쉬운 판독을 위해 pi / 180을 사전 계산할 수 double oneDegree = Math.PI / 180.0;있습니까?
Brakeroo

1
@brakeroo 답장을 보내 주셔서 감사합니다. 이것이 원래 .NET 코드이기 때문에 대답을 그대로두고 싶습니다. 물론 누구나 당신의 제안을 자유롭게 따르십시오.
마크

17

다음은 JavaScript 버전입니다.

function distanceTo(lat1, lon1, lat2, lon2, unit) {
      var rlat1 = Math.PI * lat1/180
      var rlat2 = Math.PI * lat2/180
      var rlon1 = Math.PI * lon1/180
      var rlon2 = Math.PI * lon2/180
      var theta = lon1-lon2
      var rtheta = Math.PI * theta/180
      var dist = Math.sin(rlat1) * Math.sin(rlat2) + Math.cos(rlat1) * Math.cos(rlat2) * Math.cos(rtheta);
      dist = Math.acos(dist)
      dist = dist * 180/Math.PI
      dist = dist * 60 * 1.1515
      if (unit=="K") { dist = dist * 1.609344 }
      if (unit=="N") { dist = dist * 0.8684 }
      return dist
}

10

Xamarin을 사용하고 있고 GeoCoordinate 클래스에 액세스 할 수없는 사용자는 대신 Android Location 클래스를 사용할 수 있습니다.

public static double GetDistanceBetweenCoordinates (double lat1, double lng1, double lat2, double lng2) {
            var coords1 = new Location ("");
            coords1.Latitude = lat1;
            coords1.Longitude = lng1;
            var coords2 = new Location ("");
            coords2.Latitude = lat2;
            coords2.Longitude = lng2;
            return coords1.DistanceTo (coords2);
        }

3

이 기능을 사용할 수 있습니다 :

출처 : https://www.geodatasource.com/developers/c-sharp

private double distance(double lat1, double lon1, double lat2, double lon2, char unit) {
  if ((lat1 == lat2) && (lon1 == lon2)) {
    return 0;
  }
  else {
    double theta = lon1 - lon2;
    double dist = Math.Sin(deg2rad(lat1)) * Math.Sin(deg2rad(lat2)) + Math.Cos(deg2rad(lat1)) * Math.Cos(deg2rad(lat2)) * Math.Cos(deg2rad(theta));
    dist = Math.Acos(dist);
    dist = rad2deg(dist);
    dist = dist * 60 * 1.1515;
    if (unit == 'K') {
      dist = dist * 1.609344;
    } else if (unit == 'N') {
      dist = dist * 0.8684;
    }
    return (dist);
  }
}

//:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
//::  This function converts decimal degrees to radians             :::
//:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
private double deg2rad(double deg) {
  return (deg * Math.PI / 180.0);
}

//:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
//::  This function converts radians to decimal degrees             :::
//:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
private double rad2deg(double rad) {
  return (rad / Math.PI * 180.0);
}

Console.WriteLine(distance(32.9697, -96.80322, 29.46786, -98.53506, "M"));
Console.WriteLine(distance(32.9697, -96.80322, 29.46786, -98.53506, "K"));
Console.WriteLine(distance(32.9697, -96.80322, 29.46786, -98.53506, "N"));

완벽하게 작동합니다! 감사!
Schnapz

3

이 플랫폼을위한 이 라이브러리 GeoCoordinate 가 있습니다 :

  • 단 핵증
  • .NET 4.5
  • .NET 코어
  • 윈도우 폰 8.x
  • 범용 Windows 플랫폼
  • Xamarin iOS
  • Xamarin 안드로이드

NuGet을 통해 설치가 완료됩니다.

PM> 설치 패키지 지리 좌표

용법

GeoCoordinate pin1 = new GeoCoordinate(lat, lng);
GeoCoordinate pin2 = new GeoCoordinate(lat, lng);

double distanceBetween = pin1.GetDistanceTo(pin2);

두 좌표 사이의 거리 ( 미터) 입니다.


3

Elliot Wood의 기능을 기반으로하며 C 기능에 관심이있는 사람은 작동합니다 ...

#define SIM_Degree_to_Radian(x) ((float)x * 0.017453292F)
#define SIM_PI_VALUE                         (3.14159265359)

float GPS_Distance(float lat1, float lon1, float lat2, float lon2)
{
   float theta;
   float dist;

   theta = lon1 - lon2;

   lat1 = SIM_Degree_to_Radian(lat1);
   lat2 = SIM_Degree_to_Radian(lat2);
   theta = SIM_Degree_to_Radian(theta);

   dist = (sin(lat1) * sin(lat2)) + (cos(lat1) * cos(lat2) * cos(theta));
   dist = acos(dist);

//   dist = dist * 180.0 / SIM_PI_VALUE;
//   dist = dist * 60.0 * 1.1515;
//   /* Convert to km */
//   dist = dist * 1.609344;

   dist *= 6370.693486F;

   return (dist);
}

double로 변경할 수 있습니다 . km 단위의 값을 반환합니다.


2

위도와 경도 지점 간의 거리 계산 ...

        double Lat1 = Convert.ToDouble(latitude);
        double Long1 = Convert.ToDouble(longitude);

        double Lat2 = 30.678;
        double Long2 = 45.786;
        double circumference = 40000.0; // Earth's circumference at the equator in km
        double distance = 0.0;
        double latitude1Rad = DegreesToRadians(Lat1);
        double latititude2Rad = DegreesToRadians(Lat2);
        double longitude1Rad = DegreesToRadians(Long1);
        double longitude2Rad = DegreesToRadians(Long2);
        double logitudeDiff = Math.Abs(longitude1Rad - longitude2Rad);
        if (logitudeDiff > Math.PI)
        {
            logitudeDiff = 2.0 * Math.PI - logitudeDiff;
        }
        double angleCalculation =
            Math.Acos(
              Math.Sin(latititude2Rad) * Math.Sin(latitude1Rad) +
              Math.Cos(latititude2Rad) * Math.Cos(latitude1Rad) * Math.Cos(logitudeDiff));
        distance = circumference * angleCalculation / (2.0 * Math.PI);
        return distance;

1

이것은 오래된 질문이지만 그럼에도 불구하고 답변이 성능 및 최적화와 관련하여 나를 만족시키지 못했습니다.

여기에 최적화 된 C # 변형 (변수 및 중복 계산이없는 거리 (km), Haversine Formular https://en.wikipedia.org/wiki/Haversine_formula 의 수학적 표현에 매우 가깝습니다 ).

영감을 얻은 사람 : https://rosettacode.org/wiki/Haversine_formula#C.23

public static class Haversine
{
    public static double Calculate(double lat1, double lon1, double lat2, double lon2)
    {
        double rad(double angle) => angle * 0.017453292519943295769236907684886127d; // = angle * Math.Pi / 180.0d
        double havf(double diff) => Math.Pow(Math.Sin(rad(diff) / 2d), 2); // = sin²(diff / 2)
        return 12745.6 * Math.Asin(Math.Sqrt(havf(lat2 - lat1) + Math.Cos(rad(lat1)) * Math.Cos(rad(lat2)) * havf(lon2 - lon1))); // earth radius 6.372,8‬km x 2 = 12745.6
    }
}

위키피디아의 Haversine Formular


0

이 시도:

    public double getDistance(GeoCoordinate p1, GeoCoordinate p2)
    {
        double d = p1.Latitude * 0.017453292519943295;
        double num3 = p1.Longitude * 0.017453292519943295;
        double num4 = p2.Latitude * 0.017453292519943295;
        double num5 = p2.Longitude * 0.017453292519943295;
        double num6 = num5 - num3;
        double num7 = num4 - d;
        double num8 = Math.Pow(Math.Sin(num7 / 2.0), 2.0) + ((Math.Cos(d) * Math.Cos(num4)) * Math.Pow(Math.Sin(num6 / 2.0), 2.0));
        double num9 = 2.0 * Math.Atan2(Math.Sqrt(num8), Math.Sqrt(1.0 - num8));
        return (6376500.0 * num9);
    }

0

당신은 사용할 수 있습니다 System.device.Location:

System.device.Location.GeoCoordinate gc = new System.device.Location.GeoCoordinate(){
Latitude = yourLatitudePt1,
Longitude = yourLongitudePt1
};

System.device.Location.GeoCoordinate gc2 = new System.device.Location.GeoCoordinate(){
Latitude = yourLatitudePt2,
Longitude = yourLongitudePt2
};

Double distance = gc2.getDistanceTo(gc);

행운을 빕니다


0

CPU / math 컴퓨팅 성능이 제한되는 경우 :

계산 기능이 부족한 경우 (예 : 부동 소수점 프로세서 없음, 소형 마이크로 컨트롤러 사용), 일부 삼각 함수가 엄청난 양의 CPU 시간 (예 : 3000+ 클럭주기)을 소비 할 수있는 경우가 있습니다 (예 : 3000+ 클럭주기). 근사치 만 필요합니다. 특히 CPU를 오랫동안 묶지 말아야하는 경우 CPU 오버 헤드를 최소화하기 위해 이것을 사용합니다.

/**------------------------------------------------------------------------
 * \brief  Great Circle distance approximation in km over short distances.
 *
 * Can be off by as much as 10%.
 *
 * approx_distance_in_mi = sqrt(x * x + y * y)
 *
 * where x = 69.1 * (lat2 - lat1)
 * and y = 69.1 * (lon2 - lon1) * cos(lat1/57.3)
 *//*----------------------------------------------------------------------*/
double    ApproximateDisatanceBetweenTwoLatLonsInKm(
                  double lat1, double lon1,
                  double lat2, double lon2
                  ) {
    double  ldRadians, ldCosR, x, y;

    ldRadians = (lat1 / 57.3) * 0.017453292519943295769236907684886;
    ldCosR = cos(ldRadians);
    x = 69.1 * (lat2 - lat1);
    y = 69.1 * (lon2 - lon1) * ldCosR;

    return sqrt(x * x + y * y) * 1.609344;  /* Converts mi to km. */
}

크레딧은 https://github.com/kristianmandrup/geo_vectors/blob/master/Distance%20calc%20notes.txt이동 합니다.

당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.