GVF Snake（Matlab）

有用户界面的，经典的源码。可用来学习，处理一些图像，不过有需要改进的地方，可能对有些图像就会出现问题。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 
% 		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  head.jpg
     目录           0  1999-01-06 13:49  images\
     文件       47490  1999-01-06 13:49  images\chest.pgm
     文件        1064  1998-09-30 09:50  images\heart.mat
     文件       25615  1998-09-29 10:01  images\heart.pgm
     文件        4248  1998-09-30 09:58  images\heart2.mat
     文件        4056  1998-10-01 06:18  images\heart3.mat
     文件       65594  1999-01-06 13:32  images\new.pgm
     文件        2792  1998-09-30 10:26  images\room.mat
     文件        4109  1997-08-28 12:32  images\room.pgm
     文件        2928  1998-09-30 09:43  images\room2.mat
     文件        2928  1998-09-30 09:44  images\room3.mat
     文件        3304  1998-10-01 06:20  images\u64.mat
     文件        4109  1997-08-28 12:32  images\U64.pgm
     文件        2952  1998-10-07 04:51  images\u642.mat
     文件        3800  1998-10-06 09:00  sdemo.m
     目录           0  2000-04-12 20:18  snake\
     文件         916  2000-04-12 20:06  snake\BoundMirrorEnsure.m
     文件         699  2000-04-12 20:06  snake\BoundMirrorExpand.m
     文件         488  2000-04-12 20:06  snake\BoundMirrorShrink.m
     文件        1096  1997-08-28 12:32  snake\Contents.m
     文件         608  1997-08-28 12:32  snake\dt.m
     文件         500  1997-08-28 12:32  snake\gaussianBlur.m
     文件         255  1997-08-28 12:32  snake\gaussianMask.m
     文件         895  1998-10-08 10:23  snake\GGVF.m
     文件        1640  2000-04-12 20:06  snake\GVF.m
     文件         257  1998-09-30 06:08  snake\imdisp.m
     文件        3256  1997-08-28 12:32  snake\rawread.m
     文件        1536  1998-10-07 05:15  snake\rawwrite.m
     文件        1504  1998-09-21 05:42  snake\snakedeform.m
     文件        1856  1997-08-28 12:32  snake\snakedeform2.m
............此处省略23个文件信息