Message posté par : A.Duranel (arnaud_duranel(a)yahoo.fr)
----------------------------------------
Bonjour,
J'essaye de projeter des données en Lambert 2 étendu vers le Lambert 93 avec le
package sf pour R, mais je n'obtiens pas les mêmes résultats qu'avec Circé
(v.5.4.6).
Exemple :
X: 655000 Y : 2052940 en Lambert 2 étendu
Avec Circé on obtient X : 703112.347 Y : 6486136.042 en LAMBERT-93
Avec sf :
-----------------
Code :
library(sf)
library(magrittr)
data.frame(X=655000, Y=2052940) %>%
st_as_sf(coords = c("X", "Y"), crs=st_crs("IGNF:LAMBE"))
%>%
st_transform(st_crs("IGNF:LAMB93")) %>%
st_coordinates() %>%
sprintf(fmt="%.3f")
-----------------
X :703112.908 Y : 6486136.425
On obtient exactement le même résultat si on remplace les codes IGNF par les codes EPSG
("IGNF:LAMBE" par 27572 et/ou "IGNF:LAMB93" par 2154), alors qu'on
ne devrait pas je pense ?
Merci pour votre aide.
----------------------------------------
Le message est situé
https://georezo.net/forum/viewtopic.php?pid=370636#p370636
Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net