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)
