如何在python中使用栅格数据集进行循环

import gdal

dataset = gdal.Open("path_to_raster_file")

rows = dataset.RasterYSize   # 栅格数据集的行数
cols = dataset.RasterXSize   # 栅格数据集的列数
bands = dataset.RasterCount  # 栅格数据集的波段数

for band in range(1, bands+1):
    band_data = dataset.GetRasterBand(band)
    for row in range(rows):
        scanline = band_data.ReadAsArray(0, row, cols, 1)
        for col in range(cols):
            pixel_value = scanline[0][col]
            # 在这里进行相应的操作

dataset = None