Ruimtelijk joinen en aantallen ophalen

Je gaat verder met het voorbereiden van je gegevensset voor het plotten van een geopandas-choropleth door een GeoDataFrame te maken van de bouwvergunningen die ruimtelijk zijn gejoint met de raadsdistricten. Daarna kun je het aantal bouwvergunningen per raadsdistrict bepalen.

Deze oefening maakt deel uit van de cursus

Geospatiale data visualiseren in Python

Bekijk cursus

Oefeninstructies

Maak permits_geo op basis van permits, de council_districts.crs en de geometry in permits.
Voer een ruimtelijke join uit tussen permits_geo en council_districts om bouwvergunningen te krijgen die within elk raadsdistrict liggen. Noem dit permits_by_district.
Tel vergunningen in elk district, permit_counts, door de methoden groupby() en size() te koppelen.
Maak counts_df op basis van permit_counts. Reset de index en geef de kolommen de namen district en bldg_permits.

Interactieve oefening met praktijkervaring

Probeer deze oefening door deze voorbeeldcode aan te vullen.

# Create permits_geo
permits_geo = gpd.GeoDataFrame(____, crs = ____.____, geometry = ____.____)

# Spatially join permits_geo and council_districts
permits_by_district = gpd.sjoin(____, ____, ____ = ____)
print(permits_by_district.head(2))

# Count permits in each district
permit_counts = permits_by_district.____(____).____()

# Convert permit_counts to a df with 2 columns: district and bldg_permits
counts_df = ____.to_frame()
counts_df = counts_df.____()
counts_df.____ = ['district', 'bldg_permits']
print(counts_df.head(2))

Code bewerken en uitvoeren