Ispezionare il risultato dell'overlay
Ora che abbiamo creato l'overlay degli insiemi di dati sull'uso del suolo e sui distretti, possiamo ispezionare più facilmente l'uso del suolo nei diversi distretti. Torniamo all'esempio del distretto di Muette e analizziamone l'uso del suolo.
GeoPandas e Matplotlib sono già importati. Il risultato della funzione overlay() degli esercizi precedenti è disponibile come combined.
Questo esercizio fa parte del corso
Lavorare con i dati geospaziali in Python
Istruzioni dell'esercizio
- Aggiungi una nuova colonna
'area'con l'area di ciascun poligono al GeoDataFramecombined. - Crea un sottoinsieme chiamato
land_use_muettein cui'district_name'è uguale a "Muette". - Crea un grafico di
land_use_muette, usando la colonna'class'per colorare i poligoni. - Calcola l'area totale per ciascuna
'class'diland_use_muetteusando il metodogroupby()e stampa il risultato.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# Print the first rows of the overlay result
print(combined.head())
# Add the area as a column
____ = ____
# Take a subset for the Muette district
land_use_muette = combined[____]
# Visualize the land use of the Muette district
land_use_muette.____(____)
plt.show()
# Calculate the total area for each land use class
print(land_use_muette.____(____)['area'].____() / 1000**2)