Получение RGB из спутниковых данных

Question

    from osgeo import gdal
    import rasterio as rio
    from rasterio.plot import show
    import matplotlib.pyplot as plt
    import os
    import gdal
    import numpy as np

    filepath = ('C:/Users/kelvi/Desktop/Project research/Project files/satellitedata1/MER_RR__1PRBCM20020429_081001_000003862005_00293_00848_0001.N1')

raster = gdal.Open(filepath)

rgb = np.stack((raster.GetRasterBand(b).ReadAsArray() for b in (1,2,3)))

band1 = raster.GetRasterBand(1).ReadAsArray()
band2 = raster.GetRasterBand(2).ReadAsArray()
band3 = raster.GetRasterBand(3).ReadAsArray()
band4 = raster.GetRasterBand(4).ReadAsArray()
band5 = raster.GetRasterBand(5).ReadAsArray()

Как просмотреть RGB? Это 17 полос. Но я показал здесь 5 полос. Мне нужна помощь, пожалуйста. Я больше разбираюсь в программировании.