图像预处理归一化、边缘检测等的各个matlab小程序。希望对大家有用。

				 
				[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&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)thresh);
				e((rx+1)+cx*m)=1;
				[rx,cx]=find(b(rr,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)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&abs(b(zero-1)-b(zero+1))>2*thresh);
				    e(zero(zz))=1;
				    zz=find(b(zero-1)>0&b(zero+1)2*thresh);
				    e(zero(zz))=1;
				    zz=find(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)2*thresh);
				    e(zero(zz))=1;
				end