使用GADL批量裁剪栅格数据

要使用GADL批量根据矢量边界裁剪栅格数据，可以按照以下步骤进行操作：

1. 安装GDAL库。您可以在GDAL官网上下载对应的安装包，并按照安装说明进行安装。

2. 在命令行中使用gdalwarp命令进行批量裁剪。命令格式如下：

   gdalwarp -cutline <shapefile> -crop_to_cutline <input_raster> <output_raster>
 

   其中，`<shapefile>`是您的矢量边界文件，`<input_raster>`是您要裁剪的栅格数据，`<output_raster>`是输出的裁剪后的栅格数据。

3. 如果您需要批量处理多个栅格数据，可以使用Python脚本来自动化处理。以下是一个简单的示例脚本：

   
   import os

   # 设置输入输出路径和矢量边界文件
   input_dir = "/path/to/input/rasters"
   output_dir = "/path/to/output/rasters"
   shapefile = "/path/to/shapefile.shp"

   # 遍历输入路径下的所有栅格数据
   for filename in os.listdir(input_dir):
       if filename.endswith(".tif"):
           # 构造输入输出路径
           input_raster = os.path.join(input_dir, filename)
           output_raster = os.path.join(output_dir, filename)

           # 执行裁剪
           os.system(f"gdalwarp -cutline {shapefile} -crop_to_cutline {input_raster} {output_raster}")
 

   可以将上述代码保存为一个Python脚本，并将`<shapefile>`、`<input_dir>`和`<output_dir>`替换为自己的路径。然后在命令行中运行该脚本即可自动化批量裁剪。