一 代码
import numpy as npfrom scipy import signal, miscimport matplotlib.pyplot as pltimage = misc.ascent()#二维图像数组,lena图像w = np.zeros((50,50))#全0二维数组,卷积核w[0][0]=1.0#修改参数,调整滤波器w[49][25]=1.0#可以根据需要调整image_new = signal.fftconvolve(image, w)#使用FFT算法进行卷积plt.figure()plt.imshow(image_new)#显示滤波后的图像plt.gray()plt.title('Filtered image')plt.show()
二 运行结果