如何计算python中大矩阵(平均滤波器)的每个坐标的有效局部平均值

Question

我有大矩阵,4000x4000我需要计算每个x的本地平均值11x11窗口,这个矩阵的y一般它必须是这样的

for x in range(4000)
  for y in range(4000)
    b[x,y]=mean(a[x-5:x+5,y-5:y+5]

但这会花很多时间这是一种更有效的方法吗？谢谢.

Answer 1

你基本上想要一个二维卷积.Scipy可以为你做到这一点:http://docs.scipy.org/doc/scipy/reference/generated/scipy.signal.convolve2d.html

实际上在SO上有一个类似的答案: 使用python和numpy的2d卷积