c# 根據一個給定經緯度的點和距離，搜索附近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;
        }
    }

如有更好的方法，請留言謝謝！