do we cover the newer GDAL Python Raster API introduced around GDAL 3.8–3.11
with that you can much better access band level aspects
from osgeo import gdal
ds = gdal.Open("image.tif", gdal.GA_Update)
band = ds.GetRasterBand(1)
band.SetDescription("Red")
band.SetUnitType("reflectance")
band.SetMetadataItem("wavelength", "665")
do we cover the newer GDAL Python Raster API introduced around GDAL 3.8–3.11
with that you can much better access band level aspects