如何获取 GeoDataFrame 中多边形内的所有点?

How to get all Points that are inside Polygons in GeoDataFrame?

我有两个 GeoDataFrame。一个包含点,另一个包含多边形。

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

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

我想知道是否有其他方法可以解决这个任务。

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

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)

使用经过优化的空间连接,速度应该很快。

See documentation for details.

points_within = gpd.sjoin(Points_gdf, Poly_gdf, how='inner', predicate='within')

how='inner' 可能是“可选的”(视情况而定)。