This function is a wrapper for sf::st_crs, unless
threating numeric character
strings as integers, and
accepting also UTM timezones, paths of spatial files and paths of
text files containing WKT like .prj (see details) .
st_crs2(x, ...)
x | numeric, character, or object of class sf or sfc, being:
|
---|---|
... | other parameters passed to sf::st_crs. |
An object of class crs of length 2.
See sf::st_crs for details.
License: GPL 3.0
## CRS from EPSG st_crs2(32609)#> Coordinate Reference System: #> User input: EPSG:32609 #> wkt: #> PROJCS["WGS 84 / UTM zone 9N", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",-129], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["metre",1, #> AUTHORITY["EPSG","9001"]], #> AXIS["Easting",EAST], #> AXIS["Northing",NORTH], #> AUTHORITY["EPSG","32609"]]st_crs2("EPSG:32609")#> Coordinate Reference System: #> User input: EPSG:32609 #> wkt: #> PROJCS["WGS 84 / UTM zone 9N", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",-129], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["metre",1, #> AUTHORITY["EPSG","9001"]], #> AXIS["Easting",EAST], #> AXIS["Northing",NORTH], #> AUTHORITY["EPSG","32609"]]## CRS from UTM zone st_crs2(9)#> Coordinate Reference System: #> User input: EPSG:32609 #> wkt: #> PROJCS["WGS 84 / UTM zone 9N", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",-129], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["metre",1, #> AUTHORITY["EPSG","9001"]], #> AXIS["Easting",EAST], #> AXIS["Northing",NORTH], #> AUTHORITY["EPSG","32609"]]st_crs2("09")#> Coordinate Reference System: #> User input: EPSG:32609 #> wkt: #> PROJCS["WGS 84 / UTM zone 9N", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",-129], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["metre",1, #> AUTHORITY["EPSG","9001"]], #> AXIS["Easting",EAST], #> AXIS["Northing",NORTH], #> AUTHORITY["EPSG","32609"]]st_crs2("9N")#> Coordinate Reference System: #> User input: EPSG:32609 #> wkt: #> PROJCS["WGS 84 / UTM zone 9N", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",-129], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["metre",1, #> AUTHORITY["EPSG","9001"]], #> AXIS["Easting",EAST], #> AXIS["Northing",NORTH], #> AUTHORITY["EPSG","32609"]]st_crs2("09S")#> Coordinate Reference System: #> User input: EPSG:32709 #> wkt: #> PROJCS["WGS 84 / UTM zone 9S", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",-129], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",10000000], #> UNIT["metre",1, #> AUTHORITY["EPSG","9001"]], #> AXIS["Easting",EAST], #> AXIS["Northing",NORTH], #> AUTHORITY["EPSG","32709"]]#> [1] "PROJCS[\"WGS 84 / UTM zone 9N\",GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4326\"]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"latitude_of_origin\",0],PARAMETER[\"central_meridian\",-129],PARAMETER[\"scale_factor\",0.9996],PARAMETER[\"false_easting\",500000],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH],AUTHORITY[\"EPSG\",\"32609\"]]"st_crs2(wkt_32n)#> Coordinate Reference System: #> User input: PROJCS["WGS 84 / UTM zone 9N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-129],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","32609"]] #> wkt: #> PROJCS["WGS 84 / UTM zone 9N", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",-129], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["metre",1, #> AUTHORITY["EPSG","9001"]], #> AXIS["Easting",EAST], #> AXIS["Northing",NORTH], #> AUTHORITY["EPSG","32609"]]#> Coordinate Reference System: #> User input: PROJCS["WGS 84 / UTM zone 9N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-129],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","32609"]] #> wkt: #> PROJCS["WGS 84 / UTM zone 9N", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",-129], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["metre",1, #> AUTHORITY["EPSG","9001"]], #> AXIS["Easting",EAST], #> AXIS["Northing",NORTH], #> AUTHORITY["EPSG","32609"]]## CRS from spatial file path raster_path <- system.file( "extdata/out/S2A2A_20190723_022_Barbellino_BOA_10.tif", package="sen2r" ) vector_path <- system.file( "extdata/vector/barbellino.geojson", package="sen2r" ) st_crs2(raster_path)#> Coordinate Reference System: #> User input: PROJCS["WGS 84 / UTM zone 32N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",9],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","32632"]] #> wkt: #> PROJCS["WGS 84 / UTM zone 32N", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",9], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["metre",1, #> AUTHORITY["EPSG","9001"]], #> AXIS["Easting",EAST], #> AXIS["Northing",NORTH], #> AUTHORITY["EPSG","32632"]]st_crs2(vector_path)#> Coordinate Reference System: #> User input: 32632 #> wkt: #> PROJCS["WGS 84 / UTM zone 32N", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",9], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["metre",1, #> AUTHORITY["EPSG","9001"]], #> AXIS["Easting",EAST], #> AXIS["Northing",NORTH], #> AUTHORITY["EPSG","32632"]]#> Coordinate Reference System: #> User input: PROJCS["WGS 84 / UTM zone 32N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",9],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","32632"]] #> wkt: #> PROJCS["WGS 84 / UTM zone 32N", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",9], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["metre",1, #> AUTHORITY["EPSG","9001"]], #> AXIS["Easting",EAST], #> AXIS["Northing",NORTH], #> AUTHORITY["EPSG","32632"]]#> Coordinate Reference System: #> User input: 32632 #> wkt: #> PROJCS["WGS 84 / UTM zone 32N", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",9], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["metre",1, #> AUTHORITY["EPSG","9001"]], #> AXIS["Easting",EAST], #> AXIS["Northing",NORTH], #> AUTHORITY["EPSG","32632"]]# \donttest{ ## CRS from PROJ.4 string # (avoid using this with PROJ >= 6!) st_crs2("+init=epsg:32609") # this makes use of the EPSG code#> Coordinate Reference System: #> User input: EPSG:32609 #> wkt: #> PROJCS["WGS 84 / UTM zone 9N", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",-129], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["metre",1, #> AUTHORITY["EPSG","9001"]], #> AXIS["Easting",EAST], #> AXIS["Northing",NORTH], #> AUTHORITY["EPSG","32609"]]st_crs2("+proj=utm +zone=9 +datum=WGS84 +units=m +no_defs")#> Warning: Using PROJ.4 strings is deprecated with PROJ >= 6 (see #> https://www.r-spatial.org/r/2020/03/17/wkt.html).#> Coordinate Reference System: #> User input: +proj=utm +zone=9 +datum=WGS84 +units=m +no_defs #> wkt: #> PROJCS["UTM Zone 9, Northern Hemisphere", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",-129], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["Meter",1]]#> Coordinate Reference System: #> User input: +proj=utm +zone=32 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 #> wkt: #> PROJCS["UTM Zone 32, Northern Hemisphere", #> GEOGCS["WGS 84", #> DATUM["WGS_1984", #> SPHEROID["WGS 84",6378137,298.257223563, #> AUTHORITY["EPSG","7030"]], #> AUTHORITY["EPSG","6326"]], #> PRIMEM["Greenwich",0, #> AUTHORITY["EPSG","8901"]], #> UNIT["degree",0.0174532925199433, #> AUTHORITY["EPSG","9122"]], #> AUTHORITY["EPSG","4326"]], #> PROJECTION["Transverse_Mercator"], #> PARAMETER["latitude_of_origin",0], #> PARAMETER["central_meridian",9], #> PARAMETER["scale_factor",0.9996], #> PARAMETER["false_easting",500000], #> PARAMETER["false_northing",0], #> UNIT["Meter",1]]# }