MODIS HDF格式数据处理

由于MODIS数据往往会以hdf格式数据提供，以MCD12Q1数据（https://lpdaac.usgs.gov/products/mcd12q1v006/ ） 为例。批量提取特定波段 （IGBP）分类数据

#coding=utf-8
import arcpy
import os

def file_name(file_dir):
    FileIm_list = []
    FileEx_list = []
    for root, dirs, files in os.walk(file_dir):
        for file in files:
            FileIm_list.append(file)
            fileSub = file[17:23]
            FileEx_list.append(fileSub)  
    return FileIm_list,FileEx_list
Flist = file_name("D:\\0Data_hdf")#hdf文件存放路径

for i in range(0,len(Flist[0])):
    print(i)
    print(Flist[0][i])
    #参数1：待提取的hdf数据；参数2：保存路径；参数3：提取的第几个波段
    arcpy.ExtractSubDataset_management("D:\\0Data_hdf\\%s"%(Flist[0][i]),
                                              "D:\\1HDF2IGBP\\%s"%(Flist[1][i]),"0")   
print("Done")