//添加噪聲
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)生焦點
            if rand(1,1)>=0.5 //大于0.5產(chǎn)生白點
            image2(i,j)=(255);
            else
            image2(i,j)=(0);//否則產(chǎ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)進行傅立葉變換后的頻率平面的原點在(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)為頻率平面從原點到點(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));
//如果圖像進行傅里葉轉(zhuǎn)換后立即用imshow函數(shù)顯示，則在命令行可能會顯示：Warning: Displaying real part of complex input（警告: 顯示復數(shù)輸入項的實部）。這是因為經(jīng)過傅里葉變換后的圖像矩陣大多是復數(shù)矩陣，包含實部和虛部。
//此時如要顯示圖像則需要先用abs取復數(shù)矩陣的模，再進行顯示
subplot(2,2,3);imshow(y1,[]);
//當圖像是double類型時要使用imshow（I,[]）來根據(jù)數(shù)據(jù)矩陣的數(shù)值范圍自動設(shè)置灰度圖像顯示范圍
title('高斯低通濾波器處理后');
subplot(2,2,4);imshow(y2,[]);
title('巴特沃斯低通濾波器處理后');

你是否還在尋找穩(wěn)定的海外服務器提供商？創(chuàng)新互聯(lián)www.cdcxhl.cn海外機房具備T級流量清洗系統(tǒng)配攻擊溯源，準確流量調(diào)度確保服務器高可用性，企業(yè)級服務器適合批量采購，新人活動首月15元起，快前往官網(wǎng)查看詳情吧

創(chuàng)新互聯(lián)專注于海淀企業(yè)網(wǎng)站建設(shè),成都響應式網(wǎng)站建設(shè),購物商城網(wǎng)站建設(shè)。海淀網(wǎng)站建設(shè)公司,為海淀等地區(qū)提供建站服務。全流程定制網(wǎng)站制作，專業(yè)設(shè)計，全程項目跟蹤，創(chuàng)新互聯(lián)專業(yè)和態(tài)度為您提供的服務

分享文章：Matlab低通濾波器的使用--頻域圖像增強實驗-創(chuàng)新互聯(lián)
文章起源：http://chinadenli.net/article48/jihhp.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián)，為您提供Google、品牌網(wǎng)站制作、靜態(tài)網(wǎng)站、微信小程序、網(wǎng)站設(shè)計公司、商城網(wǎng)站

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主，如果涉及侵權(quán)請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場，如需處理請聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載，或轉(zhuǎn)載時需注明來源： 創(chuàng)新互聯(lián)