## LOAD THE CUSTOM FUNCTIONS AND PACKAGES
source("myFunctions.R") ; library(reshape) ; library(rgdal)
## LOAD DATA FROM SIMULATIONS
load("Data/LandUseS1.Rda") ; load("Data/LandUseS3.Rda")
VarPS3 <- data.frame(MAILLE= LandUseS3$MAILLE,
                     (LandUseS3[, -1]- LandUseS1[, -1])* 10)
## EXTRACT THE LAND USE ABOUT ANNUAL CROPS
FOS3 <- VarPS3[, substr(names(VarPS3), 1, 2)%in% c("MA", "FO")]
PFO3 <- melt(FOS3, id= "MAILLE")
PFO3$time <- as.numeric(substr(PFO3$variable, 3, 4))
## LINEAR INTERPOLATION OF DECENIAL TRENDS
LinApprox(PFO3)
## MERGE WITH THE GEOGRAPHIC FILES
MFOS3 <- merge(DAT, CRD[, 6: 8], by= "MAILLE")
coordinates(MFOS3) <- c("X1", "X2")
## CREATE THE SPATIO-TEMPORAL ANIMATION
AnimHTML(MFOS3,
         ttle= "Variations of Forests, Scenario S3 vs S1",
         pal= "Spectral",
         brks= c(-500, -250, -100, -50, -25, -15, -5, 0.01, 105),
         t1= "Variations of Forests, Scenario S3 relative to S1 in ",
         t2= "", dir= "./output/FO3", name= "FOS3.html")