基于MYSQL获取地理位置附近的点

由于最近一年多大部分都做的是LBS相关项目，发现现在安全、社交、O2O等等的手机APP都会使用LBS技术，基于LBS位置服务的APP用户体验会更好，也更容易挖掘到用户真正的需求。下面总结了两个项目中“附近的”解决方案仅供参考，有什么问题欢迎指正。

1、近似（矩形范围内）

实现原理简单粗暴，先算出该点周围的矩形的四个点，然后使用经纬度去直接匹配数据库中的记录（参考）

<span style="color:#657b83">   //获取周围坐标
   public function returnSquarePoint($lng, $lat,$distance = 0.5){
        $earthRadius = 6378138;
        $dlng =  2 * asin(sin($distance / (2 * $earthRadius)) / cos(deg2rad($lat)));
        $dlng = rad2deg($dlng);
        $dlat = $distance/$earthRadius;
        $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)
        );
   }

   //计算两个坐标的直线距离
   public function getDistance($lat1, $lng1, $lat2, $lng2){      
        $earthRadius = 6378138; //近似地球半径米
        // 转换为弧度
        $lat1 = ($lat1 * pi()) / 180;
        $lng1 = ($lng1 * pi()) / 180;
        $lat2 = ($lat2 * pi()) / 180;
        $lng2 = ($lng2 * pi()) / 180;
       // 使用半正矢公式  用尺规来计算
       $calcLongitude = $lng2 - $lng1;
       $calcLatitude = $lat2 - $lat1;
       $stepOne = pow(sin($calcLatitude / 2), 2) + cos($lat1) * cos($lat2) * pow(sin($calcLongitude / 2), 2);  
       $stepTwo = 2 * asin(min(1, sqrt($stepOne)));
       $calculatedDistance = $earthRadius * $stepTwo;
       return round($calculatedDistance);
   }</span>

2、利用mysql 5.6+自带的空间关系函数

     后来查找资料发现mysql5.6+以上版本的空间关系函数：st_distance()，返回两点间的距离值（Mysql 5.6官方帮助文档）； 

<span style="color:#657b83">ST_Distance(g1,g2)
 
Returns the distance between g1 and g2.
 
mysql> SET @g1 = POINT(1,1), @g2 = POINT(2,2);
mysql> SELECT ST_Distance(@g1, @g2);
+-----------------------+
| ST_Distance(@g1, @g2) |
+-----------------------+
|    1.4142135623730951 |
+-----------------------+
ST_Distance() was added in MySQL 5.6.1.</span>

我项目的实现需求：

• 显示我附近的点；
• 由远及近排序；
• 计算距离；

数据库表结构：

假设我的坐标是116.397915,39.908946，得到附近100km内的点位信息:

<span style="color:#657b83">SELECT id,title,AsText(pt),(st_distance(pt, point(116.397915,39.908946))*111) AS distance FROM my_point HAVING distance<100 ORDER BY distance limit 100;</span>