Question

我在玩OpenStreetMaps数据。使用Postgis 2.2。将Osmosis加载到Postgres 9.4中，使用为ways表构建的线串几何。

我想从某些relations制作多边形（在OSM中，关系是一组方式（线串） - 一些标记为外部，一些标记为内部，而不是每种方式都必须关闭 - 多种方式一起可以表示外环）。但OSM数据并不完美 - 存在形式错误的关系，其中内部＆＃34;漏洞＆＃34;相交外圈。例如：http://www.openstreetmap.org/relation/6095901#map=17/48.93209/19.14858

我正在尝试此查询

SELECT ST_BuildArea(ST_Collect(w.linestring)) AS geom
FROM relations r
JOIN relation_members rm ON r.id = rm.relation_id AND rm.member_type = 'W'
JOIN ways w ON w.id = rm.member_id
WHERE r.id = 6095901 -- in reality here is broad condition on relation
GROUP BY r.id

，结果是错误：

ERROR:  LWGEOM_GEOS_buildArea: TopologyException: Input geom 1 is invalid: Self-intersection at or near point 19.149718131981164 48.934732947538478 at 19.149718131981164 48.934732947538478

那么，我的选择是什么？理想情况下，我想要解决这个问题＆＃34;不知何故 - 可能会忽视有问题的＆＃34;漏洞＆＃34;或者＆＃34;剪辑＆＃34;他们到外环边界等。

如果不可能，至少我希望查询到＆＃34;跳过＆＃34; （对于这样的无效关系，返回NULL几何），并且在没有结果的情况下不会因第一个错误而死亡。

编辑：原始示例案例已在OSM数据中修复。这是一个新的例子：关系ID 5636774是自相交的多面

Answer 1

我相信您可以使用ST_Union代替ST_Collect来避免这些错误。它比较慢，但它会摆脱重叠的部分。 ST_Collect的文档中有关于它的说明：

http://postgis.net/docs/ST_Collect.html

以下是我测试你的第二个例子：

create table ways ( id integer, linestring  geometry(Linestring, 4326) );
insert into ways values (1, ST_GeomFromText('LINESTRING(52.3474437 35.3635293,52.347219 35.3636168,52.3472848 35.36407,52.3475095 35.3639826,52.3474437 35.3635293 )', 4326));
insert into ways values (1, ST_GeomFromText('LINESTRING(52.3473579 35.3642849,52.347274 35.3636982,52.3468764 35.3638506,52.3468949 35.3639796, 52.3472102 35.3638587, 52.3472757 35.3643164, 52.3473579 35.3642849)', 4326));
insert into ways values (1, ST_GeomFromText('LINESTRING(52.3468455 35.3636261, 52.34673 35.3636689, 52.3468697 35.3646806, 52.3474622 35.3644612, 52.3474339 35.3642567, 52.3473579 35.3642848, 52.3472762 35.3643151, 52.346957 35.3644333,  52.3468945 35.3639806,  52.3468765 35.3638507, 52.3468455 35.3636261)', 4326));
select ST_BuildArea(ST_Collect(linestring)) from ways;
select ST_BuildArea(ST_Union(linestring)) from ways;

我还考虑使用ST_IsValid，但事实证明ST_IsValid(ST_Collect(linestring))即使对于这种糟糕的形状也是如此，所以我想这对你没有帮助。

如何从无效的线串构建区域（或防止错误）

1 个答案: