96-05-06.m

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

%96-05-06
clear all;
f=checkerboard(8);
imshow(pixeldup(f,8));

PSF=fspecial('motion',7,45);
gb=imfilter(f,PSF,'circular');
noise=imnoise(zeros(size(f)),'gaussian',0,0.001);
g=gb+noise;
g=edgetaper(g,PSF);
fr1=deconvreg(g,PSF,4);
fr2=deconvreg(g,PSF,0.4,[1e-7 1e7]);
figure;
subplot(1,2,1);
imshow(fr1,[])
subplot(1,2,2);
imshow(fr2,[])
%example 5.10
PSF2=fspecial('gaussian',7,10);
SD=0.01;
g1=imnoise(imfilter(f,PSF2),'gaussian',0,SD^2);
DAMPAR=10*SD;
LIM=ceil(size(PSF2,1)/2);
WEIGHT=zeros(size(g));
WEIGHT(LIM+1:end-LIM,LIM+1:end-LIM)=1;
NUMT=5;
fr2=deconvlucy(g1,PSF2,NUMT,DAMPAR,WEIGHT); 
NUMT1=10;
fr3=deconvlucy(g1,PSF2,NUMT1,DAMPAR,WEIGHT); 
NUMT2=20;
fr4=deconvlucy(g1,PSF2,NUMT2,DAMPAR,WEIGHT); 
NUMT3=100;
fr5=deconvlucy(g1,PSF2,NUMT3,DAMPAR,WEIGHT); 
figure;
subplot(2,3,1);
imshow(f)
subplot(2,3,2);
imshow(g1);
subplot(2,3,3);
imshow(fr2);
subplot(2,3,4);
imshow(fr3);
subplot(2,3,5);
imshow(fr4);
subplot(2,3,6);
imshow(fr5);
%example 5.11
clear all;
f=checkerboard(8);
PSF=fspecial('gaussian',7,10);
g=pixeldup(PSF,73);
SD=0.01;
gb=imnoise(imfilter(f,PSF),'gaussian',0,SD^2);
INITPSF=ones(size(PSF));
DAMPAR=10*SD;
LIM=ceil(size(PSF,1)/2);
WEIGHT=zeros(size(gb));
WEIGHT(LIM+1:end-LIM,LIM+1:end-LIM)=1;
NUMIT=5;
[fr,PSFe]= deconvblind(gb,INITPSF,NUMIT,DAMPAR,WEIGHT);
fr=pixeldup(PSFe,73);
NUMIT2=10;
[fr2,PSFe]= deconvblind(gb,INITPSF,NUMIT2,DAMPAR,WEIGHT);
fr2=pixeldup(PSFe,73);
NUMIT3=20;
[fr3,PSFe]= deconvblind(gb,INITPSF,NUMIT3,DAMPAR,WEIGHT);
fr3=pixeldup(PSFe,73);
figure;
subplot(2,2,1);
imshow(g,[])
subplot(2,2,2);
imshow(fr,[])
subplot(2,2,3);
imshow(fr2,[])
subplot(2,2,4);
imshow(fr3,[])