दो अक्षांश और देशांतर GeoCoordinates के बीच की दूरी की गणना


139

मैं दो GeoCoordinates के बीच की दूरी की गणना कर रहा हूं। मैं 3-4 अन्य ऐप्स के खिलाफ अपने ऐप का परीक्षण कर रहा हूं। जब मैं दूरी की गणना कर रहा होता हूं, तो मुझे अपनी गणना के लिए औसतन ३.३ मील की दूरी तय करनी पड़ती है, जबकि अन्य एप्लिकेशन ३.५ मील की दूरी पर होते हैं। यह उस गणना के लिए एक बड़ा अंतर है जिसे मैं प्रदर्शन करने की कोशिश कर रहा हूं। दूरी की गणना के लिए क्या कोई अच्छी श्रेणी की लाइब्रेरी हैं? मैं इसे 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;
}

मुझ से ऐसी कौनसी गलती हो जाएगी? क्या मुझे इसे पहले किमी में गणना करना चाहिए और फिर मीलों में बदलना चाहिए?


1
पृथ्वी का अर्थ त्रिज्या = 6,371 किमी = 3958.76 मील
मिच गेहूं


यह gis.stackexchange.com पर नहीं होना चाहिए
डैनियल पॉवेल

यह हो सकता है, लेकिन मेरा प्रश्न एक विंडोज फोन पर यह गणना करने के लिए अधिक संबंधित है जो थोड़ा अलग है। सूत्र समान है, लेकिन नए तरीके से कॉल करते हैं जैसे कि दूरस्थ विधि विधि आवश्यक रूप से उपलब्ध नहीं है।
जेसन एन। गेलॉर्ड

1
सुझाव दें कि आप पी / 180 स्टोर करें ताकि आपको गणना दोहराते न रहें।
बजे क्रिस कैविनेस

जवाबों:


313

GeoCoordinate वर्ग (.NET फ्रेमवर्क 4 और उच्च) पहले से ही GetDistanceToविधि।

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

return sCoord.GetDistanceTo(eCoord);

दूरी मीटरों में है।

आपको System.Device को रेफर करना होगा।


निगेल, क्या आप सुनिश्चित हैं कि फोन पर डिस्टेंटो विधि काम करेगी? मैंने सोचा कि यह WP7 के लिए जियो कॉर्डिनेट के 2.0 संस्करण का उपयोग करता है।
जेसन एन। गेलॉर्ड

1
मैंने इसकी जाँच की और डिवाइस के लिए जियोकॉर्डिनेट में एक गेटडांसटॉ विधि है जिसे आपने संदर्भित किया था (लेकिन आपके पास जो नहीं है)। कोई बड़ी बात नहीं। मैं इसे देखने के लिए परीक्षण करने जा रहा हूं कि क्या गणना में निर्मित कोई बेहतर है। धन्यवाद निगेल!
जेसन एन। गेलॉर्ड

1
मैं एक गलत प्रश्न पूछ सकता हूं, लेकिन परिणाम किस इकाई में है? क्या यह माइल्स या किलोमीटर है। मेरे द्वारा इसे कहीं भी ढूंढा जा सकता है।
सईद नेमाटी

3
@SaeedNeamati - यह भी देख रहा था, msdn.microsoft.com/en-us/library/… के अनुसार - यह मीटर में है।
एंडी बटलैंड

हां, GeoCoordinate.GetDistanceTo () मीटर में मान लौटाता है। मेरे लिए, यूएसए में, अगर यह 1610 से कम है, तो मैं इसे पैरों (मीटर * 3.28084) में बदल देता हूं अन्यथा मैं मीलों (मीटर * 0.000621371) में बदल जाता हूं। सटीकता मेरे उद्देश्यों के लिए पर्याप्त से अधिक है।
user3235770

110

GetDistance सबसे अच्छा समाधान है , लेकिन कई मामलों में हम इस विधि का उपयोग नहीं कर सकते हैं (जैसे यूनिवर्सल ऐप)

  • स्यूडोकोड एल्गोरिथ्म के लिए दूरी की गणना coorindates के बीच:

    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;
    }
  • वास्तविक विश्व सी # कार्यान्वयन , जो एक विस्तार विधियों का उपयोग करता है

    उपयोग:

    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
क्या आप इस गणना के लिए उपयोग किए गए सूत्र प्रदान कर सकते हैं या हो सकता है कि कुछ टिप्पणियां इस बारे में बताएं कि कौन सी रेखा क्या करती है? क्या मुझे परिवर्तित करने के लिए मील के बजाय सीधे मील में परिणामी दूरी को बदलना होगा?
अल्बर्टोफेडज़एम

एक अच्छे समाधान के लिए धन्यवाद, अब मैं इसे अपने डेस्कटॉप एप्लिकेशन में उपयोग कर सकता हूं।
जमशेद कामरान

मेरे UWP ऐप में बढ़िया काम किया गया जहाँ मैं जियो कॉर्डिनेट का उपयोग नहीं कर सकता।
जच ग्रीन

1
गणना 95% सत्य है। नीचे का कार्य 100% सटीक है: stackoverflow.com/a/51839058/3736063
मालेक ट्यूबैसैट

31

और यहाँ, उन लोगों के लिए जो अभी भी संतुष्ट नहीं हैं, .NET फ्रेमवर्क 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

धन्यवाद! मैं जियो कॉर्डिनेट क्लास में इस्तेमाल होने वाली पृथ्वी की वास्तविक त्रिज्या की तलाश में था।
KRoy 12

माइनर ऑप्टिमाइज़ेशन, या आसान पढ़ने के लिए, प्री-कंप्यूट पी / 180 कर सकता था double oneDegree = Math.PI / 180.0;?
ब्रेकरू

1
@brakeroo आपके उत्तर के लिए धन्यवाद। मैं इसका उत्तर छोड़ना चाहूंगा क्योंकि यह मूल .NET कोड है। किसी को भी अपने सुझाव का पालन करने के लिए स्वतंत्र महसूस कर सकते हैं, ज़ाहिर है।
मार्क

17

यहाँ जावास्क्रिप्ट संस्करण लोग और gals है

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 का उपयोग कर रहे हैं और जिओ कॉर्डिनेट वर्ग तक उनकी पहुंच नहीं है, आप उनके स्थान पर Android स्थान वर्ग का उपयोग कर सकते हैं:

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"));

अच्छी तरह से काम! धन्यवाद!
श्नापज

3

इन प्लेटफ़ॉर्म के लिए यह लाइब्रेरी GeoCoordinate है:

  • मोनो
  • .NET 4.5
  • .NET कोर
  • विंडोज फोन 8.x
  • यूनिवर्सल विंडोज प्लेटफॉर्म
  • ज़मीरिन आईओएस
  • Xamarin Android

स्थापना NuGet के माध्यम से की जाती है:

PM> इंस्टॉल-पैकेज जियो कॉर्डिनेट

प्रयोग

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

double distanceBetween = pin1.GetDistanceTo(pin2);

दो निर्देशांक के बीच की दूरी, मीटर में


3

इलियट वुड के फ़ंक्शन के आधार पर, और यदि किसी को सी फ़ंक्शन में रुचि है, तो यह काम कर रहा है ...

#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);
}

आप इसे डबल में बदल सकते हैं । यह किमी में मान लौटाता है।


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 # वैरिएंट (किमी में दूरी, बिना वैरिएबल और निरर्थक गणनाओं के, हैवेरसिन फॉर्मेलर 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
    }
}

विकिपीडिया से हैवरसाइन फॉर्मुलर


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

जब सीपीयू / गणित कंप्यूटिंग शक्ति सीमित है:

ऐसे समय होते हैं (जैसे कि मेरे काम में) जब कंप्यूटिंग शक्ति कम होती है (जैसे कोई फ्लोटिंग पॉइंट प्रोसेसर नहीं, छोटे माइक्रोकंट्रोलर के साथ काम करना) जहां कुछ ट्रिगर फ़ंक्शन सीपीयू समय (जैसे कि 3000+ घड़ी चक्र) की अत्यधिक मात्रा में ले सकते हैं, इसलिए मैं केवल एक सन्निकटन की आवश्यकता है, खासकर यदि सीपीयू को लंबे समय तक बांधना नहीं चाहिए, तो मैं सीपीयू ओवरहेड को कम करने के लिए इसका उपयोग करता हूं:

/**------------------------------------------------------------------------
 * \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 को जाता है ।

हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.