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='', ...)
RasterStack or RasterBrick
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)
Logical. Should missing values be removed?
Character. Output filename (optional)
Additional arguments as for
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)