mysql根据经纬度查询范围内数据，并根据距离排序

最近接到一个新需求，要根据经纬度查询指定范围内的数据，并且根据距离进行排序，网上找了找相关实现方法，记录一下。

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等相关函数来支持经纬度相关的计算。

SELECTs.*, (st_distance(point(lng, lat), point(-122.083235, 37.38714) ) * 111195) AS distance
FROMmarkers s
ORDER BYdistance

搭配下面根据距离计算经纬度范围代码，拼接上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;}