I'm trying to add the intersected geometries of 2 layers in a pyqgis code.
I have a grid made up by cells (so a polygon) and a line layer. I want to select the features of the grid layer where they intersect the line features (so basically the Select by Location tool used within a python script).
The code (that works) is the following:
# get the layers from the interface grid = iface.activeLayer() line = iface.activeLayer() for i in grid.getFeatures(): for l in line.getFeatures(): if i.geometry().intersects(l.geometry()): print i
It returns (correctly) 17 features.
What I'm not able to to is to add these QgsFeatures as a selection to the grid layer.