suppressPackageStartupMessages(library(stars))
if (suppressPackageStartupMessages(require(sp, quietly = TRUE)) && require(spacetime, quietly = TRUE)) {
data(air) # this loads several datasets in .GlobalEnv
dim(air)
d = st_dimensions(station = st_as_sfc(stations), time = dates)
aq = st_as_stars(list(PM10 = air), dimensions = d)
image(aperm(log(aq), 2:1), main = "NA pattern (white) in PM10 station time series")
st_geometry(aq)
}
csv = system.file("tif/ESACCI-LC-Legend.csv", package = "starsdata")
if (csv != "") {
rgb = read.csv(csv, header = TRUE, sep = ";")
luc = read_stars(system.file("tif/LUC3.tif", package = "starsdata"))
plot(luc, rgb = rgb, axes = TRUE, key.pos = NULL)
}
# https://github.com/r-spatial/stars/issues/233 :
nc = read_sf(system.file("shape/nc.shp", package = "sf"))
output = st_as_stars(st_bbox(nc), nx = 1, ny = 1, values = 1)
plot(output)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.