install.packages("cartography")

require(devtools)
devtools::install_github("Groupe-ElementR/cartography")

# Chargement de la librairie
library(cartography)

# Import de données dans la session
data(nuts2006)

ls()

# Calcul du taux de croissance annuel moyen
nuts2.df$cagr <- 100 * (((nuts2.df$pop2008/nuts2.df$pop1999)^(1/9)) - 1)

# Cartographie
choroLayer(spdf = nuts2.spdf, df = nuts2.df, var = "cagr")
title("Taux de croissance en Europe")

# Construire une palette de couleurs
cols <- carto.pal(pal1 = "green.pal", n1 = 2, 
                  pal2 = "red.pal", n2 = 4) 

# Affichage de couches d'habillage
plot(nuts0.spdf, border = NA, col = NA, bg = "#A6CAE0")
plot(world.spdf, col  = "#E3DEBF", border=NA, add=TRUE)

# Cartographie du taux de croissance annuel moyen
choroLayer(spdf = nuts2.spdf, df = nuts2.df, var = "cagr",#
           breaks = c(-2.43,-1.0,0.0,0.5,1.0,2.0,3.1), 
           col = cols,
           border = "grey40",
           lwd = 0.5, 
           legend.pos = "right",
           legend.title.txt = "taux de croissance\nannuel moyen", 
           legend.values.rnd = 2, 
           add = TRUE) 

# Affichage de couches d'habillage
plot(nuts0.spdf,border = "grey20", lwd=0.75, add=TRUE)

# Ajout des titres, légende, sources, etc.
layoutLayer(title = "Taux de croissance en Europe", 
            author = "cartography", 
            sources = "Eurostat, 2008", frame = TRUE, col = NA, 
            scale = NULL,coltitle = "black",
            south = TRUE) 

# Affichage de couches d'habillage
plot(nuts0.spdf, border = NA, col = NA, bg = "#A6CAE0")
plot(world.spdf, col  = "#E3DEBF", border=NA, add=TRUE)
plot(nuts0.spdf, col = "#D1914D",border = "grey80", add=TRUE)

# Cartographie de la population des pays en cercles proportionnels
propSymbolsLayer(spdf = nuts0.spdf, df = nuts0.df,#
                 var = "pop2008", 
                 symbols = "circle", col =  "#920000",
                 legend.pos = "right",
                 legend.title.txt = "Total\npopulation (2008)",
                 legend.style = "c")

# Ajout des titres, légende, sources, etc.
layoutLayer(title = "Countries Population in Europe",
            author = "cartography",
            sources = "Eurostat, 2008",
            scale = NULL,
            south = TRUE)

# Données sur les jumelages
head(twincities)

# Creation d'une couche de liens
twincities.spdf <- getLinkLayer(spdf = nuts2.spdf, df = twincities) 

# Affichage des liens créés
plot(twincities.spdf, lwd = 0.2)

# Affichage de couches d'habillage
plot(nuts0.spdf, border = NA, col = NA, bg = "#A6CAE0")
plot(world.spdf, col  = "#E3DEBF", border=NA, add=TRUE)
plot(nuts2.spdf, col = "#D1914D",border = "grey80", add=TRUE)

# Cartographie des liens
gradLinkLayer(spdf = twincities.spdf, df = twincities,   #
              spdfids = "i", spdfide = "j", dfids = "i", dfide = "j", 
              var = "fij", 
              breaks = c(2,5,15,20,30), 
              lwd = c(0.1,1,4,10), 
              col = "#92000090",
              legend.pos = "right", legend.frame = TRUE,
              legend.title.txt = "Number of Agreements\n(regional level)",
              add = TRUE)

# Ajout des titres, légende, sources, etc.
layoutLayer(title = "International Twinning Agreements Between Cities", 
            author = "cartography", 
            sources = "Sources: Adam Ploszaj & Wikipedia, 2011",
            scale = NULL, south = TRUE, frame = TRUE, col = NA, 
            coltitle = "black")

# Construction les polylignes des frontières inter-étatiques
nuts0.contig.spdf <- getBorders(nuts0.spdf)

plot(nuts0.spdf, bg = "#A6CAE0", col = "#D1914D",border = "grey80")
plot(nuts0.contig.spdf, col = 1:nrow(nuts0.contig.spdf), lwd = 2, add=TRUE)

head(nuts0.contig.spdf@data)

# Calcul du PIB/habitant
nuts0.df$gdpcap <- nuts0.df$gdppps2008/nuts0.df$pop2008*1000000

# Affichage de couches d'habillage
plot(nuts0.spdf, border = NA, col = NA, bg = "#A6CAE0")
plot(world.spdf, col  = "#E3DEBF", border=NA, add=TRUE)

# Cartographie du PIB/habitants
choroLayer(spdf = nuts0.spdf, df = nuts0.df, var = "gdpcap", #
           border = "grey80",
           col = carto.pal(pal1 = "kaki.pal", n1 = 6), 
           method = "quantile",
           nclass = 6, add=TRUE, 
           legend.pos = "right", 
           legend.values.rnd = -2,
           legend.title.txt = "GDP per Capita\n(in euros)")

# Plot discontinuities
discLayer(spdf = nuts0.contig.spdf, df = nuts0.df, var = "gdpcap", 
          type = "rel", 
          method = "equal", 
          nclass = 4, 
          threshold = 0.5, 
          sizemin = 0.5, 
          sizemax = 6, 
          col="red",
          legend.values.rnd = 1,
          legend.title.txt = "Discontinuities in \nGDP per Capita\n(relative)",
          legend.pos = "topright", 
          add=TRUE)

# Ajout des titres, légende, sources, etc.
layoutLayer(title = "Wealth Disparities in Europe",
            coltitle = "black", col = NA,
            sources = "Eurostat, 2011", 
            scale = NULL,
            author = "cartography", 
            frame = FALSE)

# Creation d'une grille régulière
mygrid <- getGridLayer(spdf=nuts2.spdf, cellsize = 200000)

# Affichage de la grille
plot(mygrid$spdf)

# Adaptation des données à la grille
datagrid.df <- getGridData(x = mygrid, df = nuts2.df, var = "pop2008") 
datagrid.df$densitykm <- datagrid.df$pop2008_density*1000*1000

# Affichage de couches d'habillage
plot(nuts0.spdf, border = NA, col = NA, bg = "#A6CAE0")
plot(world.spdf, col  = "#E3DEBF", border=NA, add=TRUE)

# Cartographie de la densité de population 
choroLayer(spdf = mygrid$spdf,                         #
           df = datagrid.df,
           var = "densitykm", 
           border = "grey80", 
           col = carto.pal(pal1 = "wine.pal", n1 = 6), 
           legend.pos = "topright",
           method = "q6", add = TRUE,
           legend.title.txt = "Population Density\n(inhabitant/km²)")

# Ajout des titres, légende, sources, etc.
layoutLayer(title = "Population Density", 
            coltitle = "black", col = NA,
            sources = "Eurostat, 2011", scale = NULL, 
            author = "cartography", frame = FALSE)

# Affichage de couches d'habillage
plot(nuts0.spdf, border = NA, col = NA, bg = "#A6CAE0")  #
plot(world.spdf, col = "#E3DEBF", border = NA, add = TRUE)

nuts3.df$gdppps2008 <- nuts3.df$gdppps2008 * 1e+06

# Cartographie des potentiels de PIB/Hab
smoothLayer(spdf = nuts3.spdf, df = nuts3.df, var = "gdppps2008", var2 = "pop2008", 
    span = 75000, beta = 2, mask = nuts0.spdf, legend.title.txt = "Potential\nGDP per capita", 
    add = TRUE, legend.pos = "topright", legend.values.rnd = -2)

# Ajouter un court text sur la carte
text(x = 6271272, y = 3743765, labels = "Distance function:\n- type = exponential\n- beta = 2\n- span = 75 km", 
    cex = 0.8, adj = 0, font = 3)

# Ajout des titres, légende, sources, etc.
layoutLayer(title = "Wealth Inequality in Europe", sources = "Basemap: UMS RIATE, 2015 - Data: Eurostat, 2008", 
    author = "T. Giraud, 2015")

data("nuts2006")
# extract Denmark
spdf <- nuts0.spdf[nuts0.spdf$id == "DK", ]
# Download the tiles, extent = Denmark
den <- getTiles(spdf = spdf, type = "osmtransport", crop = TRUE, zoom = 7)
class(den)

# Plot the tiles
tilesLayer(den)
# Plot countries
plot(spdf, border = "red", lwd = 2, add = TRUE)
# Map tiles sources
mtext(text = "Map data © OpenStreetMap contributors, under CC BY SA.", side = 1, 
    adj = 0, cex = 0.7, font = 3, line = -2)