## 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 ANNUAL CROPS CAS2 <- LandUseS2[, substr(names(LandUseS2), 1, 2)%in% c("MA", "CA")] PCA2 <- melt(CAS2, id= "MAILLE") PCA2$time <- as.numeric(substr(PCA2$variable, 3, 4)) ## LINEAR INTERPOLATION OF DECENIAL TRENDS LinApprox(PCA2) ## MERGE WITH THE GEOGRAPHIC FILES MCAS2 <- merge(DAT, CRD[, 6: 8], by= "MAILLE") coordinates(MCAS2) <- c("X1", "X2") ## CREATE THE SPATIO-TEMPORAL ANIMATION AnimHTML(MCAS2, ttle= "Annual crops, Scenario S2", pal= "Spectral", brks= c(-100, -30, -15, -5, 0, 5, 15, 30, 100), t1= "Scenario S2, Annual crops: (Area ", t2= "- Area 2003)* 100", dir= "./output/CA2", name= "CAS2.html")