I have these two shapefiles, one with polygon geometry and the other with point geometry, that I am working on in QGIS. I want to create a script using Python to see the intersection of points and polygons, and where there is an intersection fill in the value of a field I created.

Psuedo code for this would be something like

for polygon in Aveiro_Polygon.features():
    for point in Aveiro.features():
        if point is contained within polygon:               

Shps Point and Polygon Table Of point

  • what you are looking for is called "spatial join" – Mesa Feb 17 '19 at 9:38

Browse other questions tagged or ask your own question.