xbbyjc0.m

医学图像处理方面的 关于细胞边缘检测的

blood=imread('blood1.BMP');
[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)>225)
            I(i,y)=225;
        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:x
            if (I(i,y)>=T)
                S0=S0+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=zero(x,y);
for i=1:x
    for j=1:y
        if(I(i,j)>=T)
            Seg(i,j)=1
        end
    end
end
SI=1-Seg;
se1=strel('square',3);
SI1=imread(SI,sel);
BW=SI-SI1;
I=uint8(I);
BW1=edge(I,'sobel',0.11);
BW2=edge(I,'log',0.015);
figure(1);
imshow(I);title('Original');
figure(2);
imshow(BW2);title('Soble')
figure(3)
imshow(BW1);title('Gauss-Laplace')
figure(4)
imshow(BW);title('New algorithm')