raster = np.zeros((bands, rows, cols)) for i in range(bands): band = dataset.GetRasterBand(i+1) raster[i,:,:] = band.ReadAsArray()

这段代码是用来读取一个多波段遥感影像的数据，并将其存储在一个三维的 numpy 数组中。其中，`dataset` 是一个 GDAL 数据集对象，`bands` 是波段数，`rows` 和 `cols` 分别是影像的行数和列数。首先，我们创建一个形状为 `(bands, rows, cols)` 的全 0 数组 `raster`。然后，对于每个波段，我们获取其对应的 GDAL 波段对象 `band`，并使用 `ReadAsArray()` 方法读取其数据，将其存储在 `raster` 数组的对应位置。最终，`raster` 数组就包含了整个遥感影像的数据。