Использование scipy.weave.inline для быстрой 2D-медианной фильтрации

У меня есть узкое место в 2D-медианном фильтре (окно 3x3), который я использую на очень большом наборе изображений, и я бы хотел попробовать и оптимизировать его. Я протестировал scipy.ndimagemedian_filter, а также PIL, scipy.signalи scikits-image. Однако, просматривая SO, я узнал, что в C существует быстрый медианный фильтр O(n) (медианная фильтрация в постоянное время, см. Алгоритм скользящей медианы в C), и мне стало интересно, могу ли я реализовать это в Python, используя scipy.weave.inline? Любые предложения по альтернативному маршруту?