如何解决计算postgresql缓冲区内多边形的百分比
我有一个名为 Operation 的表,其中有多个点作为 geom,我创建了一个 100m 的缓冲区并作为新列添加到同一个表中。我有另一个名为 Residential 的表,它有多个多边形,目标是找到每个缓冲区内多边形的百分比,并将其添加为操作表中的一列。我不知道如何继续。
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));
解决方法
使用ST_Area
获得多边形的面积,使用ST_Intersection
提取它们的交集面积,然后最后使用交集面积和多边形面积计算重叠比例。
示例
给定两个重叠的多边形,p1
和 p2
,位于名为 t
的表中:
我们可以使用 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);
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。