can anyone tell me why unary_union
on a sequence of 5 or more Polygons returns an empty GeometryCollection
:
import numpy as np
from shapely.ops import unary_union
from shapely.geometry import Polygon
x = np.array([0, 0, 1, 1, 0], dtype=float)
y = np.array([0, 1, 1, 0, 0], dtype=float)
polygons = []
for i in range(5):
polygons.append(Polygon(list(zip(x, y))))
x += 1
un = unary_union(polygons)
while 4 or less Polygons returns a Polygon as expected?
polygons = []
for i in range(4):
polygons.append(Polygon(list(zip(x, y))))
x += 1
un2 = unary_union(polygons)