利用arcpy在arcgis中批量将hdf数据转tif数据

最新推荐文章于 2024-04-02 23:17:25 发布

weixin_38674806

最新推荐文章于 2024-04-02 23:17:25 发布

阅读量1k

点赞数 2

文章标签： arcgis

本文链接：https://blog.csdn.net/weixin_38674806/article/details/128216907

版权

用到的软件：Arcgis

用到的数据：MODIS数据【MOD15A2H.A2021001.h25v03.061.2021020165603.hdf】

程序内容

import arcpy
import os
from arcpy import env


arcpy.env.overwriteOutput = 1
arcpy.CheckOutExtension("Spatial")
inPath='F:/FWR/test/'
outPath='F:/FWR/test/'
arcpy.env.workspace = inPath
arcpy.env.scratchWorkspace = inPath
hdfList = arcpy.ListRasters('*','HDF')
for hdf in hdfList:
    for number in range(1,3): #注意角标。hdf文件角标由0开始，相应波段与ArcGIS中显示一致。
        eviName=hdf[1:43] + "_" + str(number) + ".tif"
        data1=arcpy.ExtractSubDataset_management(hdf,outPath + eviName,number)
        print(eviName)
print "完成"

部分代码介绍：

arcpy.env.overwriteOutput = 1 #输出文件夹里面已经有内容的,就覆盖掉arcpy.CheckOutExtension("Spatial")#打开arcgis扩展模块

核心代码

arcpy.ExtractSubDataset_ma