%Threshoid为分割阈值
[filename,pathname]=...
　　uigetfile({'*.*';'*.bmp';'*.tif';'*.png'},'select picture');　%选择图片路径
str=[pathname filename];　%合成路径+文件名
I=imread(str);
I=double(I);
[m,n,k]=size(I);%图像大小
Threshold=50/255.0;%归一化，并设置阈值为50
ColorRed=I(:,:,1);%分成3个单色图片
ColorGreen=I(:,:,2);
ColorBlue =I(:,:,3);
OutputImgRed=zeros(m,n);%3个输出单色图片的存储矩阵
OutputImgGreen=zeros(m,n);
OutputImgBlue=zeros(m,n);
Checkseed=zeros(m,n);%种子矩阵
firstlocation=1;%堆栈指针指向顶点1
for j=1:n