资源简介

使用python脚本,调用MRT批量拼接MODIS三级和四级产品。MODIS数据可以缺失某一天的数据

资源截图

代码片段和文件信息

# -*- coding:utf-8 -*-

import os

#设置MRT的data目录环境变量(请设置为自己的MRT安装路径下的data目录)
os.environ[‘MRT_DATA_DIR‘] = “D:/softs/MRT/MRT_Win/data“
#设置MRT的bin路径(请设置为自己的MRT安装路径下的bin目录)
mrtpath = r“D:\softs\MRT\MRT_Win\bin“ #请设置为自己的MRT安装路径下的bin目录
mrtmosaic = os.path.join(mrtpath “mrtmosaic.exe“) #指定MRT的mrtmosaic.exe模块的全路径
mrtmosaic = mrtmosaic.replace(“\\“ “/“) #MRT的mrtmosaic.exe模块的全路径以“/“连接
#指定数据目录
datadir = ur‘D:\works\MODIS_DATA\MOD11B1-2013(04-10)\01DATA‘
os.chdir(datadir)
#指定获取的波段
SPECTRAL_SUBSET = “1 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0“ #MOD11B1数据中的LST_Day_6km、LST_Night_6km数据集
#指定结果输出路径
outpath = ur“D:\works\MODIS_DATA\MOD11B1-2013(04-10)\02MOSAIC“
outpath = outpath.replace(“\\“ “/“)

#获取当前文件夹下所有hdf文件
allFiles = os.listdir(datadir)

#创建allDays空列表,存储所有天数
allHDFfiles = []
allDays = []

#从文件名中获取天数信息
for eachFile in allFiles:
    if os.path.splitext(eachFile)[1] == “.hdf“: #分离文件名与扩展名
        allHDFfiles.append(eachFile)
        day = eachFile.split(‘.‘)
        allDays.append(day[1])
print(“所有hdf文件如下:“)
print(“    “ + “\n    

评论

共有 条评论