Create a proxy data cube, which selects specific time slices of a data cube. The time dimension of the resulting cube will be irregular / labeled.
Usage
select_time(cube, t)
Arguments
Argument
Description
cube
source data cube
t
character vector with date/time
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-04", t1="2018-07"),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 =select_time(L8.rgb, c("2018-04", "2018-07"))L8.rgb
A data cube proxy object
Dimensions:
low high count pixel_size chunk_size values
t 2018-04-01 2018-07-01 2 P0D 1 2018-04-01,2018-07-01
y 4345299 4744931 526 759.756653992395 384
x 388941.2 766552.4 497 759.781086519115 384
Bands:
name offset scale nodata unit
1 B02 0 1 NaN
2 B03 0 1 NaN
3 B04 0 1 NaN
plot(L8.rgb, rgb=3:1)
Source Code
# select_timeSelect time slices of a data cube```{r include=FALSE}library(gdalcubes)```## DescriptionCreate a proxy data cube, which selects specific time slices of a data cube. The time dimension of the resulting cubewill be irregular / labeled.## Usage```rselect_time(cube, t)```## Arguments| Argument | Description ||:------------|:----------------------------------|| cube | source data cube || t | character vector with date/time |## 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-04", t1="2018-07"),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 =select_time(L8.rgb, c("2018-04", "2018-07"))L8.rgbplot(L8.rgb, rgb=3:1)```
