11.m

matlab的边缘检测程序

blook=imread('D:work\work.jpg');
blood=rgb2gray(blook);
[x,y]=size(blood);
b=double(blood);
N=sqrt(100)*randn(x,y);
I=b+N;
for i=1:x
    for j=1:y
        if(I(i,j)>255)I(i,j)=255;
        end
        if(I(i,j)<0)I(i,j)=0;
        end
    end
end
z0=max(max(I));
z1=min(min(I));
T=(z0+z1)/2;
TT=0;
S0=0;
n0=0;
S1=0;
n1=0;
allow=0.5;
d=abs(T-TT);
count=0;
while(d>=allow);
    count=count+1;
    for i=1:x;
        for j=1:y;
            if(I(i,j)>=T)S1=S1+I(i,j);
                n0=n0+1;
            end
            if (I(i,j)<T) S1=S1+I(i,j);
                n1=n1+1;
            end
        end
    end
    T0=S0/n0;
    T1=S1/n1;
    TT=(T0+T1)/2;
    d=abs(T-TT);
    T=TT;
end
Seg=zeros(x,y);
for i=1:x;
    for j=1:y;
        if(I(i,j)>=T)
            Seg(i,j);
        end
    end
end
S1=1-Seg;
sel=strel('square',3);
S11=imerode(S1,sel);
BW=S1-S11;
I=uint8(I);
BW1=edge(I,'sobel',0.09);
BW2=edge(I,'log',0.015);
figure(1);imshow(I);title('原始图像');
figure(2);imshow(BW1);title('soble算子')'
figure(3);imshow(BW2);title('gauss-laplace');
figure(4);imshow(BW);title('新算子');