lucC_pred_recur: Predicate Recur

In ammaciel/lucCalculus: Land Use Change Calculus

Description

Provide a predicate RECUR which evaluates as true when a location holds a class c_i, e.g. 'Forest', during two non-continuous distinct intervals t_i and t_j. Return a matrix with values within defined interval

Usage

lucC_pred_recur (raster_obj = NULL, raster_class = NULL,
time_interval1 = c("2001-01-01", "2001-01-01"),
time_interval2 = c("2002-01-01", "2005-01-01"),
label = NULL, timeline = NULL, remove_column = TRUE)

Arguments

raster_obj	Raster. A raster brick with classified images
raster_class	Character. Name of the class of interest, such as 'Forest', to research
time_interval1	Interval. A first time interval to verify if class is over or not
time_interval2	Interval. A second and non-overlapped time interval to verify if class is over or not
label	Character Vector. All labels of each value of pixel from classified raster
timeline	Character. A list of all dates of classified raster, timeline
remove_column	Boolean. Remove matrix values relating to the first time interval, only values of second interval are returned. Default is TRUE

Value

Matrix with all states which holds during a time interval

Author(s)

Adeline M. Maciel

Examples

## Not run: 
library(lucCalculus)

file <- c(system.file("extdata/raster/rasterSample.tif", package = "lucCalculus"))
rb_class <- raster::brick(file)
my_label <- c("Degradation", "Fallow_Cotton", "Forest", "Pasture", "Soy_Corn", "Soy_Cotton",
              "Soy_Fallow", "Soy_Millet", "Soy_Sunflower", "Sugarcane", "Urban_Area", "Water")
my_timeline <- c("2001-09-01", "2002-09-01", "2003-09-01", "2004-09-01", "2005-09-01",
                 "2006-09-01", "2007-09-01", "2008-09-01", "2009-09-01", "2010-09-01",
                 "2011-09-01", "2012-09-01", "2013-09-01", "2014-09-01", "2015-09-01",
                 "2016-09-01")

b <- lucC_pred_recur(raster_obj = rb_class, raster_class = "Forest",
                     time_interval1 = c("2001-09-01","2001-09-01"),
                     time_interval2 = c("2002-09-01","2016-09-01"),
                     label = my_label, timeline = my_timeline)

lucC_plot_raster_result(raster_obj = rb_class, data_mtx = b,
                        timeline = my_timeline, label = my_label,
                        custom_palette = FALSE)


## End(Not run)

ammaciel/lucCalculus documentation built on June 13, 2020, 4:57 a.m.