
Phenology extraction in Threshold method (TRS)
PhenoTrs(
x,
t = NULL,
approach = c("White", "Trs"),
trs = 0.5,
asymmetric = TRUE,
IsPlot = TRUE,
...
)# S3 method for fFIT
PhenoTrs(x, t = NULL, ...)
# S3 method for default
PhenoTrs(
x,
t = NULL,
approach = c("White", "Trs"),
trs = 0.5,
asymmetric = TRUE,
IsPlot = TRUE,
...
)
numeric vector, or fFIT
object returned by curvefit()
.
doy
vector, corresponding doy of vegetation index.
to be used to calculate phenology metrics. 'White' (White et al. 1997) or 'Trs' for simple threshold.
threshold to be used for approach "Trs", in (0, 1).
If true, background value in spring season and autumn season is regarded as different.
whether to plot?
other parameters to PhenoPlot
PhenoDeriv()
, PhenoGu()
, PhenoKl()
# `doubleLog.Beck` simulate vegetation time-series
t <- seq(1, 365, 8)
tout <- seq(1, 365, 1)
par = c( mn = 0.1 , mx = 0.7 , sos = 50 , rsp = 0.1 , eos = 250, rau = 0.1)
y <- doubleLog.Beck(par, t)
methods <- c("AG", "Beck", "Elmore", "Gu", "Zhang") # "Klos" too slow
fit <- curvefit(y, t, tout, methods)
x <- fit$model$AG # one model
par(mfrow = c(2, 2))
PhenoTrs(x)
PhenoDeriv(x)
PhenoGu(x)
PhenoKl(x)
Run the code above in your browser using DataLab