根据两点经纬度坐标计算距离

在开发中用到定位和地图功能的都会知道位置信息会有多个坐标系来表示，而不同地图厂商使用的坐标系也可能是不一样的，web、Android以及iOS可能使用了不同厂家的地图，这就出现了在不同坐标系间转换坐标的需求。 然而一些地图SDK给出的地图坐标转换接口的实现是在服务器进行的，那么这些接口的调用不但会发起网络请求而效率上也会非常的慢，无论是对于客户端离线操作还是在自己的服务后台批量转换坐标数据，都需要本地实现坐标转换的算法来满足快速运算的需求。

GPS坐标比较

坐标系	解析	使用地图
WGS84	地球坐标系，国际上通用的坐标系。设备一般包含GPS芯片或者北斗芯片获取的经纬度为WGS84地理坐标系,最基础的坐标，谷歌地图在非中国地区使用的坐标系	GPS/谷歌地图卫星
GCJ02	火星坐标系，是由中国国家测绘局制订的地理信息系统的坐标系统。并要求在中国使用的地图产品使用的都必须是加密后的坐标，而这套WGS84加密后的坐标就是gcj02。	腾讯(搜搜)地图，阿里云地图，高德地图，谷歌国内地图
BD09	百度坐标系，百度在GCJ02的基础上进行了二次加密，官方解释是为了进一步保护用户隐私	百度地图
小众坐标系	类似于百度地图，在GCJ02基础上使用自己的加密算法进行二次加密的坐标系	搜狗地图、图吧地图 等

各种坐标互转计算函数

未 测 试 根据反余弦计算直线距离

反余弦分析: https://blog.csdn.net/jk940438163/article/details/83147557#commentsedit

public class DistanceUtil {

    private static final double EARTH_RADIUS = 6371393; // 平均半径,单位：m；不是赤道半径。赤道为6378左右
 
    /**
     * @描述 反余弦进行计算
     * @参数  [lat1, lng1, lat2, lng2]
     * @返回值  double
     **/
    public static double getDistance(Double lat1,Double lng1,Double lat2,Double lng2) {
        // 经纬度（角度）转弧度。弧度用作参数，以调用Math.cos和Math.sin
        double radiansAX = Math.toRadians(lng1); // A经弧度
        double radiansAY = Math.toRadians(lat1); // A纬弧度
        double radiansBX = Math.toRadians(lng2); // B经弧度
        double radiansBY = Math.toRadians(lat2); // B纬弧度
 
        // 公式中“cosβ1cosβ2cos（α1-α2）+sinβ1sinβ2”的部分，得到∠AOB的cos值
        double cos = Math.cos(radiansAY) * Math.cos(radiansBY) * Math.cos(radiansAX - radiansBX)
            + Math.sin(radiansAY) * Math.sin(radiansBY);
//        System.out.println("cos = " + cos); // 值域[-1,1]
        double acos = Math.acos(cos); // 反余弦值
//        System.out.println("acos = " + acos); // 值域[0,π]
//        System.out.println("∠AOB = " + Math.toDegrees(acos)); // 球心角 值域[0,180]
        return EARTH_RADIUS * acos; // 最终结果
    }
    public static void main(String[] args) {
        // 121.717594,31.12055    121.817629,31.090867
        double distance = getDistance(31.12055, 121.717594, 31.090867, 121.817629);
        System.out.println("距离" + distance  + "米");
   	}
}

推荐方式 国内Gcj02坐标系两点之间的直线距离（高德地图）

根据sin和cosin计算距离误差与高度基本一样，参考博客: https://blog.csdn.net/weixin_41804049/article/details/107537320

import java.awt.*;
import java.math.BigDecimal;
public class DistanceUtil {

    private final static double EARTH_RADIUS = 6378.137;//地球半径

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


    /**
     * 计算国内坐标系两点间距离 （高德地图）
     *
     * @return double 距离 单位公里,精确到米
     */
    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 = new BigDecimal(s).setScale(3, BigDecimal.ROUND_HALF_UP).doubleValue();
        return s;
    }

    //(116.368904, 39.923423)   (116.387271, 39.772501)
    public static void main(String[] args) {
        double lon1 = 121.464636;
        double lat1 = 29.792173;
        double lon2 = 122.462831;
        double lat2 = 29.792624;

        double dist = GetDistance(lat1, lon1, lat2, lon2);
        System.out.println("输出结果：" + dist * 1000);
    }
}

在线验证API(高德)

存在误差 使用第三方jar包计算距离

import org.gavaghan.geodesy.Ellipsoid;
import org.gavaghan.geodesy.GeodeticCalculator;
import org.gavaghan.geodesy.GeodeticCurve;
import org.gavaghan.geodesy.GlobalCoordinates;
public class DistanceUtil {

    public static  double getDistanceMeter(GlobalCoordinates gpsFrom, GlobalCoordinates gpsTo, Ellipsoid ellipsoid){
	    // 创建GeodeticCalculator,调用计算方法,传入坐标系,经纬度用于计算距离
	    GeodeticCurve geoCurve = new GeodeticCalculator().calculateGeodeticCurve(ellipsoid, gpsFrom, gpsTo);
    	return geoCurve.getEllipsoidalDistance();
	}

    // (116.368904, 39.923423)   (116.387271, 39.772501)
    public static void main(String[] args) {
        // 给定两个坐标系,计算两点相差距离
	    GlobalCoordinates source = new GlobalCoordinates(29.490295, 106.486654);
	    GlobalCoordinates target = new GlobalCoordinates(29.615467, 106.581515);
	    // Sphere坐标的计算结果
	    double meter1 =getDistanceMeter(source,target,Ellipsoid.Sphere);
	    // WGS84坐标系计算结果
	    double meter2 = getDistanceMeter(source,target,Ellipsoid.WGS84);
	
	    // 计算结果Sphere 坐标系的计算结果与 WGS84坐标系的计算结果存在几十米的误差，不同的坐标系精度不同,
	    System.out.println("Sphere 坐标系的计算结果是" + meter1 + "M");
	    System.out.println("wGS84  坐标系的计算结果是" + meter2 + "M");
    }
}

参考博客

java 根据经纬度计算实际距离
计算国内Gcj02坐标系两点之间的直线距离（高德地图）
WGS84、GCJ02、BD09地图坐标系间的坐标转换及坐标距离计算