MODIS_Mosaic.py

# -*- 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