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

In phenomap: Projecting Satellite-Derived Phenology in Space

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)

phenomap documentation built on Oct. 13, 2023, 1:12 a.m.