gdal - 带读取数组如何只返回 0，我该怎么办？

Question

我是GDAL的新手。我在分类中使用 landsat 图像，但我应用这些代码行带数据只给我 0，我不明白：

raster_dataset = gdal.Open(‘LC81850552017065LGN00_B3.tif’, gdal.GA_ReadOnly) 
geo_transform = raster_dataset.GetGeoTransform() 
proj = raster_dataset.GetProjectionRef() 
bands_data = [] 

for b in range(1, raster_dataset.RasterCount+1): 
    band = raster_dataset.GetRasterBand(b) 
    bands_data.append(band.ReadAsArray())
    bands_data = np.dstack(bands_data)

Answer 1

您正在bands_data从一个列表重新分配给一个 numpy 数组，这给您带来了麻烦。np.dstack在读入所有波段后（循环完成后）执行操作。

请参阅numpy 文档以更好地了解其dstack工作原理。