Running examples in ‘gdalcubes-Ex.R’ failed
The error most likely occurred in:

> ### Name: zonal_statistics
> ### Title: Query summary statistics of data cube values over polygons
> ### Aliases: zonal_statistics
> 
> ### ** Examples
> 
> # if not already done in other examples
> if (!file.exists(file.path(tempdir(), "L8.db"))) {
+   L8_files <- list.files(system.file("L8NY18", package = "gdalcubes"),
+                          ".TIF", recursive = TRUE, full.names = TRUE)
+   create_image_collection(L8_files, "L8_L1TP", file.path(tempdir(), "L8.db"))
+ }
> L8.col = image_collection(file.path(tempdir(), "L8.db"))
> v = cube_view(srs="EPSG:32618", dy=300, dx=300, dt="P1M", 
+               aggregation = "median", resampling = "bilinear",
+               extent=list(left=388941.2, right=766552.4,
+                           bottom=4345299, top=4744931, 
+                           t0="2018-01-01", t1="2018-04-30"))
## Size of the cube in x direction does not align with dx, extent will be enlarged by 44.400000 at both sides.
## Size of the cube in y direction does not align with dy, extent will be enlarged by 134.000000 at both sides.
> L8.cube = raster_cube(L8.col, v) 
> L8.cube = select_bands(L8.cube, c("B04", "B05")) 
> L8.ndvi = apply_pixel(L8.cube, "(B05-B04)/(B05+B04)", "NDVI") 
> L8.ndvi
A GDAL data cube proxy object

Dimensions:
       low     high count pixel_size chunk_size
t  2018-01  2018-04     4        P1M          1
y  4345165  4745065  1333        300        256
x 388896.8 766596.8  1259        300        256

Bands:
  name offset scale nodata unit
1 NDVI      0     1    NaN     

> 
> # toy example: overlay NDVI data with NYC districts
> x = zonal_statistics(L8.ndvi, system.file("nycd.gpkg", package = "gdalcubes"),
+                      expr = "median(NDVI)")
Error: Missing GEOS support in GDAL installation
Error in libgdalcubes_zonal_statistics(x, geom, agg_funcs, agg_bands,  : 
  Missing GEOS support in GDAL installation
Calls: zonal_statistics -> libgdalcubes_zonal_statistics
Execution halted