Intersecar un GeoDataFrame con un Polígono
Combinando el conjunto de datos de uso del suelo y el de distritos, ahora podemos investigar cuál es el uso del suelo en un distrito concreto.
Para ello, primero necesitamos determinar la intersección del conjunto de datos de uso del suelo con un distrito dado. Tomemos de nuevo el distrito de Muette como caso de ejemplo.
Los conjuntos de datos de uso del suelo y de distritos ya se han importado como land_use y districts, y el distrito Muette se ha extraído en el polígono shapely muette. Además, GeoPandas y matplotlib ya están importados.
Este ejercicio forma parte del curso
Trabajo con datos geoespaciales en Python
Instrucciones del ejercicio
- Calcula la intersección de los polígonos de
land_usecon el polígono únicomuette. Llama al resultadoland_use_muette. - Haz un gráfico rápido de esta intersección y pasa
edgecolor='black'para ver más claramente los límites de los distintos polígonos. - Imprime las primeras cinco filas de
land_use_muette.
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# Print the land use datset and Notre-Dame district polygon
print(land_use.head())
print(type(muette))
# Calculate the intersection of the land use polygons with Notre Dame
land_use_muette = ____
# Plot the intersection
land_use_muette.____
plt.show()
# Print the first five rows of the intersection
print(land_use_muette.____)