數(shù)字圖像處理維納濾波處理運動模糊

1、數(shù)字圖像處理維納濾波處理運動模糊一、實驗九的程序:a=imread('01.jpg'fxy=a(:,:,1;fxy=double(fxy;ra,ca=size(fxy;Fuv=fftshift(fft2(fxy;subplot(1,2,1,imshow(fxy,;title('原始圖像 ' subplot(1,2,2,imshow(abs(Fuv,0,250000;u=1:ra;v=1:ca;u,v=meshgrid(u,v;nentaxy=20.*imnoise(zeros(ra,ca,'gaussian',0,0.004;figuresubp

2、lot(1,2,1,imshow(nentaxy,title('附加上的噪聲 'a=0.15;%x方向的最大移動量為 ra 的 0.2倍,可調(diào)b=0.15;%y方向的最大移動量為 ca 的 0.2倍,可調(diào)T=0.85;%移動到最大所需的時間默認(rèn)為 0.85Huv=ones(ra,ca;uv=(u-ra/2-1.*a+(v-ca/2-1.*b+eps;Huv=T.*sin(pi.*uv.*exp(-j.*pi.*uv./(pi.*uv;Guv=Huv.*Fuv;fxyyp=ifft2(Guv;fxyyp=abs(fxyyp+nentaxy;subplot(1,2,2,imshow

3、(fxyyp,;title('a=0.15,b=0.15時運動模糊的結(jié)果 ' %下面計算原始圖像的傅立葉變換估計值Guv=fftshift(fft2(fxyyp;Nuv=fftshift(fft2(nentaxy;K=Nuv.*conj(Nuv./(Fuv.*conj(Fuv;%計算 K 值Fuvyp=(Huv.*conj(Huv.*Guv./(Huv.*(Huv.*conj(Huv+K;Rtuxy=abs(ifft2(Fuvyp;figuresubplot(2,2,1,imshow(Rtuxy,title('K可知時維納濾波的結(jié)果 'K=0.02;%特殊常數(shù),

4、一般要用交互的方式確定Fuvyp=(Huv.*conj(Huv.*Guv./(Huv.*(Huv.*conj(Huv+K;Rtuxy=abs(ifft2(Fuvyp;subplot(2,2,2,imshow(Rtuxy,title('K=0.02時維納濾波的結(jié)果 ' K=0.006;%特殊常數(shù),一般要用交互的方式確定Fuvyp=(Huv.*conj(Huv.*Guv./(Huv.*(Huv.*conj(Huv+K;Rtuxy=abs(ifft2(Fuvyp;subplot(2,2,3,imshow(Rtuxy,title('K=0.006時維納濾波的結(jié)果 ' K=0.001;%特殊常數(shù),一般要用交互的方式確定Fuvyp=(Huv.*conj(Huv.*Guv./(Huv.*(Huv.*conj(Huv+K;Rtuxy=abs(ifft2(Fuvyp;subplot(2,2,4,imshow(Rtuxy,title('K=0.001時維納濾波的結(jié)果 ' 二、實驗九的運行結(jié)果:原 始 圖像 附 加 上 的 噪 聲 a=0.15,b=0.15時 運 動 模 糊 的 結(jié) 果K 可 知 時 維 納