MODIS数据通常是HDF-EOS（H4）文件，文件包含有多个数据集，如图 1。已知MOD13A3.hdf，提取其中的NDVI数据集，并以GeoTiff格式存储。参考GDAL方法，但没有找到Python+GDAL代码实例，所以网上介绍的方法没有测试成功，此外还有国外的专家制作的PyHDF软件，可以尝试。

我的代码有一个缺陷，不能自HDF文件获取其中地理参考。不过，我想到一个替代方法，MODIS产品附有行列信息（h**v**），将某一行列的文件在ENVI中另存为GeoTiff文件（GEOREF.tif），IDL读取这个文件地理参考再存储在输出的NDVI数据集中。

这算是一个折中的办法，代码中对应GeoRef参数是为包含地理参考。

附上代码（code.rar）。

图 1