Create a proxy data cube, which joins the bands of two identically shaped data cubes. The resulting cube will have bands from both input cubes.
Usage
join_bands(cube_list, cube_names =NULL)
Arguments
Argument
Description
cube_list
a list with two or more source data cubes
cube_names
list or character vector with optional name prefixes for bands in the output data cube (see Details)
Details
The number of provided cube_names must match the number of provided input cubes. If no cube_names are provided, bands of the output cube will adopt original names from the input cubes (without any prefix). If any two of the input bands have identical names, prefixes default prefixes (“X1”, “X2”, …) will be used.
Value
proxy data cube object
Note
This function returns a proxy object, i.e., it will not start any computations besides deriving the shape of the result.
Examples
# create image collection from example Landsat data only # if not already done in other examplesif (!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"), quiet =TRUE) }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-05"),srs="EPSG:32618", nx =497, ny=526, dt="P1M")L8.cube =raster_cube(L8.col, v)L8.cube.b04 =select_bands(raster_cube(L8.col, v), c("B04"))L8.cube.b05 =select_bands(raster_cube(L8.col, v), c("B05"))join_bands(list(L8.cube.b04,L8.cube.b05))
A data cube proxy object
Dimensions:
low high count pixel_size chunk_size
t 2018-01-01 2018-05-31 5 P1M 1
y 4345299 4744931 526 759.756653992395 384
x 388941.2 766552.4 497 759.781086519115 384
Bands:
name offset scale nodata unit
1 X1.B04 0 1 NaN
2 X2.B05 0 1 NaN
plot(join_bands(list(L8.cube.b04,L8.cube.b05)))
Source Code
# join_bandsJoin bands of two identically shaped data cubes```{r include=FALSE}library(gdalcubes)```## DescriptionCreate a proxy data cube, which joins the bands of two identically shaped data cubes. The resulting cubewill have bands from both input cubes.## Usage```rjoin_bands(cube_list, cube_names =NULL)```## Arguments| Argument | Description ||:------------|:----------------------------------|| cube_list | a list with two or more source data cubes || cube_names | list or character vector with optional name prefixes for bands in the output data cube (see Details) |## DetailsThe number of provided cube_names must match the number of provided input cubes.If no cube_names are provided, bands of the output cube will adopt original names from the input cubes (without any prefix). If any two of the input bands have identical names,prefixes default prefixes ("X1", "X2", ...) will be used.## Valueproxy data cube object## NoteThis function returns a proxy object, i.e., it will not start any computations besides deriving the shape of the result.## Examples```{r}# create image collection from example Landsat data only # if not already done in other examplesif (!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"), quiet =TRUE) }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-05"),srs="EPSG:32618", nx =497, ny=526, dt="P1M")L8.cube =raster_cube(L8.col, v)L8.cube.b04 =select_bands(raster_cube(L8.col, v), c("B04"))L8.cube.b05 =select_bands(raster_cube(L8.col, v), c("B05"))join_bands(list(L8.cube.b04,L8.cube.b05))plot(join_bands(list(L8.cube.b04,L8.cube.b05)))```
