PHP,Mysql-根据一个给定经纬度的点,进行附近地点查询–算法

需要计算目标周边1公里内的楼盘点

数据库几十万万条数据中查询出符合“周围1公里范围内”条件的坐标点

Mysql本身是支持空间索引的,但是在5.X版本中取消了Distance()和Related(),无法使用空间的距离函数去直接查询距离在一定范围内的点。所以,我首先想到的是,对每条数据去进行遍历,跟数据库中的每个点进行距离计算,当距离小于1公里时候,认为匹配成功。经测试,这样做确实能得到结果,但是效率极其低下,因为每条数据都得去和数据库中的几十万条数据进行比对,其耗费的时间可想而知

后来经过自己的仔细考虑,以及查询各种资料,想到一种利用正方形将方圆1公里这个圆包围起来。利用正方形的四个点,和大厦经纬度进行比较。问题转向了,如何计算正方形四个点经纬度


假设已知点的经纬度分别为lng,lat
  先实现经度范围的查询,
  在haversin公式中令φ1 = φ2
/**
     * 计算某个经纬度的周围某段距离的正方形的四个点
     *
     * @param
     *            radius 地球半径 平均6371km
     * @param
     *            lng float 经度
     * @param
     *            lat float 纬度
     * @param
     *            distance float 该点所在圆的半径,该圆与此正方形内切,默认值为1千米
     * @return array 正方形的四个点的经纬度坐标
     */
    public function returnSquarePoint($lng, $lat, $distance = 1, $radius = 6371)
    {
        $dlng = 2 * asin(sin($distance / (2 * $radius)) / cos(deg2rad($lat)));
        $dlng = rad2deg($dlng);
        
        $dlat = $distance / $radius;
        $dlat = rad2deg($dlat);
        
        return array(
            'left-top' => array(
                'lat' => $lat + $dlat,
                'lng' => $lng - $dlng
            ),
            'right-top' => array(
                'lat' => $lat + $dlat,
                'lng' => $lng + $dlng
            ),
            'left-bottom' => array(
                'lat' => $lat - $dlat,
                'lng' => $lng - $dlng
            ),
            'right-bottom' => array(
                'lat' => $lat - $dlat,
                'lng' => $lng + $dlng
            )
        );
    }


$array[0]就是起点终点坐标数组 
$room = $this->returnSquarePoint($array[0]['lng'], $array[0]['lat']);
->andwhere(['>','start_lat',$start['right-bottom']['lat']])
->andWhere(['<','start_lat',$start['left-top']['lat']])
->andWhere(['>','start_lng',$start['left-top']['lng']])
->andWhere(['<','start_lng',$start['right-bottom']['lng']]);
2.-----第二种

$range  = 180 / pi() * 1 / 6372.797; //里面的 1 就代表搜索 1km 之内,单位km
$lngR   = $range / cos($myLat * pi() / 180);
$maxLat = $myLat + $range; //最大纬度
$minLat = $myLat - $range; //最小纬度
$maxLng = $myLng + $lngR; //最大经度
$minLng = $myLng - $lngR; //最小经度
3.*----第三种
mysql两点间距,

公式如下,单位米:
第一点经纬度:lng1 lat1
第二点经纬度:lng2 lat2
round(6378.138*2*asin(sqrt(pow(sin(
(lat1*pi()/180-lat2*pi()/180)/2),2)+cos(lat1*pi()/180)*cos(lat2*pi()/180)*
pow(sin( (lng1*pi()/180-lng2*pi()/180)/2),2)))*1000)

例如:
SELECT store_id,lng,lat,
ROUND(6378.138*2*ASIN(SQRT(POW(SIN((22.299439*PI()/180-lat*PI()/180)/2),2)+COS(22.299439*PI()/180)*COS(lat*PI()/180)*POW(SIN((114.173881*PI()/180-lng*PI()/180)/2),2)))*1000)
AS
juli
FROM store_info
ORDER BY juli DESC
LIMIT 316
4.
select ROUND(6378.138*2*ASIN(SQRT(POW(SIN(($latitude*PI()/180-latitude*PI()/180)/2),2)+COS($latitude*PI()/180)*COS(latitude*PI()/180)*POW(SIN(($longitude*PI()/180-longitude*PI()/180)/2),2)))*1000) AS distance FROM shop having distance <= 5000 order by distance asc

6378.138这个是地球的直径,单位千米.latitude是用户位置的纬度,longitude是用户位置经度.latitude 为商户的纬度字段,longitude为商户的经度字段。上面一段sql计算得出根据用户经纬度计算与商户的距离。
点赞

发表评论