log_edge.m

来自「图像预处理归一化、边缘检测等的各个matlab小程序。希望对大家有用。」· M 代码 · 共 29 行

 
[m,n]=size(x);
e=repmat(logical((0)),m,n);
rr=2:m-1;cc=2:n-1;
fsize=ceil(sigma*3)*2+1;
op=fspecial('log',fsize,sigma);
op=op-sum(op(:))/prod(size(op));
b=filter2(op,a);
thresh=.75*mean2(abs(b(rr,cc)));
[rx,cx]=find(b(rr,cc)<0&b(rr,cc+1)>0&abs(b(rr,cc)-b(rr,cc+1))>thresh);
e((rx+1)+cx*m)=1;
[rx,cx]=find(b(rr,cc-1)>0&b(rr,cc)<0&abs(b(rr,cc-1)-b(rr,cc))>thresh);
e((rx+1)+cx*m)=1;
[rx,cx]=find(b(rr,cc)<0&b(rr+1,cc)>0&abs(b(rr,cc)-b(rr+1,cc))>thresh);
e((rx+1)+cx*m)=1;
[rx,cx]=find(b(rr-1,cc)>0&b(rr,cc)<0&abs(b(rr-1,cc)-b(rr,cc))>thresh);
e((rx+1)+cx*m)=1;
[rz,cz]=find(b(rr,cc)==0);
if~isempty(rz)
    zero=(rz+1)+cz*m;
    zz=find(b(zero-1)<0&b(zero+1)>0&abs(b(zero-1)-b(zero+1))>2*thresh);
    e(zero(zz))=1;
    zz=find(b(zero-1)>0&b(zero+1)<0&abs(b(zero-1)-b(zero+1))>2*thresh);
    e(zero(zz))=1;
    zz=find(b(zero-m)<0&b(zero+m)>0&abs(b(zero-m)-b(zero+m))>2*thresh);
    e(zero(zz))=1;
    zz=find(b(zero-m)>0&b(zero+m)<0&abs(b(zero-m)-b(zero+m))>2*thresh);
    e(zero(zz))=1;
end