python - python中的循环卷积

Question

有没有一种方法可以用Python在两个一维数组之间执行循环卷积，就像 Matlab 函数 cconv 一样？我尝试了 numpy.convolve 但它不一样，而且我找不到等效的。

Answer 1

您可以尝试此命令:scipy.ndimage.filters.convolve1d

您有一个名为mode的选项，您可以编写mode = wrap

这样，您就可以得到周期性边界条件作为卷积的填充

例如:

result = scipy.ndimage.convolve(image,kernel,mode='wrap')

import numpy as np
image = np.array([[0, 0, 0, 0],
                  [0, 0, 0, 1],
                  [0, 0, 0, 0]])
kernel = np.array([[1, 1, 1],
                    [1, 1, 1],
                    [1, 1, 1]])
from scipy.ndimage import convolve
convolve(image, kernel, mode='wrap')
array([[1, 0, 1, 1],   
       [1, 0, 1, 1],
       [1, 0, 1, 1]])