Maybe I am not asking the question with the appropriate English words, but I would like to basically concatenate multi bands rasters (two or more) into a single raster with all the bands intact. The rasters have the same number of bands (4), and have all the same extent.
I can't find the proper tool in QGIS or either in GDAL. I first thought merge could do the job but it ends up with 4 bands having the average value of the 3rd band of each file for instance (for the 3rd band).
I can either use QGIS, GDAL with python, or directly GDAL in the terminal - all solutions are welcome.