Discussion :

[Anne-Laure P]

Bonjour,

Je travaille actuellement sur la géolocalisation d'adresse de sujets.

La première étape est de récupérer les coordonnées de chaque sujet à partir de leur adresse. Je passe par "batchgéocodeur" qui fournit en sortie les coordonnées en degrés décimaux dans un tableau excel.

Ensuite, je les importe dans R pour transformer les coordonnées de WGS84 à Lambert 93.

J'utilise le script suivant :

Code :

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

1
2
3
4
5
6
7
8
pts<-read.delim("WGS84.txt")
ptslines<-SpatialPointsDataFrame(pts[,1:2],pts)                           
proj4string(ptslines)<-CRS("+proj=longlat +datum=WGS84 +ellps=WGS84")     
ptslinesL93<-spTransform(ptslines,CRS("+init=epsg:4171"))                 
WGS84L93<-ptslinesL93
WGS84L93<-as.data.frame(WGS84L93)
names(WGS84L93)=c("WGS84X","WGS84Y","L93X","L93Y")
write.delim(WGS84L93, "WGS84L93.txt")

Seulement lorsque je compare les coordonnées obtenues en Lambert93 à celles lues sur géoportail, cela ne correspond pas.

Est-ce que quelqu'un a une idée ?

Merci d'avance !

[ced]

Bonjour,

C'est sûrement parce que proj4, par défaut, ne prend pas en compte la grille de conversion NTF->RGF93, alors que le Géoportail oui. Et notamment le Lambert93 du code EPSG 4171.
Plus d'info sur cette page : http://geodesie.ign.fr/index.php?p=7...site_lambert93

Vous pouvez essayer avec :

Code :

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

ptslinesL93<-spTransform(ptslines,CRS("+init=LAMB93"))

Mais je doute que la version proj4 embarquée dans R ait nativement la grille de correction... Et pour l'ajouter, je ne sais pas comment faire.

[Anne-Laure P]

Bonjour,

Merci pour votre réponse. Désolée de ne pas avoir répondu plutôt, je croyais avoir basculer les réponses sur ma boite mail mais j'ai du faire une erreur.

J'ai relancé le script en utilisant "+init=LAMB93" mais cela ne fonctionne pas.
Mais merci beaucoup pour votre éclairage concernant proj4 et la grille de correction. Au moins, je comprends ce qui bloque.

En revanche, connaissez-vous un autre outil que R pour pouvoir transformer une liste de coordonnées des degrés décimaux au lambert93 ?
Qgis n'accepte pour l'instant que l'import de degrés sexagésimaux (j'ai encore beaucoup de travail pour comprendre tout ça...!!!)

Merci d'avance,
Bonne journée.

[ced]

La référence pour ça reste le logiciel Circé, fourni par l'IGN.

[Anne-Laure P]

Merci pour le tuyau. Je vais télécharger Circé et voir ce que je peux faire avec.
Comme mon voisin de bureau peut faire la manip de transformation des coordonnées très facilement sur Arcgis, j'ai cherché l'équivalent sous Qgiq. J'ai trouvé Osgeo, mais son utilisation dépasse mes compétences actuelles...

Merci de tout ces conseils en tout cas !
Bonne journée !!!