Intersecare un GeoDataFrame con un Poligono
Combinando l'insieme di dati sull'uso del suolo e quello dei distretti, possiamo ora indagare qual è l'uso del suolo in un determinato distretto.
Per farlo, dobbiamo prima determinare l'intersezione tra l'insieme di dati sull'uso del suolo e un distretto specifico. Prendiamo di nuovo il distretto della Muette come esempio.
Gli insiemi di dati sull'uso del suolo e sui distretti sono già stati importati come land_use e districts, e il distretto Muette è stato estratto nel poligono shapely muette. Inoltre, GeoPandas e matplotlib sono importati.
Questo esercizio fa parte del corso
Lavorare con i dati geospaziali in Python
Istruzioni dell'esercizio
- Calcola l'intersezione dei poligoni di
land_usecon il singolo poligonomuette. Chiama il risultatoland_use_muette. - Fai un grafico veloce di questa intersezione e passa
edgecolor='black'per vedere più chiaramente i confini dei diversi poligoni. - Stampa le prime cinque righe di
land_use_muette.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# 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.____)