import arcpy import os arcpy.env.overwriteOutput = 1 arcpy.CheckOutExtension("Spatial") inPath = 'D:\\PyCharm\\program\got\\20190101\\' ls = os.listdir(inPath) # print(ls) print(len(ls)) for i in ls: arcpy.env.workspace = inPath + i arcpy.env.scratchWorkspace = inPath + i hdfList = arcpy.ListRasters('*', 'H') if os.path.exists(r"D:\\FY_tif\\" + i): for hdf in hdfList: for number in range(0, 20): outPath = 'D:\\FY_tif\\' + str(i) + '\\' out = arcpy.ExtractSubDataset_management(hdf, outPath + str(i) + str(number) + ".tif", number) print(outPath) else: os.makedirs(r"D:\\FY_tif\\" + i) for hdf in hdfList: for number in range(0, 20): outPath = 'D:\\FY_tif\\' + str(i) + '\\' out = arcpy.ExtractSubDataset_management(hdf, outPath + str(i) + str(number) + ".tif", number) print(outPath)
HDF文件转TIF文件多波段提取
最新推荐文章于 2024-07-13 21:34:41 发布