查找相邻多边形 - 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”<b。“GEM_NR”，它们都避免了AND a。“GEM_NR”= b。“GEM_NR”的情况，即你想要避免的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 ...