获取两个经纬度距离

/*** 距离工具类 用于计算两个金纬度距离* @author zhengcj* @description:* @date 2024/11/19 10:27*/
public class CalculationLogLatDistanceUtils {/*** 地球赤道半径(km)* */public final static double EARTH_RADIUS = 6378.137;/*** 转化为弧度(rad)* */public static double rad(double d) {return d * Math.PI / 180.0;}/*** 求两经纬度距离** @param lon1*            第一点的经度* @param lat1*            第一点的纬度* @param lon2*            第二点的经度* @param lat2*            第二点的纬度* @return 两点距离，单位km* */public static double GetDistanceOne(double lon1, double lat1, double lon2,double lat2) {double r1 = rad(lat1);double r2 = rad(lon1);double a = rad(lat2);double b = rad(lon2);double s = Math.acos(Math.cos(r1) * Math.cos(a) * Math.cos(r2 - b)+ Math.sin(r1) * Math.sin(a))* EARTH_RADIUS;return s;}/*** 求两经纬度距离(google maps源码中)** @param lon1*            第一点的经度* @param lat1*            第一点的纬度* @param lon2*            第二点的经度* @param lat2*            第二点的纬度* @return 两点距离，单位km* */public static double GetDistanceTwo(double lon1, double lat1, double lon2,double lat2) {double radLat1 = rad(lat1);double radLat2 = rad(lat2);double a = radLat1 - radLat2;double b = rad(lon1) - rad(lon2);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;return s;}/*** 求两经纬度距离** @param lon1*            第一点的经度* @param lat1*            第一点的纬度* @param lon2*            第二点的经度* @param lat2*            第二点的纬度* @return 两点距离，单位km* */public static double GetDistanceThree(double lon1, double lat1,double lon2, double lat2) {double radLat1 = rad(lat1);double radLat2 = rad(lat2);double radLon1 = rad(lon1);double radLon2 = rad(lon2);if (radLat1 < 0)radLat1 = Math.PI / 2 + Math.abs(radLat1);// southif (radLat1 > 0)radLat1 = Math.PI / 2 - Math.abs(radLat1);// northif (radLon1 < 0)radLon1 = Math.PI * 2 - Math.abs(radLon1);// westif (radLat2 < 0)radLat2 = Math.PI / 2 + Math.abs(radLat2);// southif (radLat2 > 0)radLat2 = Math.PI / 2 - Math.abs(radLat2);// northif (radLon2 < 0)radLon2 = Math.PI * 2 - Math.abs(radLon2);// westdouble x1 = Math.cos(radLon1) * Math.sin(radLat1);double y1 = Math.sin(radLon1) * Math.sin(radLat1);double z1 = Math.cos(radLat1);double x2 = Math.cos(radLon2) * Math.sin(radLat2);double y2 = Math.sin(radLon2) * Math.sin(radLat2);double z2 = Math.cos(radLat2);double d = Math.pow((x1 - x2), 2) + Math.pow((y1 - y2), 2)+ Math.pow((z1 - z2), 2);// // 余弦定理求夹角// double theta = Math.acos((2 - d) / 2);
d = Math.pow(EARTH_RADIUS, 2) * d;// //余弦定理求夹角double theta = Math.acos((2 * Math.pow(EARTH_RADIUS, 2) - d)/ (2 * Math.pow(EARTH_RADIUS, 2)));double dist = theta * EARTH_RADIUS;return dist;}}