使用空间连接子查询更新 postgis 表

Question

我维护两个 postgis 表：“track_points”和“buffers”。 “track_point”表包含大量（近 10 亿个）点，“buffer”表包含大约 20 个多边形。

我想要做的是，检查包含它们的缓冲区中的所有点，并将相应的缓冲区 id 分配给点记录。在网上搜索后，我发现“空间关节”在这里可能有很大帮助。根据我在网上找到的内容，我整理了一个如下所示的查询（{schema} 只是架构名称的占位符）：

WITH join_query AS (
  SELECT
    points.id AS point_id,
    buffers.profile_id AS profile_id
  FROM {schema}.buffers AS buffers
  JOIN {schema}.track_points AS points
  ON ST_Contains(buffers.geom, points.geom)
)

UPDATE {schema}.track_points
  SET profile_id = join_query.profile_id
  FROM join_query
  WHERE id = join_query.point_id

我运行了查询，但 track_points 表中的 profile_id 值都没有改变。所以我想我的查询一定有问题？！？？

另外，是否有人建议如何更有效地实现我的目标（关于 track_points 表中的大量点）？

顺便说一下，我正在使用 Python 的 psycopg2 连接数据库。

Answer 1

如果您的点表有数十亿条记录，甚至不要尝试更新它 - 或者您可以等待几天/几周来结束此更新；）。对于这种大规模操作的完美解决方案是 CTAS（创建表作为选择）； 我假设你的多边形不相交，如果是，那么告诉我 profile_id 你想要哪个缓冲区（最大，最小....）；

create table track_points2 as
select your_columns_for_track_points(expect profile_id), b.profile_id 
  from track_points tp, buffers b
 where st_dwithin(tp.geom, b.geom,0);

接下来，删除现有表 track_points 并将其替换为新表；

drop table track_points;
alter table track_points2 alter rename to track_points;

并为您的新表创建所有需要的索引和约束。

如果您无法删除数据库中的表、更改表等，那么当然，您将不得不更新，但要准备好长时间等待。

 update track_points tp
    set profile_id=b.profile_id
   from buffers b
  where st_dwithin(tp.geom, b.geom,0);

正如我之前写的，如果您有相交的缓冲区/多边形，那么您将不得不更改更新以从许多 profile_id 选择中获得您想要的。