Intersecciones con método 'overlay' de GeoPandas

El método 'overlay' de GeoPandas permite realizar intersecciones (y otras operaciones de geoproceso) en la misma capa (auto intersección) o en capas diferentes; aunque en la misma capa puede ser difícil de interpretar. Por tanto, para comenzar, se va a realizar la operación con los shapefiles de la imagen siguiente:

El código es tan sencillo como:

1 2 3 4 5 6 7

import geopandas as gpd pol8 = gpd.GeoDataFrame.from_file("pyqgis_data/polygon8.shp") pol3 = gpd.GeoDataFrame.from_file("pyqgis_data/polygon3.shp") intersections = gpd.overlay(pol8, pol3, how='intersection')

Para graficar intersecciones con matplotlib se tiene que:

1 2 3 4

>>>from matplotlib import pyplot as plt >>>plt.ion() >>>intersections.plot() <matplotlib.axes._subplots.AxesSubplot object at 0x9c1063ec>

resultando en:

Comparando con la primera imagen se puede observar que los resultados son los esperados.