一、计算两个经纬度的距离
![](https://upload-images.jianshu.io/upload_images/1802116-a5a7bb51b0c607fb.png?imageMogr2/auto-orient/strip|imageView2/2/w/441/format/webp)
根据经纬度计算距离
- 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;
效果展示
效果展示
作者:慕凌峰
链接:https://www.jianshu.com/p/ed6ea376911e
来源:简书
著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。