建站
咨詢
售后
建站咨詢
新聞中心
這里有您想知道的互聯(lián)網(wǎng)營銷解決方案
Matlab低通濾波器的使用--頻域圖像增強(qiáng)實(shí)驗-創(chuàng)新互聯(lián)
//添加噪聲
C=imread('');
image=double(C);//轉(zhuǎn)換成double
[h,w]=size(image);
image2=uint8(zeros(h,w));
for i=1:h
for j=1:w
image2(i,j)=image(i,j);
if rand(1,1)>=0.99 //大于該值產(chǎn)生焦點(diǎn)
if rand(1,1)>=0.5 //大于0.5產(chǎn)生白點(diǎn)
image2(i,j)=(255);
else
image2(i,j)=(0);//否則產(chǎn)生黑點(diǎn)
end
end
end
end//高斯低通濾波器
f1=im2double(image2);
D0=40;
g1=fft2(f1);//傅立葉變換
g1=fftshift(g1);//轉(zhuǎn)換數(shù)據(jù)矩陣
[M1,N1]=size(g1);//如果圖像f(x,y)的尺寸為M×N，則對(-1)(x+y)f(x,y)進(jìn)行傅立葉變換后的頻率平面的原點(diǎn)在(M/2,N/2）
m1=fix(M1/2);
n1=fix(N1/2);
s1=zeros(M1,N1);
for i=1:M1
for j=1:N1
d=sqrt((i-m1)^2+(j-n1)^2);//D(u,v)為頻率平面從原點(diǎn)到點(diǎn)(u,v)的距離
h1=exp((-d^2)/(2*(D0^2)));//計算高斯低通濾波器傳遞函數(shù)
s1(i,j)=h1*g1(i,j);
end
end
y1=ifft2(ifftshift(s1));
y1=log(1+abs(y1));//巴特沃斯低通濾波器
f2=im2double(image2);
g2=fft2(f2);
g2=fftshift(g2);
n0=2;
[M2,N2]=size(g2);
m2=fix(M2/2);
n2=fix(N2/2);
s2=zeros(M2,N2);
for i=1:M2
for j=1:N2
d=sqrt((i-m2)^2+(j-n2)^2);
h2=1/(1+(d/D0)^(2*n0));//計算巴特沃斯低通濾波器傳遞函數(shù)
s2(i,j)=h2*g2(i,j);
end
end//繪畫
figure;
subplot(2,2,1);imshow(C);
title('原圖');
subplot(2,2,2);imshow(imamge2);
title('椒鹽噪聲后');
y2=ifft2(ifftshift(s2));
y2=log(1+abs(y2));
//如果圖像進(jìn)行傅里葉轉(zhuǎn)換后立即用imshow函數(shù)顯示，則在命令行可能會顯示：Warning: Displaying real part of complex input（警告: 顯示復(fù)數(shù)輸入項的實(shí)部）。這是因為經(jīng)過傅里葉變換后的圖像矩陣大多是復(fù)數(shù)矩陣，包含實(shí)部和虛部。
//此時如要顯示圖像則需要先用abs取復(fù)數(shù)矩陣的模，再進(jìn)行顯示
subplot(2,2,3);imshow(y1,[]);
//當(dāng)圖像是double類型時要使用imshow（I,[]）來根據(jù)數(shù)據(jù)矩陣的數(shù)值范圍自動設(shè)置灰度圖像顯示范圍
title('高斯低通濾波器處理后');
subplot(2,2,4);imshow(y2,[]);
title('巴特沃斯低通濾波器處理后');你是否還在尋找穩(wěn)定的海外服務(wù)器提供商？創(chuàng)新互聯(lián)www.cdcxhl.cn海外機(jī)房具備T級流量清洗系統(tǒng)配攻擊溯源，準(zhǔn)確流量調(diào)度確保服務(wù)器高可用性，企業(yè)級服務(wù)器適合批量采購，新人活動首月15元起，快前往官網(wǎng)查看詳情吧
網(wǎng)頁名稱：Matlab低通濾波器的使用--頻域圖像增強(qiáng)實(shí)驗-創(chuàng)新互聯(lián)
當(dāng)前網(wǎng)址：http://www.dlmjj.cn/article/jihhp.html
