Geo SQL查找位置附近的点

时间:2015-04-26 00:41:09

标签: mysql sql geo gtfs

我使用下面的sql:

基于http://www.scribd.com/doc/2569355/Geo-Distance-Search-with-MySQL

SELECT * , 3956 *2 * ASIN( SQRT( POWER( SIN( (
- 36.8812573 - abs( stop_lat ) ) * pi( ) /180 /2 ) , 2 ) + COS( - 36.8812573 * pi( ) /180 ) * COS( abs( stop_lat ) * pi( ) /180 ) * POWER( SIN( ( 174.63832160000004 - stop_lon ) * pi( ) /180 /2 ) , 2 ) )
) AS distance
FROM stops
HAVING distance <100
LIMIT 0 , 30

查找当前纬度和经度半径范围内的止损。

我正在使用GTFS数据,架构如下。 (纬度和经度以及最后两列)

01  7106    210 Victoria St     -36.8481990     174.7544900
0002    7108    220 Victoria St     -36.8481710     174.7523800
0003    7110    36 College Hill     -36.8485220     174.7485400
0004    7112    68 College Hill     -36.8479760     174.7466700
0005    7114    2 Jervois Rd    -36.8471670     174.7437800
0006    7116    90 Jervois Rd   -36.8453760     174.7408300
0007    7118    160 Jervois Rd  -36.8453250     174.7373100
0008    7120    206 Jervois Rd  -36.8454990     174.7354500
0009    7122    270 Jervois Rd  -36.8474010     174.7321400
0012    7121    203 Jervois Rd  -36.8473770     174.7319200
0013    7119    165 Jervois Rd  -36.8454100     174.7348400

表结构如下:

enter image description here

据我所知,查询应返回一些结果。但事实并非如此。谁能告诉我可能会缺少什么?

***更新***********

我注意到我错过了临时输出表名,所以我现在得到了一些结果。但是,它们距离我的位置5000起。有什么想法吗?

更新了SQL:

SELECT * , 3956 *2 * ASIN( SQRT( POWER( SIN( (
- 36.8812573 - abs( stop_lat ) ) * pi( ) /180 /2 ) , 2 ) + COS( - 36.8812573 * pi( ) /180 ) * COS( abs( stop_lat ) * pi( ) /180 ) * POWER( SIN( ( 174.63832160000004 - stop_lon ) * pi( ) /180 /2 ) , 2 ) )
) AS distance
FROM stops dest
HAVING distance <5100
ORDER BY distance
LIMIT 10 

谢谢,

1 个答案:

答案 0 :(得分:2)

我怀疑这个公式是罪魁祸首。即使您从可靠的来源获取它并且在实现中没有任何错误,也可以例如是一个单位混乱。

  • 出现错误:永远不需要absthe ‘haversine’ formula中没有这样的事情。在第二种情况下,它不会改变任何东西,因为cos是偶数函数,但在第一种情况下,它确实是。