ラスターレイヤーから情報を抽出する
次は、ベクターファイルに基づいてラスターレイヤーから情報を抽出します。この機能は rasterstats パッケージで提供されています。今回の演習では、採掘サイトデータセットの各地点に最も近いラスターピクセル値を取得し、すべての採掘サイトの植生タイプを判定します。
採掘サイトデータセットのサブセット（mining_sites）はすでに読み込まれており、GeoPandas と matplotlib もインポート済みです。
この演習はコースの一部です
Pythonで扱う地理空間データ
演習の手順
rasterstatsパッケージをインポートします。rasterstats.point_query()関数を使って、採掘サイトの位置（mining_sites）に対する最も近いラスタ値を抽出し、新しい列'vegetation'に代入します。- 採掘サイトデータを
'vegetation'列で色分けしてプロットします。凡例も必ず追加してください。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Import the rasterstats package
____
# Extract the nearest value in the raster for all mining sites
vegetation_raster = "central_africa_vegetation_map_foraf.tif"
mining_sites['vegetation'] = ____(____, ____, interpolate='nearest')
print(mining_sites.head())
# Replace numeric vegation types codes with description
mining_sites['vegetation'] = mining_sites['vegetation'].replace(vegetation_types)
# Make a plot indicating the vegetation type
____
plt.show()