function?finalsdf?=?lsCV（I?initsdfdelta_tmunvlambda1lambda2epsilonnum）
%?initsdf?初始化符號距離函數(shù)
%?delta_t?為步長??mu/nv為周長/面積權重系數(shù)
%?lambda1/lambda2為能量項權重系數(shù)
%?epsilon?為dirac函數(shù)參數(shù)

%?check?to?see?if?the?image?is?a?color?image...
d=ndims（I）;%d?=?length（size（in_image））;
if?d==3
????I=double（rgb2gray（I））;
elseif?d==2
????I=double（I）;
end
%delta_t*狄力克函數(shù)
dirac_delta_t?=?delta_t?*??（?1?./?pi?）?.*?（?epsilon?./?（?epsilon.^2?+?initsdf.^2?）?）;

%?計算輪廓線內(nèi)外能量
if?num>0
????H_phi?=?0.5?.*?（?1?+?（?2?./?pi?）?.*?atan（?initsdf?./?epsilon?）?）;
????H_phi_minus?=?1?-?H_phi;
????c1?=?sum（?sum（?I?.*?H_phi?）?）?/??sum（?sum（?H_phi?）?）;
????c2?=?sum（?sum（?I?.*?H_phi_minus?）?）?/?sum（?sum（?H_phi_minus?）?）;
else
????c1=220;
????c2=5;
%?????c1=0;
%?????c2=130;
end
inside?=?（?I?-?c1?）.^2;
outside?=?（?I?-?c2?）.^2;
energy_term?=?-nv?-?lambda1?.*?inside?+?lambda2?.*?outside;%能量項

%?中心差分
dx_central?=?（?circshift（?initsdf?[?0?-1?]?）?-?circshift（?initsdf?[?0?1?]?）?）?/?2;
dy_central?=?（?circshift（?initsdf?[?-1?0?]?）?-?circshift（?initsdf?[?1?0?]?）?）?/?2;

%?div（delt_phi/|delta_phi|）
abs_grad_phi?=?（?sqrt（?dx_central.^2?+?dy_central.^2?）?+?0.00001?）;
x?=?dx_central?./?abs_grad_phi;
y?=?dy_central?./?abs_grad_phi;
grad_term?=?mu?.*?divergence（?x?y?）;%梯度能量
%?
%?tt=dx_central;
%?[ab]=size（tt）;
%?c=1:a;d=1:b;
%?[CD]=meshgrid（dc）;
%?figuremesh（CDtt）;

%?phi（n+1）
finalsdf?=?initsdf?+?dirac_delta_t?.*?（?grad_term?+?energy_term?）;
%???finalsdf?=?initsdf;%保存每次迭代結(jié)果

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1570??2013-09-01?21:20??lsCV.m
?????文件????????2452??2013-11-21?21:45??lsCVTest.m