Examples for 'raster::stackApply'


Apply a function on subsets of a RasterStack or RasterBrick

Aliases: stackApply

Keywords: methods spatial

### ** Examples

r <- raster(ncol=10, nrow=10)
values(r) <- 1:ncell(r)
s <- brick(r,r,r,r,r,r)
s <- s * 1:6
b1 <- stackApply(s, indices=c(1,1,1,2,2,2), fun=sum)
b1
class      : RasterBrick 
dimensions : 10, 10, 100, 2  (nrow, ncol, ncell, nlayers)
resolution : 36, 18  (x, y)
extent     : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 
source     : memory
names      : index_1, index_2 
min values :       6,      15 
max values :     600,    1500 
b2 <- stackApply(s, indices=c(1,2,3,1,2,3), fun=sum)
b2
class      : RasterBrick 
dimensions : 10, 10, 100, 3  (nrow, ncol, ncell, nlayers)
resolution : 36, 18  (x, y)
extent     : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 
source     : memory
names      : index_1, index_2, index_3 
min values :       5,       7,       9 
max values :     500,     700,     900 

[Package raster version 3.5-15 Index]