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