基于边缘的图像分割方法

GAC模型是基于边缘的图像分割方法，当图像中的对象和背景的边界处灰度值存在较大差异，那么对象的边缘就将形成明显的轮廓。

function newim = createimage（ imuFlag ）
% CREATEIMAGE returns a color image （ 3 channel uint8 image）
% in which the current zero_level_set is in red overlaied
% onto the original image （im.）.And then re_initializes u.
[nnynnx]=size（im）;
curvIndex=zeros（5*nnx*nny2）;
curvImag=zeros（[nnynnx]）;
num=0;
for i = 2 : nnx - 1
         for j = 2 : nny - 1
            if u（ij）<0  & （u（i+1j）>=0 | u（i-1j）>=0 | u（ij+1）>=0 | u（ij-1）>=0）
               num=num+1; curvIndex（num1）=i;curvIndex（num2）=j;
               curvImag（ij）=255;
            end
         end
      end
tempim = im;
tempim（curvImag>0） = 255;
newim（ : : 1 ） = tempim;
tempim（curvImag>0） = 0;
newim（ : : 2 ） = tempim;
newim（ : : 3 ） = tempim;
% convert to uint8
newim = uint8（ newim ）;
%If Flag==1re_initialize u
if Flag
   new_u = zeros（nnxnny）;
   dist=zeros（1num）;
   for j=1:nny
     for i=1:nnx
       for k=1:num
          dist（k）=sqrt（（i-curvIndex（k1））.^2+（j-curvIndex（k2））.^2）;
       end
       new_u（ij）=min（dist）;
       if u（ij）<0
          new_u（ij）=-new_u（ij）;
       end
     end
   end
   u=new_u;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-05-18 15:39  GAC\
     文件       11274  2012-03-12 15:49  GAC\03.jpg
     文件      120054  2006-04-22 13:07  GAC\3.bmp
     文件        4063  2012-05-17 19:28  GAC\brain.jpg
     文件        1103  2007-11-20 12:07  GAC\createimage.m
     文件        7554  2012-05-17 19:29  GAC\flower.jpg
     文件        3514  2012-05-18 14:16  GAC\gacgaijinf.m
     文件        1073  2006-07-13 18:21  GAC\gauss.m
     文件        4282  2012-05-17 19:04  GAC\images （1）.jpg
     文件        8134  2004-07-05 18:59  GAC\twoObj.bmp
     文件       16434  2004-07-23 20:36  GAC\twocells.bmp