
R version 3.5.3 (2019-03-11) -- "Great Truth"
Copyright (C) 2019 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

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.

> library(spacetime)
> suppressPackageStartupMessages(library(stars))
> data(air) # this loads several datasets in .GlobalEnv
> dim(air)
space  time 
   70  4383 
> 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
bbox:           xmin: 6.28107 ymin: 47.80847 xmax: 14.78617 ymax: 54.92497
epsg (SRID):    4326
proj4string:    +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)
+ }
> 
> proc.time()
   user  system elapsed 
  1.770   0.082   1.843 
