基于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 |
要创建的数据集的名称、位置和格式。 以文件格式存储栅格数据集时,需要指定文件扩展名,具体如下:
将栅格数据集存储到地理数据库时,请勿向栅格数据集的名称添加文件扩展名。 |
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包可以详细的获取到某产品包含的波段的详细信息。
转载自:https://blog.csdn.net/qq_33657870/article/details/83345495