Een Landgebruiksgegevensset verkennen
Voor de volgende oefeningen introduceren we eerst een nieuwe gegevensset: een gegevensset over het landgebruik van Parijs (een vereenvoudigde versie op basis van de open Europese Urban Atlas). Landgebruik geeft aan voor welk type activiteit een bepaald gebied wordt gebruikt, zoals wonen of recreatie. Het is een polygoongegevensset, met een label dat de landgebruiksklasse weergeeft voor verschillende gebieden in Parijs.
In deze oefening lezen we de data in, verkennen we die visueel en berekenen we de totale oppervlakte van de verschillende landgebruiksklassen in het gebied van Parijs.
GeoPandas en matplotlib zijn al geïmporteerd.
Deze oefening maakt deel uit van de cursus
Werken met georuimtelijke data in Python
Oefeninstructies
- Lees het bestand
'paris_land_use.shp'in en ken het resultaat toe aan de variabeleland_use. - Maak een plot van
land_use, en gebruik de kolom'class'om de polygonen te kleuren. Voeg ook een legenda toe. Let op: het kan een paar seconden duren voordat de plot klaar is, omdat er veel polygonen zijn. - Voeg een nieuwe kolom
'area'toe met de oppervlakte van elke polygoon. - Bereken de totale oppervlakte in km² voor elke
'class'met de methodegroupby()en print het resultaat.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Import the land use dataset
land_use = geopandas.____
print(land_use.head())
# Make a plot of the land use with 'class' as the color
land_use.plot(____, legend=True, figsize=(15, 10))
plt.show()
# Add the area as a new column
land_use['area'] = ____
# Calculate the total area for each land use class
total_area = land_use.groupby(____)['area'].____() / 1000**2
print(total_area)