Rumslig join och räkning av tillstånd

Du fortsätter att förbereda din datamängd för att rita ett geopandas-choropleth genom att skapa en GeoDataFrame med bygglov som rumsligt joinats till kommundelsdistrikten. Därefter kan du räkna hur många bygglov som utfärdats i varje distrikt.

Den här övningen är en del av kursen

Visualisering av geospatial data i Python

Visa kurs

Övningsinstruktioner

Skapa permits_geo från permits, council_districts.crs och geometrin i permits.
Gör en rumslig join av permits_geo och council_districts för att hämta bygglov within varje kommundistrikt. Kalla resultatet permits_by_district.
Räkna tillstånden i varje distrikt, permit_counts, genom att kedja metoderna groupby() och size().
Skapa counts_df från permit_counts. Återställ indexet och namnge kolumnerna district och bldg_permits.

Interaktiv övning med praktiskt arbete

Testa den här övningen genom att slutföra den här exempelkoden.

# 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))

Redigera och kör kod