Obtenir des décomptes à partir d’une jointure spatiale

Vous allez utiliser un jeu de données sur les permis de construire délivrés à Nashville en 2017. Ce DataFrame, appelé permits, a été préchargé pour vous, ainsi que le GeoDataFrame council_districts.

Cet exercice fait partie du cours

Visualiser des données géospatiales en Python

Afficher le cours

Instructions

Créez une colonne geometry dans permits à partir de lat et lng.
Créez permits_geo, un GeoDataFrame, en utilisant permits, le CRS de council_districts, et la geometry de permits.
Utilisez une jointure spatiale pour trouver les permis délivrés within chaque district municipal. Affichez les 2 premières lignes.
Créez permit_counts pour afficher le nombre de permis dans chaque district, avec groupby() et .size(). Affichez permit_counts.

Exercice interactif pratique

Essayez cet exercice en complétant cet exemple de code.

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)

Modifier et exécuter le code