12波段TIFF图像转Numpy数组出现Reshape维度错误的解决问询

错误根源 #

报错的核心逻辑是reshape操作要求前后数组的总元素数完全一致：
你原代码中imgxyb[:,:,:7]只取了影像的前7个波段，总元素数为 高*宽*7 = 921825 *7 = 6452775，但你指定reshape后的形状为(高*宽, 12)，需要总元素数为 921825 *12 = 11061900，两者数量不匹配就触发了ValueError。
你将切片参数改为:12后仍然报错，大概率是文件夹内存在波段数不等于12的TIFF文件，导致imgxyb.shape[2]数值和12不符，出现新的长度不匹配。

修复代码 #

方案1：使用全部波段生成一维数组（推荐） #

直接去掉不必要的波段切片，自动适配每个影像的实际波段数：

# 替换原构造img1d的代码行即可
img1d = imgxyb.reshape((imgxyb.shape[0] * imgxyb.shape[1], imgxyb.shape[2]))

方案2：固定使用前12个波段 #

如果你的业务逻辑要求必须用前12个波段，显式指定切片长度和reshape最后一维长度一致，同时加波段校验跳过不符合要求的文件：

for filepath in glob.iglob('./dengue3/*.tiff'):
    elhas_raster = rio.open(filepath)
    # 新增波段校验逻辑
    if elhas_raster.count < 12:
        print(f"跳过波段数不足12的文件：{filepath}")
        continue
    elhas_arr = elhas_raster.read()
    vmin, vmax = np.nanpercentile(elhas_arr, (5,95))
    imgxyb = np.empty((elhas_raster.height, elhas_raster.width, elhas_raster.count), elhas_raster.meta['dtype'])
    for band in range(imgxyb.shape[2]):
        imgxyb[:,:,band] = elhas_raster.read(band+1)
    # 显式指定波段数，保证前后长度匹配
    use_band_num = 12
    img1d = imgxyb[:,:,:use_band_num].reshape((imgxyb.shape[0] * imgxyb.shape[1], use_band_num))

内容的提问来源于stack exchange，提问作者user15653028