phenofit (version 0.2.5-2)

get_pheno: get_pheno

Description

Get yearly vegetation phenological metrics of a curve fitting method

Usage

get_pheno(fits, method, TRS = c(0.2, 0.5), analytical = TRUE,
  smoothed.spline = FALSE, IsPlot = FALSE, showName_fitting = TRUE,
  ...)
get_pheno.fFITs(fFITs, method, TRS = c(0.2, 0.5), analytical = TRUE,
  smoothed.spline = FALSE, IsPlot = FALSE, title_left = "",
  showName_pheno = TRUE)

Arguments

fits

A list of fFITs() object, for a single curve fitting method.

method

Which fine curve fitting method to be extracted?

TRS

Threshold for PhenoTrs.

analytical

If true, numDeriv package grad and hess will be used; if false, D1 and D2 will be used.

smoothed.spline

Whether apply smooth.spline first?

IsPlot

Boolean. Whether to plot figure?

showName_fitting

Whether to show the name of fine curve fitting method in top title?

...

ignored.

fFITs

fFITs object returned by curvefit().

title_left

String of growing season flag.

showName_pheno

Whether to show names of phenological methods in top title? Generally, only show top title in the first row.

Value

List of every year phenology metrics

Examples

# NOT RUN {
library(phenofit)
# simulate vegetation time-series
fFUN = doubleLog.Beck
par  = c(
    mn  = 0.1,
    mx  = 0.7,
    sos = 50,
    rsp = 0.1,
    eos = 250,
    rau = 0.1)
t    <- seq(1, 365, 8)
tout <- seq(1, 365, 1)
y <- fFUN(par, t)
methods <- c("AG", "Beck", "Elmore", "Gu", "Zhang") # "Klos" too slow
fFITs <- curvefit(y, t, tout, methods)
# multiple years
fits  <- list(`2001` = fFITs, `2002` = fFITs)

l_param   <- get_param(fits)
d_GOF     <- get_GOF(fits)
d_fitting <- get_fitting(fits)
l_pheno   <- get_pheno(fits, "AG", IsPlot=TRUE)
# }