lucC_pred_convert: Predicate Convert

In ammaciel/lucCalculus: Land Use Change Calculus

Description

Provide a predicate CONVERT which evaluates as true when a location holds the class c_i, e.g. 'Forest', during the interval t_i, class c_j, e.g. 'Soybean', during the interval t_j and t_j is sequential of t_i. Return a matrix with values within defined interval

Usage

lucC_pred_convert (raster_obj = NULL, raster_class1 = NULL,
time_interval1 = c("2001-01-01", "2001-01-01"), relation_interval1 = "equals",
raster_class2 = NULL, time_interval2 = c("2002-01-01", "2005-01-01"),
relation_interval2 = "equals", label = NULL, timeline = NULL,
remove_column = TRUE)

Arguments

raster_obj	Raster. A raster brick with classified images
raster_class1	Character. Name of the first class of interest, such as 'Forest', to research
time_interval1	Interval. A first interval to verify if class is over or not
relation_interval1	Character. If a location HOLDS during all time interval 'equals' or can be appear in any times 'contains'. Default is 'equals'
raster_class2	Character. Name of the second class of interest, such as 'Pasture', to research
time_interval2	Interval. A second interval to verify if class is over or not
relation_interval2	Character. If a location HOLDS during all time interval 'equals' or can be appear in any times 'contains'. Default is 'equals'
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")

d <- lucC_pred_convert(raster_obj = rb_class, raster_class1 = "Forest",
                       time_interval1 = c("2012-09-01","2012-09-01"),
                       relation_interval1 = "equals",
                       raster_class2 = "Degradation",
                       time_interval2 = c("2013-09-01","2013-09-01"),
                       relation_interval2 = "equals",
                       label = my_label, timeline = my_timeline)

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


## End(Not run)

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