如何在地理数据框中获取多边形内的所有点?

Question

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

我需要从地理数据帧获取任何多边形内的所有点。

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

我想知道是否有另一种方法可以解决此任务。

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

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)

Answer 1

使用空间联接，这是经过优化的，应该很快。

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

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

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