Package check result: OK Changes to worse in reverse depends: Package: gstat Check: tests Old result: NOTE Running ‘allier.R’ [1s/1s] Comparing ‘allier.Rout’ to ‘allier.Rout.save’ ... OK Running ‘blockkr.R’ [1s/1s] Comparing ‘blockkr.Rout’ to ‘blockkr.Rout.save’ ... OK Running ‘covtable.R’ [1s/1s] Comparing ‘covtable.Rout’ to ‘covtable.Rout.save’ ... OK Running ‘cv.R’ [1s/1s] Comparing ‘cv.Rout’ to ‘cv.Rout.save’ ...68c68 < Min. :4.644 Min. :0.002371 Min. :4.727 Min. :-0.4190976 --- > Min. :4.644 Min. :0.002371 Min. :4.727 Min. :-0.4190977 Running ‘cv3d.R’ [1s/1s] Comparing ‘cv3d.Rout’ to ‘cv3d.Rout.save’ ...14c14 < 3rd Qu.: 0.161 3rd Qu.:0.240 3rd Qu.: 0.501 3rd Qu.: 0.8508 --- > 3rd Qu.: 0.161 3rd Qu.:0.240 3rd Qu.: 0.500 3rd Qu.: 0.8508 17,22c17,22 < Min. :-5.1678 Min. : 1.0 Min. :0.00678 Min. :0.0034 < 1st Qu.:-1.8749 1st Qu.:13.2 1st Qu.:0.23966 1st Qu.:0.2466 < Median : 0.2453 Median :25.5 Median :0.48668 Median :0.4525 < Mean : 0.0167 Mean :25.5 Mean :0.49966 Mean :0.4969 < 3rd Qu.: 2.0201 3rd Qu.:37.8 3rd Qu.:0.74730 3rd Qu.:0.7394 < Max. : 7.3541 Max. :50.0 Max. :0.98754 Max. :0.9872 --- > Min. :-5.168 Min. : 1.0 Min. :0.0068 Min. :0.0034 > 1st Qu.:-1.875 1st Qu.:13.2 1st Qu.:0.2397 1st Qu.:0.2466 > Median : 0.245 Median :25.5 Median :0.4867 Median :0.4525 > Mean : 0.017 Mean :25.5 Mean :0.4997 Mean :0.4969 > 3rd Qu.: 2.020 3rd Qu.:37.8 3rd Qu.:0.7473 3rd Qu.:0.7394 > Max. : 7.354 Max. :50.0 Max. :0.9875 Max. :0.9872 24,29c24,29 < Min. :0.00164 < 1st Qu.:0.18646 < Median :0.44850 < Mean :0.47142 < 3rd Qu.:0.72403 < Max. :0.99420 --- > Min. :0.0016 > 1st Qu.:0.1865 > Median :0.4485 > Mean :0.4714 > 3rd Qu.:0.7240 > Max. :0.9942 Running ‘fit.R’ [1s/1s] Comparing ‘fit.Rout’ to ‘fit.Rout.save’ ... OK Running ‘krige0.R’ [2s/2s] Comparing ‘krige0.Rout’ to ‘krige0.Rout.save’ ... OK Running ‘line.R’ [1s/2s] Comparing ‘line.Rout’ to ‘line.Rout.save’ ... OK Running ‘merge.R’ [1s/1s] Comparing ‘merge.Rout’ to ‘merge.Rout.save’ ... OK Running ‘na.action.R’ [1s/1s] Comparing ‘na.action.Rout’ to ‘na.action.Rout.save’ ... OK Running ‘rings.R’ [1s/2s] Comparing ‘rings.Rout’ to ‘rings.Rout.save’ ... OK Running ‘sim.R’ [1s/2s] Comparing ‘sim.Rout’ to ‘sim.Rout.save’ ... OK Running ‘stars.R’ [9s/10s] Comparing ‘stars.Rout’ to ‘stars.Rout.save’ ...196,201c196,201 < Min. :4.777 Min. :0.0855 (0,5]: 316 < 1st Qu.:5.238 1st Qu.:0.1373 (5,6]:1778 < Median :5.573 Median :0.1622 (6,7]: 962 < Mean :5.707 Mean :0.1853 (7,8]: 47 < 3rd Qu.:6.172 3rd Qu.:0.2116 (8,9]: 0 < Max. :7.440 Max. :0.5003 NA's :5009 --- > Min. :4.777 Min. :0.085 (0,5]: 316 > 1st Qu.:5.238 1st Qu.:0.137 (5,6]:1778 > Median :5.573 Median :0.162 (6,7]: 962 > Mean :5.707 Mean :0.185 (7,8]: 47 > 3rd Qu.:6.172 3rd Qu.:0.212 (8,9]: 0 > Max. :7.440 Max. :0.500 NA's :5009 Running ‘unproj.R’ [3s/3s] Comparing ‘unproj.Rout’ to ‘unproj.Rout.save’ ... OK Running ‘variogram.R’ [1s/1s] Comparing ‘variogram.Rout’ to ‘variogram.Rout.save’ ... OK Running ‘vdist.R’ [1s/1s] Comparing ‘vdist.Rout’ to ‘vdist.Rout.save’ ... OK Running ‘windst.R’ [6s/6s] Comparing ‘windst.Rout’ to ‘windst.Rout.save’ ... OK New result: ERROR Running ‘allier.R’ [1s/1s] Comparing ‘allier.Rout’ to ‘allier.Rout.save’ ... OK Running ‘blockkr.R’ [1s/1s] Comparing ‘blockkr.Rout’ to ‘blockkr.Rout.save’ ... OK Running ‘covtable.R’ [1s/1s] Comparing ‘covtable.Rout’ to ‘covtable.Rout.save’ ... OK Running ‘cv.R’ [2s/2s] Comparing ‘cv.Rout’ to ‘cv.Rout.save’ ...68c68 < Min. :4.644 Min. :0.002371 Min. :4.727 Min. :-0.4190976 --- > Min. :4.644 Min. :0.002371 Min. :4.727 Min. :-0.4190977 Running ‘cv3d.R’ [2s/2s] Comparing ‘cv3d.Rout’ to ‘cv3d.Rout.save’ ...14c14 < 3rd Qu.: 0.161 3rd Qu.:0.240 3rd Qu.: 0.501 3rd Qu.: 0.8508 --- > 3rd Qu.: 0.161 3rd Qu.:0.240 3rd Qu.: 0.500 3rd Qu.: 0.8508 17,22c17,22 < Min. :-5.1678 Min. : 1.0 Min. :0.00678 Min. :0.0034 < 1st Qu.:-1.8749 1st Qu.:13.2 1st Qu.:0.23966 1st Qu.:0.2466 < Median : 0.2453 Median :25.5 Median :0.48668 Median :0.4525 < Mean : 0.0167 Mean :25.5 Mean :0.49966 Mean :0.4969 < 3rd Qu.: 2.0201 3rd Qu.:37.8 3rd Qu.:0.74730 3rd Qu.:0.7394 < Max. : 7.3541 Max. :50.0 Max. :0.98754 Max. :0.9872 --- > Min. :-5.168 Min. : 1.0 Min. :0.0068 Min. :0.0034 > 1st Qu.:-1.875 1st Qu.:13.2 1st Qu.:0.2397 1st Qu.:0.2466 > Median : 0.245 Median :25.5 Median :0.4867 Median :0.4525 > Mean : 0.017 Mean :25.5 Mean :0.4997 Mean :0.4969 > 3rd Qu.: 2.020 3rd Qu.:37.8 3rd Qu.:0.7473 3rd Qu.:0.7394 > Max. : 7.354 Max. :50.0 Max. :0.9875 Max. :0.9872 24,29c24,29 < Min. :0.00164 < 1st Qu.:0.18646 < Median :0.44850 < Mean :0.47142 < 3rd Qu.:0.72403 < Max. :0.99420 --- > Min. :0.0016 > 1st Qu.:0.1865 > Median :0.4485 > Mean :0.4714 > 3rd Qu.:0.7240 > Max. :0.9942 Running ‘fit.R’ [1s/1s] Comparing ‘fit.Rout’ to ‘fit.Rout.save’ ... OK Running ‘krige0.R’ [3s/3s] Comparing ‘krige0.Rout’ to ‘krige0.Rout.save’ ... OK Running ‘line.R’ [2s/2s] Comparing ‘line.Rout’ to ‘line.Rout.save’ ... OK Running ‘merge.R’ [1s/1s] Comparing ‘merge.Rout’ to ‘merge.Rout.save’ ... OK Running ‘na.action.R’ [1s/1s] Comparing ‘na.action.Rout’ to ‘na.action.Rout.save’ ... OK Running ‘rings.R’ [1s/1s] Comparing ‘rings.Rout’ to ‘rings.Rout.save’ ... OK Running ‘sim.R’ [1s/1s] Comparing ‘sim.Rout’ to ‘sim.Rout.save’ ... OK Running ‘stars.R’ [11s/11s] Running ‘unproj.R’ [4s/4s] Comparing ‘unproj.Rout’ to ‘unproj.Rout.save’ ... OK Running ‘variogram.R’ [1s/1s] Comparing ‘variogram.Rout’ to ‘variogram.Rout.save’ ... OK Running ‘vdist.R’ [1s/1s] Comparing ‘vdist.Rout’ to ‘vdist.Rout.save’ ... OK Running ‘windst.R’ [4s/4s] Running the tests in ‘tests/stars.R’ failed. Complete output: > Sys.setenv(TZ = "UTC") > > # 0. using sp: > > suppressPackageStartupMessages(library(sp)) > demo(meuse, ask = FALSE) demo(meuse) ---- ~~~~~ > require(sp) > crs = CRS("EPSG:28992") > data("meuse") > coordinates(meuse) <- ~x+y > proj4string(meuse) <- crs > data("meuse.grid") > coordinates(meuse.grid) <- ~x+y > gridded(meuse.grid) <- TRUE > proj4string(meuse.grid) <- crs > data("meuse.riv") > meuse.riv <- SpatialPolygons(list(Polygons(list(Polygon(meuse.riv)),"meuse.riv"))) > proj4string(meuse.riv) <- crs > data("meuse.area") > meuse.area = SpatialPolygons(list(Polygons(list(Polygon(meuse.area)), "area"))) > proj4string(meuse.area) <- crs > suppressPackageStartupMessages(library(gstat)) > v = variogram(log(zinc)~1, meuse) > (v.fit = fit.variogram(v, vgm(1, "Sph", 900, 1))) model psill range 1 Nug 0.05066243 0.0000 2 Sph 0.59060780 897.0209 > k_sp = krige(log(zinc)~1, meuse[-(1:5),], meuse[1:5,], v.fit) [using ordinary kriging] > k_sp_grd = krige(log(zinc)~1, meuse, meuse.grid, v.fit) [using ordinary kriging] > > # 1. using sf: > suppressPackageStartupMessages(library(sf)) > demo(meuse_sf, ask = FALSE, echo = FALSE) > # reloads meuse as data.frame, so > demo(meuse, ask = FALSE) demo(meuse) ---- ~~~~~ > require(sp) > crs = CRS("EPSG:28992") > data("meuse") > coordinates(meuse) <- ~x+y > proj4string(meuse) <- crs > data("meuse.grid") > coordinates(meuse.grid) <- ~x+y > gridded(meuse.grid) <- TRUE > proj4string(meuse.grid) <- crs > data("meuse.riv") > meuse.riv <- SpatialPolygons(list(Polygons(list(Polygon(meuse.riv)),"meuse.riv"))) > proj4string(meuse.riv) <- crs > data("meuse.area") > meuse.area = SpatialPolygons(list(Polygons(list(Polygon(meuse.area)), "area"))) > proj4string(meuse.area) <- crs > > v = variogram(log(zinc)~1, meuse_sf) > (v.fit = fit.variogram(v, vgm(1, "Sph", 900, 1))) model psill range 1 Nug 0.05066243 0.0000 2 Sph 0.59060780 897.0209 > k_sf = krige(log(zinc)~1, meuse_sf[-(1:5),], meuse_sf[1:5,], v.fit) [using ordinary kriging] > > all.equal(k_sp, as(k_sf, "Spatial"), check.attributes = FALSE) [1] TRUE > all.equal(k_sp, as(k_sf, "Spatial"), check.attributes = TRUE) [1] "Attributes: < Component \"bbox\": Attributes: < Component \"dimnames\": Component 1: 2 string mismatches > >" [2] "Attributes: < Component \"coords\": Attributes: < Component \"dimnames\": Component 2: 2 string mismatches > >" [3] "Attributes: < Component \"coords.nrs\": Numeric: lengths (2, 0) differ >" > > # 2. using stars for grid: > > suppressPackageStartupMessages(library(stars)) > st = st_as_stars(meuse.grid) > st_crs(st) Coordinate Reference System: User input: Amersfoort / RD New wkt: PROJCRS["Amersfoort / RD New", BASEGEOGCRS["Amersfoort", DATUM["Amersfoort", ELLIPSOID["Bessel 1841",6377397.155,299.1528128, LENGTHUNIT["metre",1]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], ID["EPSG",4289]], CONVERSION["RD New", METHOD["Oblique Stereographic", ID["EPSG",9809]], PARAMETER["Latitude of natural origin",52.1561605555556, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8801]], PARAMETER["Longitude of natural origin",5.38763888888889, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8802]], PARAMETER["Scale factor at natural origin",0.9999079, SCALEUNIT["unity",1], ID["EPSG",8805]], PARAMETER["False easting",155000, LENGTHUNIT["metre",1], ID["EPSG",8806]], PARAMETER["False northing",463000, LENGTHUNIT["metre",1], ID["EPSG",8807]]], CS[Cartesian,2], AXIS["easting (X)",east, ORDER[1], LENGTHUNIT["metre",1]], AXIS["northing (Y)",north, ORDER[2], LENGTHUNIT["metre",1]], USAGE[ SCOPE["Engineering survey, topographic mapping."], AREA["Netherlands - onshore, including Waddenzee, Dutch Wadden Islands and 12-mile offshore coastal zone."], BBOX[50.75,3.2,53.7,7.22]], ID["EPSG",28992]] > > # compare inputs: > sp = as(st, "Spatial") > fullgrid(meuse.grid) = TRUE > all.equal(sp, meuse.grid["dist"], check.attributes = FALSE) [1] "Names: Lengths (5, 1) differ (string compare on first 1)" [2] "Names: 1 string mismatch" > all.equal(sp, meuse.grid["dist"], check.attributes = TRUE, use.names = FALSE) [1] "Names: Lengths (5, 1) differ (string compare on first 1)" [2] "Names: 1 string mismatch" [3] "Attributes: < Component 3: Names: 1 string mismatch >" [4] "Attributes: < Component 3: Length mismatch: comparison on first 1 components >" [5] "Attributes: < Component 3: Component 1: Mean relative difference: 1.08298 >" [6] "Attributes: < Component 4: Attributes: < Component 2: names for current but not for target > >" [7] "Attributes: < Component 4: Attributes: < Component 3: names for current but not for target > >" > > # kriging: > st_crs(st) = st_crs(meuse_sf) = NA # GDAL roundtrip messes them up! > k_st = if (Sys.getenv("USER") == "travis") { + try(krige(log(zinc)~1, meuse_sf, st, v.fit)) + } else { + krige(log(zinc)~1, meuse_sf, st, v.fit) + } [using ordinary kriging] > k_st stars object with 2 dimensions and 2 attributes attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. NA's var1.pred 4.7765547 5.2376293 5.5728839 5.7072287 6.1717619 7.4399911 5009 var1.var 0.0854949 0.1372864 0.1621838 0.1853319 0.2116152 0.5002756 5009 dimension(s): from to offset delta x/y x 1 78 178440 40 [x] y 1 104 333760 -40 [y] > > # handle factors, when going to stars? > k_sp_grd$cls = cut(k_sp_grd$var1.pred, c(0, 5, 6, 7, 8, 9)) > st_as_stars(k_sp_grd) stars object with 2 dimensions and 3 attributes attribute(s): var1.pred var1.var cls Min. :4.777 Min. :0.0855 (0,5]: 316 1st Qu.:5.238 1st Qu.:0.1373 (5,6]:1778 Median :5.573 Median :0.1622 (6,7]: 962 Mean :5.707 Mean :0.1853 (7,8]: 47 3rd Qu.:6.172 3rd Qu.:0.2116 (8,9]: 0 Max. :7.440 Max. :0.5003 NA's :5009 NA's :5009 NA's :5009 dimension(s): from to offset delta refsys x/y x 1 78 178440 40 Amersfoort / RD New [x] y 1 104 333760 -40 Amersfoort / RD New [y] > if (require(raster, quietly = TRUE)) { + print(st_as_stars(raster::stack(k_sp_grd))) # check + print(all.equal(st_redimension(st_as_stars(k_sp_grd)), st_as_stars(raster::stack(k_sp_grd)), check.attributes=FALSE)) + } stars object with 3 dimensions and 1 attribute attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. NA's var1.pred 0.0854949 0.2116778 2 2.710347 5.237542 7.439991 15027 dimension(s): from to offset delta refsys values x 1 78 178440 40 Amersfoort / RD New NULL y 1 104 333760 -40 Amersfoort / RD New NULL band 1 3 NA NA NA var1.pred, var1.var , cls x/y x [x] y [y] band [1] TRUE > > suppressPackageStartupMessages(library(spacetime)) > > tm = as.POSIXct("2019-02-25 15:37:24 CET") > n = 4 > s = stars:::st_stars(list(foo = array(1:(n^3), rep(n,3))), + stars:::create_dimensions(list( + x = stars:::create_dimension(from = 1, to = n, offset = 10, delta = 0.5), + y = stars:::create_dimension(from = 1, to = n, offset = 0, delta = -0.7), + time = stars:::create_dimension(values = tm + 1:n)), + raster = stars:::get_raster(dimensions = c("x", "y"))) + ) Error in `/.difftime`(diff(range(ud)), mean(ud)) : second argument of / cannot be a "difftime" object Calls: ... -> regular_intervals -> isTRUE -> /.difftime Execution halted Running the tests in ‘tests/windst.R’ failed. Complete output: > suppressPackageStartupMessages(library(sp)) > suppressPackageStartupMessages(library(spacetime)) > suppressPackageStartupMessages(library(gstat)) > suppressPackageStartupMessages(library(stars)) > > data(wind) > wind.loc$y = as.numeric(char2dms(as.character(wind.loc[["Latitude"]]))) > wind.loc$x = as.numeric(char2dms(as.character(wind.loc[["Longitude"]]))) > coordinates(wind.loc) = ~x+y > proj4string(wind.loc) = "+proj=longlat +datum=WGS84 +ellps=WGS84" > > wind$time = ISOdate(wind$year+1900, wind$month, wind$day) > wind$jday = as.numeric(format(wind$time, '%j')) > stations = 4:15 > windsqrt = sqrt(0.5148 * wind[stations]) # knots -> m/s > Jday = 1:366 > daymeans = colMeans( + sapply(split(windsqrt - colMeans(windsqrt), wind$jday), colMeans)) > meanwind = lowess(daymeans ~ Jday, f = 0.1)$y[wind$jday] > velocities = apply(windsqrt, 2, function(x) { x - meanwind }) > # match order of columns in wind to Code in wind.loc; > # convert to utm zone 29, to be able to do interpolation in > # proper Euclidian (projected) space: > pts = coordinates(wind.loc[match(names(wind[4:15]), wind.loc$Code),]) > pts = SpatialPoints(pts) > if (require(sp, quietly = TRUE) && require(maps, quietly = TRUE)) { + proj4string(pts) = "+proj=longlat +datum=WGS84 +ellps=WGS84" + utm29 = "+proj=utm +zone=29 +datum=WGS84 +ellps=WGS84" + pts = as(st_transform(st_as_sfc(pts), utm29), "Spatial") + # note the t() in: + w = STFDF(pts, wind$time, data.frame(values = as.vector(t(velocities)))) + + library(mapdata) + mp = map("worldHires", xlim = c(-11,-5.4), ylim = c(51,55.5), plot=FALSE) + sf = st_transform(st_as_sf(mp, fill = FALSE), utm29) + m = as(sf, "Spatial") + + # setup grid + grd = SpatialPixels(SpatialPoints(makegrid(m, n = 300)), + proj4string = m@proj4string) + # grd$t = rep(1, nrow(grd)) + #coordinates(grd) = ~x1+x2 + #gridded(grd)=TRUE + + # select april 1961: + w = w[, "1961-04"] + + covfn = function(x, y = x) { + du = spDists(coordinates(x), coordinates(y)) + t1 = as.numeric(index(x)) # time in seconds + t2 = as.numeric(index(y)) # time in seconds + dt = abs(outer(t1, t2, "-")) + # separable, product covariance model: + 0.6 * exp(-du/750000) * exp(-dt / (1.5 * 3600 * 24)) + } + + n = 10 + tgrd = seq(min(index(w)), max(index(w)), length=n) + pred = krige0(sqrt(values)~1, w, STF(grd, tgrd), covfn) + layout = list(list("sp.points", pts, first=F, cex=.5), + list("sp.lines", m, col='grey')) + wind.pr0 = STFDF(grd, tgrd, data.frame(var1.pred = pred)) + + v = vgmST("separable", + space = vgm(1, "Exp", 750000), + time = vgm(1, "Exp", 1.5 * 3600 * 24), + sill = 0.6) + wind.ST = krigeST(sqrt(values)~1, w, STF(grd, tgrd), v) + + all.equal(wind.pr0, wind.ST) + + # stars: + df = data.frame(a = rep(NA, 324*10)) + s = STF(grd, tgrd) + newd = addAttrToGeom(s, df) + wind.sta = krigeST(sqrt(values)~1, st_as_stars(w), st_as_stars(newd), v) + # 1 + plot(stars::st_as_stars(wind.ST), breaks = "equal", col = sf.colors()) + # 2 + stplot(wind.ST) + # 3 + plot(wind.sta, breaks = "equal", col = sf.colors()) + st_as_stars(wind.ST)[[1]][1:3,1:3,1] + (wind.sta)[[1]][1:3,1:3,1] + st_bbox(wind.sta) + bbox(wind.ST) + all.equal(wind.sta, stars::st_as_stars(wind.ST), check.attributes = FALSE) + + # 4: roundtrip wind.sta->STFDF->stars + rt = stars::st_as_stars(as(wind.sta, "STFDF")) + plot(rt, breaks = "equal", col = sf.colors()) + # 5: + stplot(as(wind.sta, "STFDF")) + st_bbox(rt) + + # 6: + stplot(as(st_as_stars(wind.ST), "STFDF")) + } Error in `/.difftime`(diff(range(ud)), mean(ud)) : second argument of / cannot be a "difftime" object Calls: krigeST ... create_dimension -> regular_intervals -> isTRUE -> /.difftime In addition: Warning message: In krigeST(sqrt(values) ~ 1, w, STF(grd, tgrd), v) : The spatio-temporal variogram model does not carry the strongly recommended attribute 'temporal unit'. The unit 'secs' has been assumed. krigeST could not check whether the temporal distances between locations and in the variogram coincide. Execution halted Package: sftime Check: examples New result: ERROR Running examples in ‘sftime-Ex.R’ failed The error most likely occurred in: > base::assign(".ptime", proc.time(), pos = "CheckExEnv") > ### Name: st_as_sftime > ### Title: Convert a foreign object to an 'sftime' object > ### Aliases: st_as_sftime st_as_sftime.ST st_as_sftime.Track > ### st_as_sftime.Tracks st_as_sftime.TracksCollection st_as_sftime.sftime > ### st_as_sftime.sf st_as_sftime.stars st_as_sftime.data.frame > ### st_as_sftime.ppp st_as_sftime.psp st_as_sftime.lpp > ### st_as_sftime.sftrack st_as_sftime.sftraj st_as_sftime.cubble_df > > ### ** Examples > > # modified from spacetime: > library(sp) > library(spacetime) > > sp <- cbind(x = c(0,0,1), y = c(0,1,1)) > row.names(sp) <- paste("point", 1:nrow(sp), sep="") > sp <- SpatialPoints(sp) > time <- as.POSIXct("2010-08-05") + 3600 * (10:12) > x <- STI(sp, time) > > st_as_sftime(x) Spatiotemporal feature collection with 3 features and 0 fields Geometry type: POINT Dimension: XY Bounding box: xmin: 0 ymin: 0 xmax: 1 ymax: 1 CRS: NA Time column with classes: 'POSIXct', 'POSIXt'. Ranging from 2010-08-05 10:00:00 to 2010-08-05 12:00:00. st_as_sfc.x.sp. time 1 POINT (0 0) 2010-08-05 10:00:00 2 POINT (0 1) 2010-08-05 11:00:00 3 POINT (1 1) 2010-08-05 12:00:00 > > # convert a Track object from package trajectories to an sftime object > library(trajectories) > x1_Track <- trajectories::rTrack(n = 100) > x1_Track@data$speed <- sort(rnorm(length(x1_Track))) > x1_sftime <- st_as_sftime(x1_Track) > > # convert a Tracks object from package trajectories to an sftime object > x2_Tracks <- trajectories::rTracks(m = 6) > x2_sftime <- st_as_sftime(x2_Tracks) > > # convert a TracksCollection object from package trajectories to an sftime object > x3_TracksCollection <- trajectories::rTracksCollection(p = 2, m = 3, n = 50) > x3_sftime <- st_as_sftime(x3_TracksCollection) > > # convert an sftime object to an sftime object > st_as_sftime(x3_sftime) Spatiotemporal feature collection with 300 features and 3 fields Geometry type: POINT Dimension: XY Bounding box: xmin: -28.35517 ymin: -24.66491 xmax: 38.77599 ymax: 52.90631 CRS: NA Time column with classes: 'POSIXct', 'POSIXt'. Ranging from 1970-01-01 to 1970-01-01 04:49:00. First 10 features: ones tracks_name track_name st_as_sfc.x.sp. 1 1 Tracks1 Track1 POINT (-1.539007 1.932557) 2 1 Tracks1 Track1 POINT (-0.2003776 -0.1679942) 3 1 Tracks1 Track1 POINT (-2.752113 3.166629) 4 1 Tracks1 Track1 POINT (-0.2337229 -0.3806908) 5 1 Tracks1 Track1 POINT (-3.585336 3.267737) 6 1 Tracks1 Track1 POINT (-0.2504837 0.3498938) 7 1 Tracks1 Track1 POINT (-2.653028 1.426368) 8 1 Tracks1 Track1 POINT (0.6410695 -0.5117434) 9 1 Tracks1 Track1 POINT (-2.018482 1.572974) 10 1 Tracks1 Track1 POINT (1.737529 -1.27253) time 1 1970-01-01 00:00:00 2 1970-01-01 00:00:00 3 1970-01-01 00:01:00 4 1970-01-01 00:01:00 5 1970-01-01 00:02:00 6 1970-01-01 00:02:00 7 1970-01-01 00:03:00 8 1970-01-01 00:03:00 9 1970-01-01 00:04:00 10 1970-01-01 00:04:00 > > # convert an sf object to an sftime object > g <- st_sfc(st_point(c(1, 2)), st_point(c(1, 3)), st_point(c(2, 3)), + st_point(c(2, 1)), st_point(c(3, 1))) > x4_sf <- st_sf(a = 1:5, g, time = Sys.time() + 1:5) > x4_sftime <- st_as_sftime(x4_sf) > > # convert a Tracks object from package trajectories to an sftime object > x5_stars <- stars::read_stars(system.file("nc/bcsd_obs_1999.nc", package = "stars")) pr, Error in `/.difftime`(diff(range(ud)), mean(ud)) : second argument of / cannot be a "difftime" object Calls: ... create_dimension -> regular_intervals -> isTRUE -> /.difftime Execution halted Package: stars Check: examples New result: ERROR Running examples in ‘stars-Ex.R’ failed The error most likely occurred in: > base::assign(".ptime", proc.time(), pos = "CheckExEnv") > ### Name: aggregate.stars > ### Title: spatially or temporally aggregate stars object > ### Aliases: aggregate.stars aggregate > > ### ** Examples > > # aggregate time dimension in format Date > tif = system.file("tif/L7_ETMs.tif", package = "stars") > t1 = as.Date("2018-07-31") > x = read_stars(c(tif, tif, tif, tif), along = list(time = c(t1, t1+1, t1+2, t1+3)))[,1:30,1:30] Error in `/.difftime`(diff(range(ud)), mean(ud)) : second argument of / cannot be a "difftime" object Calls: read_stars ... create_dimension -> regular_intervals -> isTRUE -> /.difftime Execution halted Package: stars Check: tests New result: ERROR Running ‘aggregate.R’ [2s/2s] Running ‘align.R’ [1s/1s] Comparing ‘align.Rout’ to ‘align.Rout.save’ ... OK Running ‘area.R’ [1s/1s] Comparing ‘area.Rout’ to ‘area.Rout.save’ ... OK Running ‘crop.R’ [4s/4s] Comparing ‘crop.Rout’ to ‘crop.Rout.save’ ... OK Running ‘curvilinear.R’ [1s/1s] Comparing ‘curvilinear.Rout’ to ‘curvilinear.Rout.save’ ... OK Running ‘datasets.R’ [1s/1s] Comparing ‘datasets.Rout’ to ‘datasets.Rout.save’ ... OK Running ‘dimensions.R’ [1s/1s] Running ‘downsample.R’ [1s/1s] Comparing ‘downsample.Rout’ to ‘downsample.Rout.save’ ... OK Running ‘ee.R’ [1s/1s] Comparing ‘ee.Rout’ to ‘ee.Rout.save’ ... OK Running ‘extract.R’ [1s/1s] Comparing ‘extract.Rout’ to ‘extract.Rout.save’ ... OK Running ‘gridtypes.R’ [1s/1s] Running ‘mdim.R’ [1s/1s] Running ‘nc.R’ [2s/2s] Comparing ‘nc.Rout’ to ‘nc.Rout.save’ ... OK Running ‘plot.R’ [2s/2s] Running ‘predict.R’ [4s/4s] Comparing ‘predict.Rout’ to ‘predict.Rout.save’ ... OK Running ‘proxy.R’ [2s/2s] Comparing ‘proxy.Rout’ to ‘proxy.Rout.save’ ... OK Running ‘raster.R’ [9s/9s] Comparing ‘raster.Rout’ to ‘raster.Rout.save’ ... OK Running ‘rasterize.R’ [2s/2s] Comparing ‘rasterize.Rout’ to ‘rasterize.Rout.save’ ... OK Running ‘rectilinear.R’ [1s/1s] Comparing ‘rectilinear.Rout’ to ‘rectilinear.Rout.save’ ... OK Running ‘redimension.R’ [1s/1s] Running ‘sf.R’ [1s/1s] Running ‘sp.R’ [2s/2s] Comparing ‘sp.Rout’ to ‘sp.Rout.save’ ... OK Running ‘spacetime.R’ [1s/1s] Running ‘spatstat.R’ [2s/2s] Comparing ‘spatstat.Rout’ to ‘spatstat.Rout.save’ ... OK Running ‘stars.R’ [8s/8s] Comparing ‘stars.Rout’ to ‘stars.Rout.save’ ... OK Running ‘subset.R’ [1s/1s] Comparing ‘subset.Rout’ to ‘subset.Rout.save’ ... OK Running ‘testthat.R’ [17s/23s] Running ‘tidy.R’ [3s/4s] Comparing ‘tidy.Rout’ to ‘tidy.Rout.save’ ... OK Running ‘transform.R’ [2s/3s] Comparing ‘transform.Rout’ to ‘transform.Rout.save’ ... OK Running ‘warp.R’ [3s/4s] Comparing ‘warp.Rout’ to ‘warp.Rout.save’ ... OK Running ‘write.R’ [2s/2s] Comparing ‘write.Rout’ to ‘write.Rout.save’ ... OK Running the tests in ‘tests/aggregate.R’ failed. Complete output: > ## IGNORE_RDIFF_BEGIN > library(stars) Loading required package: abind Loading required package: sf Linking to GEOS 3.13.0, GDAL 3.10.1, PROJ 9.5.1; sf_use_s2() is TRUE > ## IGNORE_RDIFF_END > > 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_set_crs(st_as_sfc(stations), 4326), time = dates) + + blocks = st_make_grid(st_as_sfc("POLYGON ((5.871619 47.26986, 15.03811 47.26986, 15.03811 55.05653, 5.871619 55.05653, 5.871619 47.26986))", crs = 4326), + n = c(3,3)) + + print(aq <- st_as_stars(list(PM10 = air), dimensions = d)) + print(a <- aggregate(aq, blocks, mean, na.rm = TRUE)) + } Error in `/.difftime`(diff(range(ud)), mean(ud)) : second argument of / cannot be a "difftime" object Calls: st_dimensions ... -> regular_intervals -> isTRUE -> /.difftime Execution halted Running the tests in ‘tests/dimensions.R’ failed. Complete output: > # https://github.com/r-spatial/stars/issues/295 > suppressPackageStartupMessages(library(stars)) > raw <- read_stars(system.file("nc/bcsd_obs_1999.nc", package = "stars")) pr, Error in `/.difftime`(diff(range(ud)), mean(ud)) : second argument of / cannot be a "difftime" object Calls: read_stars ... create_dimension -> regular_intervals -> isTRUE -> /.difftime Execution halted Running the tests in ‘tests/gridtypes.R’ failed. Complete output: > suppressPackageStartupMessages(library(stars)) > > # regular, but not spatial: > d = st_dimensions(a = 1:3, b = 1:3, band = c("foo", "bar")) > (st = st_as_stars(array(1:18, c(3,3,2)), dimension = d)) stars object with 3 dimensions and 1 attribute attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. X 1 5.25 9.5 9.5 13.75 18 dimension(s): from to offset delta point values a 1 3 1 1 FALSE NULL b 1 3 1 1 FALSE NULL band 1 2 NA NA FALSE foo, bar > as.data.frame(st, add_max=FALSE)[1:4,] a b band X 1 1 1 foo 1 2 2 1 foo 2 3 3 1 foo 3 4 1 2 foo 4 > try(x <- st_bbox(st)) # error Error in st_bbox.dimensions(st_dimensions(obj), ...) : dimensions table does not have x & y, nor an sfc dimension > > # regular, geotransform: > d = st_dimensions(x = 1:3, y = 1:3, band = c("foo", "bar"), .raster = c("x", "y")) > (st = st_as_stars(array(1:18, c(3,3,2)), dimension = d)) stars object with 3 dimensions and 1 attribute attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. X 1 5.25 9.5 9.5 13.75 18 dimension(s): from to offset delta point values x/y x 1 3 1 1 FALSE NULL [x] y 1 3 1 1 FALSE NULL [y] band 1 2 NA NA FALSE foo, bar > as.data.frame(st, add_max=FALSE)[1:4,] x y band X 1 1.5 1.5 foo 1 2 2.5 1.5 foo 2 3 3.5 1.5 foo 3 4 1.5 2.5 foo 4 > as.data.frame(st, add_max=TRUE)[1:4,] x y band x_max y_max X 1 1 1 foo 2 2 1 2 2 1 foo 3 2 2 3 3 1 foo 4 2 3 4 1 2 foo 2 3 4 > st_bbox(st) xmin ymin xmax ymax 1 1 4 4 > > # rectilinear with offset given: > xd = c(1, 2, 4) > d = st_dimensions(x = xd, y = 1:3, band = c("foo", "bar"), .raster = c("x", "y")) > (st = st_as_stars(array(1:18, c(3,3,2)), dimension = d)) stars object with 3 dimensions and 1 attribute attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. X 1 5.25 9.5 9.5 13.75 18 dimension(s): from to offset delta point values x/y x 1 3 NA NA FALSE [1,2), [2,4), [4,6) [x] y 1 3 1 1 FALSE NULL [y] band 1 2 NA NA FALSE foo, bar > as.data.frame(st, add_max=FALSE)[1:4,] x y band X 1 1.5 1.5 foo 1 2 3.0 1.5 foo 2 3 5.0 1.5 foo 3 4 1.5 2.5 foo 4 > as.data.frame(st, add_max=TRUE)[1:4,] x y band x_max y_max X 1 1 1 foo 2 2 1 2 2 1 foo 4 2 2 3 4 1 foo 6 2 3 4 1 2 foo 2 3 4 > st_bbox(st) xmin ymin xmax ymax 1 1 6 4 > > # rectilinear with midpoints given: > xd = c(1, 2, 4) > d = st_dimensions(x = xd, y = 1:3, band = c("foo", "bar"), .raster = c("x", "y"), cell_midpoints = TRUE) > (st = st_as_stars(array(1:18, c(3,3,2)), dimension = d)) stars object with 3 dimensions and 1 attribute attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. X 1 5.25 9.5 9.5 13.75 18 dimension(s): from to offset delta point values x/y x 1 3 NA NA FALSE [0.5,1.5), [1.5,3.0), [3.0,5.0) [x] y 1 3 0.5 1 FALSE NULL [y] band 1 2 NA NA FALSE foo, bar > as.data.frame(st, add_max=FALSE)[1:4,] x y band X 1 1.00 1 foo 1 2 2.25 1 foo 2 3 4.00 1 foo 3 4 1.00 2 foo 4 > as.data.frame(st, add_max=TRUE)[1:4,] x y band x_max y_max X 1 0.5 0.5 foo 1.5 1.5 1 2 1.5 0.5 foo 3.0 1.5 2 3 3.0 0.5 foo 5.0 1.5 3 4 0.5 1.5 foo 1.5 2.5 4 > st_bbox(st) xmin ymin xmax ymax 0.5 0.5 5.0 3.5 > > # rectilinear with midpoints given, point support > xd = c(1, 2, 4) > d = st_dimensions(x = xd, y = 1:3, band = c("foo", "bar"), .raster = c("x", "y"), cell_midpoints = TRUE, point = TRUE) > (st = st_as_stars(array(1:18, c(3,3,2)), dimension = d)) stars object with 3 dimensions and 1 attribute attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. X 1 5.25 9.5 9.5 13.75 18 dimension(s): from to offset delta point values x/y x 1 3 NA NA TRUE 1, 2, 4 [x] y 1 3 0.5 1 TRUE NULL [y] band 1 2 NA NA TRUE foo, bar > as.data.frame(st, add_max=FALSE)[1:4,] x y band X 1 1 1 foo 1 2 2 1 foo 2 3 4 1 foo 3 4 1 2 foo 4 > as.data.frame(st, add_max=TRUE)[1:4,] x y band x_max y_max X 1 1 0.5 foo 1 1.5 1 2 2 0.5 foo 2 1.5 2 3 4 0.5 foo 4 1.5 3 4 1 1.5 foo 1 2.5 4 > st_bbox(st) xmin ymin xmax ymax 1.0 0.5 4.0 3.5 > > # rectilinear with start/end given: > #xd = stars:::make_intervals(start = c(1,2,4), end = c(2, 4, 8)) > xd = c(1, 2, 4, 8) # one more than dim of the data array > d = st_dimensions(x = xd, y = 1:3, band = c("foo", "bar"), .raster = c("x", "y")) > (st = st_as_stars(array(1:18, c(3,3,2)), dimension = d)) stars object with 3 dimensions and 1 attribute attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. X 1 5.25 9.5 9.5 13.75 18 dimension(s): from to offset delta point values x/y x 1 3 NA NA FALSE [1,2), [2,4), [4,8) [x] y 1 3 1 1 FALSE NULL [y] band 1 2 NA NA FALSE foo, bar > as.data.frame(st, add_max=FALSE)[1:4,] x y band X 1 1.5 1.5 foo 1 2 3.0 1.5 foo 2 3 6.0 1.5 foo 3 4 1.5 2.5 foo 4 > as.data.frame(st, add_max=TRUE)[1:4,] x y band x_max y_max X 1 1 1 foo 2 2 1 2 2 1 foo 4 2 2 3 4 1 foo 8 2 3 4 1 2 foo 2 3 4 > st_bbox(st) xmin ymin xmax ymax 1 1 8 4 > > # with sfc: > sfc = st_sfc(st_point(0:1), st_point(2:1), st_point(4:3)) > d = st_dimensions(x = sfc, y = 1:3, band = c("foo", "bar")) > (st = st_as_stars(array(1:18, c(3,3,2)), dimension = d)) stars object with 3 dimensions and 1 attribute attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. X 1 5.25 9.5 9.5 13.75 18 dimension(s): from to offset delta point values x 1 3 NA NA TRUE POINT (0 1),...,POINT (4 3) y 1 3 1 1 FALSE NULL band 1 2 NA NA FALSE foo, bar > as.data.frame(st, add_max=FALSE)[1:4,] x y band X 1 POINT (0 1) 1 foo 1 2 POINT (2 1) 1 foo 2 3 POINT (4 3) 1 foo 3 4 POINT (0 1) 2 foo 4 > st_as_sf(st, long = TRUE)[1:4,] Simple feature collection with 4 features and 3 fields Geometry type: POINT Dimension: XY Bounding box: xmin: 0 ymin: 1 xmax: 4 ymax: 3 CRS: NA y band X x 1 1 foo 1 POINT (0 1) 2 1 foo 2 POINT (2 1) 3 1 foo 3 POINT (4 3) 4 2 foo 4 POINT (0 1) > st_bbox(st) xmin ymin xmax ymax 0 1 4 3 > > # rotated/sheared: > d = st_dimensions(x = 1:3, y = 1:3, band = c("foo", "bar"), .raster = c("x", "y"), affine = c(0.2, -0.2)) > (st = st_as_stars(array(1:18, c(3,3,2)), dimension = d)) stars object with 3 dimensions and 1 attribute attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. X 1 5.25 9.5 9.5 13.75 18 dimension(s): from to offset delta point values x/y x 1 3 1 1 FALSE NULL [x] y 1 3 1 1 FALSE NULL [y] band 1 2 NA NA FALSE foo, bar sheared raster with parameters: 0.2 -0.2 > as.data.frame(st, add_max=FALSE)[1:4,] x y X 1 1.6 1.4 1 2 2.6 1.2 2 3 3.6 1.0 3 4 1.8 2.4 4 > try(as.data.frame(st, add_max=TRUE)[1:4,]) # errors Error in st_coordinates.stars(x, add_max = add_max, center = center, ...) : add_max will not work for rotated/shared rasters > st_bbox(st) xmin ymin xmax ymax 1.0 0.4 4.6 4.0 > > # curvilinear: > set.seed(13531) > lon = st_as_stars(matrix(signif(runif(9), 2), 3, 3)) > lat = st_as_stars(matrix(signif(runif(9), 2), 3, 3)) > ll = setNames(c(lon, lat), c("lon", "lat")) > d = st_dimensions(lon = 1:3, lat = 1:3) > a = st_as_stars(list(X = array(1:9, c(3,3))), dimensions = d) > (st = st_as_stars(a, curvilinear = ll)) stars object with 2 dimensions and 1 attribute attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. X 1 3 5 5 7 9 dimension(s): from to refsys point values x/y lon 1 3 WGS 84 (CRS84) FALSE [3x3] 0.03,...,0.89 [x] lat 1 3 WGS 84 (CRS84) FALSE [3x3] 0.021,...,0.89 [y] curvilinear grid > as.data.frame(st) lon lat X 1 0.89 0.250 1 2 0.64 0.510 2 3 0.30 0.670 3 4 0.69 0.890 4 5 0.47 0.890 5 6 0.26 0.190 6 7 0.57 0.021 7 8 0.27 0.220 8 9 0.03 0.340 9 > st_bbox(st) xmin ymin xmax ymax 0.030 0.021 0.890 0.890 > > # time regular: > t = as.Date("2019-01-20") + 0:3 > d = st_dimensions(t, sfc) Error in `/.difftime`(diff(range(ud)), mean(ud)) : second argument of / cannot be a "difftime" object Calls: st_dimensions ... -> regular_intervals -> isTRUE -> /.difftime Execution halted Running the tests in ‘tests/mdim.R’ failed. Complete output: > Sys.setenv(TZ="") > suppressPackageStartupMessages(library(stars)) > matrix(1:4,2) [,1] [,2] [1,] 1 3 [2,] 2 4 > m = matrix(1:4,2) > s = st_as_stars(m) > > s = st_set_dimensions(s, st_dimensions(s), xy = NULL) > sfc = st_sfc(st_point(0:1), st_point(3:2)) > s = st_set_dimensions(s, 1, sfc) |> st_set_crs(4326) > tm = as.POSIXct("2023-03-04 12:35") + c(0, 3600) > s = st_set_dimensions(s, 2, tm) Error in `/.difftime`(diff(range(ud)), mean(ud)) : second argument of / cannot be a "difftime" object Calls: st_set_dimensions ... create_dimension -> regular_intervals -> isTRUE -> /.difftime Execution halted Running the tests in ‘tests/plot.R’ failed. Complete output: > 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) + } Error in `/.difftime`(diff(range(ud)), mean(ud)) : second argument of / cannot be a "difftime" object Calls: st_dimensions ... -> regular_intervals -> isTRUE -> /.difftime Execution halted Running the tests in ‘tests/redimension.R’ failed. Complete output: > suppressPackageStartupMessages(library(stars)) > tif = system.file("tif/L7_ETMs.tif", package = "stars") > r = read_stars(list(a = c(tif,tif), b = c(tif, tif)), proxy = TRUE) > (xx = st_redimension(r, along = list(foo = 1:4))) stars_proxy object with 1 attribute in 4 file(s): $a.L7_ETMs.tif.a.L7_ETMs.tif.b.L7_ETMs.tif.b.L7_ETMs.tif [1] "[...]/L7_ETMs.tif" "[...]/L7_ETMs.tif" "[...]/L7_ETMs.tif" [4] "[...]/L7_ETMs.tif" dimension(s): from to offset delta refsys point x/y x 1 349 288776 28.5 SIRGAS 2000 / UTM zone 25S FALSE [x] y 1 352 9120761 -28.5 SIRGAS 2000 / UTM zone 25S FALSE [y] band 1 6 NA NA NA NA foo 1 4 1 1 NA NA > (rr = c(xx, xx)) stars_proxy object with 2 attributes in 8 file(s): $a.L7_ETMs.tif.a.L7_ETMs.tif.b.L7_ETMs.tif.b.L7_ETMs.tif [1] "[...]/L7_ETMs.tif" "[...]/L7_ETMs.tif" "[...]/L7_ETMs.tif" [4] "[...]/L7_ETMs.tif" $a.L7_ETMs.tif.a.L7_ETMs.tif.b.L7_ETMs.tif.b.L7_ETMs.tif [1] "[...]/L7_ETMs.tif" "[...]/L7_ETMs.tif" "[...]/L7_ETMs.tif" [4] "[...]/L7_ETMs.tif" dimension(s): from to offset delta refsys point x/y x 1 349 288776 28.5 SIRGAS 2000 / UTM zone 25S FALSE [x] y 1 352 9120761 -28.5 SIRGAS 2000 / UTM zone 25S FALSE [y] band 1 6 NA NA NA NA foo 1 4 1 1 NA NA > # st_as_stars(rr) > (rrr = st_redimension(rr, along = list(bar = as.Date(c("2001-01-01", "2002-01-01"))))) Error in `/.difftime`(diff(range(ud)), mean(ud)) : second argument of / cannot be a "difftime" object Calls: st_redimension ... create_dimension -> regular_intervals -> isTRUE -> /.difftime Execution halted Running the tests in ‘tests/sf.R’ failed. Complete output: > 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]) stars object with 3 dimensions and 1 attribute attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. L7_ETMs.tif 17 44 58 58.19667 70 130 dimension(s): from to offset delta refsys point x/y x 1 20 288776 28.5 SIRGAS 2000 / UTM zone 25S FALSE [x] y 1 20 9120761 -28.5 SIRGAS 2000 / UTM zone 25S FALSE [y] band 1 6 NA NA NA NA > # names(x) = "L7_ETM" > > # library(abind) > # x = adrop(x) > image(x) > (sfc = st_as_sfc(x, as_points = FALSE)) Geometry set for 400 features Geometry type: POLYGON Dimension: XY Bounding box: xmin: 288776.3 ymin: 9120191 xmax: 289346.3 ymax: 9120761 Projected CRS: SIRGAS 2000 / UTM zone 25S First 5 geometries: POLYGON ((288776.3 9120761, 288804.8 9120761, 2... POLYGON ((288804.8 9120761, 288833.3 9120761, 2... POLYGON ((288833.3 9120761, 288861.8 9120761, 2... POLYGON ((288861.8 9120761, 288890.3 9120761, 2... POLYGON ((288890.3 9120761, 288918.8 9120761, 2... > plot(sfc, add =TRUE) > (sfc = st_as_sfc(x, as_points = TRUE)) Geometry set for 400 features Geometry type: POINT Dimension: XY Bounding box: xmin: 288790.5 ymin: 9120205 xmax: 289332 ymax: 9120747 Projected CRS: SIRGAS 2000 / UTM zone 25S First 5 geometries: POINT (288790.5 9120747) POINT (288819 9120747) POINT (288847.5 9120747) POINT (288876 9120747) POINT (288904.5 9120747) > 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) [1] TRUE > > 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) [1] TRUE > > # 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) stars object with 1 dimensions and 3 attributes attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. a 1 1.5 2 2 2.5 3 b 4 4.5 5 5 5.5 6 c 7 7.5 8 8 8.5 9 dimension(s): from to point values geom 1 3 TRUE POINT (0 1),...,POINT (2 3) > st_redimension(x) stars object with 2 dimensions and 1 attribute attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. a.b.c 1 3 5 5 7 9 dimension(s): from to point values geom 1 3 TRUE POINT (0 1),...,POINT (2 3) new_dim 1 3 NA a, b, c > c(x, x) stars object with 1 dimensions and 6 attributes attribute(s): Min. 1st Qu. Median Mean 3rd Qu. Max. a 1 1.5 2 2 2.5 3 b 4 4.5 5 5 5.5 6 c 7 7.5 8 8 8.5 9 a.1 1 1.5 2 2 2.5 3 b.1 4 4.5 5 5 5.5 6 c.1 7 7.5 8 8 8.5 9 dimension(s): from to point values geom 1 3 TRUE POINT (0 1),...,POINT (2 3) > (y = st_redimension(x, along = list(time = as.Date("2017-11-27") + 0:2))) Error in `/.difftime`(diff(range(ud)), mean(ud)) : second argument of / cannot be a "difftime" object Calls: st_redimension ... create_dimension -> regular_intervals -> isTRUE -> /.difftime Execution halted Running the tests in ‘tests/spacetime.R’ failed. Complete output: > suppressPackageStartupMessages(library(stars)) > if (suppressPackageStartupMessages(require(sp, quietly = TRUE)) && require(spacetime, quietly = TRUE)) { + set.seed(1331) + Sys.setenv(TZ="UTC") + example(STFDF) + + x = st_as_stars(stfdf) + xx = as(x, "STFDF") + + # as.xts.stars: + library(xts) + xts = as.xts(x) + st_as_stars(xts) + st_as_stars(xts, dimensions = st_dimensions(x)) + + library(sp) + pts = SpatialPoints(rbind(c(0,0), c(0,1), c(1,0), c(1,1))) + gridded(pts) = TRUE + + tm = xts(1:10, as.Date("2010-01-01")+0:9) + stfdf = STFDF(pts, tm, data.frame(foo = seq_len(length(pts) * length(tm)))) + x = st_as_stars(stfdf) + xx = as(x, "STFDF") + + x = st_as_stars(as(stfdf, "STSDF")) + } STFDF> sp = cbind(x = c(0,0,1), y = c(0,1,1)) STFDF> row.names(sp) = paste("point", 1:nrow(sp), sep="") STFDF> library(sp) STFDF> sp = SpatialPoints(sp) STFDF> time = as.POSIXct("2010-08-05")+3600*(10:13) STFDF> m = c(10,20,30) # means for each of the 3 point locations STFDF> mydata = rnorm(length(sp)*length(time),mean=rep(m, 4)) STFDF> IDs = paste("ID",1:length(mydata)) STFDF> mydata = data.frame(values = signif(mydata,3), ID=IDs) STFDF> stfdf = STFDF(sp, time, mydata) STFDF> stfdf An object of class "STFDF" Slot "data": values ID 1 9.38 ID 1 2 21.10 ID 2 3 29.80 ID 3 4 11.50 ID 4 5 19.80 ID 5 6 31.20 ID 6 7 10.40 ID 7 8 20.50 ID 8 9 29.80 ID 9 10 10.30 ID 10 11 21.40 ID 11 12 29.40 ID 12 Slot "sp": SpatialPoints: x y point1 0 0 point2 0 1 point3 1 1 Coordinate Reference System (CRS) arguments: NA Slot "time": timeIndex 2010-08-05 10:00:00 1 2010-08-05 11:00:00 2 2010-08-05 12:00:00 3 2010-08-05 13:00:00 4 Slot "endTime": [1] "2010-08-05 11:00:00 UTC" "2010-08-05 12:00:00 UTC" [3] "2010-08-05 13:00:00 UTC" "2010-08-05 14:00:00 UTC" STFDF> stfdf[1:2,] An object of class "STFDF" Slot "data": values ID 1 9.38 ID 1 2 21.10 ID 2 3 11.50 ID 4 4 19.80 ID 5 5 10.40 ID 7 6 20.50 ID 8 7 10.30 ID 10 8 21.40 ID 11 Slot "sp": SpatialPoints: x y point1 0 0 point2 0 1 Coordinate Reference System (CRS) arguments: NA Slot "time": timeIndex 2010-08-05 10:00:00 1 2010-08-05 11:00:00 2 2010-08-05 12:00:00 3 2010-08-05 13:00:00 4 Slot "endTime": [1] "2010-08-05 11:00:00 UTC" "2010-08-05 12:00:00 UTC" [3] "2010-08-05 13:00:00 UTC" "2010-08-05 14:00:00 UTC" STFDF> stfdf[,1:2] An object of class "STFDF" Slot "data": values ID 1 9.38 ID 1 2 21.10 ID 2 3 29.80 ID 3 4 11.50 ID 4 5 19.80 ID 5 6 31.20 ID 6 Slot "sp": SpatialPoints: x y point1 0 0 point2 0 1 point3 1 1 Coordinate Reference System (CRS) arguments: NA Slot "time": timeIndex 2010-08-05 10:00:00 1 2010-08-05 11:00:00 2 Slot "endTime": [1] "2010-08-05 11:00:00 UTC" "2010-08-05 12:00:00 UTC" STFDF> stfdf[,,2] An object of class "STFDF" Slot "data": ID 1 ID 1 2 ID 2 3 ID 3 4 ID 4 5 ID 5 6 ID 6 7 ID 7 8 ID 8 9 ID 9 10 ID 10 11 ID 11 12 ID 12 Slot "sp": SpatialPoints: x y point1 0 0 point2 0 1 point3 1 1 Coordinate Reference System (CRS) arguments: NA Slot "time": timeIndex 2010-08-05 10:00:00 1 2010-08-05 11:00:00 2 2010-08-05 12:00:00 3 2010-08-05 13:00:00 4 Slot "endTime": [1] "2010-08-05 11:00:00 UTC" "2010-08-05 12:00:00 UTC" [3] "2010-08-05 13:00:00 UTC" "2010-08-05 14:00:00 UTC" STFDF> stfdf[,,"values"] An object of class "STFDF" Slot "data": values 1 9.38 2 21.10 3 29.80 4 11.50 5 19.80 6 31.20 7 10.40 8 20.50 9 29.80 10 10.30 11 21.40 12 29.40 Slot "sp": SpatialPoints: x y point1 0 0 point2 0 1 point3 1 1 Coordinate Reference System (CRS) arguments: NA Slot "time": timeIndex 2010-08-05 10:00:00 1 2010-08-05 11:00:00 2 2010-08-05 12:00:00 3 2010-08-05 13:00:00 4 Slot "endTime": [1] "2010-08-05 11:00:00 UTC" "2010-08-05 12:00:00 UTC" [3] "2010-08-05 13:00:00 UTC" "2010-08-05 14:00:00 UTC" STFDF> stfdf[1,] values ID timeIndex 2010-08-05 10:00:00 9.38 ID 1 1 2010-08-05 11:00:00 11.50 ID 4 2 2010-08-05 12:00:00 10.40 ID 7 3 2010-08-05 13:00:00 10.30 ID 10 4 STFDF> stfdf[,2] coordinates values ID 1 (0, 0) 11.5 ID 4 2 (0, 1) 19.8 ID 5 3 (1, 1) 31.2 ID 6 STFDF> as(stfdf[,,1], "xts") point1 point2 point3 2010-08-05 10:00:00 9.38 21.1 29.8 2010-08-05 11:00:00 11.50 19.8 31.2 2010-08-05 12:00:00 10.40 20.5 29.8 2010-08-05 13:00:00 10.30 21.4 29.4 STFDF> as(stfdf[,,2], "xts") point1 point2 point3 2010-08-05 10:00:00 ID 1 ID 2 ID 3 2010-08-05 11:00:00 ID 4 ID 5 ID 6 2010-08-05 12:00:00 ID 7 ID 8 ID 9 2010-08-05 13:00:00 ID 10 ID 11 ID 12 STFDF> # examples for [[, [[<-, $ and $<- STFDF> stfdf[[1]] [1] 9.38 21.10 29.80 11.50 19.80 31.20 10.40 20.50 29.80 10.30 21.40 29.40 STFDF> stfdf[["values"]] [1] 9.38 21.10 29.80 11.50 19.80 31.20 10.40 20.50 29.80 10.30 21.40 29.40 STFDF> stfdf[["newVal"]] <- rnorm(12) STFDF> stfdf$ID [1] "ID 1" "ID 2" "ID 3" "ID 4" "ID 5" "ID 6" "ID 7" "ID 8" "ID 9" [10] "ID 10" "ID 11" "ID 12" STFDF> stfdf$ID = paste("OldIDs", 1:12, sep="") STFDF> stfdf$NewID = paste("NewIDs", 12:1, sep="") STFDF> stfdf An object of class "STFDF" Slot "data": values ID newVal NewID 1 9.38 OldIDs1 -0.3125824 NewIDs12 2 21.10 OldIDs2 -1.0536246 NewIDs11 3 29.80 OldIDs3 -0.7188880 NewIDs10 4 11.50 OldIDs4 -1.0525404 NewIDs9 5 19.80 OldIDs5 -0.8441413 NewIDs8 6 31.20 OldIDs6 1.1420848 NewIDs7 7 10.40 OldIDs7 1.9479070 NewIDs6 8 20.50 OldIDs8 -0.4512836 NewIDs5 9 29.80 OldIDs9 0.6695716 NewIDs4 10 10.30 OldIDs10 0.1269003 NewIDs3 11 21.40 OldIDs11 -0.7273297 NewIDs2 12 29.40 OldIDs12 -0.8151803 NewIDs1 Slot "sp": SpatialPoints: x y point1 0 0 point2 0 1 point3 1 1 Coordinate Reference System (CRS) arguments: NA Slot "time": timeIndex 2010-08-05 10:00:00 1 2010-08-05 11:00:00 2 2010-08-05 12:00:00 3 2010-08-05 13:00:00 4 Slot "endTime": [1] "2010-08-05 11:00:00 UTC" "2010-08-05 12:00:00 UTC" [3] "2010-08-05 13:00:00 UTC" "2010-08-05 14:00:00 UTC" STFDF> x = stfdf[stfdf[1:2,],] STFDF> all.equal(x, stfdf[1:2,]) [1] TRUE STFDF> all.equal(stfdf, stfdf[stfdf,]) # converts character to factor... [1] TRUE Error in `/.difftime`(diff(range(ud)), mean(ud)) : second argument of / cannot be a "difftime" object Calls: st_as_stars ... create_dimension -> regular_intervals -> isTRUE -> /.difftime Execution halted Running the tests in ‘tests/testthat.R’ failed. Complete output: > ## load dependencies > library(testthat) > suppressPackageStartupMessages(library(stars)) > > ## test package > test_check("stars") [ FAIL 2 | WARN 1 | SKIP 4 | PASS 58 ] ══ Skipped tests (4) ═══════════════════════════════════════════════════════════ • On CRAN (4): 'test-OpenStreetMap.R:2:3', 'test-ncdf.R:1:1', 'test-ncproxy.R:1:1', 'test-rectilinear.R:89:3' ══ Failed tests ════════════════════════════════════════════════════════════════ ── Error ('test-cubble.R:30:3'): cubble ──────────────────────────────────────── Error in ``/.difftime`(diff(range(ud)), mean(ud))`: second argument of / cannot be a "difftime" object Backtrace: ▆ 1. ├─stars::st_as_stars(a.cb) at test-cubble.R:30:3 2. └─stars:::st_as_stars.cubble_df(a.cb) 3. ├─stars:::create_dimensions(...) 4. └─stars:::create_dimension(values = times) 5. └─stars:::regular_intervals(values) 6. ├─base::isTRUE(as.numeric(abs(diff(range(ud))/mean(ud))) < epsilon) 7. └─base::`/.difftime`(diff(range(ud)), mean(ud)) ── Error ('test-stars.R:14:3'): basic st_as_stars ────────────────────────────── Error in ``/.difftime`(diff(range(ud)), mean(ud))`: second argument of / cannot be a "difftime" object Backtrace: ▆ 1. ├─stars::st_as_stars(test_list) at test-stars.R:14:3 2. └─stars:::st_as_stars.ncdfgeom(test_list) 3. └─stars:::create_dimension(...) 4. └─stars:::regular_intervals(values) 5. ├─base::isTRUE(as.numeric(abs(diff(range(ud))/mean(ud))) < epsilon) 6. └─base::`/.difftime`(diff(range(ud)), mean(ud)) [ FAIL 2 | WARN 1 | SKIP 4 | PASS 58 ] Error: Test failures Execution halted Package: stxplore Check: examples New result: ERROR Running examples in ‘stxplore-Ex.R’ failed The error most likely occurred in: > base::assign(".ptime", proc.time(), pos = "CheckExEnv") > ### Name: cancor_eof > ### Title: Performs CCA using Empirical Orthogonal Functions (EOFs) from a > ### lagged dataset > ### Aliases: cancor_eof cancor_eof.data.frame cancor_eof.stars > ### autoplot.cancoreof > > ### ** Examples > > # Dataframe example > data(SSTlonlatshort) > data(SSTdatashort) > cancor_df <- cancor_eof(x = SSTlonlatshort, + lag = 7, + n_eof = 8, + values_df = SSTdatashort) Canonical Correlation Analysis is carried out using a lagged dataset. > autoplot(cancor_df) > > # Stars example > library(dplyr) Attaching package: ‘dplyr’ The following objects are masked from ‘package:stats’: filter, lag The following objects are masked from ‘package:base’: intersect, setdiff, setequal, union > library(stars) Loading required package: abind Loading required package: sf Linking to GEOS 3.13.0, GDAL 3.10.1, PROJ 9.5.1; sf_use_s2() is TRUE > # Create a stars object from a data frame > precip_df <- NOAA_df_1990[NOAA_df_1990$proc == 'Precip', ] %>% + filter(date >= "1992-02-01" & date <= "1992-02-28") > precip <- precip_df[ ,c('lat', 'lon', 'date', 'z')] > st_precip <- st_as_stars(precip, dims = c("lon", "lat", "date")) Error in `/.difftime`(diff(range(ud)), mean(ud)) : second argument of / cannot be a "difftime" object Calls: st_as_stars ... create_dimension -> regular_intervals -> isTRUE -> /.difftime Execution halted Package: stxplore Check: re-building of vignette outputs New result: ERROR Error(s) in re-building vignettes: ... --- re-building ‘stxplore.Rmd’ using rmarkdown ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/ssnap1-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 22312 bytes Input file size = 22414 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 18988 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 18988 Output IDAT size = 18988 bytes (3324 bytes decrease) Output file size = 19066 bytes (3348 bytes = 14.94% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/ssnap2-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 36500 bytes Input file size = 36626 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 33160 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 33160 Output IDAT size = 33160 bytes (3340 bytes decrease) Output file size = 33238 bytes (3388 bytes = 9.25% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/tsnap1-1.png 288x288 pixels, 8 bits/pixel, 256 colors in palette Reducing image to 8 bits/pixel, grayscale Input IDAT size = 16809 bytes Input file size = 17691 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 16449 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 16449 Output IDAT size = 16449 bytes (360 bytes decrease) Output file size = 16527 bytes (1164 bytes = 6.58% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/sem1-1.png 288x288 pixels, 8 bits/pixel, 256 colors in palette Reducing image to 8 bits/pixel, grayscale Input IDAT size = 13607 bytes Input file size = 14477 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 12756 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 12756 Output IDAT size = 12756 bytes (851 bytes decrease) Output file size = 12834 bytes (1643 bytes = 11.35% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/tem1-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 32297 bytes Input file size = 32411 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 IDAT size = 29510 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 5 IDAT size = 29510 Output IDAT size = 29510 bytes (2787 bytes decrease) Output file size = 29588 bytes (2823 bytes = 8.71% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/hovm1-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 22100 bytes Input file size = 22202 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 17483 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 17483 Output IDAT size = 17483 bytes (4617 bytes decrease) Output file size = 17561 bytes (4641 bytes = 20.90% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/hovm2-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 22686 bytes Input file size = 22788 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 17798 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 17798 Output IDAT size = 17798 bytes (4888 bytes decrease) Output file size = 17876 bytes (4912 bytes = 21.56% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/ridgeline-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 25661 bytes Input file size = 25775 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 20695 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 20695 Output IDAT size = 20695 bytes (4966 bytes decrease) Output file size = 20773 bytes (5002 bytes = 19.41% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/stecov1-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 14844 bytes Input file size = 14934 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 11122 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 11122 Output IDAT size = 11122 bytes (3722 bytes decrease) Output file size = 11200 bytes (3734 bytes = 25.00% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/stecov2-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 15178 bytes Input file size = 15268 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 11537 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 11537 Output IDAT size = 11537 bytes (3641 bytes decrease) Output file size = 11615 bytes (3653 bytes = 23.93% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/stecov3-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 14554 bytes Input file size = 14644 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 10961 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 10961 Output IDAT size = 10961 bytes (3593 bytes decrease) Output file size = 11039 bytes (3605 bytes = 24.62% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/stecov3-2.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 14830 bytes Input file size = 14920 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 11400 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 11400 Output IDAT size = 11400 bytes (3430 bytes decrease) Output file size = 11478 bytes (3442 bytes = 23.07% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/stsemiv3-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 12242 bytes Input file size = 12332 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 9671 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 9671 Output IDAT size = 9671 bytes (2571 bytes decrease) Output file size = 9749 bytes (2583 bytes = 20.95% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/stsemiv4-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 12020 bytes Input file size = 12110 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 9423 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 9423 Output IDAT size = 9423 bytes (2597 bytes decrease) Output file size = 9501 bytes (2609 bytes = 21.54% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/stsemiv5-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 11939 bytes Input file size = 12029 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 9449 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 9449 Output IDAT size = 9449 bytes (2490 bytes decrease) Output file size = 9527 bytes (2502 bytes = 20.80% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/eof1-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 31740 bytes Input file size = 31854 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 24750 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 24750 Output IDAT size = 24750 bytes (6990 bytes decrease) Output file size = 24828 bytes (7026 bytes = 22.06% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/eof2-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 30633 bytes Input file size = 30747 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 24216 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 24216 Output IDAT size = 24216 bytes (6417 bytes decrease) Output file size = 24294 bytes (6453 bytes = 20.99% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_files/figure-html/cancor-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 32482 bytes Input file size = 32596 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 25761 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 25761 Output IDAT size = 25761 bytes (6721 bytes decrease) Output file size = 25839 bytes (6757 bytes = 20.73% decrease) --- finished re-building ‘stxplore.Rmd’ --- re-building ‘stxplore_stars.Rmd’ using rmarkdown ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_stars_files/figure-html/ssnap1-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 55938 bytes Input file size = 56088 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 32827 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 32827 Output IDAT size = 32827 bytes (23111 bytes decrease) Output file size = 32905 bytes (23183 bytes = 41.33% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_stars_files/figure-html/tempsnap-1.png 288x288 pixels, 8 bits/pixel, 256 colors in palette Reducing image to 8 bits/pixel, grayscale Input IDAT size = 10090 bytes Input file size = 10960 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 9732 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 9732 Output IDAT size = 9732 bytes (358 bytes decrease) Output file size = 9810 bytes (1150 bytes = 10.49% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_stars_files/figure-html/spmeans-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 33385 bytes Input file size = 33511 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 26545 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 26545 Output IDAT size = 26545 bytes (6840 bytes decrease) Output file size = 26623 bytes (6888 bytes = 20.55% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_stars_files/figure-html/tempmeans-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 27896 bytes Input file size = 28010 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 IDAT size = 25255 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 5 IDAT size = 25255 Output IDAT size = 25255 bytes (2641 bytes decrease) Output file size = 25333 bytes (2677 bytes = 9.56% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_stars_files/figure-html/hovmoller1-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 12873 bytes Input file size = 12963 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 9816 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 9816 Output IDAT size = 9816 bytes (3057 bytes decrease) Output file size = 9894 bytes (3069 bytes = 23.68% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_stars_files/figure-html/hovmoller2-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 12259 bytes Input file size = 12349 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 9338 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 9338 Output IDAT size = 9338 bytes (2921 bytes decrease) Output file size = 9416 bytes (2933 bytes = 23.75% decrease) ** Processing: /home/hornik/tmp/CRAN_recheck/stxplore.Rcheck/vign_test/stxplore/vignettes/stxplore_stars_files/figure-html/ridgeline-1.png 288x288 pixels, 3x8 bits/pixel, RGB Input IDAT size = 23876 bytes Input file size = 23978 bytes Trying: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 19763 zc = 9 zm = 8 zs = 1 f = 0 zc = 1 zm = 8 zs = 2 f = 0 zc = 9 zm = 8 zs = 3 f = 0 zc = 9 zm = 8 zs = 0 f = 5 zc = 9 zm = 8 zs = 1 f = 5 zc = 1 zm = 8 zs = 2 f = 5 zc = 9 zm = 8 zs = 3 f = 5 Selecting parameters: zc = 9 zm = 8 zs = 0 f = 0 IDAT size = 19763 Output IDAT size = 19763 bytes (4113 bytes decrease) Output file size = 19841 bytes (4137 bytes = 17.25% decrease) Quitting from lines 131-137 [empcov1] (stxplore_stars.Rmd) Error: processing vignette 'stxplore_stars.Rmd' failed with diagnostics: second argument of / cannot be a "difftime" object --- failed re-building ‘stxplore_stars.Rmd’ SUMMARY: processing the following file failed: ‘stxplore_stars.Rmd’ Error: Vignette re-building failed. Execution halted