查找相邻多边形 - postgis查询

Question

问题： 我有一张满是多边形的桌子（一个国家的区域）。这些多边形中的一些具有可能或可能不被检查的特定属性。在这种情况下，该属性称为“spread”，“checked”值为1.

现在我想运行一个查找所有“已检查”多边形的查询。 类似的东西：

SELECT * FROM gemstat WHERE spread = 1

然后我想在每个尚未检查的相邻多边形中将“spread”属性设置为“1”。 （我也想设置第二个属性，但这只是一个小的补充）

首先，让我们从一个查询开始，该查询选择所有多边形的多边形，并使用spread-value = 1

SELECT (b."GEM_NR")
FROM gemstat_simple5 as a
JOIN gemstat_simple5 as b
ON ST_Touches((a.the_geom),b.the_geom)
where a.spread =1;

此查询返回所有多边形的多边形，其中spread = 1

现在我想根据该子查询的结果更新表。 这是通过John Powell aka Barca提供的这段代码实现的（参见下面的回答和评论）：

Update gemstat_simple5 gem set spread=1, time=2
FROM (
   SELECT (b."GEM_NR")
   FROM gemstat_simple5 as a,
   gemstat_simple5 as b
   WHERE ST_Touches(a.the_geom, b.the_geom) 
   AND a."GEM_NR" != b."GEM_NR"
   AND a.spread = 1
) as subquery
WHERE gem."GEM_NR" = subquery."GEM_NR"

运行此查询，它将设置扩展为1的属性和相邻多边形的2个时间，同时不触及spread = 1的原始多边形。 因此，它为我的问题提供了完美的答案。

Answer 1

如果您询问如何基于子查询更新表，只查找那些具有邻居的多边形（即，它们接触另一个多边形），那么以下内容应该适合您。

Update gemstat_simple5 gem set spread=1, time=2
  FROM (
     SELECT (b."GEM_NR")
       FROM gemstat_simple5 as a,
            gemstat_simple5 as b
       WHERE ST_Touches(a.the_geom, b.the_geom) 
       AND a."GEM_NR" != b."GEM_NR"
       AND a.spread = 1
     ) as subquery
 WHERE gem."GEM_NR" = subquery."GEM_NR"

注意，我已经把AND a。“GEM_NR”＆lt; b。“GEM_NR”，它们都避免了AND a。“GEM_NR”= b。“GEM_NR”，即你想要避免的本身的情况，并且还将成对比较减少了一半。我也使用a，b where方法而不是st_touches上的连接b，这是相同的，但我发现空间连接更令人困惑。最后，您只需将要更新的表GEM_NR与在子查询中找到的表等同起来。

Answer 2

直接答案是

SELECT a.*
FROM polygon1 as a
JOIN polygon1 as b
ON st_intersects((st_buffer(a.the_geom,0.00001)),b.the_geom) 
where b.id = 561334;

在我的机器中测试过它还会给出你在哪个类中给出的id。 我的多边形表的空间参考 - EPSG：4326 ......