Get or set a coordinate reference system (projection)

Get or set the coordinate reference system (CRS) of a Raster* object.

# S4 method for ANY
crs(x, asText=FALSE, ...)
# S4 method for Raster
wkt(obj)
crs(x, ...) <- value

projection(x, asText=TRUE)
projection(x) <- value

Arguments

x

Raster* or Spatial object

obj

Raster*, Spatial, or CRS object

asText

logical. If TRUE, the projection is returned as text. Otherwise a CRS-class object is returned

...

additional arguments. None implemented

value

CRS object or a character string describing a projection and datum in the PROJ.4 format

Value

Raster*, Spatial*, or character object

Note

crs replaces earlier function projection. For compatibility with sp you can use proj4string instead of crs. wkt returns the "well-known-text" representation of the crs.

See also

projectRaster, spTransform

Details

projections are done by with the PROJ library

Examples

r <- raster()
crs(r)
#> Coordinate Reference System:
#> Deprecated Proj.4 representation: +proj=longlat +datum=WGS84 
crs(r) <- "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100 +datum=WGS84"
crs(r)
#> Coordinate Reference System:
#> Deprecated Proj.4 representation:
#>  +proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100 +datum=WGS84 
w <- wkt(r)
w
#> [1] ""
cat(w, "\n")
#>