信号滤波具有径向对称带通的图像

Question

我的任务是对图像执行“任意”傅里叶滤波-巴特沃斯、贝塞尔和其他类型的低/高/带通和抑制滤波器。天真的实现：

# b, a are the result of scipy.signal.bessel or .butter, etc
# lfilter could be filtfilt too
filt_y = signal.lfilter(b, a, img, axis=0)
filt_both = signal.lfilter(b, a, filt_y, axis=1)

将沿穿过k空间原点的十字传递效果，但保持拐角区域不变。有没有办法仍然使用scipy.signal来实现这一点，但是带内和带外区域都有环形空间？我不想要像cv2那样的DoG或其他类型的基于内核的过滤器，而是一个在scipy.signal中实现的因果过滤器

原则上，所需要的只是能够为每个样本传入频率值np.sqrt（fxx**2+fyy**2），但我不知道如何在信号API中做到这一点