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,
  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

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)

phenofit documentation built on April 2, 2020, 5:07 p.m.