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 }