# \donttest{
## Prepare file paths
td <- tempdir()
in_raster <- file.path(td, "europe.tif")
out_raster <- file.path(td, "europe_small.tif")
file.copy(system.file("extdata/europe.tif", package = "gdalUtilities"),
to = td)
## Shrink a tiff by 50% in both x and y dimensions
gdal_translate(in_raster, out_raster, outsize = c("50%","50%"))
## Check that it worked
if(require(terra)) {
r1 <- rast(in_raster)
r1[is.na(r1)] <- 0
r1 <- as.factor(r1)
rat <- levels(r1)[[1]]
rat[["landcover"]] <- c("water", "land")
levels(r1) <- rat
r2 <- rast(out_raster)
r2[is.na(r2)] <- 0
r2 <- as.factor(r2)
rat <- levels(r2)[[1]]
rat[["landcover"]] <- c("water", "land")
levels(r2) <- rat
op <- par(mfcol = c(1, 2))
plot(r1, col = c("lightblue", "brown"), legend = FALSE)
plot(r2, col = c("lightblue", "brown"), legend = FALSE)
par(op) ## Reset pre-existing parameters
}
# }
Run the code above in your browser using DataLab