[1054] Select only the records from one GeoSeries that intersect with the polygons from another GeoSeries

To select only the records from one GeoSeries that intersect with the polygons from another GeoSeries in GeoPandas, you can use the intersects method along with boolean indexing. Here’s a step-by-step guide:

1. Import GeoPandas:

   import geopandas as gpd

2. Load your GeoSeries:

   # Assuming you have two GeoDataFrames, gdf1 and gdf2
   gdf1 = gpd.read_file('path_to_your_first_shapefile.shp')
   gdf2 = gpd.read_file('path_to_your_second_shapefile.shp')

3. Check for intersections:

   # Create a boolean series indicating whether each geometry in gdf1 intersects with any geometry in gdf2
   intersects = gdf1.geometry.apply(lambda x: gdf2.geometry.intersects(x).any())

4. Filter the GeoSeries:

   # Select only the records from gdf1 that intersect with gdf2
   intersecting_records = gdf1[intersects]

Here’s the complete code snippet:

import geopandas as gpd
 
# Load your GeoDataFrames
gdf1 = gpd.read_file('path_to_your_first_shapefile.shp')
gdf2 = gpd.read_file('path_to_your_second_shapefile.shp')
 
# Create a boolean series indicating intersections
intersects = gdf1.geometry.apply(lambda x: gdf2.geometry.intersects(x).any())
 
# Filter the GeoDataFrame
intersecting_records = gdf1[intersects]
 
# Display the result
print(intersecting_records)

This code will give you a GeoDataFrame containing only the records from gdf1 that intersect with any polygon in gdf2.

Would you like more details on any specific part of this process?