Sys.setenv("AWS_NO_SIGN_REQUEST" = "YES")
dsn <- "/vsis3/us-west-2.opendata.source.coop/tge-labs/aef/v1/annual/2023/17N/xfbc3s1k0nrye0ub4-0000008192-0000000000.tiff"
r <- rast(sprintf("vrt://%s?ovr=3", dsn))[[1]]
par(mfrow = c(2, 1))
plot(r)
plot(r)
imdsn <- "WMTS:https://services.arcgisonline.com/arcgis/rest/services/World_Imagery/MapServer/WMTS/1.0.0/WMTSCapabilities.xml,layer=World_Imagery"
im <- vapour::gdal_raster_nara(imdsn, target_ext = c(336160, 418080, 4177920, 4259840), target_crs = "EPSG:32617", target_dim = c(1024, 0))
ximage::ximage(im, add = T)
Last active
April 14, 2026 23:15
-
-
Save mdsumner/6f592e6bb444aa974c2365ac435bb965 to your computer and use it in GitHub Desktop.
mdsumner
commented
Apr 14, 2026
Author
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment