基于arcpy处理 modis hdf格式转geotif

1：modis转tif.有一下几种方法

MRT(MODIS REPROJECTION TOOL)

HEG(The HDF-EOS To GeoTIFF Conversion Tool)

IDL调用MRT

2：利用arcgis中集成的arcpy工具包，实现hdf转tif格式

ExtractSubDataset_management (in_raster, out_raster, {subdataset_index})

参数	说明	数据类型
in_raster	包含要提取的图层的 HDF 或 NITF 数据集。	Raster Layer
out_raster	要创建的数据集的名称、位置和格式。 以文件格式存储栅格数据集时，需要指定文件扩展名，具体如下： .bil - Esri BIL .bip - Esri BIP .bmp - BMP .bsq - Esri BSQ .dat - ENVI DAT .gif - GIF .img - ERDAS IMAGINE 文件 .jpg - JPEG .jp2 - JPEG 2000 .png - PNG .tif - TIFF 无扩展名 - Esri GRID 将栅格数据集存储到地理数据库时，请勿向栅格数据集的名称添加文件扩展名。 将栅格数据集存储到 JPEG 文件、JPEG 2000 文件或地理数据库时，可在“环境设置”中指定压缩类型和压缩质量。	Raster Dataset
subdataset_index [subdataset_index,...] (可选)	想要提取的子数据集。	Value Table

import arcpy
import os
for dirpath,dirnames,filenames in os.walk(r'R:\study\wanghao\python\data_process\da'):
    for file in filenames:
        if file.endswith('.hdf'):
            arcpy.ExtractSubDataset_management("R:\study\wanghao\python\data_process\da\\"+file,
                                               "R:\study\wanghao\python\data_process\da_arcpy\\"+file[:-4]+".tif", "2;3;4")

上述代码提取modis某一产品的2，3，4数据子集。可以基于gdal或arcpy进行投影的添加，进行后续分析。

pyhdf包可以详细的获取到某产品包含的波段的详细信息。