java - 当您在java中知道经度和纬度时，以米为单位计算距离

Question

可能重复：
在 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 个十进制数字。

这是我的坐标给出的第三种方式，所以这个值的代码将是首选的:)

Answer 1

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

Answer 2

您可以使用Java Geodesy Library for GPS，它使用考虑了地球表面曲率的文森蒂公式。

实现是这样的：

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

结果距离以米为单位。

Answer 3

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