基于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包可以详细的获取到某产品包含的波段的详细信息。

转载自:https://blog.csdn.net/qq_33657870/article/details/83345495

You may also like...