Discussion :

[Tristan8776]

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) :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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