thermal image

If any channel left, a black and white raster is generated.
thermal<-ToRGB(my.raster$band$tirs1, my.raster$band$tirs2)