private static final double DEF_PI = 3.14159265359;

private static final double DEF_2PI = 6.28318530712;

private static final double DEF_PI180 = 0.01745329252;

private static final double DEF_R = 6370693.5; // 地球半径

/**
 * 利用勾股定理计算，适用两点距离很近的情况
 *
 * @param lon1 第一个坐标经度
 * @param lat1 第一个坐标纬度
 * @param lon2 第二个坐标经度
 * @param lat2 第二个坐标纬度
 * @return
 */
public static double getShortDistance(double lon1, double lat1, double lon2, double lat2) {
    double ew1;
    double ns1;
    double ew2;
    double ns2;

    double dx;
    double dy;
    double dew;

    double distance;
    // 角度转换为弧度
    ew1 = lon1 * DEF_PI180;
    ns1 = lat1 * DEF_PI180;
    ew2 = lon2 * DEF_PI180;
    ns2 = lat2 * DEF_PI180;
    // 经度差
    dew = ew1 - ew2;
    // 若跨东经和西经180度，进行调整
    if (dew > DEF_PI) {
        dew = DEF_2PI - dew;
    } else if (dew < -DEF_PI) {
        dew = DEF_2PI + dew;
    }
    dx = DEF_R * Math.cos(ns1) * dew; // 东西方向长度（在纬度圆上的投影长度）
    dy = DEF_R * (ns1 - ns2); // 南北方向长度（在经度圆上的投影长度）
    // 勾股定理求你斜边长
    distance = Math.sqrt(dx * dx + dy * dy);
    return distance;
}

/**
 * 按标准的球面大圆劣弧长度计算，适用距离较远的情况
 *
 * @param lon1 第一个坐标经度
 * @param lat1 第一个坐标纬度
 * @param lon2 第二个坐标经度
 * @param lat2 第二个坐标纬度
 * @return
 */
public static double getLongDistance(double lon1, double lat1, double lon2, double lat2) {
    double ew1;
    double ns1;
    double ew2;
    double ns2;
    double distance;
    // 角度转换为弧度
    ew1 = lon1 * DEF_PI180;
    ns1 = lat1 * DEF_PI180;
    ew2 = lon2 * DEF_PI180;
    ns2 = lat2 * DEF_PI180;
    // 求大圆劣弧与球心所夹的角（弧度）
    distance = Math.sin(ns1) * Math.sin(ns2) + Math.cos(ns1) * Math.cos(ns2) * Math.cos(ew1 - ew2);
    // 调整到[-1..1]范围内，避免溢出
    if (distance > 1.0) {
        distance = 1.0;
    } else if (distance < -1.0) {
        distance = -1.0;
    }
    // 求大圆劣弧长度
    distance = DEF_R * Math.acos(distance);
    return distance;
}

关注公众号“大模型全栈程序员”回复“小程序”获取1000个小程序打包源码。更多免费资源在http://www.gitweixin.com/?p=2627