I followed the instructions here. My input layers are 7 .tif images, each representing a spectral band of a Landsat image. I'm attempting to extract these bands using GDAL:
raster_dataset = gdal.Open(raster_data_path, gdal.GA_ReadOnly) bands_data =  for b in range(1, raster_dataset.RasterCount+1): band = raster_dataset.GetRasterBand(b) bands_data.append(band.ReadAsArray())
Issue is, raster_dataset.RasterCount == 1 and I'm getting arrays full of zeros.