python调用MRT批处理MODIS数据.rar

# encoding:utf-8 import os os.environ['MRT_DATA_DIR'] = r'C:\ModisTools\MRT\data' # MRT安装路径下 mrtpath = r'C:\ModisTools\MRT\bin' mrtmosaic = os.path.join(mrtpath, 'mrtmosaic.exe') mrtmosaic = mrtmosaic.replace('\\', '/') datadir = r'。。。\hdf' # 原始HDF文件所在路径 datadir = datadir.replace('\\', '/') """ """ # spectural_subset = '1' # 波段 0代表不处理，1代表处理，比如HDF有3个波段，处理第2个，则参数应该为：'0 1 0'。具体第几个可以先用MRT打开一个HDF文件看一下 outpath = r'。。。\hdf_mosaic' # 拼接后输出HDF路径 outpath = outpath.replace('\\', '/') allfiles = os.listdir(datadir) # 列出所有HDF文件 # ------按照每天分组------------ allHDFfiles = [] allDays = [] for eachfile in allfiles: if os.path.splitext(eachfile)[1] == '.hdf': allHDFfiles.append(eachfile) day = eachfile.split('.') allDays.append(day[2]) # [2]这个参数根据文件名确定，就是代表日期的那个字段所在位置 allDays = set(allDays) for eachday in allDays: f = open(outpath + '/' + eachday + '.txt', 'a') for i in range(len(allHDFfiles)): if allHDFfiles[i].split('.')[2] == eachday: # [2]这个参数根据文件名确定，就是代表日期的那个字段所在位置 f.write(datadir+'/' + allHDFfiles[i]+'\n') f.close() # ------------------------------ # ---------执行拼接操作---------- for eachday in allDays: try: if os.path.exists(outpath + '/' + eachday + '.hdf'): continue mosaicFiles = '{0} -i '.format(mrtmosaic) + outpath + '/' + eachday + '.txt' + ' -s "{0}" -o {1}.hdf'.format(spectural_subset, eachday) os.system(mosaicFiles) moveFile = "move {0}.hdf {1}".format(eachday, outpath) os.system(moveFile) except: print(eachday + "was wrong") print("all has done")