GAC算法matlab实现

可以用，注释是我自己写的，已经够全了我感觉

function newim = createimage（ imuFlag ）
%%% 子函数：createimage（）实现当前零水平集（演化曲线）在原图上的显示
%%%         和嵌入函数u的重新初始化
%%% 参数说明：
%%% im     - 原始图像
%%% u      - 当前水平集函数
%%% Flag   - 是否进行重新初始化的标志，0：否；1：是
%%% newim  - 返回的当前零水平集（红色曲线）叠加在原图上的新图像

[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;
newim = uint8（ newim ）;

%%- 如果Flag=1，则重新初始化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

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

     文件     120054  2014-04-01 16:39  GAC\3.bmp

     文件       1259  2014-04-01 16:39  GAC\createimage.m

     文件       4076  2014-04-03 18:07  GAC\GAC_upwind.asv

     文件       4076  2014-04-04 16:02  GAC\GAC_upwind.m

     文件       1300  2014-04-04 13:03  GAC\gauss.asv

     文件       1152  2014-04-01 16:39  GAC\gauss.m

     目录          0  2014-04-09 19:48  GAC

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

               131917                    7