Commit 94cb5e01 authored by Gaetano Raffaele's avatar Gaetano Raffaele
Browse files

Temporal smoothing (Savitzky Golay) script (pre).

parents
import glob
import datetime
import subprocess
def temporalSmoothing(fld,template='MOD13Q1.006__250m_16_days_NDVI_',date_pos=34,date_len=7,date_jul=True,ext='tif',sgdeg=4,sgrad=6):
lst = sorted(glob.glob(fld + '/' + template + '*.' + ext))
cmd = ['otbcli_TemporalSmoothing','-il']
with open(fld + '/dates.txt','wb') as df:
for f in lst:
dt = f[date_pos:date_pos+date_len]
if date_jul:
y,jd = int(f[0:4]),int(f[4:7])
dt = (datetime.datetime(y,1,1) + datetime.timedelta(jd)).strftime('%Y%m%d')
df.write(dt+'\n')
cmd += [fld + '/' + f]
cmd += ['-dates',fld+'/dates','-out',fld+'/SmoothedSeries.tif','uint16','-interp','sg','-interp.sg.deg',str(sgdeg),'-interp.sg.rad',str(sgrad)]
subprocess.call(cmd)
\ No newline at end of file
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment