R: Merge SpatRaster or SpatExtent objects, or a SpatVector with...

merge {terra} | R Documentation

## Merge SpatRaster or SpatExtent objects, or a SpatVector with a data.frame

### Description

Merge SpatRasters to form a new SpatRaster object with a larger spatial extent. If objects overlap, the values get priority in the same order as the arguments. See `classify` to merge a `SpatRaster` and a `data.frame`. You can also merge SpatExtent objects.

There is a also a method for merging SpatVector with a data.frame; that is, to join the data.frame to the attribute table of the SpatVector..

### Usage

    ## S4 method for signature 'SpatRaster,SpatRaster'
    merge(x, y, ..., filename="", overwrite=FALSE, wopt=list())

    ## S4 method for signature 'SpatExtent,SpatExtent'
    merge(x, y, ...)

    ## S4 method for signature 'SpatVector,data.frame'
    merge(x, y, ...)

### Arguments

`x` | SpatRaster or SpatExtent
`y` | object of same class as `x`
`...` | if `x` is a SpatRaster or SpatVector: additional objects of the same class as `x`. If `x` is a SpatVector, the same arguments as in `merge`
`filename` | character. Output filename
`overwrite` | logical. If `TRUE`, `filename` is overwritten
`wopt` | list with named options for writing files as in `writeRaster`

### Details

The SpatRaster objects must have the same origin and spatial resolution. In areas where the SpatRaster objects overlap, the values of the SpatRaster that is last in the sequence of arguments will be retained.

### Value

SpatRaster or SpatExtent

### Note

You can use `merge` with `do.call` to merge a list of SpatRasters (see example). But note that if the list is named, these names are used by `merge`. So if all elements are named, there should be one element with a SpatRaster called `x` and another one called `y`. For example with `names(x)[1:2] <- c("x"m "y")`. You can also removed the names of the the first two elements (assuming these are SpatRasters) with `names(x)[1:2] <- ""`.

### See Also

`mosaic`

### Examples

    x <- rast(xmin=-110, xmax=-50, ymin=40, ymax=70, ncols=60, nrows=30)
    y <- rast(xmin=-80, xmax=-20, ymax=60, ymin=30)
    res(y) <- res(x)
    values(x) <- 1:ncell(x)
    values(y) <- 1:ncell(y)
    mr <- merge(x, y)
    plot(mr)
    mr <- merge(y, x)

    # if you have many SpatRaster objects in a list
    # you can use do.call:
    s <- list(x, y)
    # add arguments such as filename
    s$filename <- ""
    m <- do.call(merge, s)

    ##
    # SpatVector with data.frame
    f <- system.file("ex/lux.shp", package="terra")
    p <- vect(f)
    dfr <- data.frame(District=p$NAME_1, Canton=p$NAME_2, Value=round(runif(length(p), 100, 1000)))
    dfr <- dfr[1:5, ]
    pm <- merge(p, dfr, all.x=TRUE, by.x=c('NAME_1', 'NAME_2'), by.y=c('District', 'Canton'))
    pm
    values(pm)

[Package terra version 1.2-10 Index]