Skip to content

Instantly share code, notes, and snippets.

@mdsumner
Created September 21, 2021 12:51
Show Gist options
  • Star 1 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save mdsumner/39e7901ee1d107e98c393d9332004d01 to your computer and use it in GitHub Desktop.
Save mdsumner/39e7901ee1d107e98c393d9332004d01 to your computer and use it in GitHub Desktop.
alos <- "/vsicurl/https://opentopography.s3.sdsc.edu/raster/AW3D30/AW3D30_global.vrt"
##eng_wales in Hugh's example (but I reduce the dimensions, for now)
g <- list(extent = c(82667.22, 655646.4, 5337.574, 657533.7), dimension = c(28575, 29006), projection = "EPSG:27700")
library(gdalio)
gdalio_set_default_grid(g)
source(system.file("raster_format/raster_format.codeR", package = "gdalio", mustWork = TRUE))
r <- gdalio_raster(alos)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment