一次从栅格中提取所有波段的最快方法（python / gdal）

Question

我有一些带有大量波段的高光谱图像，我想对其进行分析。我的脚本需要能够一次访问所有波段。

目前，我正在通过以下方式实现这一目标：

bands = np.asarray([dataset.GetRasterBand(n+1) for n in range(dataset.RasterCount)])

这很好用，但似乎这一步在我的处理工作流程中占用了大量时间，我怀疑有更好的方法。另外，我认为以这种方式使用numpy列表推导是不好的做法（？）。

numpy或gdal是否有任何内置方法可以加快速度？

Answer 1

在GDAL中，波段和波段数据之间存在差异。假设你想要后者，只需使用：

data = dataset.ReadAsArray()