get_pheno: get_pheno

In phenofit: Extract Remote Sensing Vegetation Phenology

Description

Get yearly vegetation phenological metrics of a curve fitting method

Usage

get_pheno(
  fits,
  method,
  TRS = c(0.2, 0.5, 0.6),
  analytical = TRUE,
  smoothed.spline = FALSE,
  IsPlot = FALSE,
  show_title = 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?
show_title	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 phenological methods names in the top panel?

Value

List of every year phenology metrics

Note

Please note that only a single fine curve fitting method allowed here!

Examples

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)

Example output

phenofit documentation built on Oct. 15, 2021, 5:09 p.m.