根據一個給定經緯度的點和距離,搜索附近5公裡範圍的地點 –合理利用演算法 功能需求:需要查詢一個站點(已知該站點經緯度)5公裡範圍內的其它站點。 方 法 一 :對每條記錄,去進行遍歷,跟資料庫中的每一個點進行距離計算,當距離小於5公裡時,認為匹配(效率極其低下,耗時長)。 方 法 二 :先過濾出大概的 ...
根據一個給定經緯度的點和距離,搜索附近5公裡範圍的地點 –合理利用演算法
功能需求:需要查詢一個站點(已知該站點經緯度)5公裡範圍內的其它站點。
方 法 一 :對每條記錄,去進行遍歷,跟資料庫中的每一個點進行距離計算,當距離小於5公裡時,認為匹配(效率極其低下,耗時長)。
方 法 二 :先過濾出大概的經緯度範圍再進行計算。我們求取的雖然是個圓,但我們可以先求出該圓的外接正方形(比方說正方形的四個點),先算出該點周圍的矩形的四個點的經緯度,然後使用正方形的經緯度去直接匹配資料庫中的記錄,然後在用方法一進行篩選(相對於只使用方法一效率高了非常多 –合理利用演算法)。
public static class DistanceHelper { /// <summary> /// 根據一個給定經緯度的點和距離,進行附近地點查詢 /// </summary> /// <param name="longitude">經度</param> /// <param name="latitude">緯度</param> /// <param name="distance">距離(單位:公裡或千米)</param> /// <returns>返回一個範圍的4個點,最小緯度和緯度,最大經度和緯度</returns> public static PositionModel FindNeighPosition(double longitude, double latitude, double distance) { //先計算查詢點的經緯度範圍 double r = 6378.137;//地球半徑千米 double dis = distance;//千米距離 double dlng = 2 * Math.Asin(Math.Sin(dis / (2 * r)) / Math.Cos(latitude * Math.PI / 180)); dlng = dlng * 180 / Math.PI;//角度轉為弧度 double dlat = dis / r; dlat = dlat * 180 / Math.PI; double minlat = latitude - dlat; double maxlat = latitude + dlat; double minlng = longitude - dlng; double maxlng = longitude + dlng; return new PositionModel { MinLat = minlat, MaxLat = maxlat, MinLng = minlng, MaxLng = maxlng }; } /// <summary> /// 計算兩點位置的距離,返回兩點的距離,單位:公裡或千米 /// 該公式為GOOGLE提供,誤差小於0.2米 /// </summary> /// <param name="lat1">第一點緯度</param> /// <param name="lng1">第一點經度</param> /// <param name="lat2">第二點緯度</param> /// <param name="lng2">第二點經度</param> /// <returns>返回兩點的距離,單位:公裡或千米</returns> public static double GetDistance(double lat1, double lng1, double lat2, double lng2) { //地球半徑,單位米 double EARTH_RADIUS = 6378137; double radLat1 = Rad(lat1); double radLng1 = Rad(lng1); double radLat2 = Rad(lat2); double radLng2 = Rad(lng2); double a = radLat1 - radLat2; double b = radLng1 - radLng2; double result = 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))) * EARTH_RADIUS; return result / 1000; } /// <summary> /// 經緯度轉化成弧度 /// </summary> /// <param name="d"></param> /// <returns></returns> private static double Rad(double d) { return (double)d * Math.PI / 180d; } }
如有更好的方法,請留言謝謝!