Question

我试图创建一个numpy meshgrid并将其转换为Shapely多边形。我可以用一种非常强力的方法解决这个问题，但感觉必须要有一个很好的技巧才能实现这个目标，但我还没有想到它。

这让我获得了积分点（假设在Jupyter中运行） -

import numpy as np
from matplotlib import pyplot

fig = pyplot.figure(figsize=(10, 10))
ax = fig.add_subplot(111, aspect='equal')

x,y = np.mgrid[-5:-1:8j, 1:5:8j]
ax.plot(x,y, 'o', color='#000000')
pyplot.show()

现在需要水平和垂直连接所有这些点以形成Shapely多边形。我对此的第一次尝试是生成一个Shapely MultiLineString来绘制垂直和水平线，然后对其执行多边形操作。这导致仅创建主外部多边形 - 这是由于MultiLineString仅包含外部多边形上的顶点。

我知道使用栅格和GDAL可能更合理，但我的情况要求最终结果是Shapely多边形。

感谢您追踪解决方案的任何帮助！

Answer 1

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

import numpy as np
from shapely.geometry import MultiLineString
from shapely.ops import polygonize

x = 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)))

numpy meshgrid到Shapely多边形

1 个答案: