一起添加多个单波段图像

Question

我想将31张单波段图像加在一起以制作一张图像。我设法使用rasterio在python中打开它们。图像添加步骤无效。

我试图将这些图像添加为numpy数组。任何wat numpy.add都有两个参数，我不确定如何在for循环中使用它。

for i in range (214,220):
    data = rasterio.open("et"+str(i)+'.img')
    #show(data, cmap='PiYG')
    print 'image {val} opened'.format(val=i)
    band = data.read(1)
    et = np.add(band, band)

我在q gis栅格计算器中尝试过此操作。最终结果是不同的。