## RasterBrick of structural Essential Biodiversity Variables
## covering the extent of a location in the northern Amazon basin
## (Colombia) is imported:
path. <- system.file('amazon.grd',package = 'ecochange')
amazon <- brick(path.)
## Changes in layers of tree-canopy cover (TC) are computed by
## processing the 'amazon' brick:
def <- echanges(amazon, eco = 'TC',
change = 'lossyear',
eco_range = c(1,80),
get_unaffected = TRUE,
binary_output = FALSE,
mc.cores = 2)
## Function 'EBVstats' is used to compute ecosystem statistics
st_amazon <- EBVstats(def)
## A plot of the 'st_amazon' object
plot.EBVstats(st_amazon,
cex = 1.5,
xlab = 'Year',
ylab = 'Canopy cover (%)',
main = 'Ecosystem changes',
sub = 'Municipality: Cartagena del Chaira',
fill = 'Layer')
Run the code above in your browser using DataLab