create_events: Create phenological events (onset/diapause/mortality)

In barrks: Calculate Bark Beetle Phenology Using Different Models

Create phenological events (onset/diapause/mortality)

Description

Generate onset, diapause or mortality manually to be able to run phenology() with observed or arbitrary inputs.

Usage

create_onset(
  template,
  doys = NULL,
  stations = NULL,
  dates = NULL,
  .quiet = FALSE
)

create_diapause(
  template,
  doys = NULL,
  stations = NULL,
  dates = NULL,
  .quiet = FALSE
)

create_mortality(
  template,
  doys = NULL,
  stations = NULL,
  dates = NULL,
  .quiet = FALSE
)

Arguments

template	SpatRaster or data frame that determines the spatial and temporal extent of the result. If a single-layer SpatRaster was passed, the temporal extent should be defined by using dates.
doys	Numeric vector, (multi-layer) SpatRaster or data frame that specifies the days of year when the event is triggered. Vectors will define the events globally whereas SpatRasters allow spatially explict definitions. For the creation of events based on stations, data frames are used. In that case, the field station specifies the station name and doy indicates the respective day of year.
stations	If template is a SpatRaster and doys is a data frame, stations should be passed to define which cells are affected. See stations_create() for details.
dates	Dates to define the temporal extent of the output if template is a single-layer SpatRaster.
.quiet	If TRUE, messages are suppressed.

Value

A logical multi-layer SpatRaster. Each layer represents one date.

Functions

• create_onset(): Create a onset.

• create_diapause(): Create a diapause.

• create_mortality(): Create mortality events.

Examples

# load sample data
d <- barrks_data('stations')

# create onset, diapause, mortality
on <- create_onset(d, lubridate::yday('2015-04-15'))
dia <- create_diapause(d, lubridate::yday('2015-08-15'))
mort <- create_mortality(d, lubridate::yday('2015-11-15'))

# claculate phenologe
p <- phenology('phenips-clim', d, .quiet = TRUE,
               .onset = on, .diapause = dia, .mortality = mort)

# plot development
plot_development_diagram(p, .lwd = 4)

barrks documentation built on April 3, 2025, 9:47 p.m.