data-raw/DATASET.R

library(terra)

x <- rast("~/Downloads/N53W118.hgt")
x <- project(x, "epsg:3402", method = "bilinear")
x <- crop(x, ext(c(310000, 350000, 5880000, 5920000)))

srtm <- resample(x, rast(ext(x), resolution = c(100, 100)), method = "bilinear")
writeRaster(srtm, "inst/extdata/srtm.tif", overwrite = TRUE)
stevenpawley/sagacmd documentation built on Sept. 12, 2024, 2:18 p.m.