有没有相当于非常好的st_make_grid的方法sf来自 r-spatial 的包在Python中?该方法在多边形的边界框上创建矩形网格几何体。
我想做与中提出的解决方案完全相同的事情this问题,例如将一个多边形分成若干个面积相同的正方形I choose。感谢您的帮助。
或者,我可以使用rpy2在 r 中运行一个脚本来执行st_make_grid
方法,以形状多边形作为输入并输出正方形多边形,以便用形状读取。这对要处理的许多多边形有效吗?
这对要处理的许多多边形有效吗?
当然不是。没有内置的 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)])
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)