%%%?“C-V”模型（Active?contour?without?edge）

clear?all;
close?all;
clc;

%Img=imread（‘synthetic.png‘）;
%Img?=?imread（‘europe_night_ceshi.jpg‘）;?%%%%%%好用
%Img?=?imread（‘noisyStar_SNR_20_to_10.bmp‘）;?%%%%%%好用，效果很好
%Img?=?imread（‘2.bmp‘）;?%%%%%%好用，效果很好，選取合適的系數nu
%Img?=?imread（‘synthetic1.png‘）;??%%%%%%好用，效果很好

%?I?=?imread（‘01.tif‘）;

Img?=?imread（‘2.bmp‘）;

%?
%?imshow（I）
%?rect?=?getrect（gcf）;
%?Img?=?I（rect（2）:rect（2）+rect（4）rect（1）:rect（1）+rect（3））;

%?imshowMy（IM）

%?Img=imread（‘s:\b.tif‘）;
%?Img=double（Img）;
Img=double（rgb2gray（Img））;
Img=imresize（Img[8080]）;
figure（1）;?imshow（uint8（Img））;

[nxny]=size（Img）;
%%-?將初始曲線C設置為圓

%ic=50.5;?????????
%jc=20.5;

ic=floor（nx/2）;?????????%?計算初始圓形曲線的圓心
jc=floor（ny/2）;
r=25;
%r=ic/1.2;?????????????????%?圓形曲線的半徑

%%-?初始化u為距離函數
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‘）;


%%-?初始化參數
epsilon=1.0;????????????%?Heaviside函數參數設置
?
%%-?對于不同的圖片，nu的差異很大，數值太小時，輪廓線容易散。
nu=0.1*255*255;??%%%%?europe_night_ceshi.jpg
%nu=0.01*255*255;???%%%%??noisyStar_SNR_20_to_10.bmp
%nu=0.1*255*255;????%%%%??2.bmp
%nu=250;???