基于matlab二种子大脑图像分割代码

图像分割算法，区域生长法，基于二个种子点分割脑部MRI图像，程序完整，希望跟大家分享。

image=imread（‘mri1.jpg‘）;
I=rgb2gray（image）;
figureimshow（I）title（‘原始图像‘）
I=double（I）/255;
[MN]=size（I）;
x1=65;y1=56;x2=72;y2=113;
seed1=I（x1y1）;             %将生长起始点灰度值存入seed中
seed2=I（x2y2）; 
Y=zeros（MN）;             %作一个全零与原图像等大的图像矩阵Y，作为输出图像矩阵
Z=zeros（MN）; 
Y（x1y1）=I（x1y1）;          %将Y中与所取点相对应位置的点设置为与原图像相同的灰度
Z（x2y2）=I（x2y2）;          %将Z中与所取点相对应位置的点设置为与原图像相同的灰度

sum1=seed1;               %储存符合区域生长条件的点的灰度值的和
suit1=1;                   %储存符合区域生长条件的点的个数
count1=1;                 %记录每次判断一点周围八点符合条件的新点的数目
threshold1=0.05555;         %判断域值
while count1>0
 s1=0;                    %记录判断一点周围八点时，符合条件的新点的灰度值之和
 count1=0;
 for i1=1:M
   for j1=1:ceil（N/2）
     if Y（i1j1）~=0
      if （i1-1）>0 & （i1+1）<（M+1） & （j1-1）>1 & （j1+1）<（（ceil（N/2））+1）
 %判断此点是否为图像边界上的点
       for u=-1:1                                %判断点周围八点是否合和域值条件
        for v=-1:1                               %uv为偏移量
          if  Y（i1+uj1+v）==0 & abs（I（i1+uj1+v）-seed1）<=threshold1    
%判断是

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3209  2008-08-05 14:07  基于matlab的有关医学中人脑部的分割\seg_twoseeds.m

     文件      47468  2005-05-26 09:10  基于matlab的有关医学中人脑部的分割\mri1.jpg

    ..A.SH.      9216  2008-08-05 14:10  基于matlab的有关医学中人脑部的分割\Thumbs.db

     目录          0  2008-08-05 14:10  基于matlab的有关医学中人脑部的分割

----------- ---------  ---------- -----  ----

                59893                    4