Pythonを使用したラスターデータのゾナル統計分析方法

1. GDALとNumPyをインストールする:

   pip install gdal numpy

2. ラスターデータを読み込む:

   from osgeo import gdal
   raster_path = "path/to/your/raster.tif"
   raster_dataset = gdal.Open(raster_path)

3. ゾーンベクターデータを読み込む:

   import geopandas as gpd
   vector_path = "path/to/your/vector.shp"
   vector_data = gpd.read_file(vector_path)

4. ゾナル統計を計算する:

   from rasterstats import zonal_stats
   stats = zonal_stats(vector_data, raster_path, stats=["mean", "max", "min"])

   この例では、平均値、最大値、最小値の統計が計算されます。他の統計メトリックスも利用可能です。

5. 結果を確認する:

   for zone_stat in stats:
      print(zone_stat)

   各ゾーンの統計結果が表示されます。

これらの手順を実行することで、Pythonを使用してラスターデータのゾナル統計分析を行うことができます。適宜、必要なライブラリやデータのパスを修正してください。