时间戳

首页 > 解决方案 > 如何获取 GeoDataFrame 中多边形内的所有点？

python - 如何获取 GeoDataFrame 中多边形内的所有点？

问题描述

我有两个 GeoDataFrames。一个包含点，另一个包含多边形。

我需要从 GeoDataFrame 获取任何多边形内的所有点。

我尝试遍历所有多边形并检查 Point 是否在该多边形内。此解决方案有效，但速度很慢。

我想知道是否有另一种方法来解决这个任务。

编辑：我的解决方案如下所示：

for i in range(len(Poly_gdf.index)):
    inter = Points_gdf[Points_gdf.intersects(Poly_gdf.loc[i,'geometry'])]
    if not inter.empty:
        for i in inter['geometry'].values:
            points.append(i)

标签： pythongeometrygeopandasshapely

解决方案

使用空间连接，它经过优化并且应该很快。

有关详细信息，请参阅文档。

points_within = gpd.sjoin(Points_gdf, Poly_gdf, op='within')

推荐阅读