根据给定点的经纬度和特定距离取范围类的数据

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

根据经纬度计算距离
  • 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;

效果展示

效果展示

推荐阅读更多精彩内容

  • 这些经纬线是怎样定出来的呢?地球是在不停地绕地轴旋转(地轴是一根通过地球南北两极和地球中心的假想线),在地球中腰画...
    Akitas阅读 305评论 0 2
  • 这些经纬线是怎样定出来的呢?地球是在不停地绕地轴旋转(地轴是一根通过地球南北两极和地球中心的假想线),在地球中腰画...
    DragonersLi阅读 1,089评论 0 7
  • 1 序: 很多新接触GIS的人员对地图投影以及坐标系统很难理解,甚至做GIS开发做了好几年的人也有这方面的疑惑,地...
    三维GIS那点事_王跃军阅读 8,018评论 4 38
  • 知识科普 这些经纬线是怎样定出来的呢?地球是在不停地绕地轴旋转(地轴是一根通过地球南北两极和地球中心的 假想线),...
    Yumazhiyao阅读 1,052评论 1 3
  • 虽然是个很热的秋天,但也是个丰收的季节. 画点水果……是的……画完之后我都给吃了,不浪费 #Hello October#
    艾米札记阅读 28评论 2 0