## LOAD THE CUSTOM FUNCTIONS AND PACKAGES
source("myFunctions.R") ; library(reshape) ; library(rgdal)
## LOAD DATA FROM SIMULATIONS
load("Data/LandUseS2.Rda") ; load("Data/LandUseS4.Rda")
VarPS4 <- data.frame(MAILLE= LandUseS4$MAILLE,
                     (LandUseS4[, -1]- LandUseS2[, -1])* 10)
## EXTRACT THE LAND USE ABOUT ANNUAL CROPS
PRS4 <- VarPS4[, substr(names(VarPS4), 1, 2)%in% c("MA", "PR")]
PPR4 <- melt(PRS4, id= "MAILLE")
PPR4$time <- as.numeric(substr(PPR4$variable, 3, 4))
## LINEAR INTERPOLATION OF DECENIAL TRENDS
LinApprox(PPR4)
## MERGE WITH THE GEOGRAPHIC FILES
MPRS4 <- merge(DAT, CRD[, 6: 8], by= "MAILLE")
coordinates(MPRS4) <- c("X1", "X2")
## CREATE THE SPATIO-TEMPORAL ANIMATION
AnimHTML(MPRS4,
         ttle= "Variations of Pastures, Scenario S4 vs S2",
         pal= "Spectral",
         brks= c(-105, -.01, 25, 50, 75, 100, 125, 150, 500),
         t1= "Variations of Pastures, Scenario S4 relative to S2 in ",
         t2= "", dir= "./output/PR4", name= "PRS4.html")