Afficher des shapefiles ensemble

Bonjour à tous,

Novice sur R, je commence une étude ou je vais être amené à programmer dans ce langage. Dans un premier temps, j'essaie de prendre en main les packages permettant de faire du SIG.

Je rencontre d'ores et déjà un petit problème (qui n'en est sans doute pas vraiment un). Je m'explique :

J'ai 4 shapefiles de base.
Je les ouvre sur R avec à l'aide de la fonction read.OGR de la library(sp) :
Pièce jointe 235659

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18

#Telechargement des packages de cartographie   lib<-c( "maptools", "rgdal", "rgeos", "sp") for (i in 1:length(lib))install.packages(lib[i],dependencies = T) library(rgdal) library(maptools) library(sp)   # Espace de travail path<-'D:/Users/bourgetr/Desktop/STEP_StationsQualité/Distance_STEP_STQUAL'     #Importation des couches   UH<-readOGR(dsn='D:/Users/bourgetr/Desktop/STEP_StationsQualité/Distance_STEP_STQUAL', layer="UH_SAV") TRONCON<-readOGR(dsn = path,layer = 'TRONCON_HYDRO_DSAV') STEP<-readOGR(dsn = path,layer = 'STEP_2012') STQUAL<-readOGR(dsn = path,layer = 'STATIONS QUALITE DSAV')

---

Avant de projeter ces couches je vérifie si elles ont le même SCR :

Code:

---

1 2	proj4string(STQUAL)==proj4string(TRONCON) [1] FALSE

---

J'ai donc 2 couches dans le même SCR : "UH" et "TRONCON" et les 2 autres sont dans des systèmes de projection différents.
Je décide de définir un nouveau système de coordonnées pour les 4 couches :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14

#Definition du système de coordonnées EPSG<-make_EPSG()# Liste des SCR EPSG_Lambert93<- EPSG[EPSG$code==2154,][1,]#Récuperation des coordonnées du SCR ciblé   proj4string(UH)<-"+proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs" proj4string(STEP)<-"+proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs" proj4string(TRONCON)<-"+proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs" proj4string(STQUAL)<-"+proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"   #Reprojection des couches dans le SCR UH<-spTransform(UH,CRS(EPSG_Lambert93$prj4)) STEP<-spTransform(STEP,CRS(EPSG_Lambert93$prj4)) STQUAL<-spTransform(STQUAL,CRS(EPSG_Lambert93$prj4)) TRONCON<-spTransform(TRONCON,CRS(EPSG_Lambert93$prj4))

---

Cependant, je ne peux toujours pas afficher mes 4 couches simultanément après cette manip qui doit comporter une ou plusieurs erreurs.
J'ai cherché mais je n'ai rien trouvé.
Est-ce que quelqu'un pourrait m'aider ?!

Merci d'avance.

Tristan