lucC_relation_in: Allen Relation In

In ammaciel/lucCalculus: Land Use Change Calculus

Description

Provide an Allen's interval relation IN which aggregates the relations DURING, STARTS and FINISHES. See more at (ALLEN, J. F. "Maintaining knowledge about temporal intervals". Communications of the ACM, v(26), 11, 1983, 832-843. DOI: http://dx.doi.org/10.1145/182.358434)

Usage

lucC_relation_in(first_raster = NULL, second_raster = NULL)

Arguments

first_raster	Matrix. An interval between two dates.
second_raster	Matrix. An interval between two dates.

Value

Data set with merge of two data sets

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")

lucC_plot_raster(raster_obj = rb_class, timeline = my_timeline, label = my_label,
                 custom_palette = FALSE, plot_ncol = 4)

a <- lucC_pred_holds(raster_obj = rb_class, raster_class = "Pasture",
                     time_interval = c("2005-09-01","2009-09-01"),
                     relation_interval = "equals", label = my_label,
                     timeline = my_timeline)

b <- lucC_pred_holds(raster_obj = rb_class, raster_class = "Degradation",
                     time_interval = c("2005-09-01","2011-09-01"),
                     relation_interval = "equals", label = my_label,
                     timeline = my_timeline)

# in
c <- lucC_relation_in(first_raster = a, second_raster = b)

lucC_plot_raster_result(raster_obj = rb_class, data_mtx = c, timeline = my_timeline,
                        label = my_label, custom_palette = FALSE, plot_ncol = 4)


## End(Not run)

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