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