crop: Cut out a geographic subset

Description Usage Arguments Value See Also Examples

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

1
2
3
4
5
6
7
8
## S4 method for signature 'SpatRaster'
crop(x, y, snap="near", filename="", ...)

## S4 method for signature 'SpatVector,ANY'
crop(x, y)

## S4 method for signature '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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
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")

Example output

terra version 0.9.8

terra documentation built on Aug. 20, 2021, 9:07 a.m.