图像分割算法 GVF snake matlab编写程序

snake模型的图像分割算法，用matlab编写而成。你可以访问那上面的网站，那里还有作者的关于GVF snake的论文

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 
% 		SNAKE DEMO ver 1.0 beta
% 
%   Interface by
%
%	 Dejan Tomazevic	（e-mail dejan.tomazevic@kiss.uni-lj.si）	10/2/98
%   Copyright （c） 1998 by Dejan Tomazevic
%	 BIPROG Faculty of Electrical Engineering University of Ljubljana
%
%   Snake and GVF functions by
%   Chenyang Xu and Jerry L. Prince 6/17/97
%   Copyright （c） 1996-97 by Chenyang Xu and Jerry L. Prince
%   Image Analysis and Communications Lab Johns Hopkins University
%
%
%	 This program was written during Dejan Tomazevic visit Johns Hopkins University in
% 	 purpose to show the use of traditional snakes and GVF snakes M fuctions.
%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%% PROGRAM BEGINS %%%%%%%%%%


%%%%%%% INICIALIZATIONA %%%%%%%%%%%%%%
clear;

% Add the tools to the current path
p = path;
path（p ‘./snake;./subf‘）;

%%%%%%%%%%%%% define global variables %%%%%%%%%%%%%%%%%%%
%%%% image processing variables
global Image1;							% original image
global sigma;
global mu;
global alpha beta gamma kappa dmin dmax;		% parameters for the snake
global NoGVFIterations;				% number of GVF iterations
global NoSnakeIterations;		   % number of Snake iterations
global SchangeInFieldType;
global VectorFieldButt;				% define the type of vector filed
global SnakeON;						% indicate if snake is visible
global IncSnakeRadius;				% inicializaton snake radius
global CircleOn;						% inicalization snake will be circle
global SnakeDotsON;				   % 1 if snake dots should be displeyed
global GradientOn;					% 1 if gradient is applayed with blur
global XSnake YSnake;				% conture of the snake
global XSnakeInc YSnakeInc;   	% incicialization conture of the snake


sigma=0;
mu=0.1;
alpha=0.05;
beta=0;
gamma=1;
kappa=0.6;			
dmin=0.5;
dmax=2;
SchangeInFieldType=1;

NoGVFIterations=80;
NoSnakeIterations=40;

VectorFieldButt（1）=1;				% standard field
VectorFieldButt（2）=0;				% GVF filed
VectorFieldButt（3）=1;				% normalized GVF
SnakeON=0;								% snake is not drown on the picture at the begining
IncSnakeRadius=0.5;

CircleOn=1;
SnakeDotsON=1;
GradientOn=0;

%%%% interface variables
global VectorOfLocalMenuHD;		%vector of local objects that can be arased
VectorOfLocalMenuHD=[];
global HDmainf;			 			%main figure handle
global HDorigPic;						%original picture axes handle
global HDbluredPic;				   %blured picture axes handle
global HDvectorFPic;					%handle of vector field picture
global HDhelpAxes;					%handle of help window
global HDhelpSlider1 HDhelpSlider2;	%handle of the slider for help text 
global HDvectorOfTexthd;			% vector of help text handels
global xsize ysize;					%size of the picture
global HDMenuAxes;					%Axes for the help
global HDSnakeLine;					%vector of Handles of Snake lines on the picture
global FileName PathName;			% name of the file
global ExampleNo;						% number of examp

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

     文件      31456  1998-10-02 08:33  snake_demo\head.jpg

     文件      47490  1999-01-06 13:49  snake_demo\images\chest.pgm

     文件       1064  1998-09-30 09:50  snake_demo\images\heart.mat

     文件      25615  1998-09-29 10:01  snake_demo\images\heart.pgm

     文件       4248  1998-09-30 09:58  snake_demo\images\heart2.mat

     文件       4056  1998-10-01 06:18  snake_demo\images\heart3.mat

     文件      65594  1999-01-06 13:32  snake_demo\images\new.pgm

     文件       2792  1998-09-30 10:26  snake_demo\images\room.mat

     文件       4109  1997-08-28 12:32  snake_demo\images\room.pgm

     文件       2928  1998-09-30 09:43  snake_demo\images\room2.mat

     文件       2928  1998-09-30 09:44  snake_demo\images\room3.mat

     文件       3304  1998-10-01 06:20  snake_demo\images\u64.mat

     文件       4109  1997-08-28 12:32  snake_demo\images\U64.pgm

     文件       2952  1998-10-07 04:51  snake_demo\images\u642.mat

     文件       3800  1998-10-06 09:00  snake_demo\sdemo.m

     文件        916  2000-04-12 20:06  snake_demo\snake\BoundMirrorEnsure.m

     文件        699  2000-04-12 20:06  snake_demo\snake\BoundMirrorExpand.m

     文件        488  2000-04-12 20:06  snake_demo\snake\BoundMirrorShrink.m

     文件       1096  1997-08-28 12:32  snake_demo\snake\Contents.m

     文件        599  2008-04-20 23:31  snake_demo\snake\dt.asv

     文件        608  1997-08-28 12:32  snake_demo\snake\dt.m

     文件        500  1997-08-28 12:32  snake_demo\snake\gaussianBlur.m

     文件        255  1997-08-28 12:32  snake_demo\snake\gaussianMask.m

     文件        895  1998-10-08 10:23  snake_demo\snake\GGVF.m

     文件       1640  2000-04-12 20:06  snake_demo\snake\GVF.m

     文件        257  1998-09-30 06:08  snake_demo\snake\imdisp.m

     文件       3256  1997-08-28 12:32  snake_demo\snake\rawread.m

     文件       1536  1998-10-07 05:15  snake_demo\snake\rawwrite.m

     文件       1504  1998-09-21 05:42  snake_demo\snake\snakedeform.m

     文件       1856  1997-08-28 12:32  snake_demo\snake\snakedeform2.m

............此处省略33个文件信息