数学形态学的边缘检测程序

blood = imread（‘jx1.jpg‘）;
[xyz]=size（blood）;                % 求出图象大小
b=double（blood）;                  
N =sqrt（100） * randn（xyz）;        %  生成方差为10的白噪声
I=b+N;                            %  噪声干扰图象
for i=1:x                         % 实际图象的灰度为0～255
    for j=1:y
        if （I（ij）>255）
            I（ij）=255;
        end
        if （I（ij）<0）
            I（ij）=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（ij）>=T）
                S0=S0+I（ij）;
                n0=n0+1;
            end
            if （I