numpy meshgrid到Shapely多边形

numpy meshgrid到Shapely多边形

在构造MultiLineString之前，基本上必须定义每行。

import numpy as npfrom shapely.geometry import MultiLineStringfrom shapely.ops import polygonizex = np.linspace(-5, -1, 8)y = np.linspace(1, 5, 8)hlines = [((x1, yi), (x2, yi)) for x1, x2 in zip(x[:-1], x[1:]) for yi in y]vlines = [((xi, y1), (xi, y2)) for y1, y2 in zip(y[:-1], y[1:]) for xi in x]grids = list(polygonize(MultiLineString(hlines + vlines)))