# A défaut de pouvoir installer le package "sf", on installe "sp".
# Pour plus d'information cf.doc officielle avec listes des packages dispo : https://cran.r-project.org/web/packages/available_packages_by_name.html
# Ensuite, on installe les librairies spatiales "gdal" et "rgeos" 
install.packages(pkgs="sp", dependencies=T)
install.packages(pkgs=c("rgdal", "rgeos"), dependencies=T)
	c("rgdal", "rgeos")

# i)Import des données spatiales vectorielles avec le package readOGR pour GDAL
# cf.https://cran.r-project.org/web/packages/rgdal/rgdal.pdf
# Pour les données de l'atelier : il s'agit de shapefiles des communes belges en libre accès sur: http://www.atlas-belgique.be/cms2/index.php?page=cartodata_fr

library(rgdal)
Commune <- readOGR(dsn="le chemin des données")
 # rem: on importe le shapefile des communes de Belgique stocké sur notre PC avec readOGR (car données vectorielles) en l'assigant à une variable nommée ici "Communes"

# ii) Visualisation des données avec les fonctions de base et des packages spécifiques dans sp
library(sp)
Commune
class(Commune)
head(Commune@data)
	#rem: on accéde aux données attributaires et visualise les 5 1eres lignes de la table
plot(Commune)
	#rem: on affiche graphiquement ces données

# iii)cartographier les données via intégration de données aléatoires
 dim(Commune)
	#rem: on affiche la dimension du jeux de données (=n)
 xUnif <- runif(n=589, min=5, max=10)
	#rem: on génère nos valeurs (en utilisant loi uniforme) cf.fonction "?runif"
  summary(xUnif)
  X11(); hist(xUnif, col="pink")
 	#rem: on génère l'histogramme de notre variable xUnif en rose
 xNorm <- rnorm(589, 7, 3)
  summary(xNorm)
  X11(); hist(xNorm, col="green")

#Maintenant les insère ds les données spatiales
 head(Commune@data)
 Commune@data <- data.frame(Commune@data, xUnif)
	#rem: on transforme notre vecteur ligne en vecteur collone et l'assigne à notre objet existant

?spplot
 spplot(obj=Commune, zcol = "xUnif")
	# rem: on cartographie grâce à la fonction spplot cf.doc officielle : https://cran.r-project.org/web/packages/sp/index.html


# Pour aller plus loin cf.données du Tuto: spGallery sur https://edzer.github.io/sp/
browseURL("https://edzer.github.io/sp/")
class(meuse)
x11(); plot(meuse)