mirror of
https://github.com/appelmar/gdalcubes.git
synced 2025-02-23 07:54:15 +01:00
60 lines
1.6 KiB
R
60 lines
1.6 KiB
R
library(gdalcubes)
|
|
v = cube_view(srs = "EPSG:4326", extent = list(left = 5, right = 10, bottom = 48, top = 53,
|
|
t0 = "2021-01-01", t1 = "2021-12-31"), dt = "P1M",
|
|
dx = 1, dy = 1)
|
|
v
|
|
|
|
|
|
gdalcubes:::.raster_cube_dummy(v, 1, 1.0) |>
|
|
aggregate_space(dx = 2, dy = 2, method = "sum") |>
|
|
as_array() -> x
|
|
|
|
expect_true(all(x[1,1,,] == matrix(c(1,2,2,2,4,4,2,4,4), nrow = 3, byrow = TRUE)))
|
|
expect_true(all(x[1,12,,] == matrix(c(1,2,2,2,4,4,2,4,4), nrow = 3, byrow = TRUE)))
|
|
|
|
|
|
|
|
gdalcubes:::.raster_cube_dummy(v, 1, 1.0, chunking = c(1,4,4)) |>
|
|
aggregate_space(dx = 2, dy = 2, method = "count") |>
|
|
as_array() -> y
|
|
expect_equal(x,y)
|
|
|
|
|
|
gdalcubes:::.raster_cube_dummy(v, 1, 1.0, chunking = c(1,2,10)) |>
|
|
aggregate_space(dx = 2, dy = 2, method = "count") |>
|
|
as_array() -> z
|
|
expect_equal(x,z)
|
|
|
|
|
|
gdalcubes:::.raster_cube_dummy(v, 3, 2.0, chunking = c(1,2,10)) |>
|
|
aggregate_space(dx = 3, dy = 3, method = "mean") |>
|
|
as_array() -> z
|
|
expect_true(all(z == 2))
|
|
|
|
|
|
gdalcubes:::.raster_cube_dummy(v, 3, 2.0, chunking = c(3,4,4)) |>
|
|
aggregate_space(dx = 3, dy = 3, method = "mean") |>
|
|
as_array() -> z
|
|
expect_true(all(z == 2))
|
|
|
|
|
|
gdalcubes:::.raster_cube_dummy(v, 3, 2.0, chunking = c(3,4,4)) |>
|
|
aggregate_space(dx = 3, dy = 3, method = "var") |>
|
|
as_array() -> z
|
|
expect_true(all(z == 0))
|
|
|
|
|
|
gdalcubes:::.raster_cube_dummy(v, 3, 2.0, chunking = c(1,3,3)) |>
|
|
aggregate_space(dx = 3, dy = 3, method = "sd") |>
|
|
as_array() -> z
|
|
expect_true(all(z == 0))
|
|
|
|
|
|
gdalcubes:::.raster_cube_dummy(v, 2, 1.0, chunking = c(1,3,3)) |>
|
|
aggregate_space(dx = 3, dy = 3, method = "prod") |>
|
|
as_array() -> z
|
|
expect_true(all(z == 1))
|
|
|
|
|
|
|
|
|