如何就地应用卷积（在 Python 中）

Question

我需要就地对图像执行卷积，就地我的意思是当结构元素应用于不同的像素时，我希望前面步骤的结果覆盖图像。 把它放在上下文中，这在 Gauss-Seidel 迭代中很有用。

我目前正在使用scipy.ndimage.convolve1d ，它显然没有像我描述的那样进行就地卷积。 我知道如何编写一个 kernel 来做到这一点（基本上是一个 for 循环，您在迭代像素时覆盖现有元素）。

有没有办法使用 python 库或其他东西来实现这一点？

Answer 1

如果我理解正确，您正在尝试运行 2D 无限滤波器响应滤波器。

如果您的过滤器是可分离的，您可以通过应用lfilter来做到这一点，每个轴一次。 此外，此方法将转移 output 可能更适合使用filtfilt 。