## 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 URS4 <- VarPS4[, substr(names(VarPS4), 1, 2)%in% c("MA", "UR")] PUR4 <- melt(URS4, id= "MAILLE") PUR4$time <- as.numeric(substr(PUR4$variable, 3, 4)) ## LINEAR INTERPOLATION OF DECENIAL TRENDS LinApprox(PUR4) ## MERGE WITH THE GEOGRAPHIC FILES MURS4 <- merge(DAT, CRD[, 6: 8], by= "MAILLE") coordinates(MURS4) <- c("X1", "X2") ## CREATE THE SPATIO-TEMPORAL ANIMATION AnimHTML(MURS4, ttle= "Variations of Urban Area, Scenario S4 vs S2", pal= "Spectral", brks= c(-500, -250, -100, -50, -25, -15, -5, 10, 105), t1= "Variations of Urban Area, Scenario S4 relative to S2 in ", t2= "", dir= "./output/UR4", name= "URS4.html")