开发中经常会遇到计算两个点(经纬度)之间的距离或者计算最近门店的场景,下面简单实现一下如何计算两个经纬度之间相隔的距离。
方法一:工具类实现
1、导入geodesy的maven依赖 或者到阿里云maven仓库下载jar包
<dependency>
<groupId>org.gavaghan</groupId>
<artifactId>geodesy</artifactId>
<version>1.1.3</version>
</dependency>
2、实现计算
import org.gavaghan.geodesy.Ellipsoid;
import org.gavaghan.geodesy.GeodeticCalculator;
import org.gavaghan.geodesy.GeodeticCurve;
import org.gavaghan.geodesy.GlobalCoordinates;
public class CaculateDistanceTest
{
public static void main(String[] args)
{
GlobalCoordinates source =new GlobalCoordinates(29.490295, 106.486654);
GlobalCoordinates target =new GlobalCoordinates(29.615467, 106.581515);
double meter1 = getDistanceMeter(source, target, Ellipsoid.Sphere);
double meter2 = getDistanceMeter(source, target, Ellipsoid.WGS84);
System.out.println("Sphere坐标系计算结果:"+meter1 +"米");
System.out.println("WGS84坐标系计算结果:"+meter2 +"米");
}
public static double getDistanceMeter(GlobalCoordinates gpsFrom, GlobalCoordinates gpsTo, Ellipsoid ellipsoid)
{
//创建GeodeticCalculator,调用计算方法,传入坐标系、经纬度用于计算距离
GeodeticCurve geoCurve =new GeodeticCalculator().calculateGeodeticCurve(ellipsoid, gpsFrom, gpsTo);
return geoCurve.getEllipsoidalDistance();
}
}
3、输出结果:
对比百度地图,计算结果和Sphere坐标系计算结果一致,表明计算结果正确,WGS84坐标系的计算结果存在几十米的误差。不同的坐标系精度不同,计算结果不一样。大家根据实际情况自己选择。
方法二:自己实现
/**
* 地图工具类
*
* @author 1000112682
* @version MapUtil 1.0 2021/4/26 10:51
*/
public class MapUtil {
/**
* 地球半径
*/
private static double EARTH_RADIUS = 6378138.0;
private static double rad(double d) {
return d * Math.PI / 180.0;
}
/**
* 判断一个点是否在圆形区域内
*
* @param circleLng 圆心经度
* @param circleLat 圆心纬度
* @param coordLng 坐标经度
* @param coordLat 坐标纬度
* @param radius 半径
*/
public static boolean isInCircle(double circleLng, double circleLat, double coordLng,
double coordLat, double radius) {
return getDistance(circleLng, circleLat, coordLng, coordLat) <= radius;
}
/**
* 通过经纬度获取距离(单位:米)
*
* @param lng1 经度1
* @param lat1 纬度1
* @param lng2 经度2
* @param lat2 纬度2
* @return 距离
*/
public static double getDistance(double lng1, double lat1, double lng2, double lat2){
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;
return s;
}
public static void main(String[] args) {
System.out.println(getDistance(106.486654,29.490295,106.581515,29.615467));
}
}