timing_test_california_wa.md

dsn <- "/vsicurl/https://e84-earth-search-sentinel-data.s3.us-west-2.amazonaws.com/sentinel-2-c1-l2a/1/L/AC/2023/1/S2A_T01LAC_20230103T221938_L2A/TCI.tif"
library(vapour)
system.time(d <- gdal_raster_data(dsn, bands = 1:3, target_dim = c(1373, 1373)))

0.5 - 2 seconds in California

10-15 seconds in Western Australia or Tasmania

with terra, analogously

dsn <- "/vsicurl/https://e84-earth-search-sentinel-data.s3.us-west-2.amazonaws.com/sentinel-2-c1-l2a/1/L/AC/2023/1/S2A_T01LAC_20230103T221938_L2A/TCI.tif"
library(terra)
r <- rast(dsn)
system.time(d <- project(r, rast(ext(r), ncols = 1373, nrows = 1373, crs = crs(r)), by_util = TRUE))
#user  system elapsed 
#0.230   0.160   0.659

with gdalraster

dsn <- "/vsicurl/https://e84-earth-search-sentinel-data.s3.us-west-2.amazonaws.com/sentinel-2-c1-l2a/1/L/AC/2023/1/S2A_T01LAC_20230103T221938_L2A/TCI.tif"

            
library(gdalraster)
system.time({
warp(dsn, tf <- tempfile(fileext = ".tif"), t_srs = "", 
     cl_arg = c("-ts", 1373, 1373))

ds <- new(GDALRaster, tf)
dat <- read_ds(ds)
unlink(tf)
})
#   user  system elapsed 
#  0.212   0.042   0.620