数字图像处理维纳滤波处理运动模糊

数字图像处理维纳滤波处理运动模糊一、实验九的程序：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);figuresubplot(1,2,1),imshow(nentaxy,[]),title('附加上的噪声')a=0.15;%x方向的最大移动量为ra的0.2倍，可调b=0.15;%y方向的最大移动量为ca的0.2倍，可调T=0.85;%移动到最大所需的时间默认为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(fxyyp,[]);title('a=0.15,b=0.15时运动模糊的结果')%下面计算原始图像的傅立叶变换估计值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可知时维纳滤波的结果')K=0.02;%特殊常数，一般要用交互的方式确定Fuvyp=(Huv.*conj(Huv)).*Guv./(Huv.*(Huv.*conj(Huv)+K));Rtuxy=abs(ifft2(Fuvyp));subplot(2,2,2),imshow(Rtuxy,[]),title('K=0.02时维纳滤波的结果')K=0.006;%特殊常数，一般要用交互的方式确定Fuvyp=(Huv.*conj(Huv)).*Guv./(Huv.*(Huv.*conj(Huv)+K));Rtuxy=abs(ifft2(Fuvyp));subplot(2,2,3),imshow(Rtuxy,[]),title('K=0.006时维纳滤波的结果')K=0.001;%特殊常数，一般要用交互的方式确定Fuvyp=(Huv.*conj(Huv)).*Guv./(Huv.*(Huv.*conj(Huv)+K));Rtuxy=abs(ifft2(Fuvyp));subplot(2,2,4),imshow(Rtuxy,[]),title('K=0.001时维纳滤波的结果')二、实验九的运行结果：原始图像附加上的噪声a=0.15,b=0.15时运动模糊的结果K可知时维纳滤波的结果K=0.02时维纳滤波的结果K=0.006时维纳滤波的结果K=0.001时维纳滤波的结果