ポリゴンと GeoDataFrame の交差をとる
土地利用データセットと地区データセットを組み合わせることで、特定の地区での土地利用を調べられるようになります。
そのために、まずは土地利用データセットと指定した地区との交差を求めます。例として、再び「Muette」地区を使いましょう。
土地利用と地区の各データセットはそれぞれ land_use と districts としてすでに読み込まれており、Muette 地区は muette という shapely のポリゴンに抽出済みです。さらに、GeoPandas と matplotlib もインポートされています。
この演習はコースの一部です
Pythonで扱う地理空間データ
演習の手順
land_useの各ポリゴンと単一のmuetteポリゴンとの交差を計算し、結果をland_use_muetteと名付けてください。- この交差結果を手早くプロットし、異なるポリゴンの境界が見やすいように
edgecolor='black'を指定してください。 land_use_muetteの先頭5行を出力してください。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Print the land use datset and Notre-Dame district polygon
print(land_use.head())
print(type(muette))
# Calculate the intersection of the land use polygons with Notre Dame
land_use_muette = ____
# Plot the intersection
land_use_muette.____
plt.show()
# Print the first five rows of the intersection
print(land_use_muette.____)