Weighted mean of rasters — weighted.mean • raster

Computes the weighted mean for each cell of a number or raster layers. The weights can be spatially variable or not.

# S4 method for RasterStackBrick,vector
weighted.mean(x, w, na.rm=FALSE, filename='', ...)

# S4 method for RasterStackBrick,RasterStackBrick
weighted.mean(x, w, na.rm=FALSE,filename='', ...)

Arguments

x: RasterStack or RasterBrick
w: A vector of weights (one number for each layer), or for spatially variable weights, a RasterStack or RasterBrick with weights (should have the same extent, resolution and number of layers as x)
na.rm: Logical. Should missing values be removed?
filename: Character. Output filename (optional)
...: Additional arguments as for writeRaster

Value

RasterLayer

See also

Summary-methods, weighted.mean

Examples

b <- brick(system.file("external/rlogo.grd", package="raster"))

# give least weight to first layer, most to last layer
wm1 <- weighted.mean(b, w=1:3)

# spatially varying weights
# weigh by column number
w1 <- init(b, v='col')

# weigh by row number
w2 <- init(b, v='row')
w <- stack(w1, w2, w2)

wm2 <- weighted.mean(b, w=w)