当你知道java中的经度和纬度时,计算以米为单位的距离

时间:2009-05-08 01:38:08

标签: java distance latitude-longitude

  

可能重复:
  Working with latitude/longitude values in 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 个答案:

答案 0 :(得分:177)

基于another question on 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;
    }

答案 1 :(得分:24)

您可以使用Java Geodesy Library for GPS,它使用考虑了地球表面曲率的Vincenty's formulae

实施如下:

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

结果距离以米为单位。

答案 2 :(得分:4)

在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;
}