Commit 2a653407 authored by Gaetano Raffaele's avatar Gaetano Raffaele
Browse files

Temporal smoothing (Savitzky Golay) script (pre).

parent e63d1747
import glob
import datetime
import os
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):
......@@ -10,7 +11,7 @@ def temporalSmoothing(fld,template='MOD13Q1.006__250m_16_days_NDVI_',date_pos=34
with open(fld + '/dates.txt','wb') as df:
for f in lst:
dt = f[date_pos:date_pos+date_len]
dt = os.path.basename(f)[date_pos:date_pos+date_len]
if date_jul:
y,jd = int(dt[0:4]),int(dt[4:7])
dt = (datetime.datetime(y,1,1) + datetime.timedelta(jd)).strftime('%Y%m%d')
......
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