I have a collection of multi-band Sentinel-2 rasters. From these, I apply the Extract Band Raster Function. On the resulting mosaic dataset, I need to retain the initial band names. For example, if I extract B4 (red) and B8 (ir) bands, in the resulting mosaic dataset, I will have B1 and B2 as names, but they can be easily misinterpreted as B1 (coastal aerosol) and B2 (blue) bands.
Is there a way in ArcMap or Python to retain the original band names?