Question

我想要交叉两张桌子。第一个表代表大约5000万个点，第二个表是世界上所有国家的多边形层。我想获得与该多边形相交的所有点。

SELECT d.id, d.geom 
FROM export d, world_boundaries b 
WHERE (b.cntry_name = 'UK') 
  AND d.date_inserted >= '2012-06-01' 
  AND d.geom && b.wkb_geometry 
  AND intersects(d.geom, b.wkb_geometry);

此查询非常简单，但运行时间超过4小时。我在每个表的几何列上构建了GIST索引，并对它们进行了VACUUM ANALYZE。仍然没有性能提升。我使用Postgres 8.4和PostGIS 1.5运行CENTOS 6。有人能说明如何加快速度吗？将查询限制为1000到10000条记录时，我会很快得到结果。当我试图抓住完整的结果集时，它会拖延。想法？

更新：我现在看到我必须优化我的查询作为此过程的第一步。我得到了这样的信封

select astext(st_envelope(wkb_geometry)) as e 
from world_borders 
where cntry_name = 'UK'

现在，作为整个查询的一部分包含/执行此操作的最有效方法是什么？

Answer 1

尝试使用EXPLAIN（和LIMIT）运行它以查看索引是否正在使用。

由于真正的交集检查是那里最慢的操作，可能是针对子查询的ST_Collect运行它（除了ST_Intersects检查之外的所有内容）都会有帮助。那样只有一个调用，如果多几何结构足够快，最终结果可能会更好。

EDIT1：好吧，事实证明它不是那么理想，因为除非你强制坐标到3d（也保持id），需要额外的查找来获得几何id：

SELECT d.id, d.geom
FROM
(
    SELECT *
    FROM
    ( 
        SELECT ST_Collect(d.geom)
        FROM export d, world_boundaries b 
        WHERE (b.cntry_name = 'UK') 
        AND d.date_inserted >= '2012-06-01' 
        AND d.geom && b.wkb_geometry
    ) as c, world_boundaries b 
    WHERE (b.cntry_name = 'UK')
    AND ST_Intersection(c.geom, b.wkb_geometry);
) as e, export d
WHERE (ST_Dump(e.geom)).geom = d.geom

Postgresql Spatial Query太慢了

1 个答案: