TSGFspline: Temporal smoothing and gap filling using splines

In greenbrown: Land Surface Phenology and Trend Analysis

Description

This function fills gaps in a time series by using na.spline and smoothes the time series by using smooth.spline

Usage

TSGFspline(Yt, interpolate = FALSE, 
    ...)

Arguments

Yt	univariate time series of class ts.
interpolate	Should the smoothed and gap filled time series be interpolated to daily values by using na.spline?
...	further arguments (currently not used)

Value

The function returns a gap-filled and smoothed version of the time series.

Author(s)

Matthias Forkel <[e-mail]> [aut, cre]

See Also

TsPP

Examples

# introduce random gaps 
gaps <- ndvi
gaps[runif(100, 1, length(ndvi))] <- NA
plot(gaps)

# do smoothing and gap filling
tsgf <- TSGFspline(gaps)
plot(gaps)
lines(tsgf, col="red")

# compare original data with gap-filled data
plot(ndvi[is.na(gaps)], tsgf[is.na(gaps)], xlab="original", ylab="gap filled")
abline(0,1)
r <- cor(ndvi[is.na(gaps)], tsgf[is.na(gaps)])
legend("topleft", paste("Cor =", round(r, 3)))

greenbrown documentation built on Dec. 18, 2020, 3:02 p.m.