Commit 5580b801 authored by Gaetano Raffaele's avatar Gaetano Raffaele
Browse files

Temporal smoothing (Savitzky Golay) script (pre).

parent deab96a9
......@@ -3,7 +3,16 @@ 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,clip_valid=23,clip_both_sides=True):
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,
clip_valid=23,
clip_both_sides=True):
lst = sorted(glob.glob(fld + '/' + template + '*.' + ext))
......@@ -35,4 +44,4 @@ def temporalSmoothing(fld,template='MOD13Q1.006__250m_16_days_NDVI_',date_pos=34
cdf.write(dts[n])
cmd += ['{' + ';'.join(expr) + '}']
print cmd
#subprocess.call(cmd)
\ No newline at end of file
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