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

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

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

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

/**
     * 计算某个经纬度的周围某段距离的正方形的四个点
     *
     * @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计算得出根据用户经纬度计算与商户的距离。