从 Shapely 中的多边形中删除重复点

Question

要从列表中删除重复列表，Python 中有几种不错的方法 - for example:

a = [[ 9.1514622, 47.1166004 ], [ 9.1513045, 47.1164599 ], [ 9.1516278, 47.1163001 ], [ 9.1517832, 47.1164408 ], [ 9.1514622, 47.1166004 ] ] 

print len(a) # 5
b_set = set(map(tuple,a))
b = map(list,b_set)
print len(b) # 4

但不幸的是，我不得不将我的列表转换为 Shapely Polygon 对象，我需要在其中简化几何图形并执行一些其他地理功能。

from shapely.geometry import Polygon
a = [[[ 9.1514622, 47.1166004 ], [ 9.1513045, 47.1164599 ], [ 9.1516278, 47.1163001 ], [ 9.1517832, 47.1164408 ], [ 9.1514622, 47.1166004 ] ] ]
polys = [Polygon(item) for item in a] # convert list to polygon
print len(polys) # prints 5

展示了如何从多边形列表中删除重复的多边形，但是如何从点列表中删除重复的点，作为形状多边形？

我想可以将其转换回列表，删除重复项，然后重新转换为多边形。

但这似乎过于复杂了。关于如何做到这一点有什么想法吗？

Answer 1

让我们以您问题中的数据为例。您有一个坐标列表：

L = [[ 9.1514622, 47.1166004 ], [ 9.1513045, 47.1164599 ], [ 9.1516278, 47.1163001 ], [ 9.1517832, 47.1164408 ], [ 9.1514622, 47.1166004 ]]

然后转换成 Polygon:

P = Polygon(L)

现在，似乎 L 是多余的，因为最后一点与第一点相同。但这实际上不是问题，因为 Shapely 无论如何都会复制第一个点（以关闭多边形的边界）。您可以通过以下方式查看：

P = Polygon(L)
print(list(P.exterior.coords))
#[(9.1514622, 47.1166004), (9.1513045, 47.1164599), (9.1516278, 47.1163001), (9.1517832, 47.1164408), (9.1514622, 47.1166004)]

#now skip the last point
P = Polygon(L[:-1])
print(list(P.exterior.coords))
#[(9.1514622, 47.1166004), (9.1513045, 47.1164599), (9.1516278, 47.1163001), (9.1517832, 47.1164408), (9.1514622, 47.1166004)]

以防出现重复点"inside" L，例如：

L = [[ 9.1514622, 47.1166004 ], [ 9.1513045, 47.1164599 ], [ 9.1513045, 47.1164599 ], [ 9.1516278, 47.1163001 ], [ 9.1517832, 47.1164408 ], [9.1514622, 47.1166004 ]]

然后可以使用零容忍的simplify方法消除它（为了不引入side-effects）：

print(list(Polygon(L).simplify(0).exterior.coords))
#[(9.1514622, 47.1166004), (9.1513045, 47.1164599), (9.1516278, 47.1163001), (9.1517832, 47.1164408), (9.1514622, 47.1166004)]

从 Shapely 中的多边形中删除重复点

Remove a duplicate point from polygon in Shapely

python

polygon

geospatial

shapely