Granska överlappningsresultatet
Nu när vi har skapat överlappningen mellan markanvändnings- och distriktsdata kan vi enklare granska markanvändningen för de olika distrikten. Låt oss återgå till exempeldistriktet Muette och undersöka markanvändningen där.
GeoPandas och Matplotlib är redan importerade. Resultatet av funktionen overlay() från de föregående övningarna finns tillgängligt som combined.
Den här övningen är en del av kursen
Arbeta med geospatial data i Python
Övningsinstruktioner
- Lägg till en ny kolumn
'area'med arean för varje polygon i GeoDataFramecombined. - Skapa ett urval kallat
land_use_muettedär'district_name'är lika med "Muette". - Skapa ett diagram över
land_use_muetteoch använd kolumnen'class'för att färgsätta polygonerna. - Beräkna den totala arean för varje
'class'iland_use_muettemed hjälp av metodengroupby(), och skriv ut resultatet.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# 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)