tuxhfu4.asv

来自「包括MATLAB7.0实用指南（下）中所有例子的源代码」· ASV 代码 · 共 31 行

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)