一线召、計算兩個經(jīng)緯度的距離

根據(jù)經(jīng)緯度計算距離

• R為地球半徑，可取平均值 6371km多矮；
• φ1, φ2 表示兩點的緯度缓淹；
• Δλ 表示兩點經(jīng)度的差值。

    /**   
     * 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;    
        }    
        
        /**   
         * 通過經(jīng)緯度獲取距離(單位：米)   
         * @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;    
        }    
    }    

二塔逃、根據(jù)某一經(jīng)緯度和距離計算附近的經(jīng)緯度

• R為地球半徑讯壶，可取平均值 6371km；
• φ1, φ2 表示兩點的緯度患雏；
• Δλ 表示兩點經(jīng)度的差值鹏溯。

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);

• 確定四個點坐標（即兩個經(jīng)度和兩個緯度）

    // 坐標 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);

    // 通過計算可以得到上下左右四個點的經(jīng)緯度，即肤舞，兩個經(jīng)度紫新，兩個緯度
    // 坐標 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;

效果展示

效果展示