I have two shapefiles: Italian postcodes; and Italian municipalities. Both are polygon format and both contains the name of each relevant shape. I'd like to add to the postcode shapefile a field for the postcode's municipality. I believe this is a spatial query. How do I do this using GeoPandas?
1 Answers
Such a spatial query can be done with the spatial join function sjoin of GeoPandas. It will look something like:
res = geopandas.sjoin(postcodes, municipalities, op='contains')
where the op will depend on the specific spatial relationship of your data.
See https://geopandas.readthedocs.io/en/latest/mergingdata.html#spatial-joins for more.
Related