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

Description Usage Arguments Value Examples

View source: R/mapTrend.R

Description

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

Usage

1
2
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

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

Related to mapTrend in phenomap...