オーバーレイ結果を確認する

土地利用データと地区データのオーバーレイを作成したので、各地区の土地利用をより簡単に確認できるようになりました。例として使っていた Muette 地区に戻り、その地区の土地利用を見てみましょう。

GeoPandas と Matplotlib はすでにインポート済みです。前の演習で overlay() 関数から得た結果は combined として利用できます。

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

Pythonで扱う地理空間データ

コースを見る

演習の手順

各ポリゴンの面積を格納する新しい列 'area' を combined GeoDataFrame に追加します。
'district_name' が "Muette" に等しいサブセット land_use_muette を作成します。
'class' 列でポリゴンに色を付けて、land_use_muette をプロットします。
groupby() メソッドを使って、land_use_muette の各 'class' ごとの合計面積を計算し、結果を出力します。

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

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

# Print the first rows of the overlay result
print(combined.head())

# Add the area as a column
____ = ____

# Take a subset for the Muette district
land_use_muette = combined[____]

# Visualize the land use of the Muette district
land_use_muette.____(____)
plt.show()

# Calculate the total area for each land use class
print(land_use_muette.____(____)['area'].____() / 1000**2)

コードを編集して実行