Data set of available spatial references (CRS) in tibble format.
Format
A tibble with 9,608 rows and fields:
- projtype
Projection type (
"ProjectedCoordinateSystems","GeographicCoordinateSystems"or"VerticalCoordinateSystems").- wkid
Well-Known ID.
- latestWkid
Most recent
wkid, ifwkidis deprecated.- authority
wkidauthority (Esri or EPSG).- deprecated
Logical indicating whether
wkidis deprecated.- description
Human-readable description of the
wkid.- areaname
Use area of the
wkid.- wkt
Representation of
wkidin Well-Known Text (WKT). Useful when working with sf or terra.
Details
This data set is useful when using the outsr argument.
Some projection IDs have changed over time. For example, Web Mercator
wkid = 102100 is deprecated and is currently wkid = 3857. However, both
values work and return similar results.
See also
Other datasets:
arc_categories
Examples
# \donttest{
# Get all possible values.
data("arc_spatial_references")
arc_spatial_references
#> # A tibble: 9,608 × 8
#> projtype wkid latestWkid authority deprecated description areaname wkt
#> <chr> <int> <int> <chr> <lgl> <chr> <chr> <chr>
#> 1 ProjectedCo… 2000 2000 EPSG FALSE Anguilla 1… Anguill… "PRO…
#> 2 ProjectedCo… 2001 2001 EPSG FALSE Antigua 19… Antigua… "PRO…
#> 3 ProjectedCo… 2002 2002 EPSG FALSE Dominica 1… Dominic… "PRO…
#> 4 ProjectedCo… 2003 2003 EPSG FALSE Grenada 19… Grenada… "PRO…
#> 5 ProjectedCo… 2004 2004 EPSG FALSE Montserrat… Montser… "PRO…
#> 6 ProjectedCo… 2005 2005 EPSG FALSE St. Kitts … St Kitt… "PRO…
#> 7 ProjectedCo… 2006 2006 EPSG FALSE St. Lucia … St Luci… "PRO…
#> 8 ProjectedCo… 2007 2007 EPSG FALSE St. Vincen… St Vinc… "PRO…
#> 9 ProjectedCo… 2008 2008 EPSG TRUE NAD 1927 S… Canada … "PRO…
#> 10 ProjectedCo… 2009 2009 EPSG FALSE NAD 1927 C… Canada … "PRO…
#> # ℹ 9,598 more rows
# Request with deprecated Web Mercator.
library(dplyr)
wkid <- arc_spatial_references |>
filter(latestWkid == 3857 & deprecated) |>
slice(1)
glimpse(wkid)
#> Rows: 1
#> Columns: 8
#> $ projtype <chr> "ProjectedCoordinateSystems"
#> $ wkid <int> 102100
#> $ latestWkid <int> 3857
#> $ authority <chr> "Esri"
#> $ deprecated <lgl> TRUE
#> $ description <chr> "WGS 1984 Web Mercator (auxiliary sphere)"
#> $ areaname <chr> "World - 85~S to 85~N"
#> $ wkt <chr> "PROJCS[\"WGS_1984_Web_Mercator_Auxiliary_Sphere\",GEOGCS[…
add <- arc_geo("London, United Kingdom", outsr = wkid$wkid)
# Note values for `lat`, `lon` and `wkid`. `latestWkid` gives the current
# valid `wkid`.
add |>
select(lat, lon, wkid, latestWkid) |>
glimpse()
#> Rows: 1
#> Columns: 4
#> $ lat <dbl> 6711544
#> $ lon <dbl> -14215.35
#> $ wkid <int> 102100
#> $ latestWkid <int> 3857
# Try with `sf`.
try(sf::st_crs(wkid$wkid))
#> Warning: GDAL Error 1: PROJ: proj_create_from_database: crs not found
#> Coordinate Reference System: NA
# Try the latest WKID.
try(sf::st_crs(wkid$latestWkid))
#> Coordinate Reference System:
#> User input: EPSG:3857
#> wkt:
#> PROJCRS["WGS 84 / Pseudo-Mercator",
#> BASEGEOGCRS["WGS 84",
#> ENSEMBLE["World Geodetic System 1984 ensemble",
#> MEMBER["World Geodetic System 1984 (Transit)"],
#> MEMBER["World Geodetic System 1984 (G730)"],
#> MEMBER["World Geodetic System 1984 (G873)"],
#> MEMBER["World Geodetic System 1984 (G1150)"],
#> MEMBER["World Geodetic System 1984 (G1674)"],
#> MEMBER["World Geodetic System 1984 (G1762)"],
#> MEMBER["World Geodetic System 1984 (G2139)"],
#> ELLIPSOID["WGS 84",6378137,298.257223563,
#> LENGTHUNIT["metre",1]],
#> ENSEMBLEACCURACY[2.0]],
#> PRIMEM["Greenwich",0,
#> ANGLEUNIT["degree",0.0174532925199433]],
#> ID["EPSG",4326]],
#> CONVERSION["Popular Visualisation Pseudo-Mercator",
#> METHOD["Popular Visualisation Pseudo Mercator",
#> ID["EPSG",1024]],
#> PARAMETER["Latitude of natural origin",0,
#> ANGLEUNIT["degree",0.0174532925199433],
#> ID["EPSG",8801]],
#> PARAMETER["Longitude of natural origin",0,
#> ANGLEUNIT["degree",0.0174532925199433],
#> ID["EPSG",8802]],
#> PARAMETER["False easting",0,
#> LENGTHUNIT["metre",1],
#> ID["EPSG",8806]],
#> PARAMETER["False northing",0,
#> LENGTHUNIT["metre",1],
#> ID["EPSG",8807]]],
#> CS[Cartesian,2],
#> AXIS["easting (X)",east,
#> ORDER[1],
#> LENGTHUNIT["metre",1]],
#> AXIS["northing (Y)",north,
#> ORDER[2],
#> LENGTHUNIT["metre",1]],
#> USAGE[
#> SCOPE["Web mapping and visualisation."],
#> AREA["World between 85.06°S and 85.06°N."],
#> BBOX[-85.06,-180,85.06,180]],
#> ID["EPSG",3857]]
# Or use WKT.
try(sf::st_crs(wkid$wkt))
#> Coordinate Reference System:
#> User input: PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator_Auxiliary_Sphere"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Standard_Parallel_1",0.0],PARAMETER["Auxiliary_Sphere_Type",0.0],UNIT["Meter",1.0]]
#> wkt:
#> PROJCRS["WGS 84 / Pseudo-Mercator",
#> BASEGEOGCRS["WGS 84",
#> DATUM["World Geodetic System 1984",
#> ELLIPSOID["WGS 84",6378137,298.257223563,
#> LENGTHUNIT["metre",1]],
#> ID["EPSG",6326]],
#> PRIMEM["Greenwich",0,
#> ANGLEUNIT["Degree",0.0174532925199433]]],
#> CONVERSION["unnamed",
#> METHOD["Popular Visualisation Pseudo Mercator",
#> ID["EPSG",1024]],
#> PARAMETER["Latitude of natural origin",0,
#> ANGLEUNIT["Degree",0.0174532925199433],
#> ID["EPSG",8801]],
#> PARAMETER["Longitude of natural origin",0,
#> ANGLEUNIT["Degree",0.0174532925199433],
#> ID["EPSG",8802]],
#> PARAMETER["False easting",0,
#> LENGTHUNIT["metre",1],
#> ID["EPSG",8806]],
#> PARAMETER["False northing",0,
#> LENGTHUNIT["metre",1],
#> ID["EPSG",8807]]],
#> CS[Cartesian,2],
#> AXIS["(E)",east,
#> ORDER[1],
#> LENGTHUNIT["metre",1,
#> ID["EPSG",9001]]],
#> AXIS["(N)",north,
#> ORDER[2],
#> LENGTHUNIT["metre",1,
#> ID["EPSG",9001]]]]
# }