LBS附近地点搜索算法资料集合

方法一：基于球面距离搜索附近地点

• 附近地点搜索初探

    点评：需要使用2个字段进行查询，查询效率差，适用于数据量较少的小型应用。

• Haversine formula - Wikipedia

    维基百科推荐使用 Haversine 公式计算球面距离

• Great-circle distance - Wikipedia

 

方法二：基于Geohash算法搜索附近地点

• geohash：用字符串实现附近地点搜索

    点评：精度可控，使用单字段字符串前缀查询，查询效率较高，适用于大数据的应用。

• Geohash - Wikipedia
• geohash演示

 

其它方法：

• MySQL中的空间扩展

 

PHP源代码

 1 //
 2 // 使用 Haversine 公式计算两个地理坐标点之间的球面距离。
 3 //
 4 // 参数说明:
 5 // $lat1        点1的纬度值，单位：degree
 6 // $lng1        点1的经度值，单位：degree
 7 // $lat2        点2的纬度值，单位：degree
 8 // $lng2        点2的经度值，单位：degree
 9 //
10 // 返回值：两点之间的球面距离，单位: km
11 //
12 function haversine_distance($lat1, $lng1, $lat2, $lng2) {
13   $EARTH_RADIUS = 6371.00;   // 地球平均半径，6371km
14   
15   // 角度转换成弧度
16   $rlat1 = deg2rad($lat1);
17   $rlng1 = deg2rad($lng1);
18   $rlat2 = deg2rad($lat2);
19   $rlng2 = deg2rad($lng2);
20   
21   $rlat_diff = abs($rlat1 - $rlat2);
22   $rlng_diff = abs($rlng1 - $rlng2);
23   
24   $h = pow(sin($rlat_diff/2), 2) + cos($rlat1) * cos($rlat2) * pow(sin($rlng_diff/2), 2);
25   return 2 * $EARTH_RADIUS * asin(sqrt($h));
26 }
27 
28 //
29 // 根据到中心点的距离计算东西两侧的经度边界范围，采用 Haversine 公式计算。
30 //
31 // 参数说明：
32 // $lat         中心点的纬度，单位：degree
33 // $lng         中心点的经度，单位：degree
34 // $distance    到中心点的距离，单位：km
35 // 
36 // 返回值: array($lat_diff, $lng_diff)
37 // $lat_diff    到中心点的纬度边界范围，单位: degree
38 // $lng_diff    到中心点的纬度边界范围，单位: degree
39 //
40 function lbs_haversine_diff($lat, $lng, $distance) {
41   $EARTH_RADIUS = 6371.00;   // 地球平均半径，6371km
42   
43   // 纬度边界
44   $rlat_diff = $distance / $EARTH_RADIUS;
45   $dlat_diff = rad2deg($rlat_diff);
46   
47   // 经度边界
48   $rlat = deg2rad($lat);
49   $rlng_diff = 2 * asin( sin($distance / (2 * $EARTH_RADIUS)) / cos($rlat) );
50   $dlng_diff = rad2deg($rlng_diff);
51   
52   return array($dlat_diff, $dlng_diff);
53 }