Question

我试图在php应用程序中实现geospacial搜索。目前，我使用以下查询查找给定纬度和经度的10公里范围内的点：

SELECT * FROM (
      SELECT *,
             (6378.1 * ACOS(
                  COS(RADIANS(48.856614)) * COS(RADIANS(latitude))
                * COS(RADIANS(2.3522219) - RADIANS(longitude))
                + SIN(RADIANS(48.856614))
                * SIN(RADIANS(latitude)))
             ) AS distance
      FROM `destinations_info`
      WHERE latitude BETWEEN 48.856614  - (10 / 69)
               AND 48.856614  + (10 / 69)
        AND longitude BETWEEN 2.3522219 - (10 / (69 * COS(RADIANS(48.856614))))
               AND 2.3522219 + (10 / (69 * COS(RADIANS(48.856614))))
      ) d
WHERE distance <= 10
LIMIT 1

只要我不搜索存储在POI表中的确切纬度和经度，这就可以正常工作。

例如，如果我在我的poi表中有以下条目

id     name     latitude    longitude
1      Paris    48.856614   2.3522219000000177

我用lat = 48.856614和long = 2.3522219000000177来调用该查询，但我没有得到任何结果。据我所知，这是因为以下操作在mysql中返回NULL而不是0：

SELECT (6378.1 * acos(cos(radians(48.856614)) * cos(radians(48.856614)) * cos( radians(2.3522219) - radians(2.3522219)) + sin(radians(48.856614)) * sin(radians(48.856614))))

在Chrome计算器I get 0 rad中运行相同内容。

查询中是否有问题，或者我是否需要包含NULL结果以及使用＆＃34; OR IS NULL＆＃34;

Answer 1

IFNULL（...），0）AS Distance也可以帮助您解决错误。

Answer 2

问题是调用acos的值不在-1和1的范围内，因此它返回null。

SELECT cos(radians(48.856614)) * cos(radians(48.856614)) * cos( radians(2.3522219) - radians(2.3522219)) + sin(radians(48.856614)) * sin(radians(48.856614))

1.0000000000000002

这似乎是某种舍入问题。

因为，arc cos只定义在-1和1之间，你可能应该预先检查它的参数，或者如果计算出错了就可以返回null。

使用hasrsine公式的MySQL地理空间搜索在同一点上返回null

2 个答案: