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

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

Cette activité fait partie du cours

Visualisation de données géospatiales avec Python

Voir le cours

Instructions de l’exercice

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 émis within chaque district de conseil. Affichez les 2 premières lignes.
Créez permit_counts pour afficher le nombre de permis dans chaque district, en utilisant groupby() et .size(). Affichez permit_counts.

Exercice interactif pratique

Essayez cet exercice en complétant ce code d’exemple.

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