使用scipy.weave.inline进行快速2D中值滤波

Question

我在2D中值滤波器(3x3窗口)中存在瓶颈,我在一组非常大的图像上使用,我想尝试优化它.我测试过scipy.ndimagemedian_filter,以及PIL,scipy.signal和scikits-image.然而,浏览SO我已经知道C中有一个快速的O(n)中值滤波器(恒定时间中的中值滤波,参见C中的滚动中值算法),我想知道我是否可以使用scipy在Python中实现它. weave.inline？有关替代路线的任何建议吗？

Answer 1

我不知道底层算法，但 scikits-image 有一个滚动中值滤波器。

否则，我建议用Cython（C/Python pidgin 语言）编写它。请务必查看有关使用 numpy 数组的卷积示例/教程。