数字图像处理-图像分割：Snake主动轮廓模型 Matlab代码及运行结果

% =========================================================================
%                   Snakes：Active Contour Models
% =========================================================================
% By gujinjin 2012/12/10-12  Sunny
% 基于KASS等的论文思想
% 参考文献：
% [1] KASS etc.Snakes：Active Contour Models
% [2] CSDN 博客 - Author：乐不思蜀Tone
% [3] Ritwik Kumar(Harvard University),D.Kroon（Twente University)的程序
% [4] 《数学建模与数学实验》
% ------
clc;clf;clear all;

% =========================================================================
%                      获取手动取点坐标
% =========================================================================
% 读取显示图像
%I = imread('Coronary_MPR.jpg');
I =  imread('test1.png');
% 转化为双精度型
%I = im2double(I); 
% 若为彩色，转化为灰度
%if(size(I,3)==3), I=rgb2gray(I); end
figure(1),imshow(I);
%---------------------------
%        高斯滤波
%---------------------------
sigma=1;
% 创建特定形式的二维高斯滤波器H
H = fspecial('gaussian',ceil(3*sigma), sigma);
% 对图像进行高斯滤波,返回和I等大小矩阵
Igs = filter2(H,I,'same');
%---------------------------
%      获取Snake的点坐标
%---------------------------
figure(2),imshow(Igs);
x=[];y=[];c=1;N=100; %定义取点个数c,上限N
% 获取User手动取点的坐标
% [x,y]=getpts
while c