Ottenere i conteggi da un join spaziale

Userai un insieme di dati dei permessi edilizi rilasciati a Nashville nel 2017. Questo DataFrame, chiamato permits, è già stato caricato per te, insieme al GeoDataFrame council_districts.

Questo esercizio fa parte del corso

Visualizzare dati geospaziali in Python

Visualizza corso

Istruzioni dell'esercizio

Crea una colonna geometry in permits a partire da lat e lng.
Crea permits_geo, un GeoDataFrame, usando permits, il CRS di council_districts e la geometry da permits.
Usa un join spaziale per trovare i permessi rilasciati within ciascun distretto comunale. Stampa le prime 2 righe.
Crea permit_counts per mostrare il conteggio dei permessi all'interno di ogni distretto, usando groupby() e .size(). Stampa permit_counts.

esercizio interattivo pratico

Prova questo esercizio completando questo codice di esempio.

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)

Modifica ed esegui il codice