Aantallen vinden via een ruimtelijke join

Je gaat een gegevensset gebruiken met de bouwvergunningen die in 2017 in Nashville zijn afgegeven. Deze DataFrame, permits, is alvast voor je ingeladen, samen met de GeoDataFrame council_districts.

Deze oefening maakt deel uit van de cursus

Geospatiale data visualiseren in Python

Bekijk cursus

Oefeninstructies

Maak een geometry-kolom in permits op basis van lat en lng.
Maak permits_geo, een GeoDataFrame, met permits, de crs van council_districts en de geometry uit permits.
Gebruik een ruimtelijke join om vergunningen te vinden die within elk district vallen. Print de eerste 2 rijen.
Maak permit_counts om het aantal vergunningen per district te tonen, met groupby() en .size(). Print permit_counts.

Interactieve oefening met praktijkervaring

Probeer deze oefening door deze voorbeeldcode aan te vullen.

import geopandas as gpd

# Create a geometry column in permits from lat and lng
permits[____] = gpd.points_from_xy(____.____ , ____.____)

# Build a GeoDataFrame: permits_geo
permits_geo = gpd.GeoDataFrame(____, crs = ____.crs, geometry = ____.geometry)

# Spatial join of permits_geo and council_districts
permits_by_district = gpd.____(permits_geo, council_districts, predicate = ____)
print(permits_by_district.head(2))

# Create permit_counts
permit_counts = permits_by_district.groupby([____]).____()
print(permit_counts)

Code bewerken en uitvoeren