%%%?本程序采用半隱式方案實(shí)現(xiàn)變分水平集圖像分割方法中的
%%%?“C-V”模型（Active?contour?without?edge）
clear?all;
close?all;
clc;

Img=imread（‘3-1.bmp‘）;
Img=double（rgb2gray（Img））;
%?Img=imresize（Img[8080]）;
figure（1）;?imshow（uint8（Img））;

[nxny]=size（Img）;
%%-?將初始曲線C設(shè)置為圓
ic=floor（nx/2）;?????????%?計(jì)算初始圓形曲線的圓心
jc=floor（ny/2）;
r=ic/3;?????????????????%?圓形曲線的半徑

%%-?初始化u為距離函數(shù)
u?=?zeros（[nxny]）;?????
for?i=1:nx
????for?j=1:ny
????????u（ij）=?r-sqrt（（i-ic）.^2+（j-jc）.^2）;
????end
end
%%-?將初始圓形曲線疊加在原始圖片上
figure（2）;??????????????
imshow（uint8（Img））;
hold?on;
[ch]?=?contour（u[0?0]‘r‘）;


%%-?初始化參數(shù)
epsilon=1.0;????????????%?Heaviside函數(shù)參數(shù)設(shè)置
nu=250;?????????????????
delta_t=0.1;

nn=0;
%%-?迭代計(jì)算開(kāi)始
for?n=1:1000
????%%-?計(jì)算正則化的Heavside函數(shù)
????H_u?=?0.5*（1+（2/pi）*atan（u/epsilon））;

????%%-?由當(dāng)前u計(jì)算參數(shù)c1和c2
????c1=sum（sum（H_u.*Img））/sum（sum（H_u））;
????c2=sum（sum（（1-H_u）.*Img））/sum（sum（1-H_u））;

????%%-?由當(dāng)前c1和c2更新u
????delta_H?=?（1/p