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

import rasterio
from rasterio.plot import show

# 打开栅格文件
with rasterio.open('path_to_your_raster_file.tif') as src:
    # 获取栅格的宽度和高度
    width = src.width
    height = src.height

    # 遍历每个像素
    for y in range(height):
        for x in range(width):
            # 读取每个波段的像素值
            pixel_values = []
            for band in src.read_masks():
                pixel_value = src.read(band, window=((y, y+1), (x, x+1))).item()
                pixel_values.append(pixel_value)
            
            # 处理像素值
            # 在这里添加你的处理逻辑
            print(f"Pixel at ({x}, {y}): {pixel_values}")

# 可选：显示栅格图像（仅用于调试）
# show(src)

pip install rasterio