2つのポリゴンの交差
この演習では、2つの個別のポリゴンを使います。1つは districts データセットから抽出したムエット地区（Muette）、もう1つは land_use データセットから抽出したブローニュの都市緑地（パリ西部の大規模な公園）です。これら2つのポリゴンは、すでに muette と park_boulogne という変数に代入されています。
まず2つのポリゴンを可視化します。重なってはいますが、公園の全体がムエット地区内にあるわけではないことが分かります。重なっている部分を特定しましょう。
GeoPandas と matplotlib はすでにインポートされています。
この演習はコースの一部です
Pythonで扱う地理空間データ
演習の手順
- 2つのポリゴンを1つの地図にプロットし、重なり具合を視覚的に確認します。
park_boulogneとmuetteのポリゴンの交差を計算します。- 公園が地区の面積に占める割合を出力します。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Plot the two polygons
geopandas.GeoSeries([park_boulogne, muette]).plot(alpha=0.5, color=['green', 'blue'])
plt.show()
# Calculate the intersection of both polygons
intersection = ____
# Plot the intersection
geopandas.GeoSeries([intersection]).plot()
plt.show()
# Print proportion of district area that occupied park
print(____ / ____)