在Python中循环遍历多波段栅格中的每个像素

pip install gdal

import gdal

# 打开栅格文件
dataset = gdal.Open('path/to/raster.tif')

# 获取栅格的宽度和高度
width = dataset.RasterXSize
height = dataset.RasterYSize

# 获取栅格的波段数
band_count = dataset.RasterCount

# 循环遍历每个像素
for y in range(height):
    for x in range(width):
        # 读取每个像素的值
        pixel_values = []
        for band_index in range(band_count):
            band = dataset.GetRasterBand(band_index + 1)
            pixel_value = band.ReadAsArray(x, y, 1, 1)[0, 0]
            pixel_values.append(pixel_value)
        
        # 在这里进行对每个像素的处理
        # ...

# 关闭栅格文件
dataset = None