带读数据如何只返回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 docs以更好地了解dstack的工作原理。