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

1、数字图像处理维纳滤波处理运动模糊一、实验九的程序: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倍,可调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

3、(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;%特殊常数,

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时维纳滤波的结果 ' 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 可 知 时 维 纳