tests/align.R

In stars: Spatiotemporal Arrays, Raster and Vector Data Cubes

suppressPackageStartupMessages(library(stars))
set.seed(1331)
st_sample(st_as_sfc(st_bbox(L7_ETMs)), 1)
b = st_buffer(st_sample(st_as_sfc(st_bbox(L7_ETMs)), 1), 335)
bba = stars:::st_align(st_bbox(b), st_dimensions(L7_ETMs))

plot(L7_ETMs[,,,1], extent = st_buffer(b, 20), reset = FALSE, key.pos = NULL)
plot(b, add = TRUE, col = NA, border = 'red')
plot(st_as_sfc(st_bbox(b)), col = NA, border = 'red', add = TRUE)
plot(st_as_sfc(bba), col = NA, border = 'blue', add = TRUE)
st_bbox(b)
bba

plot(L7_ETMs[,,,1], reset = FALSE, key.pos = NULL)
plot(b, col = 'red', add = TRUE)
plot(st_as_sfc(bba), col = NA, border = 'blue', add = TRUE)

b0 = st_buffer(b, -100)
L7 = st_as_stars(L7_ETMs)[b0][,,,1] |> st_normalize()
bbb = stars:::st_align(st_bbox(b), st_dimensions(L7))
plot(L7, extent = st_buffer(b, 120), reset = FALSE, key.pos = NULL)
plot(st_as_sfc(st_bbox(L7)), col = NA, border = 'black', add = TRUE)
plot(st_as_sfc(bbb), col = NA, border = 'blue', add = TRUE)
plot(b, col = NA, border = 'red', add = TRUE)

all.equal(bbb, bba)