維納濾波實現(xiàn)圖像恢復(fù).doc

《維納濾波實現(xiàn)圖像恢復(fù).doc》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、第六次MATLAB大作業(yè)一．實驗題目用式5.6.3對圖像進(jìn)行模糊處理,用逆濾波和維納濾波恢復(fù)圖像二．程序清單clear;I=imread('rice.tif');figure;subplot(2,2,1);imshow(I);title('原圖像');[m,n]=size(I);F=fftshift(fft2(I));k=0.0025;foru=1:mforv=1:nH(u,v)=exp((-k)*(((u-m/2)^2+(v-n/2)^2)^(5/6)));endendG=F.*H;I0=real(ifft2(fftshift(G)));I1=imnoise(uint8(I0),'ga

2、ussian',0,0.001)subplot(2,2,2);imshow(uint8(I1));title('模糊退化且添加高斯噪聲的圖像');F0=fftshift(fft2(I1));F1=F0./H;I2=ifft2(fftshift(F1));subplot(2,2,3);imshow(uint8(I2));title('全逆濾波復(fù)原圖');K=0.1;foru=1:mforv=1:nH(u,v)=exp(-k*(((u-m/2)^2+(v-n/2)^2)^(5/6)));H0(u,v)=(abs(H(u,v)))^2;H1(u,v)=H0(u,v)/(H(u,v)*(H0(u

3、,v)+K));endendF2=H1.*F0;I3=ifft2(fftshift(F2));subplot(2,2,4);imshow(uint8(I3));title('維納濾波復(fù)原圖');三．運行結(jié)果四．結(jié)果分析如運行結(jié)果所示，其中退化函數(shù)中常數(shù)k=0.0025（劇烈湍流）時，全逆濾波的結(jié)果不可用了，而維納濾波的結(jié)果比較接近源圖像；而且經(jīng)過模糊退化加高斯噪聲后，維納濾波能消除一些噪聲影響，復(fù)原圖的效果也比較好。改變常數(shù)k=0.00025(輕微湍流)時，全逆濾波結(jié)果有所改善。但總體來看維納濾波的效果更好些。五．作業(yè)時遇到的問題最初寫的程序中圖像(未添加噪聲)經(jīng)傅立葉變換后沒有經(jīng)過中心

4、化，直接運用所給退化函數(shù)經(jīng)過濾波。運行結(jié)果是逆濾波復(fù)原圖基本和原圖一樣，而維納濾波結(jié)果和經(jīng)模糊退化的圖像一樣，即逆濾波的結(jié)果更好，與預(yù)期效果完全相反，目前還沒想明白原因。