get_pheno: get_pheno

Description Usage Arguments Value Note Examples

View source: R/get_pheno.R

Description

Get yearly vegetation phenological metrics of a curve fitting method

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
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.