postgresql - ST_ApproximateMedialAxis 在最终近似值上留下骨架线

Question

我想近似一个封闭的多边形，但 PostGIS 给我的不是一个干净的线串，而是带有 3 个线串的多线串。都是因为直骨架上留下的那条小尾巴。 有没有合适或好的方法来处理这个问题？

Answer 1

您可以st_approximatemedialaxis()对输出进行后处理，以通过某些标准（例如点数）过滤掉这些尾部：

select st_astext(geom)
from (select (st_dump(x.geom)).geom
    from (values (st_geomfromtext('MULTILINESTRING( (27 80, 18 68, 29 48, 55 58, 53 76, 27 80), (55 58, 58 57))')),
    (st_geomfromtext('MULTILINESTRING( (46 34, 32 21, 44 7, 67 9, 67 29, 46 34), (46 34, 46 36))')),
    (st_geomfromtext('MULTILINESTRING( (69 66, 61 48, 75 40, 94 54, 88 68, 69 66), (61 48, 58 47))'))) as x(geom)) y
where st_npoints(geom)>2

返回

LINESTRING(27 80,18 68,29 48,55 58,53 76,27 80)
LINESTRING(46 34,32 21,44 7,67 9,67 29,46 34)
LINESTRING(69 66,61 48,75 40,94 54,88 68,69 66)