最近接到一个新需求,要根据经纬度查询指定范围内的数据,并且根据距离进行排序,网上找了找相关实现方法,记录一下。
java代码如下
public final class DistanceUtils {
/**
* 地球半径,单位 km
*/
private static final double EARTH_RADIUS = 6378.137;
/**
* 根据经纬度,计算两点间的距离
*
* @param longitude1 第一个点的经度
* @param latitude1 第一个点的纬度
* @param longitude2 第二个点的经度
* @param latitude2 第二个点的纬度
* @return 返回距离 单位千米
*/
public static double getDistance(double longitude1, double latitude1, double longitude2, double latitude2) {
// 纬度
double lat1 = Math.toRadians(latitude1);
double lat2 = Math.toRadians(latitude2);
// 经度
double lng1 = Math.toRadians(longitude1);
double lng2 = Math.toRadians(longitude2);
// 纬度之差
double a = lat1 - lat2;
// 经度之差
double b = lng1 - lng2;
// 计算两点距离的公式
double s = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(a / 2), 2) +
Math.cos(lat1) * Math.cos(lat2) * Math.pow(Math.sin(b / 2), 2)));
// 弧长乘地球半径, 返回单位: 千米
s = s * EARTH_RADIUS;
return s;
}
public static void main(String[] args) {
double d = getDistance(116.308479, 39.983171, 116.353454, 39.996059);
System.out.println(d);
}
}
返回的单位是千米,如果想返回米,可以修改地球半径的单位从千米到米,并且由于该结果是double类型的,所以还可以借助Math.round方法进行四舍五入为long类型,然后精确到米:
// ......
// 弧长乘地球半径(6378137)
s = s * EARTH_RADIUS;
// 返回类型: long,单位: 米
return Math.round(s * 10000) / 10000;
MySQL实现,我是采用的这种方法,实现起来更简单
测试表结构如下
CREATE TABLE `markers` (
`id` INT NOT NULL AUTO_INCREMENT PRIMARY KEY ,
`name` VARCHAR( 60 ) NOT NULL ,
`address` VARCHAR( 80 ) NOT NULL ,
`lat` FLOAT( 10, 6 ) NOT NULL ,
`lng` FLOAT( 10, 6 ) NOT NULL
) ENGINE = MYISAM ;
当然存储引擎可以是InnoDB。
需要用到st_distance函数
MySQL其实在很早就提供了这种存储经纬度及相关运算的功能,这种数据类型叫做空间数据类型,而对应的索引被称为空间索引,但由于MySQL之前的版本对InnoDB支持的并不是太好,所以使用的并不多。不过MySQL5.6和MySQL5.7对此进行了优化,添加了st_distance等相关函数来支持经纬度相关的计算。
SELECT
s.*,
(st_distance(point(lng, lat), point(-122.083235, 37.38714) ) * 111195) AS distance
FROM
markers s
ORDER BY
distance
搭配下面根据距离计算经纬度范围代码,拼接上sql条件即可
/**
* 根据范围计算经纬度
*
* @param longitude 纬度
* @param latitude 经度
* @param dis 范围 单位为千米
* @return 结果
*/
public static Map<String, Double> calcLatitudeAndLongitude(double longitude, double latitude, double dis) {
double dlng = 2 * Math.asin(Math.sin(dis / (2 * EARTH_RADIUS)) / Math.cos(latitude * Math.PI / 180));
//角度转为弧度
dlng = dlng * 180 / Math.PI;
double dlat = dis / EARTH_RADIUS;
dlat = dlat * 180 / Math.PI;
//最小纬度
double minlat = latitude - dlat;
//最大纬度
double maxlat = latitude + dlat;
//最小经度
double minlng = longitude - dlng;
//最大经度
double maxlng = longitude + dlng;
HashMap<String, Double> hashMap = new HashMap<>();
hashMap.put("minlng", minlng);
hashMap.put("maxlng", maxlng);
hashMap.put("minlat", minlat);
hashMap.put("maxlat", maxlat);
return hashMap;
}
参考链接:
https://blog.csdn.net/qq_35071164/article/details/90602501
点击直达