根据与另一个 table 的几何交集将列数据添加到 table - PostgreSQL

Add column data to a table based on geometry intersection with another table - PostgreSQL

我的城市定义了一些 Urban Revitalization Perimeters 使用 PostGIS 我已经把城市分开了 Blocks within Perimeters 和 Buildings within Perimeters。我现在想做的是用两个新列更新建筑物 table,将它们与它们所在的街区相关联。

以下是Blocks的示例table,qd_setor和qd_fiscal这两列是我需要传递给Buildings的信息table

 qd_setor | qd_fiscal |  wkb_geometry 
----------+-----------+---------------
 041      | 038       | 010600...
 086      | 302       | 010600...
 041      | 076       | 010600...
 075      | 054       | 010600...
 019      | 055       | 010600...

这是建筑示例 table:

  ed_id  | ed_area | ed_altura |  ed_distrito  | wkb_geometry
---------+---------+-----------+---------------+--------------
 2739627 | 187.211 |     6.922 | SANTA_CECILIA | 010600...
 2739628 |  240.39 |    43.382 | SANTA_CECILIA | 010600...
 2739781 |  51.758 |     4.668 | SANTA_CECILIA | 010600...
 2739629 | 388.264 |    52.907 | SANTA_CECILIA | 010600...

这是我用来创建街区和建筑物的查询 tables:

CREATE TABLE quadras_opurb AS
SELECT q.qd_setor, q.qd_fiscal, q.wkb_geometry
FROM quadra_fiscal q
WHERE EXISTS (
SELECT ou_nome
FROM op_urbana_valid o
WHERE
ST_Intersects(q.wkb_geometry, o.wkb_geometry)
);

这是预期的结果,所以一旦完成,我就不必测试所需街区和建筑物之间的交集,我可以按 qd_setor 和 qd_fiscal 进行过滤:

  ed_id  | ed_area | ed_altura |  ed_distrito  | wkb_geometry | qd_setor | qd_fiscal 
---------+---------+-----------+---------------+--------------+----------+------------
 2739627 | 187.211 |     6.922 | SANTA_CECILIA | 010600...    | 041      | 038       
 2739628 |  240.39 |    43.382 | SANTA_CECILIA | 010600...    | 041      | 038       
 2739781 |  51.758 |     4.668 | SANTA_CECILIA | 010600...    | 041      | 038       
 2739629 | 388.264 |    52.907 | SANTA_CECILIA | 010600...    | 086      | 302
 2739631 | 187.244 |    22.865 | SANTA_CECILIA | 010600...    | 086      | 302  

请注意,qd_setor 和 qd_fiscal 将重复多次,因为一个街区中有多个建筑物。

感谢您的帮助!

所以我弄明白了,我使用了以下查询:

CREATE TABLE ed_quadra AS
SELECT
q.qd_setor AS setor,
q.qd_fiscal AS quadra,
e.*
FROM ed_opurb AS e
JOIN quadras_opurb AS q
ON ST_Intersects(e.wkb_geometry, q.wkb_geometry)
ORDER BY q.qd_setor ASC, q.qd_fiscal ASC;