一、计算两个经纬度的距离

根据经纬度计算距离

• R为地球半径，可取平均值 6371km；
• φ1, φ2 表示两点的纬度；
• Δλ 表示两点经度的差值。

    /**   
     * Created by yuliang on 2015/3/20.   
     */    
    public class LocationUtils {    
        private static double EARTH_RADIUS = 6378.137;    
        
        private static double rad(double d) {    
            return d * Math.PI / 180.0;    
        }    
        
        /**   
         * 通过经纬度获取距离(单位：米)   
         * @param lat1   
         * @param lng1   
         * @param lat2   
         * @param lng2   
         * @return   
         */    
        public static double getDistance(double lat1, double lng1, double lat2,    
                                         double lng2) {    
            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;    
            s = s*1000;    
            return s;    
        }    
    }    

二、根据某一经纬度和距离计算附近的经纬度

• R为地球半径，可取平均值 6371km；
• φ1, φ2 表示两点的纬度；
• Δλ 表示两点经度的差值。

1、计算步骤：

• 先来求东西两侧的的范围边界。在haversin公式中令φ1 = φ2，可得：

计算公式

double dlng = 2 * Math.asin(Math.sin(gpsAround.getDistance() / (2 * EARTH_RADIUS)) / Math.cos(gpsAround.getLatitude().doubleValue()));
// 弧度转换成角度
dlng = Math.toRadians(dlng);

• 然后求南北两侧的范围边界，在haversin公式中令 Δλ = 0，可得:

公式

    // 然后求南北两侧的范围边界，在haversin公式中令 Δλ = 0
    double dlat = gpsAround.getDistance() / EARTH_RADIUS;
    // 弧度转换成角度
    dlng = Math.toDegrees(dlat);

• 确定四个点坐标（即两个经度和两个纬度）

    // 坐标 10
    GPS left_bottom = new ShopAddress().new GPS();
    left_bottom.setLatitude(gpsAround.getLatitude().subtract(new BigDecimal(dlat)));
    left_bottom.setLongitude(gpsAround.getLongitude().subtract(new BigDecimal(dlng)));

    // 坐标 01
    GPS right_top = new ShopAddress().new GPS();
    right_top.setLatitude(gpsAround.getLatitude().add(new BigDecimal(dlat)));
    right_top.setLongitude(gpsAround.getLongitude().add(new BigDecimal(dlng)));

2、用例

private final static double EARTH_RADIUS = 6371.393;// 地球半径

  private static Map<String, GPS> getNearbyGCJ(GPSAround gpsAround) {
    if (null == gpsAround) {
      return null;
    }

    // 求东西两侧的的范围边界。在haversin公式中令φ1 = φ2(维度相同),传入的距离为 米，转换为千米
    double dlng = 2 * Math.asin(Math.sin((gpsAround.getDistance() / 1000) / (2 * EARTH_RADIUS))
        / Math.cos(gpsAround.getLatitude().doubleValue() * Math.PI / 180));
    // 弧度转换成角度
    // dlng = Math.toRadians(dlng);
    dlng = Math.toDegrees(dlng);

    // 然后求南北两侧的范围边界，在haversin公式中令 Δλ = 0
    double dlat = (gpsAround.getDistance() / 1000) / EARTH_RADIUS;
    // 弧度转换成角度
    dlat = Math.toDegrees(dlat);

    // 通过计算可以得到上下左右四个点的经纬度，即，两个经度，两个纬度
    // 坐标 10
    GPS left_bottom = new GPS(
        gpsAround.getLongitude().subtract(new BigDecimal(dlng)).setScale(8,
            BigDecimal.ROUND_HALF_UP),
        gpsAround.getLatitude().subtract(new BigDecimal(dlat)).setScale(8,
            BigDecimal.ROUND_HALF_UP));

    // 坐标 01
    GPS right_top = new GPS(
        gpsAround.getLongitude().add(new BigDecimal(dlng)).setScale(8, BigDecimal.ROUND_HALF_UP),
        gpsAround.getLatitude().add(new BigDecimal(dlat)).setScale(8, BigDecimal.ROUND_HALF_UP));

    Map<String, GPS> gpsMap = new HashMap<>();
    gpsMap.put("top", right_top);
    gpsMap.put("bottom", left_bottom);

    return gpsMap;
  }

三、mysql计算逻辑（性能较差）

SELECT
    shop_id ,
    shop_name ,
    lng ,
    lat ,
    POWER(lat - 40.05748 , 2) + POWER(lng - 116.30759 , 2) * POWER(COS((lat + 40.05748) / 2) , 2) AS distance
FROM
    shop_list
HAVING
    distance < 1000
ORDER BY
    distance
LIMIT 100;

效果展示

效果展示