java根据GPS经纬度坐标计算两点的距离算法

在开发项目中用到根据两点计算距离的算法,找了很多。不过很多是存在问题的。现将其整理了一下。供大家参考:


/**
 * Created by yuliang on 2015/3/20.
 */
public class LocationUtils {
    private static double EARTH_RADIUS = 6378.137;

    private static double rad(double d) {
        return d * Math.PI / 180.0;
    }

    /**
     * 通过经纬度获取距离(单位:米)
     * @param lat1
     * @param lng1
     * @param lat2
     * @param lng2
     * @return
     */
    public static double getDistance(double lat1, double lng1, double lat2,
                                     double lng2) {
        double radLat1 = rad(lat1);
        double radLat2 = rad(lat2);
        double a = radLat1 - radLat2;
        double b = rad(lng1) - rad(lng2);
        double s = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(a / 2), 2)
                + Math.cos(radLat1) * Math.cos(radLat2)
                * Math.pow(Math.sin(b / 2), 2)));
        s = s * EARTH_RADIUS;
        s = Math.round(s * 10000d) / 10000d;
        s = s*1000;
        return s;
    }
}

这样的计算方式得到的距离并非真实的距离,可以说是逻辑距离,但其距离也已经很准确。不过毕竟是通过逻辑计算得到的距离,若要求高准确性的距离信息的话,还是借助第三方的地图api接口获取比较合适。

  • 2
    点赞
  • 16
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Vincenty算法是一种用于计算地球上两点之间的距离算法,它考虑了地球的椭球形状和旋转。以下是使用Java实现Vincenty算法计算两点之间距离的代码示例: ```java public class VincentyDistanceCalculator { private static final double EQUATORIAL_RADIUS = 6378137.0; // 地球赤道半径,单位: private static final double POLAR_RADIUS = 6356752.3; // 地球极半径,单位: public static double calculateDistance(double lat1, double lon1, double lat2, double lon2) { double a = EQUATORIAL_RADIUS; double b = POLAR_RADIUS; double f = (a - b) / a; double L = Math.toRadians(lon2 - lon1); double U1 = Math.atan((1 - f) * Math.tan(Math.toRadians(lat1))); double U2 = Math.atan((1 - f) * Math.tan(Math.toRadians(lat2))); double sinU1 = Math.sin(U1); double cosU1 = Math.cos(U1); double sinU2 = Math.sin(U2); double cosU2 = Math.cos(U2); double lambda = L; double lambdaP = 2 * Math.PI; int iterLimit = 100; while (Math.abs(lambda - lambdaP) > 1e-12 && iterLimit-- > 0) { double sinLambda = Math.sin(lambda); double cosLambda = Math.cos(lambda); double sinSigma = Math.sqrt((cosU2 * sinLambda) * (cosU2 * sinLambda) + (cosU1 * sinU2 - sinU1 * cosU2 * cosLambda) * (cosU1 * sinU2 - sinU1 * cosU2 * cosLambda)); if (sinSigma == 0) { return 0.0; // 两点重合 } double cosSigma = sinU1 * sinU2 + cosU1 * cosU2 * cosLambda; double sigma = Math.atan2(sinSigma, cosSigma); double sinAlpha = cosU1 * cosU2 * sinLambda / sinSigma; double cosSqAlpha = 1 - sinAlpha * sinAlpha; double cos2SigmaM = cosSigma - 2 * sinU1 * sinU2 / cosSqAlpha; double C = f / 16 * cosSqAlpha * (4 + f * (4 - 3 * cosSqAlpha)); lambdaP = lambda; lambda = L + (1 - C) * f * sinAlpha * (sigma + C * sinSigma * (cos2SigmaM + C * cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM))); } if (iterLimit == 0) { throw new RuntimeException("过程未收敛"); } double uSq = cosSqAlpha * (a * a - b * b) / (b * b); double A = 1 + uSq / 16384 * (4096 + uSq * (-768 + uSq * (320 - 175 * uSq))); double B = uSq / 1024 * (256 + uSq * (-128 + uSq * (74 - 47 * uSq))); double deltaSigma = B * sinSigma * (cos2SigmaM + B / 4 * (cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM) - B / 6 * cos2SigmaM * (-3 + 4 * sinSigma * sinSigma) * (-3 + 4 * cos2SigmaM * cos2SigmaM))); double s = b * A * (sigma - deltaSigma); return s; } } ``` 该算法中的参数lat1、lon1、lat2和lon2分别表示两点的纬度和经度,单位为度。函数返回值表示两点之间的距离,单位为。需要注意的是,该算法中的距离是基于地球表面的,不考虑地球内部的变化。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值