24、栅格数据处理：属性表与虚拟栅格格式

栅格数据处理：属性表与虚拟栅格格式

1. 栅格属性表

在处理栅格数据时，属性表是一个重要的概念。整数型栅格数据集可以拥有属性表，但与矢量属性表相比，其字段数量通常较少。栅格属性表中的每条记录对应一个特定的像素值，而非单个要素。例如，所有值为 56 的像素在属性表中共享同一条记录，因为每个像素并不代表单个要素，而具有相同值的多个像素应代表相同的事物，如特定颜色、高程、土地利用分类等。

对于许多栅格数据而言，属性表可能并无实际意义。例如，在航空照片中，很难为不同的像素值附加有意义的属性。实际上，栅格属性表在处理分类数据（如土地覆盖或土壤类型）时最为有用，因为此时需要包含每个类别的相关信息。

以下是一个为分类高程栅格添加属性表的示例：

import os
from osgeo import gdal

os.chdir(r'D:\osgeopy-data\Switzerland')
ds = gdal.Open('dem_class2.tif')
band = ds.GetRasterBand(1)
band.SetNoDataValue(-1)
rat = gdal.RasterAttributeTable()
rat.CreateColumn(
    'Value', gdal.GFT_Integer, gdal.GFU_Name)
rat.CreateColumn(
    'Count', gdal.GFT_Integer, gdal.GFU_PixelCount)
rat.CreateColumn(
    'Elevation', gdal.GFT_String, gdal.GFU_Gene