## LOAD THE CUSTOM FUNCTIONS source("myFunctions.R") ## LOAD DATA FROM SIMULATIONS AND PACKAGES load("Data/LandUseS1.Rda") ; library(reshape) ; library(rgdal) ## EXTRACT THE LAND USE ABOUT PERENNIAL CROPS PES1 <- LandUseS1[, substr(names(LandUseS1), 1, 2)%in% c("MA", "PE")] PPE1 <- melt(PES1, id= "MAILLE") PPE1$time <- as.numeric(substr(PPE1$variable, 3, 4)) ## LINEAR INTERPOLATION OF DECENIAL TRENDS LinApprox(PPE1) ## MERGE WITH THE GEOGRAPHIC FILES MPES1 <- merge(DAT, CRD[, 6: 8], by= "MAILLE") coordinates(MPES1) <- c("X1", "X2") ## CREATE THE SPATIO-TEMPORAL ANIMATION AnimHTML(MPES1, ttle= "Perennial Crops, Scenario S1", pal= "Spectral", brks= c(-100, -30, -15, -5, 0, 5, 15, 30, 100), t1= "Scenario S1, Perenial Crops: (Area ", t2= "- Area 2003)* 100", dir= "./output/PE1/", name= "PES1.html")