空間結合と件数の集計

geopandas のコロプレス図を作成する準備として、建築許可データを議会地区に空間結合した GeoDataFrame を作成します。これができれば、各議会地区で発行された建築許可の件数を集計できるようになります。

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

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

演習の手順

permits、council_districts.crs、および permits の geometry を使って permits_geo を作成します。
permits_geo と council_districts を空間結合し、各議会地区の内部（within）にある建築許可を取得します。結果は permits_by_district とします。
groupby() と size() をメソッドチェーンして、各地区の許可件数 permit_counts を数えます。
permit_counts から counts_df を作成します。インデックスをリセットし、列名を district と bldg_permits にします。

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

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

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

コードを編集して実行