在 MySQL 中对空间表执行子查询

Question

我有一个polygons 的空间表。我正在尝试有效地在 MySQL 空间表中找到 polygon 和 lat/lng point 的内容。

我知道ST_Within 可用于MySQL 5.6，但是我的版本早于这个版本，所以我将我可以使用的两个功能组合成比单独使用更有效/更准确的东西。

这些是原生MBRWITHIN 函数和自定义GISWithin from a MySQL forum post。 MBRWITHIN 通常给我大约 2 行的点，应该只有一个但很快，GISWithin 给出正确的行，但是很慢。

所以我正在尝试编写一个查询：

最初使用 快速 MBRWITHIN 功能。

然后使用较慢的方法将此列表细化为单个正确的行 GISWithin 现在在大约 2 行而不是 9000 行上运行。

我尝试过各种子查询，比如

SET @x = -0.3578;
SET @y = 51.477;
SET @point = CONCAT('POINT(',@x,' ',@y,')');


SELECT * FROM 
(SELECT `geometry` FROM world_borders WHERE MBRWITHIN( @point , `geometry` )) AS TT 
WHERE GISWithin(GeomFromText(@point), `geometry`); 

和

SELECT * FROM world_borders WHERE GISWithin(GeomFromText(@point), `geometry`) IN (SELECT ogr_fid FROM world_borders WHERE MBRWITHIN( @point , `geometry` ));

只为了得到错误信息

#1416 - Cannot get geometry object from data you send to the GEOMETRY field 

有谁知道我如何编写这样的查询或函数？

我表中的多边形字段称为geometry，类型为Geometry - 从ogr2ogr 导入，供知道的人使用。

Answer 1

错误来自尝试在需要几何对象的地方使用点的 WKT（文本）表示，因为 @point 包含文本字符串，而您并未始终使用 GeomFromText() 来转换它。

但是您根本不需要连接字符串并将其转换为几何，因为 MySQL 有一个内置的POINT() 函数可以直接生成几何对象。

SET @x = -0.3578;
SET @y = 51.477;
SET @point = POINT(@x,@y);

@point 现在包含一个适当的几何对象，查询应该可以通过简单地将两个测试与AND 组合在一个查询中来实现。

由于两个条件都必须为真，优化器应该首先尝试做最简单的事情——意识到MBRWithin()可以使用空间索引来解析，它应该尝试从索引中找到匹配的行，然后执行GISWithin() 函数仅在它找到的行上运行，这将进一步消除结果集中任何不使用存储函数计算为真的行。

SELECT * 
  FROM world_borders 
 WHERE MBRWithin(@point,`geometry`)
   AND GISWithin(@point,`geometry`);