mapPheno {phenomap}

Convert a series of raster files to a single phenology raster.

Description

Convert a series of raster files to a single phenology raster.

Usage

mapPheno(
  File_List = NA,
  PhenoFactor = NA,
  phase = NA,
  threshold = NA,
  year = NA,
  NDVI = NA,
  VIQ = NA,
  DOY = NA,
  PR = NA,
  SnowExtent = NA,
  verbose = FALSE
)

Arguments

File_List

List of raster files

PhenoFactor

Character string; type of dataset to analyze (e.g., "VI", "Snow")

phase

Character string; name of phenophase to be measured (e.g., "greenup", "snowmelt", "senescence" or other arguments passed to phenex::phenophase())

threshold

Float threshold GWI value to be projected. Use only for VI option.

year

Integer Year (YYYY)

NDVI

Integer Band number of NDVI band in raster files

VIQ

Integer Band number of VI Quality layer in raster files

DOY

Integer Band number of Composite Day of Year layer in raster files

PR

Integer Band Number of PR layer in raster files

SnowExtent

Integer Band number of Maximum_Snow_Extent in raster files

verbose

TRUE or FALSE (Default = FALSE)

Value

Raster object with extent=extent(terra::rast(File_List)[1]) and CRS = crs(terra::rast(File_List)[1]). Digital numbers are expressed as Day of Year.

Examples

## Not run: 
fpath <- system.file("extdata", package="phenomap")
File_List <- paste(fpath, list.files(path = fpath, pattern=c("TinyCrop_")), sep="/")
File_List

PhenoFactor = "VI"
phase = "greenup"
threshold = 0.5
year = 2016
NDVI = 1
VIQ = 3
DOY = 4
PR = 5
verbose = TRUE

Sample.Greenup <- mapPheno(File_List = File_List, PhenoFactor = PhenoFactor,
                           phase = phase, threshold = threshold, year = year,
                           NDVI = NDVI, VIQ = VIQ, DOY = DOY, PR = PR,
                           SnowExtent=SnowExtent,
                           verbose = verbose)

## End(Not run)