3D 局部平均值和使用 3D 卷积

Question

我是 python 新手，对 Matlab 更熟悉。如果我的问题不适合本论坛，请随时指出。

我正在尝试以非常快的速度制作局部平均值。就像我试图通过为每个新像素平均多个像素来减少图像中的像素数量一样，除了我是在 3D 中进行的。

想象一个 1000x1000x6 的阵列。我将这个数组分成多个 10x10x3 的小数组。然后我想计算所有这些小数组的平均值并将它们重新组合在一起以构建我的数组。

我在 Matlab 上做的方式是用convn(array,seed,'valid')，这是一个多维卷积函数。

在 python 中最简单的方法是什么？ 谢谢 手机版

Answer 1

我认为你能找到的最接近 convn 的是 SciPy 的 convolve。下面是例子

import numpy as np
from scipy.ndimage import convolve
M = np.random.random((1000, 1000, 6))
seed = np.ones((3, 3, 3)) * 0.1 / 27.
N = convolve(M, seed, mode='constant', cval=0)

mode='constant', cval=0 只是零填充。

不确定这是否是您需要的，但这是一个开始

文档：https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.ndimage.filters.convolve.html