python中等效的st_make_grid方法

Question

是否有与 r-spatial in python 中的sf 包的非常好的st_make_grid 方法等效？该方法在多边形的边界框上创建矩形网格几何体。

我想做与this 问题中提出的解决方案完全相同的操作，例如将一个多边形分成几个面积相同的正方形，我选择。感谢您的帮助。

---

或者，我可以使用rpy2 在 r 中运行一个脚本，该脚本执行st_make_grid 方法，该方法将一个匀称的多边形作为输入并输出方形多边形，以便用匀称的方式读取。这对许多要处理的多边形有效吗？

Answer 1

这对许多要处理的多边形有效吗？

当然不是。没有内置的 Python 版本，但下面的函数可以解决问题。如果您需要性能，请确保在您的环境中安装了pygeos。

def make_grid(polygon, edge_size):
    """
    polygon : shapely.geometry
    edge_size : length of the grid cell
    """
    from itertools import product
    import numpy as np
    import geopandas as gpd
    
    bounds = polygon.bounds
    x_coords = np.arange(bounds[0] + edge_size/2, bounds[2], edge_size)
    y_coords = np.arange(bounds[1] + edge_size/2, bounds[3], edge_size)
    combinations = np.array(list(product(x_coords, y_coords)))
    squares = gpd.points_from_xy(combinations[:, 0], combinations[:, 1]).buffer(edge_size / 2, cap_style=3)
    return gpd.GeoSeries(squares[squares.intersects(polygon)])