scipy的二维图像卷积运算

一 代码

1. import numpy as np
2. from scipy import signal, misc
3. import matplotlib.pyplot as plt
5. image = misc.ascent()#二维图像数组，lena图像
6. w = np.zeros((50,50))#全0二维数组，卷积核
7. w[0][0]=1.0#修改参数，调整滤波器
8. w[49][25]=1.0#可以根据需要调整
9. image_new = signal.fftconvolve(image, w)#使用FFT算法进行卷积
11. plt.figure()
12. plt.imshow(image_new)#显示滤波后的图像
13. plt.gray()
14. plt.title('Filtered image')
15. plt.show()

二 运行结果