Dear Gabor, Paul and Ryan,
Many thanks for all your guidance concerning the application of the
filter/smoother. I need to temporally smooth 32 years 5 km resolution
bi-monthly global NDVI dataset. On average each pixel have got up to 768
cases/images with data missing at random. I consider the dat
On Fri, Sep 26, 2014 at 3:32 AM, Erick Okuto wrote:
> Dear Paul and Henrik,
> I have a time series with some missing data points that i need smoothed
> using Savitzky-Golay filter. Related question was asked here
> http://thr3ads.net/r-help/2012/11/2121748-Savitzky-Golay-filtering-with-missing-da
You can perform the equivalent filtering by hand using something like:
k = (window_length - 1) / 2
for i in k to len(x) - k
xh,yh = x[points-k:points+k], y[points-k:points+k]
idx = isfinite(xh) and isfinite(yh)
xh,yh = xh[idx],yh[idx]
p = polyfit(xh,yh, degr
Dear Paul and Henrik,
I have a time series with some missing data points that i need smoothed
using Savitzky-Golay filter. Related question was asked here
http://thr3ads.net/r-help/2012/11/2121748-Savitzky-Golay-filtering-with-missing-data
but no straight forward answer was posted. However, Henrik
4 matches
Mail list logo
