Intersección de dos polígonos
En este ejercicio, vamos a usar 2 polígonos individuales: el distrito de Muette extraído del conjunto de datos districts, y la zona verde urbana de Boulogne, un gran parque público al oeste de París, extraída del conjunto de datos land_use. Los dos polígonos ya se han asignado a las variables muette y park_boulogne.
Primero vamos a visualizar ambos polígonos. Verás que se superponen, pero el parque no está completamente dentro del distrito de Muette. Vamos a determinar la parte que se solapa.
GeoPandas y matplotlib ya están importados.
Este ejercicio forma parte del curso
Trabajo con datos geoespaciales en Python
Instrucciones del ejercicio
- Representa los dos polígonos en un único mapa para examinar visualmente el grado de solapamiento.
- Calcula la intersección de los polígonos
park_boulogneymuette. - Imprime la proporción del área del distrito que ocupa el parque.
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# Plot the two polygons
geopandas.GeoSeries([park_boulogne, muette]).plot(alpha=0.5, color=['green', 'blue'])
plt.show()
# Calculate the intersection of both polygons
intersection = ____
# Plot the intersection
geopandas.GeoSeries([intersection]).plot()
plt.show()
# Print proportion of district area that occupied park
print(____ / ____)