码农之家

专注优质代码开发,为软件行业发展贡献力量

根据一个给定经纬度的点和距离,搜索附近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;

        }

    }

0 评论数

添加新评论

带(*)的是必填项,不能输入html代码