给你一个完整的例子,希望对你有帮助(频域滤波基本可以参照这个例子)
I = imread('rice.png');
J1=imnoise(I,'salt & pepper'); % 叠加椒盐噪声
f=double(J1); % 数据类型转换,MATLAB不支持图像的无符号整型的计算
g=fft2(f); % 傅立叶变换
g=fftshift(g); % 转换数据矩阵
[M,N]=size(g);
nn=2; % 二阶巴特沃斯(Butterworth)低通滤波器
d0=20; %此参数对处理结果影响很大
m=fix(M/2); n=fix(N/2);
for i=1:M
for j=1:N
d=sqrt((i-m)^2+(j-n)^2);
h=1/(1+0.414*(d/d0)^(2*nn)); % 计算低通滤波器传递函数
result(i,j)=h*g(i,j);
end
end
result=ifftshift(result);
J2=ifft2(result);
J3=uint8(real(J2));
figure % 对比显示结果
subplot(131),imshow(I),title('原图')
subplot(132),imshow(J1),title('椒盐噪声干扰后')
subplot(133),imshow(J3),title('低通滤波结果')
[ 本帖最后由 花如月 于 2008-6-3 11:44 编辑 ] |