//添加噪聲
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)));//計(jì)算高斯低通濾波器傳遞函數(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));//計(jì)算巴特沃斯低通濾波器傳遞函數(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ù)顯示,則在命令行可能會(huì)顯示:Warning: Displaying real part of complex input(警告: 顯示復(fù)數(shù)輸入項(xiàng)的實(shí)部)。這是因?yàn)榻?jīng)過傅里葉變換后的圖像矩陣大多是復(fù)數(shù)矩陣,包含實(shí)部和虛部。
//此時(shí)如要顯示圖像則需要先用abs取復(fù)數(shù)矩陣的模,再進(jìn)行顯示
subplot(2,2,3);imshow(y1,[]);
//當(dāng)圖像是double類型時(shí)要使用imshow(I,[])來根據(jù)數(shù)據(jù)矩陣的數(shù)值范圍自動(dòng)設(shè)置灰度圖像顯示范圍
title('高斯低通濾波器處理后');
subplot(2,2,4);imshow(y2,[]);
title('巴特沃斯低通濾波器處理后');
你是否還在尋找穩(wěn)定的海外服務(wù)器提供商?創(chuàng)新互聯(lián)www.cdcxhl.cn海外機(jī)房具備T級(jí)流量清洗系統(tǒng)配攻擊溯源,準(zhǔn)確流量調(diào)度確保服務(wù)器高可用性,企業(yè)級(jí)服務(wù)器適合批量采購,新人活動(dòng)首月15元起,快前往官網(wǎng)查看詳情吧
分享文章:Matlab低通濾波器的使用--頻域圖像增強(qiáng)實(shí)驗(yàn)-創(chuàng)新互聯(lián)
文章起源:http://jinyejixie.com/article48/jihhp.html
成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供Google、品牌網(wǎng)站制作、靜態(tài)網(wǎng)站、微信小程序、網(wǎng)站設(shè)計(jì)公司、商城網(wǎng)站
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請盡快告知,我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場,如需處理請聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時(shí)需注明來源: 創(chuàng)新互聯(lián)
猜你還喜歡下面的內(nèi)容