terra (version 1.1-17)

crs: Get or set a coordinate reference system

Description

Get or set the coordinate reference system (CRS), also referred to as a "projection" of a SpatRaster or SpatVector object.

Setting a new CRS does not change the data itself, it just changes the label. So you should only set the CRS of a dataset (if it does not come with one) to what it *is*, not to what you would *like it to be*. See project to *transform* spatial from one CRS to another.

Usage

# S4 method for SpatRaster
crs(x, proj4=FALSE, describe=FALSE)
# S4 method for SpatVector
crs(x, proj4=FALSE, describe=FALSE)
# S4 method for SpatRaster
crs(x)<-value
# S4 method for SpatVector
crs(x)<-value

Arguments

x

SpatRaster or SpatVector

proj4

logical. If TRUE the crs is returned in PROJ.4 notation. But note that GDAL now only supports the WGS84 and NAD83 datums with this notation

describe

logical. If TRUE the name, EPSG code, and the name and extent of the area of use are returned if known

value

character string describing a coordinate reference system. This can be in a WKT format, as a EPSG code, or a PROJ.4 "+" format (but see Note)

Value

character or modified SpatRaster/Vector

Examples

# NOT RUN {
r <- rast()
crs(r)
crs(r, describe=TRUE, proj4=TRUE)

crs(r) <- "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100 +ellps=WGS84"
crs(r)

# You can also use epsg codes
crs(r)  <- "epsg:25831"
crs(r, describe=TRUE)$area
# }