R: Cut out a geographic subset

crop {terra} | R Documentation

## 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 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

    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")

[Package terra version 1.2-10 Index]