将多个只有单一波段的栅格数据融合成一个多波段的栅格数据，最终想要得到如下所示的栅格数据。尝试了用stack,brick和merge都不符合。


#class      : RasterLayer
#band       : 1  (of  10  bands)
#dimensions : 1000, 1000, 1e+06  (nrow, ncol, ncell)
#resolution : 10, 10  (x, y)
#extent     : 356780, 366780, 341500, 351500  (xmin, xmax, ymin, ymax)
#crs        : +proj=utm +zone=33 +datum=WGS84 +units=m +no_defs
#source     : S2A_T33NUD_20180104_Subset
#names      : S2A_T33NUD_20180104_Subset
#values     : -32768, 32767  (min, max)