resample: Transfer values of a SpatRaster to another one with a...

In rspatial/terra: Spatial Data Analysis

Transfer values of a SpatRaster to another one with a different geometry

Description

resample transfers values between SpatRaster objects that do not align (have a different origin and/or resolution). See project to change the coordinate reference system (crs).

If the origin and extent of the input and output are the same, you should consider using these other functions instead: aggregate, disagg, extend or crop.

Usage

## S4 method for signature 'SpatRaster,SpatRaster'
resample(x, y, method, threads=FALSE, by_util=FALSE, filename="", ...)

Arguments

x	SpatRaster to be resampled
y	SpatRaster with the geometry that x should be resampled to
method	character. Method used for estimating the new cell values. One of: bilinear: bilinear interpolation (3x3 cell window). This is used by default if the first layer of x is not categorical average: This can be a good choice with continuous variables if the output cells overlap with multiple input cells. near: nearest neighbor. This is used by default if the first layer of x is categorical. This method is not a good choice for continuous values. mode: The modal value. This can be a good choice for categrical rasters, if the output cells overlap with multiple input cells. cubic: cubic interpolation (5x5 cell window). cubicspline: cubic B-spline interpolation. (5x5 cell window). lanczos: Lanczos windowed sinc resampling. (7x7 cell window). sum: the weighted sum of all non-NA contributing grid cells. min, q1, median, q3, max: the minimum, first quartile, median, third quartile, or maximum value. rms: the root-mean-square value of all non-NA contributing grid cells.
threads	logical. If TRUE multiple threads are used (faster for large files)
by_util	logical. If TRUE the GDAL warp utility is used
filename	character. Output filename
...	additional arguments for writing files as in writeRaster

Value

SpatRaster

See Also

aggregate, disagg, crop, project

Examples

r <- rast(nrows=3, ncols=3, xmin=0, xmax=10, ymin=0, ymax=10)
values(r) <- 1:ncell(r)
s <- rast(nrows=25, ncols=30, xmin=1, xmax=11, ymin=-1, ymax=11)
x <- resample(r, s, method="bilinear")

opar <- par(no.readonly =TRUE)
par(mfrow=c(1,2))
plot(r)
plot(x)
par(opar)

rspatial/terra documentation built on Feb. 9, 2025, 4:06 p.m.