需求

将MODIS的hdf文件，需要转换为tif文件。

简介

ExtractSubDataset(in_raster, out_raster, {subdataset_index})

• in_raster:包含要提取的图层的 HDF 或 NITF 数据集。
• out_raster:要创建的数据集的名称、位置和格式。
• subdataset_index:想要提取的子数据集。

代码

import os
import arcpy
arcpy.env.overwriteOutput = 1
arcpy.CheckOutExtension("Spatial")
inPath='D:\\GPP\\ecology\\mcd13q1\\raw\\'
outPath='D:\\GPP\\ecology\\mcd13q1\\test\\'
arcpy.env.workspace = inPath
hdfList = arcpy.ListRasters('*','HDF')
for hdf in hdfList:
    Name = hdf[0:13] + ".tif"
        # 提取波段0，对应波段与ArcGIS中显示的一致
    data = arcpy.ExtractSubDataset_management(hdf,outPath + Name, "0")
print "OK！"

参考博文

Python批处理多波段HDF文件为TIF文件
利用Python（ArcGIS环境下的ArcPy）批量将hdf数据特定波段输出为tif格式