terra (version 0.6-9)

project: Change the coordinate reference system

Description

Change the coordinate reference system ("project") of a SpatVector or SpatRaster.

Usage

# S4 method for SpatVector
project(x, y, ...)
# S4 method for SpatRaster
project(x, y, method="bilinear", filename="", overwrite=FALSE, wopt=list(), ...)

Arguments

x

SpatVector

y

SpatRaster or character if x is a SpatRaster; character if x is a SpatVector. The character should be a PROJ.4 description of the output coordinate reference system (crs). If x is a SpatVector, you can also provide another object from which a crs can be extracted with crs

method

character. Method used for estimating the new cell values. One of:

near: nearest neighbour. This method is fast, and it can be the preferred method if the cell values represent classes. It is not a good choice for continuous values. bilinear: bilinear interpolation. Default. cubic: cubic interpolation. cubicspline: cubic spline interpolation. lanczos: Lanczos windowed sinc resampling. average: the weighted average of all contributing cells excluding NA. med: the median value of of all contributing cells excluding NA. mode: the most frequent value. max: the maximum value of all contributing cells excluding NA. min: the minimum value of all contributing cells excluding NA. q1: the first quartile of all contributing cells excluding NA. q3: the third quartile of all contributing cells excluding NA.

filename

character. Output filename. Optional

overwrite

logical. If TRUE, filename is overwritten

wopt

list. Options for writing files as in writeRaster

...

additional arguments. None implemented

Value

SpatVector or SpatRaster

See Also

resample

Examples

# NOT RUN {
# SpatRaster
a <- rast(ncol=40, nrow=40, xmin=-110, xmax=-90, ymin=40, ymax=60, 
          crs="+proj=longlat +datum=WGS84")
values(a) = 1:ncell(a)
newcrs="+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100 +datum=WGS84"
b <- rast(ncol=94, nrow=124, xmin=-944881, xmax=935118, ymin=4664377, ymax=7144377, crs=newcrs)
w <- project(a, b)

# y can also be a crs
if  (gdal_version() >= "3.0.0") {
	ww <- project(a, newcrs)
}

# SpatVector
f <- system.file("ex/lux.shp", package="terra")
v <- vect(f)
crs <- "+proj=moll +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84"
p <- project(v, crs)
p
# }