【R】stack raster layers

using R and library "raster", "rgdal"

library(raster)

library(rgdal)

# read layers
band1 <- stack("E:/tmp_data/ALOS/bbb1.tif")
band2 <- stack("E:/tmp_data/ALOS/bbb2.tif")
band3 <- stack("E:/tmp_data/ALOS/bbb3.tif")
band4 <- stack("E:/tmp_data/ALOS/bbb4.tif")

# stack layers
stack.band <- stack(band1, band2, band3, band4)

# output the stacked image
writeRaster(stack.band, "E:/tmp_data/ALOS/tsushima.tif", overwrite=TRUE, datatype="INT2U")

# plot
plotRGB(bands_stack, r=1, g=2, b=3)