tests/sf.R

In stars: Spatiotemporal Arrays, Raster and Vector Data Cubes

suppressPackageStartupMessages(library(stars))
library(sf)
#jp2 = system.file("jp2/B01.jp2", package = "stars")
#(x = read_stars(jp2, options = c("OVERVIEW_LEVEL=3")))
tif = system.file("tif/L7_ETMs.tif", package = "stars")
#(x = read_stars(tif, options = c("OVERVIEW_LEVEL=3"))) # FIXME: breaks on windows
(x = read_stars(tif)[,1:20,1:20])
# names(x) = "L7_ETM"

# library(abind)
# x = adrop(x)
image(x)
(sfc = st_as_sfc(x, as_points = FALSE))
plot(sfc, add  =TRUE)
(sfc = st_as_sfc(x, as_points = TRUE))
plot(sfc, add = TRUE)

sf = st_as_sf(x, as_points = FALSE)
plot(sf, border = NA)
sf = st_as_sf(x[,,,1], as_points = FALSE, merge = TRUE)
plot(sf, border = NA)

sfc1 <- st_as_sfc(x, as_points = TRUE, use_cpp = TRUE, na.rm = FALSE)
sfc2 <- st_as_sfc(x, as_points = TRUE, use_cpp = FALSE)
identical(sfc1, sfc2)

sfc1 <- st_as_sfc(x, as_points = FALSE, use_cpp = TRUE, na.rm = FALSE)
sfc2 <- st_as_sfc(x, as_points = FALSE, use_cpp = FALSE)
identical(sfc1, sfc2)

# sf -> stars -> sf
x = st_sfc(st_point(0:1), st_point(1:2), st_point(2:3))
m = matrix(1:9,3)
colnames(m) = c("a", "b", "c")
foo = st_sf(m, geom = x)
x = st_as_stars(foo)
c(x)
st_redimension(x)
c(x, x)
(y = st_redimension(x, along = list(time = as.Date("2017-11-27") + 0:2)))
st_as_sf(y)
st_as_sf(st_as_stars(foo))

c(x,x,along=1) # NEEDS FIXING!
# st_as_stars.sfc
st_as_stars(st_geometry(st_as_sf(y)), nx = 10, ny = 10, as_points = FALSE)