MOD13 NDVI数据处理教程
适用于用MODIS数据做植被长时序分析,并且用气候数据做残差分析分析人类活动的
前言
从官网下载的MODIS数据
提示:以下是本篇文章正文内容,下面案例可供参考
一、MODIS NDVI 重命名
下载后得MODIS元数据名称特别长,并且是按照天作为时间单位,重命名成几月的方面后续处理
代码如下(示例):
#-*- coding:UTF-8 -*-
import datetime
import os
startTime = datetime.datetime.now()
def ydToymd(yd):
'''Change the date style
'2019060' --> '20190301'
'''
year = int(yd[:4])
day = int(yd[-3:])
ymd = datetime.date(year,1,1)+datetime.timedelta(day-1)
ymd_str = ymd.strftime('%Y%m%d')
return ymd_str
# data path
dataPath =r"F:\MODISNDVI\2014"
dataFormat=".tif"
files = os.listdir(dataPath)
for f in files:
fileFullPath = dataPath+'\\'+f
fnameList = f.split('.')
yearDay = fnameList[1][1:]
ymd = ydToymd(yearDay)
newFileFullPath = dataPath + '\\' + 'MODIS.' + ymd + dataFormat
os.rename(fileFullPath,newFileFullPath)
print "Start time : "+startTime.strftime('%Y-%m-%d-%H:%M:%S')
print "End time : "+datetime.datetime.now().strftime('%y-%m-%d-%H:%M:%S')