inst/doc/oceanis.R

In oceanis: Cartography for Statistical Analysis

## ----echo=FALSE---------------------------------------------------------------
knitr::opts_chunk$set(collapse = TRUE)

## -----------------------------------------------------------------------------
donnees_monoloc <- rio::import(file = system.file("data/donnees_monoloc.rda",
                                                  package = "oceanis"))

## ---- fig.width = 5-----------------------------------------------------------
library(oceanis)

# chargement des donnees
donnees_monoloc <- rio::import(file = system.file("data/donnees_monoloc.rda",
                                                  package = "oceanis"))
# visualisation de la distribution de la variable VAR_AN_MOY
distrib_variable(data = donnees_monoloc,
                 varRatio = "VAR_AN_MOY",
                 nbClasses = 4)

## ----fondCarte, fig.height = 3, fig.width = 3---------------------------------
# chemin du fond de carte .shp
path_to_shp <- system.file("extdata",
                           "dep_francemetro_2018.shp",
                           package = "oceanis")
# import de l'objet sf
depm <- sf::st_read(dsn = path_to_shp, quiet = TRUE)

# visualisation de la geometrie
# modification des marges
par(mai = c(0,0,0,0))
# contour des departements de France metropolitaine
plot(sf::st_geometry(depm))

## ----zonageAFacon, fig.height = 5, fig.width = 5------------------------------
library(oceanis)

# chargement des donnees :
# donnees a facon avec variable de regroupement (ZE2010)
donnees_a_facon <- rio::import(file = system.file("data/donnees_a_facon.rda",
                                                  package = "oceanis"))
# fond communal des departements 13, 30, 83 et 84
com_dep_13_30_83_84 <- sf::st_as_sf(rio::import(file = system.file("data/com_dep_13_30_83_84.rda",
                                                      package = "oceanis")))

# chemin du fond de carte .shp
path_to_shp <- system.file("extdata",
                           "dep_francemetro_2018.shp",
                           package = "oceanis")
# import en objet sf
depm <- sf::st_read(dsn = path_to_shp, quiet = TRUE, stringsAsFactors = FALSE)

# creation du zonage des zones d'emploi des Bouches-du-Rhone (partie entiere des ze)
ze13etplus <- zonage_a_facon(fondMaille = com_dep_13_30_83_84,
                             groupe = donnees_a_facon,
                             idMaille = "DEPCOM",
                             idGroupe = "ZE2010",
                             libGroupe = "LIB_ZE2010",
                             fondContour = NULL)
# creation du zonage des zones d'emploi des Bouches-du-Rhone (partie tronquee des ze au contour du departement)
ze13 <- zonage_a_facon(fondMaille = com_dep_13_30_83_84,
                       groupe = donnees_a_facon,
                       idMaille = "DEPCOM",
                       idGroupe = "ZE2010",
                       libGroupe = "LIB_ZE2010",
                       fondContour = depm[depm$CODE=="13",])

# visualisation de la geometrie
# modification des marges
par(mai = c(0,0,0,0))
# fond des ze des Bouches-du-Rhone en 1er pour fixer le niveau de zoom, en bleu
plot(sf::st_geometry(ze13etplus),
     col = "powderblue",
     border = "transparent")
# fond de la partie tronquee des ze des Bouches-du-Rhone, en rouge
plot(sf::st_geometry(ze13),
     col = "lightsalmon",
     border = "transparent",
     add = TRUE)
# contour des communes
plot(sf::st_geometry(com_dep_13_30_83_84),
     col = "transparent",
     border = "lavender",
     add = TRUE)
# contour de la partie tronquee des ze des Bouches-du-Rhone, en rouge
plot(sf::st_geometry(ze13),
     col = "transparent",
     border = "indianred",
     lwd = 3,
     add = TRUE)
# contour des ze des Bouches-du-Rhone, en bleu
plot(sf::st_geometry(ze13etplus),
     col = "transparent",
     border = "steelblue",
     lwd = 3,
     add = TRUE)

## ----shinyRondsAnalyseClasses, eval = FALSE-----------------------------------
#  library(oceanis)
#  library(shiny)
#  
#  # chargement des donnees
#  donnees_monoloc <- rio::import(file = system.file("data/donnees_monoloc.rda",
#                                                    package = "oceanis"))
#  
#  # import du fond des departements
#  depm <- sf::st_read(dsn = system.file("extdata",
#                                        "dep_francemetro_2018.shp",
#                                        package = "oceanis"),
#                      quiet = TRUE,
#                      stringsAsFactors = FALSE)
#  # import du fond des regions
#  regm <- sf::st_read(dsn = system.file("extdata",
#                                        "reg_francemetro_2018.shp",
#                                        package = "oceanis"),
#                      quiet = TRUE,
#                      stringsAsFactors = FALSE)
#  # import du fond de France metropolitaine
#  fram <- sf::st_read(dsn = system.file("extdata",
#                                        "francemetro_2018.shp",
#                                        package = "oceanis"),
#                      quiet = TRUE,
#                      stringsAsFactors = FALSE)
#  
#  shiny_classes_ronds(data = donnees_monoloc,
#                      fondMaille = depm[depm$REG %in% c("93","94"),],
#                      fondMailleElargi = depm,
#                      fondContour = fram,
#                      fondSuppl = regm,
#                      idData = "COD_DEP",
#                      varVolume = "POP_2015",
#                      varRatio = "VAR_AN_MOY")

## ----shinyJoignantes, eval = FALSE--------------------------------------------
#  library(oceanis)
#  library(shiny)
#  
#  # chargement des donnees
#  donnees_biloc <- rio::import(file = system.file("data/donnees_biloc.rda",
#                                                  package = "oceanis"))
#  
#  # import du fond des regions
#  regm <- sf::st_read(dsn = system.file("extdata",
#                                        "reg_francemetro_2018.shp",
#                                        package = "oceanis"),
#                      quiet = TRUE,
#                      stringsAsFactors = FALSE)
#  # import du fond des departements
#  depm <- sf::st_read(dsn = system.file("extdata",
#                                        "dep_francemetro_2018.shp",
#                                        package = "oceanis"),
#                      quiet = TRUE,
#                      stringsAsFactors = FALSE)
#  # import du fond de France metropolitaine
#  fram <- sf::st_read(dsn = system.file("extdata",
#                                        "francemetro_2018.shp",
#                                        package = "oceanis"),
#                      quiet = TRUE,
#                      stringsAsFactors = FALSE)
#  
#  shiny_joignantes(data = donnees_biloc,
#                   fondMaille = regm,
#                   typeMaille = "REG",
#                   fondContour = fram,
#                   fondSuppl = depm,
#                   idDataDepart = "REG_DEPART",
#                   idDataArrivee = "REG_ARRIVEE",
#                   varFlux = "MIGR",
#                   decalageAllerRetour = 10,
#                   decalageCentroid = 20)

## ----leafletAnalyseClassesRonds, fig.height = 6, fig.width = 9----------------
library(oceanis)
library(leaflet)
library(sf)
library(classInt)
library(leaflet.extras)
library(stringr)

# chargement des donnees
donnees_monoloc <- rio::import(file = system.file("data/donnees_monoloc.rda",
                                                  package = "oceanis"))

# import du fond des departements
depm <- sf::st_read(dsn = system.file("extdata",
                                      "dep_francemetro_2018.shp",
                                      package = "oceanis"),
                    quiet = TRUE,
                    stringsAsFactors = FALSE)
# import du fond des regions
regm <- sf::st_read(dsn = system.file("extdata",
                                      "reg_francemetro_2018.shp",
                                      package = "oceanis"),
                    quiet = TRUE,
                    stringsAsFactors = FALSE)

# affichage de la carte
map <- leaflet_ronds_classes(data = donnees_monoloc,
                             fondMaille = depm[depm$REG=="93",],
                             fondMailleElargi = depm,
                             fondSuppl = regm,
                             idData = "COD_DEP",
                             varVolume = "POP_2015",
                             varRatio = "VAR_AN_MOY")

# affichage du rayon du rond le plus grand en metres
rayon_ronds(map)

# affichage de la carte avec des rayons de ronds plus grands
map <- leaflet_ronds_classes(data = donnees_monoloc,
                             fondMaille = depm[depm$REG == "93",],
                             fondMailleElargi = depm,
                             fondSuppl = regm,
                             idData = "COD_DEP",
                             varVolume = "POP_2015",
                             varRatio = "VAR_AN_MOY",
                             rayonRond = 29000,
                             rapportRond = NULL)

# affichage du rapport du rond le plus grand, a recuperer pour permettre la comparaison de plusieurs cartes
rapport_ronds(map)

# ajout de la legende des ronds (position par defaut)
map <- add_legende_ronds(map = map,
                         titre = "Population en 2015",
                         zoom = 6)
# ajout de la legende des classes (position par defaut)
map <- add_legende_classes(map = map,
                           titre = "Variation ann.moy. 2010-2015",
                           zoom = 6)

# modification de la position de la legende des ronds et du niveau de zoom
map <- add_legende_ronds(map = map,
                         titre = "Population en 2015",
                         lng = 8,
                         lat = 44,
                         zoom = 8)
# modification de la position de la legende des classes et du niveau de zoom
map <- add_legende_classes(map = map,
                           titre = "Variation ann.moy. 2010-2015",
                           lng = 7.5,
                           lat = 43.5,
                           zoom = 8)

# ajout d'une source a la carte
map <- add_source(map = map,
                  source = "Source : INSEE - RP2016")
# ajout d'un titre a la carte
map <- add_titre(map = map,
                 titre = "Population des departements de la region Provence-Alpes-Cote d'Azur en 2015 et son evolution depuis 2010")

# affichage de la palette par defaut
# affiche_palette(nomPalette = "defaut")
recup_palette(stylePalette = "defaut", nbPos = 3)

# affichage de la palette Insee_Rouge
recup_palette(stylePalette = "Insee_Rouge", nbPos = 3)

# modification du style de la palette et de la bordure des ronds
map <- set_couleur_classes(map = map,
                           stylePalette = "Insee_Rouge",
                           colBorder = "grey")

# modification de l'opacite de la representation elargie
map <- set_opacite_elargi(map = map,
                          opacite = 0.3)

map

## ----plotSaphirs, fig.height = 6, fig.width = 7-------------------------------
library(oceanis)

# chargement des donnees
donnees_biloc_saphirs <- rio::import(file = system.file("data/donnees_biloc_saphirs.rda",
                                                        package = "oceanis"))

# import du fond des regions
regm <- sf::st_read(dsn = system.file("extdata",
                                      "reg_francemetro_2018.shp",
                                      package = "oceanis"),
                    quiet = TRUE,
                    stringsAsFactors = FALSE)
# import du fond de France metropolitaine
fram <- sf::st_read(dsn = system.file("extdata",
                                      "francemetro_2018.shp",
                                      package = "oceanis"), 
                    quiet = TRUE,
                    stringsAsFactors = FALSE)
# import du fond des pays
paysm <- sf::st_read(dsn = system.file("extdata",
                                       "paysf_2018.shp",
                                       package = "oceanis"),
                     quiet = TRUE,
                     stringsAsFactors = FALSE)
# import du fond de mer
merm <- sf::st_read(dsn = system.file("extdata",
                                      "merf_2018.shp",
                                      package = "oceanis"),
                    quiet = TRUE,
                    stringsAsFactors = FALSE)

# affichage de la carte
fond_saphirs <- plot_saphirs(data = donnees_biloc_saphirs,
                             fondMaille = regm,
                             typeMaille = "REG",
                             idDataDepart = "REG_DEPART",
                             idDataArrivee = "REG_ARRIVEE",
                             varFlux = "MIGR",
                             direction = "Ent",
                             titreLeg = "Entrees",
                             xLeg = 1100000,
                             yLeg = 6470000,
                             titreCarte = "Migrations residentielles vers l'Ile-de-France",
                             sourceCarte = "Source : INSEE - RP2016",
                             colBorder = "transparent",
                             colBorderMaille = "grey")

# construction de la table des etiquettes
etiquettes <- coordonnees_etiquettes(fondMaille = regm,
                                     listeCode = as.character(regm$CODE))
# modification des valeurs (latitude Y, longitude X, taille, couleur et style de police)
etiquettes[etiquettes$CODE=="24","Y"] <- 6680000
etiquettes[etiquettes$CODE=="27","Y"] <- 6660000
etiquettes[etiquettes$CODE=="28","X"] <- 410000
etiquettes[etiquettes$CODE=="32","Y"] <- 7015000
etiquettes[etiquettes$CODE=="44","X"] <- 955000
etiquettes[etiquettes$CODE=="52","X"] <- 330000
etiquettes[etiquettes$CODE=="52","Y"] <- 6700000
etiquettes[etiquettes$CODE=="53","X"] <- 215000
etiquettes[etiquettes$CODE=="75","Y"] <- 6420000
etiquettes[etiquettes$CODE=="76","Y"] <- 6270000
etiquettes[etiquettes$CODE=="84","Y"] <- 6455000
etiquettes[etiquettes$CODE=="93","Y"] <- 6290000
etiquettes[etiquettes$CODE=="94","Y"] <- 6120000
etiquettes[etiquettes$CODE!="11","TAILLE"] <- 0.6
etiquettes[etiquettes$CODE=="11","COL"] <- "#002D7F"
etiquettes[etiquettes$CODE!="11","FONT"] <- 1

# affichage de la carte
fond_saphirs <- plot_saphirs(data = donnees_biloc_saphirs,
                             fondMaille = regm,
                             typeMaille = "REG",
                             idDataDepart = "REG_DEPART",
                             idDataArrivee = "REG_ARRIVEE",
                             varFlux = "MIGR",
                             direction = "Ent",
                             titreLeg = "Entrees",
                             xLeg = 1150000,
                             yLeg = 6470000,
                             titreCarte = "Migrations residentielles vers l'Ile-de-France",
                             sourceCarte = "Source : INSEE - RP2016",
                             etiquettes = etiquettes,
                             colBorder = "transparent",
                             colBorderMaille = "grey")

# ajout de colonnes dans les fonds pour modifier leur apparence
# couleur de remplissage : COL
# couleur des contours : BORDER
# epaisseur des contours : EPAISEEUR
merm$COL <- "lightsteelblue"
merm$BORDER <- "lightsteelblue"
paysm$COL <- "gray"
paysm$BORDER <- "white"
fram$BORDER <- "darkgray"
fram$EPAISSEUR <- 2

# creation des listes des fonds d'habillage, en-dessous et au-dessus de l'analyse
fondSousAnalyse <- list(merm,paysm)
fondSurAnalyse <- list(fram)

# affichage de la carte
fond_saphirs <- plot_saphirs(data = donnees_biloc_saphirs,
                             fondMaille = regm,
                             fondSousAnalyse = fondSousAnalyse,
                             fondSurAnalyse = fondSurAnalyse,
                             typeMaille = "REG",
                             idDataDepart = "REG_DEPART",
                             idDataArrivee = "REG_ARRIVEE",
                             varFlux = "MIGR",
                             direction = "Ent",
                             titreLeg = "Entrees",
                             xLeg = 1150000,
                             yLeg = 6470000,
                             titreCarte = "Migrations residentielles vers l'Ile-de-France",
                             sourceCarte = "Source : INSEE - RP2016",
                             etiquettes = etiquettes,
                             colBorder = "transparent",
                             colBorderMaille = "grey")


## ----exportImage, eval = FALSE------------------------------------------------
#  
#  library(grDevices)
#  
#  jpeg(filename = "sortie.jpg",
#       quality = 100,
#       width = 16,
#       height = 18,
#       units = "cm",
#       res = 120)
#  
#  # affichage de la carte
#  fond_saphirs <- plot_saphirs(data = donnees_biloc_saphirs,
#                               fondMaille = regm,
#                               fondSousAnalyse = fondSousAnalyse,
#                               fondSurAnalyse = fondSurAnalyse,
#                               typeMaille = "REG",
#                               idDataDepart = "REG_DEPART",
#                               idDataArrivee = "REG_ARRIVEE",
#                               varFlux = "MIGR",
#                               direction = "Ent",
#                               titreLeg = "Entrees",
#                               xLeg = 1150000,
#                               yLeg = 6470000,
#                               titreCarte = "Migrations residentielles vers l'Ile-de-France",
#                               sourceCarte = "Source : INSEE - RP2016",
#                               etiquettes = etiquettes,
#                               colBorder = "transparent",
#                               colBorderMaille = "grey")
#  
#  dev.off()
#  
#  pdf(file = "sortie.pdf",
#      width = 10,
#      height = 10)
#  
#  # affichage de la carte
#  fond_saphirs <- plot_saphirs(data = donnees_biloc_saphirs,
#                               fondMaille = regm,
#                               fondSousAnalyse = fondSousAnalyse,
#                               fondSurAnalyse = fondSurAnalyse,
#                               typeMaille = "REG",
#                               idDataDepart = "REG_DEPART",
#                               idDataArrivee = "REG_ARRIVEE",
#                               varFlux = "MIGR",
#                               direction = "Ent",
#                               titreLeg = "Entrees",
#                               xLeg = 1150000,
#                               yLeg = 6470000,
#                               titreCarte = "Migrations residentielles vers l'Ile-de-France",
#                               sourceCarte = "Source : INSEE - RP2016",
#                               etiquettes = etiquettes,
#                               colBorder = "transparent",
#                               colBorderMaille = "grey")
#  
#  dev.off()

## ----exportQgis, eval = FALSE-------------------------------------------------
#  library(oceanis)
#  library(leaflet)
#  
#  # chargement des donnees
#  donnees_monoloc <- rio::import(file = system.file("data/donnees_monoloc.rda",
#                                                    package = "oceanis"))
#  
#  # import du fond des departements
#  depm <- sf::st_read(dsn = system.file("extdata",
#                                        "dep_francemetro_2018.shp",
#                                        package = "oceanis"),
#                      quiet = TRUE,
#                      stringsAsFactors = FALSE)
#  # import du fond des regions
#  regm <- sf::st_read(dsn = system.file("extdata",
#                                        "reg_francemetro_2018.shp",
#                                        package = "oceanis"),
#                      quiet = TRUE,
#                      stringsAsFactors = FALSE)
#  
#  # affichage de la carte avec des rayons de ronds plus grands
#  map <- leaflet_ronds_classes(data = donnees_monoloc,
#                               fondMaille = depm[depm$REG == "93",],
#                               fondMailleElargi = depm,
#                               fondSuppl = regm,
#                               idData = "COD_DEP",
#                               varVolume = "POP_2015",
#                               varRatio = "VAR_AN_MOY",
#                               rayonRond = 29000,
#                               rapportRond = NULL)
#  
#  # modification de la position de la legende des ronds et du niveau de zoom
#  map <- add_legende_ronds(map = map,
#                           titre = "Population en 2015",
#                           lng = 8.5,
#                           lat = 45,
#                           zoom = 8)
#  # modification de la position de la legende des classes et du niveau de zoom
#  map <- add_legende_classes(map = map,
#                             titre = "Variation ann.moy. 2010-2015",
#                             lng = 8,
#                             lat = 44.5,
#                             zoom = 8)
#  # ajout d'une source a la carte
#  map <- add_source(map = map,
#                    source = "Source : INSEE - RP2016")
#  # ajout d'un titre a la carte
#  map <- add_titre(map = map,
#                   titre = "Population des departements de la region Provence-Alpes-Cote d'Azur en 2015 et son evolution depuis 2010")
#  # modification de la couleur de bordure des ronds
#  map <- set_bordure_ronds(map = map,
#                           colBorderPos = "grey")
#  # modification du style de la palette
#  map <- set_couleur_classes(map = map,
#                             stylePalette = "Insee_Rouge")
#  # modification de l'opacite de la representation elargie
#  map <- set_opacite_elargi(map,
#                            opacite = 0.6)
#  
#  export_qgis_ronds_classes(map,
#                            cheminDossier = getwd(),
#                            nomFichier = "export_carte_rp_ac",
#                            titre1 = "Population des departements de la region Provence-Alpes-Cote d'Azur en 2015 et son evolution depuis 2010",
#                            titre2 = "",
#                            source = "Source : INSEE - RP2016")