tuxhfu4.m

包括MATLAB7.0实用指南（下）中所有例子的源代码

i=imread('cameraman.tif');
%i=i(50+[1:256],2+[1:256],:);
%imshow(i);
psf=fspecial('motion',13,45);
figure;
imshow(psf,[],'notruesize');
b=imfilter(i,psf,'circ','conv');
%v=0.002
%bn=imnoise(b,'gaussian',0,v);
figure;
imshow(b);
%np=v*prod(size(i));
npsf=ones(size(psf));
[j p]=deconvblind(b,npsf,30);
figure;
imshow(j);
figure;
imshow(p,[],'notruesize');
w=edge(i,'sobel',28);
se1=strel('disk',1);
se2=strel('line',13,45);
w=~imdilate(w,[se1 se2]);
w=padarray(w(2:end-3,2:end-3),[2 2]);
figure;
imshow(w);
p1=p;
p1(find(p1<0.01))=0;
W=double(w);
[j2 p2]=deconvblind(b,p1,30,[],W);
figure;
imshow(j2);
figure;
imshow(p2,[],'notruesize');