Cordonnées LAMBERT à convertir

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.

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

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:

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 ?

Bonjour,

Merci pour ton aide.

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