理论基础知识,高通部分已经讲过,这里不重复了。本部分做低通滤波,保留细节。即频谱中间部分保留,与高通相反。
import cv2
import numpy as np
import matplotlib.pyplot as plt
o=cv2.imread('image\\lena256.bmp',0) #读入图片
dft=cv2.dft(np.float32(o),flags=cv2.DFT_COMPLEX_OUTPUT)
dshift=np.fft.fftshift(dft) #傅里叶变换
rows,cols=o.shape #行 列
crow,ccol=int(rows/2),int(cols/2) #中心
mask=np.zeros((rows,cols,2),np.int8) #生成掩膜,2个通道,256位
mask1=np.zeros((rows,cols,2),np.int8) #生成掩膜,2个通道,256位
mask[crow-30:crow+30,ccol-30:ccol+30]=1 #中心区域设置为1,大小30*30
mask1[crow-10:crow+10,ccol-10:ccol+10]=1 #中心区域设置为1,大小10*10
md=dshift*mask
ishift=np.fft.ifftshift(md)
io=cv2.idft(ishift)
io=cv2.magnitude(io[:,:,0],io[:,:,1]) #逆傅里叶变换
md1=dshift*mask1
ishift1=np.f