乐队如何读取数组只是返回我0，我该怎么办？

Question

我是GDAL的新手。我使用的分类TM遥感影像，但我将这些代码行带数据只是给我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的工作原理。