当你知道经纬度时,在Java中计算距离的单位是米。

77

可能是重复内容:
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个小数。

我的坐标是用第三种方式给出的,因此首选该值的代码 :)


我成功地使用了这个网站上的代码:* 经纬度计算 我需要对答案进行一些计算,以将其转换为欧洲米,但这很简单 :) - Espen Herseth Halvorsen
2
http://en.wikipedia.org/wiki/Haversine_formula - dotjoe
如果您理解数学公式,那么这个页面应该足以在任何编程语言中计算两点之间的距离。http://en.wikipedia.org/wiki/Geographical_distance - jan
只是出于好奇。由于地球不是一个完美的球体(赤道半径为6,378.137公里,极半径为6,356.752公里),误差会有多大?下面的答案中的3958.75英里是使用体积半径的近似值。 - KarlP
2
如果你关心地球不完全是一个球体,可以查看Vicenty公式--http://www.movable-type.co.uk/scripts/latlong-vincenty.html。 - mob
显示剩余3条评论
3个回答

193

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

17
为什么要将数据类型转换为浮点数,然后再转回浮点数? - Steve Kuo
3
返回 (float) (dist * meterConversion)。其中,meterConversion 是一个常数,表示从某个单位到米的转换因子,而 dist 则是需要进行转换的距离值。 - mob
6
地球的半径是3958.75英里,而不是公里。因此,您上面的代码返回英里。 - swinefeaster
5
倒数第二行有一个meterConversion常量,它可以将英里转换为米。 - Keith Irwin
9
请使用6369作为地球半径而不是英里,然后在转换为公制单位之后进行计算。 - Eric Darchis
显示剩余12条评论

27

你可以使用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

得出的距离单位为米。


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

“a” 可以是负数吗? - Xi Wei
@XiWei No. 纬度在-pi/2和pi/2之间,因此它们的余弦值为正数。而其他因素则是平方的。 - Dawood ibn Kareem

网页内容由stack overflow 提供, 点击上面的
可以查看英文原文,
原文链接