我有两个地理数据框。一个包含点,另一个包含多边形。
我需要从 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)
Run Code Online (Sandbox Code Playgroud)