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

In phenomap: Projecting Satellite-Derived Phenology in Space

Description

Convert a series of phenology raster files to a single long-term trend raster.

Usage

mapTrend(File_List, Year_List, parallel = FALSE, n.cores = NULL,
  verbose = FALSE)

Arguments

File_List	List of phenology raster files (i.e. those produced in 'mapPheno')
Year_List	Vector of Integer Year (YYYY) with length > 5
parallel	TRUE or FALSE (Default = FALSE) if TRUE, use parallel backend through plyr::aaply
n.cores	Integer number of cores to be used for parallel processing (only use if parallel = TRUE)
verbose	TRUE or FALSE (Default = FALSE)

Value

Raster object with extent=extent(raster(File_List)[1]) and CRS = crs(raster(File_List)[1]). Layer 1 is the slope estimate of the linear model relating green-up timing (Day of Year) to time (Year). Layer 2 is the p-value of the slope estimate. Layer 3 is the standard error of the slope estimate. Layer 4 is the r-squared value for the linear model.

Examples

## Not run: 

fpath <- system.file("extdata", package="phenomap")
File_List.Trend <- paste(fpath, list.files(path = fpath, pattern=c("Sample_Greenup_")), sep="/")

Year_List <- 2011:2016 # Tell it what years you're using
n.cores <- 4 # Set up parallel computing

phenotrend <- mapTrend(File_List = File_List.Trend,
                             Year_List = Year_List,
                             parallel = TRUE,
                             n.cores = n.cores,
                             verbose=TRUE)


## End(Not run)

phenomap documentation built on May 2, 2019, 11:01 a.m.