terra (version 1.1-17)

crop: Cut out a geographic subset

Description

Cut out a part of a SpatRaster with a SpatExtent, or another object from which an extent can be obtained. With a SpatRaster you can only extract rectangular areas, but see mask for setting cell values within SpatRaster to NA.

You can crop a SpatVector with a rectangle, or with another vector (if these are not polygons, the minimum convex hull is used). Unlike with intersect the geometries and attributes of y are not trasferred to the output.

Usage

# S4 method for SpatRaster
crop(x, y, snap="near", filename="", ...)
# S4 method for SpatVector,ANY
crop(x, y)
# S4 method for SpatVector,SpatVector
crop(x, y)

Arguments

x

SpatRaster or SpatVector

y

SpatExtent or other object that has a SpatExtent (ext returns a SpatExtent), or a SpatVector to crop another SpatVector

snap

character. One of "near", "in", or "out"

filename

character. Output filename

...

additional arguments for writing files as in writeRaster

Value

SpatRaster

See Also

intersect

Examples

# NOT RUN {
r <- rast(xmin=0, xmax=10, ymin=0, ymax=10, nrows=25, ncols=25)
values(r) <- 1:ncell(r)
e <- ext(-5, 5, -5, 5)
rc <- crop(r, e)

# vector
f <- system.file("ex/lux.shp", package="terra")
v <- vect(f)
e <- ext(6.15, 6.3, 49.7, 49.8)
x <- crop(v, e)
#plot(x, "NAME_1")
# }