4 Dec
2024
4 Dec
'24
11:56 a.m.
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
4 Dec
4 Dec
1:28 p.m.
Message posté par : Yves Jacolin
----------------------------------------
Bonjour,
Si justement (on est d'accord), on doit avoir des différences si on passe des codes
EPSG au code IGNF pour le lambert 2 étendue. Le code IGNF prend en compte la grille et est
donc plus précis.
Il faut comprendre comment fonctionne la reprojection : lambert 2 étendue -> WGS84
-> Lambert 93
WGS84 et Lambert 93 sont suffisament proche au niveau de l'ellipsoide pour que le
passage de l'un à l'autre se limite à un changement d'ordonnée à
l'origine.
Pour la transformation de lambert 2 étendue à WGS84, c'est complétement différent.
L'ellipsoide étant incompatible, le simple décalage de l'ordonnée à l'origine
provoque des imprécisions. Le fait de définir le code à IGNF:lambe implique une correction
via la grille NGF.
D'où vient le fait qu'on n'a pas de différence ? Est ce que R utilise la
grille pour la reprojection ?
Y.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=370641#p370641
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
2:30 p.m.
Message posté par : magr
----------------------------------------
Il faut utiliser sf_proj_pipeline()
-----------------
Code :
library(sf)
sf_proj_network(TRUE)
source <- data.frame(X = 655000, Y = 2052940) |>
st_as_sf(coords = c("X", "Y"),
crs = st_crs("EPSG:27572"))
source |>
st_transform(st_crs("EPSG:2154")) |>
st_coordinates() |>
sprintf(fmt="%.3f")
(p <- sf_proj_pipelines("EPSG:27572", "EPSG:2154"))
source |>
st_transform("EPSG:2154", pipeline = p[1, ]$definition) |>
st_coordinates() |>
sprintf(fmt="%.3f")
-----------------
On retombe sur les mêmes valeurs que Circée
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=370644#p370644
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
2:06 p.m.
Message posté par : A.Duranel
----------------------------------------
Merci pour ces précisions.
"D'où vient le fait qu'on n'a pas de différence ? Est ce que R utilise la
grille pour la reprojection ?"
J'ai l'impression que non vu que la conversion donne exactement le même résultat
qu'on parte de IGNF:LAMBE ou de EPSG:27572. Y a-t-il un moyen de faire en sorte que sf
utilise cette grille et produise le même résultat que Circé ?
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=370642#p370642
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
4:04 p.m.
Message posté par : A.Duranel
----------------------------------------
Fantastique, merci !
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=370648#p370648
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
311
days inactive
311
days old
4 comments
1 participants
participants (1)
-
geolibre＠georezo.net