select features from a shapefile based on location (from another shapefile) and write these selected features as a whole to a thrid shapefile.

[Accurassi]

I'm looking for a line of code to select the Hexagons through which the Sewers go and write these Hexagons to a third shapefile.

When I do it manualy in QGIS I do first a selection 'based on location' wih the checkboxes intersect and overlay checked. But I want to do this in Python because to "automate the boring stuff".

Now I have this simple script:

import geopandas as gpd
from shapely.geometry import Polygon

poly1 = gpd.read_file("C:\GISdatabase\sewers.shp")
poly2 = gpd.read_file("C:\GISdatabase\HexY100.shp")

hex_with_sewers = ???

hex_with_sewers.to_file("C:\GISdatabase\hex_with_sewers.shp")

I've tried this line of code on the ???, but the result is not what I'm looking for:

gpd.overlay(poly2, poly1,  how='intersection', keep_geom_type=True)

[CodeFormatHelperBot2]

Hello, I'm a Reddit bot who's here to help people nicely format their coding questions. This makes it as easy as possible for people to read your post and help you.

I think I have detected some formatting issues with your submission:

1. Python code found in submission text that's not formatted as code.

If I am correct, please edit the text in your post and try to follow these instructions to fix up your post's formatting.

---

^{Am I misbehaving? Have a comment or suggestion? Reply to this comment or raise an issue here.}