Het overlay-resultaat bekijken
Nu we de overlay van de gegevenssets met grondgebruik en arrondissementen hebben gemaakt, kunnen we het grondgebruik per arrondissement makkelijker onderzoeken. Laten we teruggaan naar het voorbeeldarrondissement Muette en het grondgebruik daar bekijken.
GeoPandas en Matplotlib zijn al geïmporteerd. Het resultaat van de functie overlay() uit de vorige oefeningen is beschikbaar als combined.
Deze oefening maakt deel uit van de cursus
Werken met georuimtelijke data in Python
Oefeninstructies
- Voeg een nieuwe kolom
'area'toe met de oppervlakte van elk polygoon aan de GeoDataFramecombined. - Maak een subset
land_use_muettewaarin'district_name'gelijk is aan "Muette". - Maak een plot van
land_use_muetteen gebruik de kolom'class'om de polygonen te kleuren. - Bereken de totale oppervlakte per
'class'vanland_use_muettemetgroupby()en print het resultaat.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# 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)