Search code examples
rrasterterrarasterstats

Set raster metadata (band names) in R on disk

I have a large raster on file. It can be read like so:

my_rast <- rast("my_large_file.tif")
names(my_rast)
[1] "value_a"

Using R or command-line GDAL I would like to change the name of the layer to something different, without totally re-writing the whole raster. How can I do this?

Note: I see there is a command in the Python rasterio library to do this, but I don't see a similar one in R.

Solution

  • You can now use the update(x, names=TRUE) to write the changed layer names to disk (for file formats that support layer names, such as GTiff).

    Example data

    library(terra)
#terra 1.6.29
s <- rast(system.file("ex/logo.tif", package="terra"))   
names(s)
#[1] "red"   "green" "blue" 
fname <- paste0(tempfile(), ".tif")
x <- writeRaster(s, fname)

    Change the names, and write them to the file.

    names(x) <- c("A", "B", "C")
update(x, names=TRUE)

    Show that it worked

    r <- rast(fname)
names(r)
#[1] "A" "B" "C"