I would like to select polygons that contain at least a point. I can use QGIS's tool called "Select by location: Select all buildings by location of point".
Is there a python alternative?
So far, I wrote a jupyter notebook and worked with GeoPandas.
I have tried
import pandas as pd polygon_layer = gpd.read_file(r'file.shp') excel = pd.read_excel('file2.xlsx') points_layer = gpd.GeoDataFrame(excel, geometry=gpd.points_from_xy(excel.X, excel.Y)) subset = gpd.sjoin(polygon_layer, points_layer, how='inner', predicate='within')
However, the spatial join does not work.