4 Les systèmes de coordonnées

4.1 Consulter le système de coordonnées d’un objet

La fonction st_crs() permet de consulter le système de coordonnées utilisé par un objet sf.

library(sf) 
st_crs(x = com)

#> Coordinate Reference System:
#>   User input: RGF93 v1 / Lambert-93 
#>   wkt:
#> PROJCRS["RGF93 v1 / Lambert-93",
#>     BASEGEOGCRS["RGF93 v1",
#>         DATUM["Reseau Geodesique Francais 1993 v1",
#>             ELLIPSOID["GRS 1980",6378137,298.257222101,
#>                 LENGTHUNIT["metre",1]]],
#>         PRIMEM["Greenwich",0,
#>             ANGLEUNIT["degree",0.0174532925199433]],
#>         ID["EPSG",4171]],
#>     CONVERSION["Lambert-93",
#>         METHOD["Lambert Conic Conformal (2SP)",
#>             ID["EPSG",9802]],
#>         PARAMETER["Latitude of false origin",46.5,
#>             ANGLEUNIT["degree",0.0174532925199433],
#>             ID["EPSG",8821]],
#>         PARAMETER["Longitude of false origin",3,
#>             ANGLEUNIT["degree",0.0174532925199433],
#>             ID["EPSG",8822]],
#>         PARAMETER["Latitude of 1st standard parallel",49,
#>             ANGLEUNIT["degree",0.0174532925199433],
#>             ID["EPSG",8823]],
#>         PARAMETER["Latitude of 2nd standard parallel",44,
#>             ANGLEUNIT["degree",0.0174532925199433],
#>             ID["EPSG",8824]],
#>         PARAMETER["Easting at false origin",700000,
#>             LENGTHUNIT["metre",1],
#>             ID["EPSG",8826]],
#>         PARAMETER["Northing at false origin",6600000,
#>             LENGTHUNIT["metre",1],
#>             ID["EPSG",8827]]],
#>     CS[Cartesian,2],
#>         AXIS["easting (X)",east,
#>             ORDER[1],
#>             LENGTHUNIT["metre",1]],
#>         AXIS["northing (Y)",north,
#>             ORDER[2],
#>             LENGTHUNIT["metre",1]],
#>     USAGE[
#>         SCOPE["Engineering survey, topographic mapping."],
#>         AREA["France - onshore and offshore, mainland and Corsica (France métropolitaine including Corsica)."],
#>         BBOX[41.15,-9.86,51.56,10.38]],
#>     ID["EPSG",2154]]

4.2 Modifier le système de coordonnées d’un objet

La fonction st_transform() permet de changer le système de coordonnées d’un objet sf, de le reprojeter.

mf_map(com, expandBB = c(0, .12, 0, 0))
mf_graticule(x = com)
mf_title("RGF93 / Lambert-93")
# changement de projection
com_reproj <- st_transform(x = com, crs = "EPSG:3035")

mf_map(com_reproj, expandBB = c(0, .12, .0, 0))
mf_graticule(x = com_reproj)
mf_title("ETRS89-extended / LAEA Europe")

Si l’on souhaite projeter un objet spatial utilisant des coordonnées géographiques (lon/lat), le package crsuggest (Walker, 2022) propose des projections adaptées à l’emprise de l’objet.

Le site CRS Explorer met à disposition les références de très nombreux systèmes de coordonnées.