資源簡介
使用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????
評論
共有 條評論