a96-04-29.m

Matlab數位影音處理 影像方面的處理

clear all;
f=checkerboard(8);
PSF=fspecial('motion',7,45);
gb=imfilter(f,PSF,'circular');
noise=imnoise(zeros(size(f)),'gaussian',0,0.001);
g=gb+noise;
fr1=deconvwnr(g,PSF);
sn=abs(fft2(noise)).^2;
na=sum(sn(:))/prod(size(noise));
sf=abs(fft2(f)).^2;
fa=sum(sf(:))/prod(size(f));
R=na/fa;
fr2=deconvwnr(g,PSF,R);
NCORR=fftshift(real(ifft2(sn)));
ICORR=fftshift(real(ifft2(sf)));
fr3=deconvwnr(g,PSF,NCORR,ICORR);

%imshow(pixeldup(f,8),[])
figure;
subplot(2,4,1);
imshow(f,[])
subplot(2,4,2);
imshow(gb,[])
subplot(2,4,3);
imshow(noise,[])
subplot(2,4,4);
imshow(g)
subplot(2,4,5);
imshow(fr1,[])
subplot(2,4,6);
imshow(fr2,[])
subplot(2,4,7);
imshow(fr3,[])
subplot(2,4,8);
imshow(g)