修复Python中gdal.Warp导出MODIS HDF为TIFF的报错问题

错误根源 #

你遇到的TypeError: Warp() got multiple values for argument 'destNameOrDestDS'，是因为调用gdal.Warp时同时混用了位置参数和关键字参数传递输出路径，或者参数顺序完全颠倒，导致destNameOrDestDS这个必填参数被重复赋值。

修正后的完整脚本 #

以下是可以直接运行的优化代码，包含HDF子数据集提取、裁剪导出的完整流程：

import gdal
import glob
import os

# 配置GDAL避免编码和警告问题
gdal.SetConfigOption('GDAL_FILENAME_IS_UTF8', 'YES')
gdal.SetConfigOption('CPL_LOG', '/dev/null')

# 自定义路径（替换为你的实际路径）
INPUT_DIR = '/your/hdf/folder/path'
OUTPUT_DIR = '/your/tiff/output/folder'
CUT_SHAPE = '/your/extent.shp'

# 创建输出目录（自动处理不存在的情况）
if not os.path.exists(OUTPUT_DIR):
    os.makedirs(OUTPUT_DIR)

# 获取所有HDF文件
hdf_list = glob.glob(os.path.join(INPUT_DIR, '*.hdf'))

for hdf_path in hdf_list:
    # 生成输出文件名
    file_name = os.path.splitext(os.path.basename(hdf_path))[0]
    output_tif = os.path.join(OUTPUT_DIR, f"{file_name}_NDVI_clipped.tif")
    
    # 打开HDF并提取索引0的NDVI子数据集
    hdf_ds = gdal.Open(hdf_path)
    # GetSubDatasets()返回列表，每个元素是(子数据集路径, 描述)，取第一个（索引0）
    ndvi_subset = hdf_ds.GetSubDatasets()[0][0]
    ndvi_ds = gdal.Open(ndvi_subset)
    
    # 配置Warp裁剪参数
    warp_opts = gdal.WarpOptions(
        cutlineDSName=CUT_SHAPE,
        cropToCutline=True,
        dstNodata=-9999,  # MODIS NDVI默认无值，可根据产品调整
        format='GTiff',
        # 如果矢量和NDVI坐标系不一致，添加下面一行（示例为WGS84）
        # dstSRS='EPSG:4326'
    )
    
    # 正确调用gdal.Warp：第一个参数是输出路径，第二个是输入数据源，options传配置
    gdal.Warp(output_tif, ndvi_ds, options=warp_opts)
    
    # 释放内存资源
    ndvi_ds = None
    hdf_ds = None
    print(f"已完成：{output_tif}")

关键优化点说明 #

1. 子数据集正确读取：MODIS HDF是多数据集格式，必须通过GetSubDatasets()获取子数据集路径，再单独打开才能作为Warp的输入源。
2. Warp参数规范：gdal.Warp的第一个位置参数固定是输出路径（即destNameOrDestDS），第二个是输入数据源，所有额外配置通过options参数传递，避免参数冲突。
3. 资源管理：手动将数据源对象设为None，释放GDAL占用的内存，避免批量处理时内存泄漏。
4. 兼容性处理：添加坐标系匹配选项，如果你的裁剪矢量和MODIS NDVI（通常是正弦投影）坐标系不同，需指定dstSRS转换坐标系。

额外注意事项 #

• 先打印hdf_ds.GetSubDatasets()确认子数据集索引：部分MODIS产品的NDVI可能不是索引0，需根据实际输出调整。
• 批量处理时建议添加异常捕获，比如用try-except包裹循环内的代码，避免单个文件出错导致整个脚本终止。

内容的提问来源于stack exchange，提问作者tazrart