chan-vese的matlab代码实现

这是chan-vese模型的matlab代码实现，对于c-v模型想了解的用户还是很有帮助的。

function indicator = checkstop（oldnewdt）
% indicate whether we should performance further iteraions or stop

% Copyright （c） 2009 
% Yue Wu @ ECE Department Tufts University
% All Rights Reserved  

layer = size（new3）;

for i = 1:layer
    old_{i} = old（::i）;
    new_{i} = new（::i）;
end

if layer
    ind = find（abs（new）<=.5）;
    M = length（ind）;
    Q = sum（abs（new（ind）-old（ind）））./M;
    if Q<=dt*.18^2
        indicator = 1;
    else
        indicator = 0;
    end
else
    ind1 = find（abs（old_{1}）<1）;
    ind2 = find（abs（old_{2}）<1）;
    M1 = length（ind1）;
    M2 = length（ind2）;
    Q1 = sum（abs（new_{1}（ind1）-old_{1}（ind1）））./M1;
    Q2 = sum（abs（new_{2}（ind2）-old_{2}（ind2）））./M2;
    if Q1<=dt*.18^2 && Q2<=dt*.18^2
        indicator = 1;
    else
        indicator = 0;
    end
end
return

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件       16663  2014-02-12 12:44  4colors.JPG
     文件       24291  2014-02-12 12:44  anti-mass.jpg
     文件        2879  2014-02-12 12:44  brain.jpg
     文件         848  2014-02-12 12:44  checkstop.m
     文件       14373  2014-02-12 12:44  chenvese.m
     文件        1897  2014-02-12 12:44  demo chenvese.m
     文件      136210  2014-02-12 12:44  flowers.jpg
     文件         389  2014-02-12 12:44  Heaviside.m
     文件      809765  2014-02-12 12:44  HELP.pdf
     文件         767  2014-02-12 12:44  kappa.m
     文件        2101  2014-02-12 12:44  maskcircle2.m
     文件        4052  2014-02-12 12:44  READ ME.txt
     文件         855  2014-02-12 12:44  reinitialization.m
     文件         706  2014-02-12 12:44  showphi.m
     文件        1306  2014-02-12 12:44  license.txt