Removing Duplicate Geometries In Shapely

I have a list of Shapely polygons. From that list I want to extract only unique polygons removing the duplicates. How to do it in a faster way? (My list contains thousands of polyg

Solution 1:

Note that you should not use intersects(), as that will identify any polygons that overlap at all as duplicates. Use equals() instead.

You can create a list that stores the unique polygons, and then for each polygon in your list, loop over the polygons stored in the outer list, and if none of them are equal to the new one, add this to the list, you can use any() function for that.

Example:

from shapely.geometry import Polygon

lists = [[(1,1),(2,2),(3,3),(4,4)], [(6,6),(7,7),(8,8),(9,9)], [(1,1),(2,2),(3,3),(4,4)]]
polys = [Polygon(item) for item in lists] ##This is given condition
uniqpolies = []
for poly in polys:

    ifnot any(p.equals(poly) for p in uniqpolies):
        uniqpolies.append(poly)

From documentation on any():

Return True if any element of the iterable is true. If the iterable is empty, return False.