I've been managing within QGIS two shapefiles: the first is a spatial square grid, and the second a city region. Now I would like to intersect the grid keeping only the cells that are overlapping the city, returning the whole cells on the boundaries. So far I used the Intersect function which returns this:
Is there a way (within QGIS or Python, for instance) to keep the whole cells overlapping the boundaries of the blue region?