lvbo.rar_lvbo_维纳lvbo_维纳滤波 chirp

pic=imread('user.jpg'); pic=rgb2gray(pic); [r,c]=size(pic); pic=double(pic); hh=zeros(r,c); figure; v1=8; v2=16; v3=32; for i=1:13 for j=1:13 hh(i,j)=exp(sqrt(i^2+j^2)/240); end end hm=hh(1:13,1:13); pic2=imfilter(pic,hm); maxp=max(max(pic2)); minp=min(min(pic2)); pic2=(pic2-minp)./(maxp-minp).*255; h=fft2(pic2)./fft2(pic); %noise and show1 noisev1=sqrt(v1).*randn(size(pic2)); pic3v1=uint8(pic2+noisev1); subplot(3,3,1);imshow(pic3v1,[]); noisev2=sqrt(v2).*randn(size(pic2)); pic3v2=uint8(pic2+noisev2); subplot(3,3,4);imshow(pic3v2,[]); noisev3=sqrt(v3).*randn(size(pic2)); pic3v3=uint8(pic2+noisev3); subplot(3,3,7);imshow(pic3v3,[]); pic3v1=double(pic3v1); fpic3v1=fft2(pic3v1); pic3v2=double(pic3v2); fpic3v2=fft2(pic3v2); pic3v3=double(pic3v3); fpic3v3=fft2(pic3v3); fp=abs(fft2(pic2)).^2; sp=sum(fp(:))/prod(size(pic)); pc = fftshift(real(ifft2(fp))); fn1=abs(fft2(noisev1)).^2; sn1=sum(fn1(:))/prod(size(noisev1)); nc1=fftshift(real(ifft2(fn1))); fn2=abs(fft2(noisev2)).^2; sn2=sum(fn2(:))/prod(size(noisev2)); nc2=fftshift(real(ifft2(fn2))); fn3=abs(fft2(noisev3)).^2; sn3=sum(fn3(:))/prod(size(noisev3)); nc3=fftshift(real(ifft2(fn3))); %nilvbo pic4v1=deconvwnr(uint8(pic3v1),hm,nc1,pc); subplot(3,3,2);imshow(pic4v1,[]); pic4v2=deconvwnr(uint8(pic3v2),hm,nc2,pc); subplot(3,3,5);imshow(pic4v2,[]); pic4v3=deconvwnr(uint8(pic3v3),hm,nc3,pc); subplot(3,3,8);imshow(pic4v3,[]); %wienner m1=(1./h).*((abs(h)).^2)./((abs(h)).^2+fn1./fp); m2=(1./h).*((abs(h)).^2)./((abs(h)).^2+fn2./fp); m3=(1./h).*((abs(h)).^2)./((abs(h)).^2+fn3./fp); fpic5v1=m1.*fpic3v1; pic5v1=abs(ifft2(fpic5v1)); fpic5v2=m2.*fpic3v2; pic5v2=abs(ifft2(fpic5v2)); fpic5v3=m3.*fpic3v2; pic5v3=abs(ifft2(fpic5v3)); %show subplot(3,3,3);imshow(pic5v1,[]); subplot(3,3,6);imshow(pic5v2,[]); subplot(3,3,9);imshow(pic5v3,[]);