空間結合から件数を求める

2017 年に Nashville で発行された建築許可のデータセットを使用します。permits という DataFrame は、council_districts GeoDataFrame とあわせて、あらかじめ読み込まれています。

この演習はコースの一部です

Pythonで可視化する地理空間データ

コースを見る

演習の手順

lat と lng から permits に geometry 列を作成します。
permits、council_districts の crs、そして permits の geometry を使って、GeoDataFrame の permits_geo を作成します。
空間結合を使って、各 council district の within に発行された許可を取得します。先頭 2 行を表示してください。
groupby() と .size() を使って、各 district 内の許可件数を示す permit_counts を作成します。permit_counts を表示してください。

実践的なインタラクティブ演習

このサンプルコードを完成させて、この演習に挑戦してみましょう。

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)

コードを編集して実行