數(shù)字圖像處理維納濾波處理運(yùn)動(dòng)模糊.doc

《數(shù)字圖像處理維納濾波處理運(yùn)動(dòng)模糊.doc》由會(huì)員上傳分享，免費(fèi)在線(xiàn)閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。

1、數(shù)字圖像處理維納濾波處理運(yùn)動(dòng)模糊一、實(shí)驗(yàn)九的程序：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);fig

2、uresubplot(1,2,1),imshow(nentaxy,[]),title('附加上的噪聲')a=0.15;%x方向的最大移動(dòng)量為ra的0.2倍，可調(diào)b=0.15;%y方向的最大移動(dòng)量為ca的0.2倍，可調(diào)T=0.85;%移動(dòng)到最大所需的時(shí)間默認(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;subplo

3、t(1,2,2),imshow(fxyyp,[]);title('a=0.15,b=0.15時(shí)運(yùn)動(dòng)模糊的結(jié)果')%下面計(jì)算原始圖像的傅立葉變換估計(jì)值Guv=fftshift(fft2(fxyyp));Nuv=fftshift(fft2(nentaxy));K=Nuv.*conj(Nuv)./(Fuv.*conj(Fuv));%計(jì)算K值Fuvyp=(Huv.*conj(Huv)).*Guv./(Huv.*(Huv.*conj(Huv)+K));Rtuxy=abs(ifft2(Fuvyp));figuresubplot(2,2,1),imshow(Rtuxy,[]),title(

4、'K可知時(shí)維納濾波的結(jié)果')K=0.02;%特殊常數(shù)，一般要用交互的方式確定Fuvyp=(Huv.*conj(Huv)).*Guv./(Huv.*(Huv.*conj(Huv)+K));Rtuxy=abs(ifft2(Fuvyp));subplot(2,2,2),imshow(Rtuxy,[]),title('K=0.02時(shí)維納濾波的結(jié)果')K=0.006;%特殊常數(shù)，一般要用交互的方式確定Fuvyp=(Huv.*conj(Huv)).*Guv./(Huv.*(Huv.*conj(Huv)+K));Rtuxy=abs(ifft2(Fuvyp));subplot(2,2,3),

5、imshow(Rtuxy,[]),title('K=0.006時(shí)維納濾波的結(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時(shí)維納濾波的結(jié)果')二、實(shí)驗(yàn)九的運(yùn)行結(jié)果：