数字图像处理与机器视觉-基于MATLAB实现 第6章 灰度变换与滤波.pptVIP

6.5.4带阻滤波 图6-29 理想带阻滤波器 图6-30 巴斯沃特带阻滤波器 利用高斯带阻滤波消除噪声示例代码如下： I = imread(lena.bmp); imshow(I) [M,N] = size(I); P = I; for i = 1:M for j = 1:N P(i,j) = P(i,j)+20*sin(20*i)+20*sin(20*j);%添加周期噪声 end end figure,imshow(P); Diff_PI = double(P)-double(I);%噪声 figure,imshow((P-I),[]); IF = fftshift(fft2(I));%对原始图像作傅里叶变换，并将原点移至 中心 IFV = log(1+abs(IF)); PF = fftshift(fft2(P)); PFV = log(1+abs(PF)); 6.5.4带阻滤波 figure,imshow(IFV,[]) % figure,imshow(PFV,[]) freq = 90; %设置带阻滤波器中心频率-可调（要 会调参） width=5; %设置带阻滤波器频带宽度 ff = ones(M,N); for i = 1:M for j = 1:N ff(i,j) = 1-exp(-0.5*((((i-