gdalcubes (version 0.2.5)

reduce: Reduce all bands of a data cube over the time dimension with a single reducer function

Description

Create a proxy data cube, which applies a single reducer function over per-band pixel time series of a data cube

Usage

reduce(cube, reducer = c("mean", "median", "min", "max"))

Arguments

cube

source data cube

reducer

reducer function, currently "min", "max", "median", "mean", "count", "sd", "var", or "sum"

Value

proxy data cube object

Examples

# NOT RUN {
# create image collection from example Landsat data only 
# if not already done in other examples
if (!file.exists(file.path(tempdir(), "L8.db"))) {
  L8_files <- list.files(system.file("L8NY18", package = "gdalcubes"),
                         ".TIF", recursive = TRUE, full.names = TRUE)
  create_image_collection(L8_files, "L8_L1TP", file.path(tempdir(), "L8.db")) 
}

L8.col = image_collection(file.path(tempdir(), "L8.db"))
v = cube_view(extent=list(left=388941.2, right=766552.4, 
              bottom=4345299, top=4744931, t0="2018-01", t1="2018-12"),
              srs="EPSG:32618", nx = 497, ny=526, dt="P1M")
L8.cube = raster_cube(L8.col, v) 
L8.rgb = select_bands(L8.cube, c("B02", "B03", "B04"))
L8.rgb.median = reduce(L8.rgb, "median")  
L8.rgb.median
# }