
R Under development (unstable) (2026-02-12 r89409 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2026 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> 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)
+ }
Geometry set for 70 features 
Geometry type: POINT
Dimension:     XY
Bounding box:  xmin: 6.28107 ymin: 47.80847 xmax: 14.78617 ymax: 54.92497
Geodetic CRS:  +proj=longlat +datum=WGS84 +no_defs
First 5 geometries:
POINT (9.585911 53.67057)
POINT (9.68503 53.52418)
POINT (9.791584 54.07312)
POINT (13.64701 52.44775)
POINT (13.29635 52.65315)
> 
> 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)
Warning message:
In plot.stars(output) :
  breaks="quantile" leads to a single class; maybe try breaks="equal" instead?
> 
> proc.time()
   user  system elapsed 
   2.73    0.28    2.98 
