python - scipy.signal.convolve2d 计算什么?

标签 python scipy

<分区>

我目前对

的输出有点困惑
#!/usr/bin/env python

import scipy.signal

image = [[1, 2, 3, 4, 5, 6, 7],
         [8, 9, 10, 11, 12, 13, 14],
         [15, 16, 17, 18, 19, 20, 21],
         [22, 23, 24, 25, 26, 27, 28],
         [29, 30, 31, 32, 33, 34, 35],
         [36, 37, 38, 39, 40, 41, 42],
         [43, 44, 45, 46, 47, 48, 49]]

filter_kernel = [[-1, 1, -1],
                 [-2, 3, 1],
                 [2, -6, 0]]

res = scipy.signal.convolve2d(image, filter_kernel,
                              mode='same', boundary='fill', fillvalue=0)
print(res)

原来是

[[  -2   -8   -7   -6   -5   -4   28]
 [   3   -7  -10  -13  -16  -19   14]
 [ -18  -28  -31  -34  -37  -40    0]
 [ -39  -49  -52  -55  -58  -61  -14]
 [ -60  -70  -73  -76  -79  -82  -28]
 [ -81  -91  -94  -97 -100 -103  -42]
 [-101  -61  -63  -65  -67  -69  -57]]

我希望左上角的元素是 3*1 + 1*2 + (-6) *8 + 0*9 = -43(省略填充的零)。

我认为这会通过在左/右和上/下添加一个 0 将矩阵图像\in R^{7x7} 扩展到 R^{9x9}。然后我认为 filter_kernel 可以通过在 image 上“滑动”来计算。在每个位置,图像中的数字与内核中的数字逐点相乘。将九个产品相加并写入res

但是,它是-2。显然,会发生一些不同的事情。

最佳答案

卷积反转它所处理的函数之一的方向。检查 definition on Wikipedia :一个函数用 τ 参数化,另一个函数用 -τ 参数化。这同样适用于二维卷积。

您需要镜像内核以获得预期的结果:

filter_kernel = [[0, -6, 2],
                 [1, 3, -2],
                 [-1, 1, -1]]

res = scipy.signal.convolve2d(image, filter_kernel,
                              mode='same', boundary='fill', fillvalue=0)
print(res[0, 0])
# -43

关于python - scipy.signal.convolve2d 计算什么?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/41613155/

相关文章:

python - Django Wsgi ModuleNotFoundError : No module named 'project_name'

python - 有没有一种简单的方法来判断文件指针所在的行号?

python - 使用 scipy.stats.hypergeom 从 SciPy 中的牌组绘制卡片

python - 对数缩放(应用函数?)稀疏矩阵

python - 使用python发送unsigned char

python - Flask-SQLAlchemy 数据类型

python - 如何组合多个 NumPy bool 数组?

python - 将密集矩阵从文件直接读取到稀疏 numpy 数组中?

pandas - 巨大的稀疏数据帧到 scipy 稀疏矩阵,无需密集变换

python - 将坐标元组信息转换为numpy数组