空间聚类 - 将聚类属性（ID）与属于聚类的几何图形相关联

Question

空间聚类 - 将聚类属性（ID）与属于聚类的几何图形相关联

3

我在将一组集群几何体与他们自己的所有权关联时遇到了一些问题。

数据

我有一个包含一组几何图形的表格，

buildings {
gid integer,
geom geometry(Multipoligon,4326)
}

我已经在“建筑物”表上运行了 ST_ClusterWithin函数，并设置了一定的阈值进行聚类分析。从该聚类分析中，我得到了一个名为“clusters”的表。

clusters {
cid Integer,
geom geometry(GeometryCollection,4326)
}

问题

我希望能够提取所有几何信息及其关联的聚类信息，并将其放入表格中。

clustered_building {
gid Integer
cid Integer
geom geometry(Multipoligon,4326)
}

 gid |    cid     |       geom            |   
-----+------------+-----------------------+
  1  |     1      | multypoligon(...)     |
  2  |     1      | multypoligon(...)     |
  3  |     1      | multypoligon(...)     |
  4  |     2      | multypoligon(...)     |
  5  |     3      | multypoligon(...)     |
  6  |     3      | multypoligon(...)     |

我做了什么（但不起作用）

我一直在尝试使用两个函数ST_GeometryN/ST_NumGeometries解析每个MultyGeometry，并使用这个查询从ST_Geometry手册页面的标准示例中提取聚类信息。

INSERT INTO  clustered_building (cid, c_item , geom)
SELECT sel.cid, n, ST_GeometryN(sel.geom, n) as singlegeom
FROM ( SELECT cid, geom, ST_NumGeometries(geom) as num
       FROM clusters") AS sel
       CROSS JOIN generate_series(1,sel.num) n
WHERE n <= ST_NumGeometries(sel.geom);

如果我强制使用一系列10个查询，它需要几秒钟的时间。

CROSS JOIN generate_series(1,10)

但是当我要求根据每个GeometryCollection中的项目数量生成一系列时，它卡住了。而且，由于我丢失了“gid”，这个查询不允许我将单个几何图形链接到其自己的建筑表中的特征。

请问有人能帮帮我吗？谢谢。

Stefano

- entalpia

1个回答

网页内容由stack overflow 提供, 点击上面的

可以查看英文原文，
原文链接

- John Powell · Accepted Answer

我没有你的数据，但是使用一些虚拟值，其中ID 1、2和3相交，4和5，您可以执行以下操作：

WITH 
  temp (id, geom) AS 
    (VALUES (1, ST_Buffer(ST_Makepoint(0, 0), 2)),
    (2, ST_Buffer(ST_MakePoint(1, 1), 2)),
    (3, ST_Buffer(ST_MakePoint(2, 2), 2)), 
    (4, ST_Buffer(ST_MakePoint(9, 9), 2)), 
    (5, ST_Buffer(ST_MakePoint(10, 10), 2))),
  clusters(geom) as 
    (SELECT 
        ST_Makevalid(
          ST_CollectionExtract(
              unnest(ST_ClusterIntersecting(geom)), 3)) 
      FROM temp
    )
 SELECT array_agg(temp.id), cl.geom 
   FROM clusters cl, temp 
  WHERE ST_Intersects(cl.geom, temp.geom) 
GROUP BY cl.geom;

如果你将最终的cl.geom用ST_AsText包装起来，你会看到类似以下的内容：

{1,2,3} | 多边形(((2.81905966523328 0.180940334766718,2.66293922460509 -0.111140466039203,2.4142135623731 -0.414213562373094,2.11114046603921 -0.662939224605089,1.81905966523328 -0.819059665233282,1.84775906502257 -0.765366864730179,1.96157056080646 -0.390180644032256,2 0,2 3.08780778723872e-16,2 0,2.39018064403226 0.0384294391935396,2.76536686473018 0.152240934977427,2.81905966523328 0.180940334766718))...... {4,5} | 多边形(((10.8190596652333 8.18094033476672,10.6629392246051 7.8888595339608,10.4142135623731 7.58578643762691,10.1111404660392 7.33706077539491,9.76536686473018 7.15224093497743,9.39018064403226 7.03842943919354,9 7,8.60981935596775 7.03842943919354,8.23463313526982 7.15224093497743,7.8888595339608 7.33706077539491,7.58578643762691 7.5857864376269,7.33706077539491 7.88885953396079,7.15224093497743 8.23463313526982

您可以看到，编号为1、2、3的多边形属于第一个多边形集合，而编号为4、5的多边形属于另一个多边形集合。

总体思路是对数据进行聚类，然后使用array_agg将返回的聚类与原始数据进行交集运算，将相同的id分组在一起，以便返回的多边形包含原始的id。使用ST_CollectionExtract和参数3，结合unnest，将由ST_ClusterIntersecting返回的几何集合拆分成行，返回每个连续的聚类作为（多）边形。ST_MakeValid用于处理当您将几何图形与其他相关几何图形（例如原始多边形与聚类多边形）进行交集运算时，可能会出现奇怪的舍入效应和GEOS错误等问题。

我最近在gis.stackexchange上回答了一个类似的问题，你可能会觉得有用。