Discussion :

[marou1991]

Bonjour,

J'ai trouvé un fichier avec des cordonnées avec un format que je connais pas dit 'LAMBERT' avez-vous une idée sur ce format et une manière de les convertir en cordonnées GPS (Lat, long) s'il vous plait ?

Voilà un extrait de mes données :

LAMBERT_X	LAMBERT_Y
880792.7	6544216
882609.4	6542611
882620.0	6542692
882596.8	6542755
882534.7	6542907
882650.2	6542927
882534.7	6542907
904975.9	6579348
904804.2	6579509
905028.5	6579380

Merci de votre aide.

[olivier.decourt]

Bonjour Marou.
En cherchant rapidement sur Internet j'ai trouvé le CRS de la projection Lambert ici. Du coup, on peut aisément créer un objet de type SpatialPoints (package {sp}) avec ce CRS, puis faire la conversion en coordonnées GPS classiques (WGS84 si je ne me trompe pas).

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
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
library(sp) # objets géographiques
library(leaflet) # visualisation par cartes interactives
 
LAMBERT_X <- c(880792.7 ,
               882609.4 ,
               882620.0 ,
               882596.8 ,
               882534.7 ,
               882650.2 ,
               882534.7 ,
               904975.9 ,
               904804.2 ,
               905028.5 )
LAMBERT_Y <- c(6544216,
               6542611,
               6542692,
               6542755,
               6542907,
               6542927,
               6542907,
               6579348,
               6579509,
               6579380)
# on transforme en objet sp avec la projection Lambert
coord <- SpatialPoints(data.frame(LAMBERT_X, LAMBERT_Y),
                       proj4string=CRS(paste("+proj=lcc +lat_1=44 +lat_2=49",
                                             "+lat_0=46.5 +lon_0=3 +x_0=700000",
                                             "+y_0=6600000 +ellps=GRS80",
                                             "+units=m +no_defs")))
# puis en coordonnées WGS84
coord2 <- sp::spTransform(coord, 
            CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +"))
 
# on vérifie avec {leaflet} où se trouvent ces points
coord2 %>% 
  leaflet() %>% 
  addTiles() %>% 
  addMarkers()

Les coordonnées que tu indiques sont dans le département de l'Ain ; est-ce bien ce qui est prévu ?

[marou1991]

Bonjour,

Merci pour ton aide.

ca marche parfaitement bien. Oui c'est le département de l'Ain.