《巴特沃斯(Butterworth)低通濾波器對受噪聲干擾的圖像進(jìn)行平滑處理》由會員分享,可在線閱讀,更多相關(guān)《巴特沃斯(Butterworth)低通濾波器對受噪聲干擾的圖像進(jìn)行平滑處理(4頁珍藏版)》請?jiān)谘b配圖網(wǎng)上搜索。
1、%頻域增強(qiáng)(巴特沃斯原型)
%利用巴特沃斯(Butterworth )低通濾波器對受噪聲干擾的圖像進(jìn)行平滑處理
clc;
clear;
J1=imread(le na.gif); subplot(3,2,1);
imshow(JI);
title(原圖);
f=double(J1);
g=fft2(f); % 傅立葉變換
g=fftshift(g); % 轉(zhuǎn)換數(shù)據(jù)矩陣
subplot(3,2,2);
x=0:1:255;
y=0:1:255;
[x,y]=meshgrid(x,y);
z=log(abs(g)); % 取幅度
mesh(z); %以三維坐標(biāo)顯示該圖
2、像頻譜圖 title( 原圖頻譜);
[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
計(jì)算低通濾波器傳遞函數(shù)
d=sqrt((i-m)A2+(j-門)八2);
h=1/(1+0.414*(d/d0)A(2* nn)); % result(i,j)=h*g(i,j);
end end
subplot(3,2,4);
x=0:1:255;
y=0:1:255;
[x,y]=meshgrid(x,y);
z=log(
3、abs(result)); % 取幅度
mesh(z); %以三維坐標(biāo)顯示該圖像頻譜圖
title( 低通濾波后的頻譜);
result=ifftshift(result);
J2=ifft2(result);
J3=ui nt8(abs(J2));
subplot(3,2,3);
imshow(J3);
title( 低通濾波后的圖像);
%利用巴特沃斯(Butterworth )高通濾波器對受噪聲干擾的圖像進(jìn)行平滑處理
nn=2; % 二階巴特沃斯(Butterworth) 高通濾波器
d0=5;
m=fix(M/2);
n=fix(N/2);
for i=1:
4、M
for j=1:N
d=sqrt((i-m)A2+(j-門)八2);
if (d==0)
h=0;
else
h=1/(1+0.414*(d0/d)A(2* nn));% 計(jì)算傳遞函數(shù)
end
result(i,j)=h*g(i,j);
end
end
subplot(3,2,7);
x=0:1:255;
y=0:1:255;
[x,y]=meshgrid(x,y);
z=log(abs(result)); % 取幅度
mesh(z); %以三維坐標(biāo)顯示該圖像頻譜圖 title( 高通濾波后的頻譜);
result=ifftshift(result);
J2=ifft2(result);
J3=ui nt8(abs(J2));
subplot(3,2,5);
imshow(J3);
title( 高通濾波后的圖像);