在postgresql中计算缓冲区内多边形的百分比

Question

我有一个名为 Operation 的 table，它有多个点作为 geom，我创建了一个 100m 的缓冲区并作为新列添加到同一个 table。我有另一个名为 Residential 的 table，它有多个多边形，目标是找到每个缓冲区内的多边形百分比，并将其添加为操作 table 中的一列。我不确定如何进行此操作。

SELECT AddGeometryColumn ('public','operations','buffer',4326,'POLYGON',2);
UPDATE operations SET buffer = ST_Buffer(geom::geography,100)::geometry;

ALTER TABLE operations ADD COLUMN pts int;

UPDATE operations o 
SET pts = (SELECT count(*) FROM supermarket s
           WHERE ST_Contains(o.buffer,s.geom));

到此为止，以下几行不适合获取百分比。如何解决这个问题。

    ALTER TABLE public."Operation" ADD COLUMN res_percent double precision;
UPDATE public."Operation"  
SELECT      
  ST_Intersection(ST_MakeValid(r.geom),o.buffer) AS intersection,   
  ST_Area(ST_Intersection(ST_MakeValid(r.geom),o.buffer))/ST_Area(r.geom)*100)) 
FROM public."Residential" r, public."Operation" o 
WHERE ST_Intersects(o.buffer,ST_MakeValid(r.geom));

dbfiddle

Answer 1

使用ST_Area to get the area of your polygons, extract the area of their intersection using ST_Intersection，最后用交集面积和多边形面积计算重叠比例

例子

给定两个重叠的多边形，p1 和 p2，在一个名为 t 的 table 中：

我们可以使用ST_Intersection得到两个多边形的交点：

SELECT ST_Intersection(p1,p2) FROM t;

现在我们可以用ST_Area来计算这个交点的面积：

SELECT ST_Area(ST_Intersection(p1,p2)) FROM t;

    st_area     
----------------
 601950.9523732
(1 row)

因此，在交叉点和多边形中使用 ST_Area，您可以计算多边形与另一个多边形重叠的百分比，例如

SELECT 
  ST_Area(ST_Intersection(p1,p2))/ST_Area(p2)*100 AS perc_p2, 
  ST_Area(ST_Intersection(p1,p2))/ST_Area(p1)*100 AS perc_p1
FROM t;

     perc_p2      |     perc_p1      
------------------+------------------
 30.0839473794556 | 37.9061430278047
(1 row)

演示：db<>fiddle

根据您的描述，您的查询应该看起来像这样：

SELECT   
  ST_Intersection(r.geom,o.buffer) AS intersection,
  ST_Area(ST_Intersection(r.geom,o.buffer))/ST_Area(r.geom)*100
FROM residential r, operations o
WHERE ST_Intersects(o.buffer,r.geom);

在postgresql中计算缓冲区内多边形的百分比

Calculate the percentage of polygon inside the buffer in postgresql

sql

postgresql

postgis

geospatial

spatial-query

例子