PostGIS如何解决多边形自相交处理PostGIS如何解决多边形自相交处理

115 阅读 0 评论 76 点赞

我是靠谱客的博主粗暴金针菇，这篇文章主要介绍PostGIS如何解决多边形自相交处理PostGIS如何解决多边形自相交处理，现在分享给大家，希望可以做个参考。

PostGIS如何解决多边形自相交处理

postgis 是一个很强大的空间数据库，最近查看很多资料和文档，为了解决空间边界自相加问题。
知乎有一篇写到处理这类函数，但是结果可能会出现自相交部分出现镂空，这并不是我想要的。
知乎直通车

复制代码

1
2
这里先感谢公司一位大佬提供的解决方案,这里为了表达敬仰,大伙一起隔空喊一声 ：帆哥,牛逼

首先我希望得到的图形是这样的：
在这里插入图片描述
我在地图上画边界时这样的：

相信看到这懂的都懂了，是不是很兴奋：“我感觉到了，一切都回来了，大师兄~”

复制代码

1
2
3
4
如图：左边的是自相交图形, 右边的是正常图形
问题：当图形格式为集合类型(FeatureCollection),该数据集中存在N个多多边形,其中又包含自相交图形,
如何得到一个叠加后正确且完整图形?

直接上答案：
复制粘贴以下SQL格式，替换GeoJSON内容即可:

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
SELECT ST_Union(t.geom) FROM (
    WITH data AS (SELECT
    #{
   这里放GeoJSON格式数据}::json AS fc)
    SELECT
      row_number() OVER () AS gid,
      ST_MakeValid(ST_GeomFromGeoJSON(feat->>'geometry')) AS geom,
      feat->'properties' AS properties
    FROM (
      SELECT json_array_elements(fc->'features') AS feat
      FROM DATA
    ) AS f
) AS t

解题思路:
上方图形对应数据(GeoJSON完整格式)
这里提供一个地图绘制工具非常好用 GeoJSON 来绘制图形

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
{
   
	"type": "FeatureCollection",
	"features": [{
   
		"type": "Feature",
		"properties": {
   },
		"geometry": {
   
			"type": "Polygon",
			"coordinates": [
				[
					[120.52757263183594, 31.495432803134843],
					[120.59211730957031, 31.31199502365151],
					[120.89012145996094, 31.35950051982242],
					[120.74729919433594, 31.283245492650792],
					[120.52757263183594, 31.495432803134843]
				]