資源描述：
《圖像復原維納濾波邊緣提取》由會員上傳分享，免費在線閱讀，更多相關內(nèi)容在教育資源-天天文庫。
1、圖像復原1).針對圖像水平運動所引起的圖像退化，給出連續(xù)函數(shù)退化模型的系統(tǒng)傳遞函數(shù)H（u,v）2).編程用維納濾波實現(xiàn)對運動模糊圖像的復原。解：首先將設計一個水平運動模糊濾波器對無噪的Lenna.bmp進行濾波，得到模糊后的圖像，然后進行維納濾波，也就是逆濾波（要消除H的零點，避免出現(xiàn)0/0）。程序如下：A=imread('Lenna.bmp');F=fft2(A);F=fftshift(F);N=512;u=1:N/2;T=1;a=0.05;%選擇水平運動參數(shù)x=pi*a*u;H=zeros(N,N);H(1,1+N/2:N)=(T./x).*sin(x).*
2、exp(-i*x);%計算濾波系數(shù)H(1,1:N/2)=fliplr(conj(H(1,1+N/2:N)));%利用共軛對稱性forn=2:NH(n,:)=H(1,:);endfigure,mesh(abs(H))%濾波器幅度特性G=F.*H;G=fftshift(G);AA=uint8(ifft2(G));figure,imshow(AA)%模糊的圖像%維納濾波，也是逆濾波H2=H+(H==0);%去除濾波器零點，將H中的0變?yōu)?，避免后面計算出現(xiàn)0/0G2=fftshift(G);F2=G2./H2;F2=fftshift(F2);A2=uint8(ifft
3、2(F2));figure,imshow(A2)實驗結(jié)果如下：產(chǎn)生模糊圖像的濾波器：模糊后的Lenna維納濾波的結(jié)果：向模糊的圖片加入均值為0、方差0.001的高斯噪聲（程序幾乎同上，不再給出）采用直接逆濾波（去除零點）和維納濾波（其中K通過嘗試，取10）的結(jié)果對比如下：如果采用半徑受限的逆濾波,得到結(jié)果如下：可見，半徑受限的逆濾波產(chǎn)生了圖像模糊，同時也濾除了部分噪聲，有振鈴效應；而維納濾波產(chǎn)生的圖像更為清晰，振鈴效應不明顯，但是噪聲比較明顯。二.圖像復原1).針對圖像水平平移a后與原圖像疊加形成的重影圖像所引起的圖像退化，給出連續(xù)函數(shù)退化模型的系統(tǒng)傳遞函數(shù)H（
4、u,v）2).編程用維納濾波實現(xiàn)對一個重影圖像的復原。解：1）設系統(tǒng)沖擊響應為h(x,y),則對上式做傅里葉變換，，即為連續(xù)函數(shù)退化模型的系統(tǒng)傳遞函數(shù)H（u,v）2）A=imread('Lenna.bmp');[M,N]=size(A);a=12;%水平位移12pixelB=[A(:,N-a+1:N)A(:,1:N-a)];B=B*0.5+A*0.5;figure,subplot(121),imshow(B),title('重影圖像（無噪）')u=1:N;H=zeros(M,N);H(1,:)=0.5*(1+exp(-1i*2*pi*(u-N/2)*a/N));
5、forn=2:MH(n,:)=H(1,:);end%k=10;H2=H+(abs(H)<1e-10);%避免除以很小的數(shù)導致誤差，去閾值為1e-10F=fftshift(fft2(B));%F_R=F./H2;%無噪時，維納濾波即逆濾波%F_R=F./H2.*(1-k./(k+(abs(H2))^2));%有噪聲時A_re=uint8(real(ifft2(fftshift(F_R))));figure,subplot(122),imshow(A_re),title('維納濾波后')向重影的圖片加入均值為0、方差0.001的高斯噪聲，選取合適的閾值（經(jīng)試驗，可取
6、0.1,即H2=H+(abs(H)<0.1)，k可取10），實驗結(jié)果如下：閾值選取過小時，噪聲明顯；閾值選取過大時，振鈴效應明顯。K增大時，噪聲變小，圖像亮度下降。三拍攝一張包含硬幣、橡皮、鑰匙等小物件的圖像獲得上述物體的邊緣二值圖像解：拍攝圖像如下：程序如下：A=imread('edge_test.jpg');%1280×960×3A=rgb2gray(A);A=imresize(A,0.25);%為了方便在屏幕上顯示，并且容易觀察到細節(jié)，長寬均縮小為原來的1/4figure,imshow(A)BW1=edge(A,'sobel',0.04);%選取合適的閾值
7、BW2=edge(A,'prewitt',0.04);BW3=edge(A,'roberts',0.04);BW4=edge(A,'log',0.003);BW5=edge(A,'zerocross',0.003);BW6=edge(A,'canny',0.08);figure,subplot(131),imshow(BW1),title('sobel,thresh=0.04')subplot(132),imshow(BW2),title('prewitt,thresh=0.04')subplot(133),imshow(BW3),title('roberts,
8、thresh=0.04'
