A function to get pixel values from a GDAL landcover raster

landcover-raster-attributes.py

import collections

from osgeo import gdal

RasterStats = collections.namedtuple('RasterStats', ['min', 'max'])

def get_raster_values_from_attr_table(filepath, band_index):
    gdal_dataset = gdal.OpenEx(filepath, gdal.OF_RASTER)
    band = gdal_dataset.GetRasterBand(1)
    attr_table = band.GetDefaultRAT()
    if attr_table:
        value_column_index = 0
        for col_index in range(attr_table.GetColumnCount()):
            if attr_table.GetNameOfCol(col_index).lower() == 'value':
                value_column_index = col_index
                break

        raster_values = set()
        for row_index in attr_table.GetRowCount():
            raster_values.add(attr_table.GetValueAsInt(
                row_index, value_column_index))
    return raster_values