I created a script to separate individual bands from two multiband .tif files in order to combine them in a "custom" multiband .tif. I also clip them into 3 smaller study areas. When
I display them in QGIS and rasterio and they are completely different. In QGIS, all dark colours are gone. I think this is because I missed the alpha band. Nevertheless, the all other raster values are like the originals.