Вычислите расстояние в метрах, когда вы знаете долготу и широту в java [дубликат]


Возможные Дубликаты:
работа со значениями широты/долготы в Java

дублирую:

Мне нужно вычислить расстояние между двумя точками, заданными двумя координатами. Проект, над которым я работаю это Java-проект, поэтому Java-код будет отличным, но псевдокод также может быть дан, тогда я могу реализовать его сам:)

Как вы, наверное, знаете, есть три способа представления координат:

  • Градусы: Минуты: Секунды (49°30 ' 00 "N, 123°30' 00 " W)
  • Градусы: десятичные минуты (49°30.0', -123°30.0'), (49d30. 0m, - 123d30.0')
  • десятичные Градусы (49.5000°, -123.5000°), как правило, с 4-6 десятичными числами.

Это третий способ мои координаты приведены, поэтому код для этого значения будет предпочтительным :)

3 73

3 ответа:

на основе еще один вопрос о stackoverflow, я получил этот код.. Это вычисляет результат в метрах, а не в километрах :)

 public static float distFrom(float lat1, float lng1, float lat2, float lng2) {
    double earthRadius = 6371000; //meters
    double dLat = Math.toRadians(lat2-lat1);
    double dLng = Math.toRadians(lng2-lng1);
    double a = Math.sin(dLat/2) * Math.sin(dLat/2) +
               Math.cos(Math.toRadians(lat1)) * Math.cos(Math.toRadians(lat2)) *
               Math.sin(dLng/2) * Math.sin(dLng/2);
    double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));
    float dist = (float) (earthRadius * c);

    return dist;
    }

можно использовать библиотека геодезии Java для GPS использует Vincenty по формулам, которая учитывает кривизну поверхности Земли.

реализация идет следующим образом:

import org.gavaghan.geodesy.*;

...

GeodeticCalculator geoCalc = new GeodeticCalculator();

Ellipsoid reference = Ellipsoid.WGS84;  

GlobalPosition pointA = new GlobalPosition(latitude, longitude, 0.0); // Point A

GlobalPosition userPos = new GlobalPosition(userLat, userLon, 0.0); // Point B

double distance = geoCalc.calculateGeodeticCurve(reference, userPos, pointA).getEllipsoidalDistance(); // Distance between Point A and Point B

полученное расстояние в метрах.

В C++ это делается так:

#define LOCAL_PI 3.1415926535897932385 

double ToRadians(double degrees) 
{
  double radians = degrees * LOCAL_PI / 180;
  return radians;
}

double DirectDistance(double lat1, double lng1, double lat2, double lng2) 
{
  double earthRadius = 3958.75;
  double dLat = ToRadians(lat2-lat1);
  double dLng = ToRadians(lng2-lng1);
  double a = sin(dLat/2) * sin(dLat/2) + 
             cos(ToRadians(lat1)) * cos(ToRadians(lat2)) * 
             sin(dLng/2) * sin(dLng/2);
  double c = 2 * atan2(sqrt(a), sqrt(1-a));
  double dist = earthRadius * c;
  double meterConversion = 1609.00;
  return dist * meterConversion;
}