水平集level set

李春明2005年关于水平集的文章，及其文章所对应的代码，可运行，希望对大家有帮助

% This Matlab file demomstrates the level set method in Li et al‘s paper
%    “Level Set Evolution Without Re-initialization: A New Variational Formulation“
%    in Proceedings of CVPR‘05 vol. 1 pp. 430-436.
% Author: Chunming Li all rights reserved.
% E-mail: li_chunming@hotmail.com
% URL:  http://www.engr.uconn.edu/~cmli/

clear all;
close all;
Img = imread（‘twocells.bmp‘）;  % The same cell image in the paper is used here
Img=double（Img（::1））;
sigma=1.5;    % scale parameter in Gaussian kernel for smoothing.
G=fspecial（‘gaussian‘15sigma）;
Img_smooth=conv2（ImgG‘same‘）;  % smooth image by Gaussiin convolution
[IxIy]=gradient（Img_smooth）;
f=Ix.^2+Iy.^2;
g=1./（1+f）;  % edge indicator function.

epsilon=1.5; % the papramater in the definition of smoothed Dirac function

timestep=5;  % time step
mu=0.2/timestep;  % coefficient of the internal （penalizing） energy term P（\phi）
          % Note: the product timestep*mu must be less than 0.25 for stability!

lambda=5; % coefficient of the weighted length term Lg（\phi）
alf=1.5;  % coefficient of the weighted area term Ag（\phi）;
          % Note: Choose a positive（negative） alf if the initial contour is outside（inside） the object.


% define initial level set function （LSF） as -c0 0 c0 at points outside on
% the boundary and inside of a region R respectively.
[nrow ncol]=size（Img）;  
c0=4;   
initialLSF=c0*ones（nrowncol）;
w=8;
initialLSF（w+1:end-w w+1:end-w）=0;  % zero level set is on the boundary of R. 
                                     % Note: this can be commented out. The intial LSF does NOT necessarily need a zero level set.
                                     
initialLSF（w+2:end-w-1 w+2: end-w-1）=-c0; % negative constant -c0 inside of R postive constant c0 outside of R.
u=initialLSF;
figure;imagesc（Img）;colormap（gray）;hold on;
[ch] = contour（u[0 0]‘r‘）;                          
title（‘Initial contour‘）;

% start level set evolution
for n=1:300
    u=EVOLUTION（u g lambda mu alf epsilon timestep 1）;     
    if mod（n20）==0
        pause（0.001）;
        imagesc（Img）;colormap（gray）;hold on;
        [ch] = contour（u[0 0]‘r‘）; 
        iterNum=[num2str（n） ‘ iterations‘];        
        title（iterNum）;
        hold off;
    end
end
imagesc（Img）;colormap（gray）;hold on;
[ch] = contour（u[0 0]‘r‘）; 
totalIterNum=[num2str（n） ‘ iterations‘];  
title（[‘Final contour ‘ totalIterNum]）;

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

     文件       2465  2007-03-20 14:26  LevelSet_ChunmingLi_v0\Demo1.m

     文件       2492  2007-03-20 14:27  LevelSet_ChunmingLi_v0\Demo1_ManualBinaryInitial.m

     文件       2464  2007-03-20 14:26  LevelSet_ChunmingLi_v0\Demo2.m

     文件       2467  2007-03-20 14:28  LevelSet_ChunmingLi_v0\Demo2_ManualBinaryInitial.m

     文件       2071  2006-10-30 02:20  LevelSet_ChunmingLi_v0\EVOLUTION.m

     文件     801246  2006-02-03 16:17  LevelSet_ChunmingLi_v0\levelset_CVPR05.pdf

    ..A.SH.      6144  2007-03-20 14:23  LevelSet_ChunmingLi_v0\Thumbs.db

     文件      16434  2004-07-23 19:36  LevelSet_ChunmingLi_v0\twocells.bmp

     文件       8134  2004-07-05 18:59  LevelSet_ChunmingLi_v0\twoObj.bmp

     目录          0  2007-03-20 14:23  LevelSet_ChunmingLi_v0

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

               843917                    10