跟据经纬度实现附近搜索及按距离排序

《跟据经纬度实现附近搜索》

网上有很多资料，mysql空间数据库、矩形算法、geohash等，当数据上了百万之后mysql空间数据库方法是最强最精确的（查询前100条数据只需5秒左右）。

接下来推出一个计算方法，查询速度是mysql空间数据库算法的2倍

$lng是你的经度，$lat是你的纬度

SELECT lng,lat,
            (POWER(MOD(ABS(lng - $lng),360),2) + POWER(ABS(lat - $lat),2)) AS distance 
            FROM `user_location`
            ORDER BY distance LIMIT 100
        

经测试，在100万数据中取前100条数据只需2.5秒左右。

注意以上查询出来的distance只是一个系数,而非两点间距离，只用在排序。

以下查询语句可以约等于单位米m（去掉*1000则可约等于单位千米km）：

SELECT lng,lat,
            (SQRT(POWER(MOD(ABS(lng - {$lng}),360),2) + POWER(ABS(lat - {$lat}),2))*80*1000) AS distance 
            FROM `user_location`
            ORDER BY distance LIMIT 100
        

####################################
另外的几种算法还是在这里展示一下：

一、距形算法

define(EARTH_RADIUS, 6371);//地球半径，平均半径为6371km
/**
 *计算某个经纬度的周围某段距离的正方形的四个点
 *
 *@param lng float 经度
 *@param lat float 纬度
 *@param distance float 该点所在圆的半径，该圆与此正方形内切，默认值为0.5千米
 *@return array 正方形的四个点的经纬度坐标
 */
function returnSquarePoint($lng, $lat,$distance = 0.5){
    $dlng =  2 * asin(sin($distance / (2 * EARTH_RADIUS)) / cos(deg2rad($lat)));
    $dlng = rad2deg($dlng);
    $dlat = $distance/EARTH_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)
    );
}

//使用此函数计算得到结果后，带入sql查询。
$squares = returnSquarePoint($lng, $lat);
$info_sql = "select id,locateinfo,lat,lng from `lbs_info` where 
    lat<>0 and lat>{$squares['right-bottom']['lat']} and 
    lat<{$squares['left-top']['lat']} and 
    lng>{$squares['left-top']['lng']} and 
    lng<{$squares['right-bottom']['lng']} ";
        

二、 空间数据库算法

以下location字段是跟据经纬度来生成的空间数据，如：
location字段的type设为point

"update feed set location=GEOMFROMTEXT('point({$lat} {$lng})') where id='{$id}'"

mysql空间数据查询

SET @center = GEOMFROMTEXT('POINT(35.801559 -10.501577)');
SET @radius = 4000;
SET @bbox = CONCAT('POLYGON((',
            X(@center) - @radius, ' ', Y(@center) - @radius, ',',
            X(@center) + @radius, ' ', Y(@center) - @radius, ',',
            X(@center) + @radius, ' ', Y(@center) + @radius, ',',
            X(@center) - @radius, ' ', Y(@center) + @radius, ',',
            X(@center) - @radius, ' ', Y(@center) - @radius, '))'
    );

SELECT id,lng,lat,
            SQRT(POW( ABS( X(location) - X(@center)), 2) + POW( ABS(Y(location) - Y(@center)), 2 )) AS distance 
            FROM `user_location` WHERE 1=1 
            AND INTERSECTS( location, GEOMFROMTEXT(@bbox) ) 
            AND SQRT(POW( ABS( X(location) - X(@center)), 2) + POW( ABS(Y(location) - Y(@center)), 2 )) < @radius 
            ORDER BY distance LIMIT 20
        

三、geo算法

http://blog.sina.com.cn/s/blog_62ba0fdd0100tul4.html