Combinación de bandas en R

1. Carga las librerías necesarias para el proceso: rgdal y raster.

• raster: para leer, manipular y escribir imágenes raster.
• rgdal: para leer y escribir imágenes en diferentes formatos vectoriales y raster.

library(rgdal)
library(raster)

2. Utiliza la función setwd() para establecer el directorio de trabajo donde se encuentran las imágenes de bandas individuales.

setwd("Ruta")

3. Utiliza la función raster() para cargar cada una de las bandas individuales en una variable distinta.

banda1 <- raster("LC08_L1TP_010063_20220614_20220617_02_T1_B5.tif")
banda2 <- raster("LC08_L1TP_010063_20220614_20220617_02_T1_B6.tif")
banda3 <- raster("LC08_L1TP_010063_20220614_20220617_02_T1_B2.tif")

4. Crea una lista con las bandas individuales cargadas en las variables anteriores (banda1, banda2, banda3).

lista_bandas <- list(banda1, banda2, banda3)

5. Utiliza la función stack() para combinar las bandas individuales en una sola imagen compuesta.

imagen_compuesta <- stack(lista_bandas)

6. Utiliza la función writeRaster() para guardar la imagen compuesta en el directorio de trabajo establecido anteriormente.

writeRaster(imagen_compuesta, "Nombre del archivo de salida.tif")