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

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;figuresubplot(1,2,1,imshow(nentaxy,title

2、(附加上的噪声 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;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=f

3、ftshift(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(

4、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 可 知 时 维 纳 滤 波 的