Forums » Support »
Time series filtering
Added by Dirk Olonscheck over 4 years ago
Hello,
I am unsure how to use the lowpass filter correctly.
I have a yearly timeseries and like to transform it do decadal, i.e. omitting all variability on timscales smaller than 10 years.
I tried with
cdo lowpass,0.1 data.nc data_decadal.nc
Is this correct? The description in the CDO Manual on how to determine fmin and fmax confused me.
Thanks,
Dirk
hi Dirk!
could you upload the data?
cheers
ralf
Hi Ralf,
here is the data, the first file without the filtering, the second one with the lowpass filter applied. Is this correct?
Afterwards I do (with both files)
cdo timstd tsc_control_yr_detrend.nc tsc_control_yr_detrend_timstd.nc
to estimate the variability, comparing interannual and decadal.
I am simply not sure if "cdo lowpass,0.1" really removes all variability on timescales smaller than 10 years.
Thanks for your support!
Dirk
Hi Dirk!
here is a plot of the original and the filtered data:
I didn't use the second upload, but called CDO instead. From my point of view it looks good. But I need to check the frequencies itself in more depth.
Forgot to mention, that I used
cdo -sellonlatbox,-6,-2,-2,0
to select a single location from your input file. This location I used to timeseries analysis.
Thanks Ralf!
This looks correct indeed. I can reconstruct your plot. I would be happy to hear whether the frequencies are as intended, but the result is promising. So "cdo lowpass,0.1" seems to be correct for omitting frequencies higher than 10 years for annual data.
Best,
Dirk
IMO 0.1 is the right value for 10-year frequency according to the documentation. '1' represents 'yearly'. Hence a 10-times smaller frequency must represent decadel.
The examples in the documentation use a value of 0.5 for 2-years variability. So, your value of 0.1 makes perfectly sense to me.
cheers
ralf
Thanks a lot! That's what I wanted to ensure. Makes totally sense for me now 
Have a nice weekend,
Dirk