import numpy as np
from scipy import signal, misc
import matplotlib.pyplot as plt
image = 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()

import numpy as np
from scipy import signal, misc
import matplotlib.pyplot as plt
image = misc.ascent()
w = signal.gaussian(50,10.0)
image_new = signal.sepfir2d(image, w, w)
plt.figure()
plt.imshow(image_new)#显示滤波后的图像
plt.gray()
plt.title('Filtered image')
plt.show()