oisst_pawsey.md

    library(raadtools)
#> Loading required package: raster
#> Loading required package: sp
#> global option 'raadfiles.data.roots' set:
#> '
#> /rdsi/PRIVATE/raad/data               
#>  /rdsi/PRIVATE/raad/data_local         
#>  /rdsi/PRIVATE/raad/data_staging       
#>  /rdsi/PRIVATE/raad/data_deprecated    
#>  /rdsi/PUBLIC/raad/data
#> '
#> Uploading raad file cache as at 2024-07-22 04:21:31 (1608044 files listed)

    library(terra)
#> terra 1.7.81

files <- sstfiles()


files$oldname <- files$fullname
files$fullname <- gsub("/rdsi/PUBLIC/raad/data", "/vsicurl/https://projects.pawsey.org.au", files$fullname)
files$fullname <- sprintf("vrt://%s?unscale=true&ot=Float32&a_srs=EPSG:4326&sd_name=sst", files$fullname)

Sys.setenv("GDAL_DISABLE_READDIR_ON_OPEN"="EMPTY_DIR")
i <- sample(1:nrow(files), 1)
## nearest for warper because we don't want resampling
system.time(r1 <- project(terra::rast(files$fullname[i]), rast(res = 0.25), by_util = TRUE, method = "near"))
#>    user  system elapsed 
#>   0.687   0.042   4.669

system.time(r2 <- readsst(files$date[i], lon180 = TRUE))
#>    user  system elapsed 
#>   2.046   0.024   2.200

## very small difference
plot(r1 - rast(r2))

^{Created on 2024-07-22 with reprex v2.0.2}