seedgrow.m

自己编的阈值分割的源程序

function seedgrow(seedx,seedy,num)
load file;
T=imread(file);
[m,n]=size(T);
T=double(T);
delta=50;
Logical(seedx,seedy)=num;
if (seedx-1)>=1&&Logical(seedx,seedy)==0&&abs(T(seedx,seedy)-T(seedx-1,seedy))<=delta
     seedgrow(seedx-1,seedy,num);
 else
     return
end
if (seedy-1)>=1&&Logical(seedx,seedy)==0&&abs(T(seedx,seedy)-T(seedx,seedy-1))<=delta
    seedgrow(seedx,seedy-1,num);
else
    return
end
if (seedx+1)<=n&&Logical(seedx,seedy)==0&&abs(T(seedx,seedy)-T(seedx+1,seedy))<=delta
   seedgrow(seedx+1,seedy,num);
else
    return
end
if (seedy+1)<=m&&Logical(seedx,seedy)==0&&abs(T(seedx,seedy)-T(seedx,seedy+1))<=delta
    seedgrow(seedx,seedy+1,num);
else
    return    
end
return