医学图像分割 matlab程序

附录B
主程序：
clc;
clear;
close all;
[fname,fpath]=uigetfile('D:\MATLAB1\work\brain_tumor detection\brain_tumor\*.bmp','Open a File',100,200);
s=strcat(fpath,fname);
if fname~=0
   II=imread(s);
end
imshow(II,[]);
I=double(II(:,:,1));
title(' 原始图像')
[OB1,NUM1]=interactive(II);
uiwait(msgbox('Please input the coordinate of the seed:','Notice message!','modal'));
[XX,YY]=ginput(1);
y=round(XX);
x=round(YY);
[m,n]=size(I);
D = zeros(m, n);  
%%%%%%%%111111%%%%%%
k=0;
flag=0;
xx=x;
yy=y;
while 1
    L=I(x,y)
    [m,g]=fenleiqi(OB1,L);
    [K,a]=max(m);
    if (L==K)|(g(a)==g(NUM1+1))

        D(x,y)=1;
        k=0;
    else
        k=k+1;
        if k==25
           break;
    end
        if (k==5)
          x=x+1;
           flag=flag+1;
           flag=mod(flag,2);
           continue;
        end
    end
  if flag==0
     y=y+1;
  end
  if flag==1
     y=y-1;
  end
end
%%%%%%%%22222%%%%%%
k=0;
flag=1;
x=xx;
y=yy;
while 1
    L=I(x,y)
    [m,g]=fenleiqi(OB1,L)
    [K,a]=max(m)