使用numpy计算栅格图层的逐像素均值

Question

由于两个栅格（raster1和raster2）重叠在一起，我想通过计算每个重叠像素的平均值来创建新的栅格；即，得到的新栅格计算如下：

new = [[mean(1,3), mean(1,3), mean(1,3), mean(1,3), mean(1,3)],[mean(2,4),mean(2,4),mean(2,4),mean(2,4),mean(2,4)]]


import numpy as np    
raster1 = np.array([[1,1,1,1,1],[2,2,2,2,2]])
raster2 = np.array([[3,3,3,3,3],[4,4,4,4,4]])

new = np.mean(raster1,raster2,axis=1)
print (new.tolist())

有什么问题吗？

- 2964502

Answer 1

也许我误解了你的意思，但是你想要什么？

raster = (raster1 + raster2) / 2

实际上，在这种情况下，您甚至不需要使用np.mean，只需使用矩阵运算即可。 np.mean用于计算特定轴上单个矩阵的平均值，因此情况不同。

- Ray

Answer 2

应该是：

它应该是

new = np.mean([raster1,raster2],axis=1)

用括号。实际上我猜应该是

。

new = np.mean([raster1,raster2],axis=0)

np.mean的第一个参数应该是整个数组，例如参考http://wiki.scipy.org/Numpy_Example_List_With_Doc#mean。

- oschoudhury