【图像分割】基于snake模型的图像分割matlab 源码

一、简介

2 原理解释

二、源代码


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',15,sigma);
Img_smooth=conv2(Img,G,'same');  % smooth image by Gaussiin convolution
[Ix,Iy]=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 functiontimestep=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 L(\phi)
alf=1.5;  % coefficient of the weighted area term A(\phi);% Note: Choose smaller value for weak object bounday, such as the cell image in this demo.% define initial level set function (LSF) as -c, 0, c at points outside, on
% the boundary, and inside of a region R, respectively.
[nrow, ncol]=size(Img);
c0=4;
initialLSF=c0*ones(nrow,ncol);
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 -c inside of R, postive constant c outside of R.
u=initialLSF;
figure;imagesc(Img);colormap(gray);hold on;
[c,h] = contour(u,[0 0],'r');
title('Initial contour');% start level set evolution
for n=1:300u=EVOLUTION(u, g ,lambda, mu, alf, epsilon, timestep, 1);  pause(0.001);if mod(n,20)==0imagesc(Img);colormap(gray);hold on;[c,h] = contour(u,[0 0],'r'); iterNum=[num2str(n), ' iterations'];        title(iterNum);hold off;end
end
% define initial level set function (LSF) as -c, 0, c at points outside, on
% the boundary, and inside of a region R, respectively.
[nrow, ncol]=size(Img);
c=4;
initialLSF=c*ones(nrow,ncol);
w=10;
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)=-c; % negative constant -c inside of R, postive constant c outside of R.
u=initialLSF;
figure;imagesc(Img, [0, 255]);colormap(gray);hold on;
[c,h] = contour(u,[0 0],'r');
title('Initial contour');% start level set evolution
for n=1:500u=EVOLUTION(u, g ,lambda, mu, alf, epsilon, timestep, 1);  pause(0.001);if mod(n,20)==0imagesc(Img, [0, 255]);colormap(gray);hold on;[c,h] = contour(u,[0 0],'r'); iterNum=[num2str(n), ' iterations'];        title(iterNum);hold off;end
end

三、运行结果

【图像分割】基于snake模型的图像分割matlab 源码相关推荐

【预测模型】基于灰度模型预测房价matlab源码
1 简介自2006年开始,房价保持较长时间持续攀升,房价过高已成为一个全社会关注的民生问题.从市场价格内在发展趋势和外部因素影响两方面建立房价构成模型,将居民消费物价指数,房地产开发综合景气指数,土 ...
【数学建模】基于matlab船舶三自由度MMG模型【含Matlab源码 1925期】
⛄一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[数学建模]基于matlab船舶三自由度MMG模型[含Matlab源码 1925期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式 ...
qr-rls算法matlab实现,【预测模型】基于RLS算法进行预测matlab源码
一.简介 1 概述递归最小二乘(RLS)算法是一种典型的数据处理方法,由著名学者高斯在1795年提出,高斯认为,根据所获得的观测数据来推断未知参数时,未知参数最可能的值是这样一个数据,即它使各项实际 ...
【数学建模】基于matlab模拟疫情SEIRS模型【含Matlab源码 2214期】
一.⛄SEIR模型简介 1 SEIR模型简介如果所研究的传染病有一定的潜伏期,与病人接触过的健康人并不马上患病,而是成为病原体的携带者,归入 E 类.此时有: 仍有守恒关系 S(t) + E(t) ...
【元胞自动机】基于matlab元胞自动机3D森林火灾模型【含Matlab源码 656期】
⛄一.元胞自动机简介 1 元胞自动机发展历程最初的元胞自动机是由冯 · 诺依曼在 1950 年代为模拟生物细胞的自我复制而提出的. 但是并未受到学术界重视. 1970 年, 剑桥大学的约翰 · 何 ...
【元胞自动机】基于matlab元胞自动机模拟SEIR传播模型和采取隔离措施的SEIR模型【含Matlab源码 2181期】
⛄一.元胞自动机简介 1 元胞自动机发展历程最初的元胞自动机是由冯 · 诺依曼在 1950 年代为模拟生物细胞的自我复制而提出的. 但是并未受到学术界重视.1970 年, 剑桥大学的约翰 · 何顿 ...
【飞行器】基于matlab GUI四旋飞行器模型【含Matlab源码 2075期】
⛄一.四旋翼飞行器简介 0 引言四旋翼飞行器由于具有可垂直起降.机动性强.操作方便等诸多优点,在军事和民用场合得到广泛应用,从而成为众多学者的研究热点.四旋翼飞行器是具有四输入.六输出的欠驱动.非 ...
【图像识别】基于HSV和RGB模型水果分类matlab源码含 GUI
一.简介 1 RGB颜色空间: RGB(red,green,blue)颜色空间最常用的用途就是显示器系统(计算机.电视机等都是采用RGB颜色空间来进行图像显示).一般来说,电脑,电视机等是利用三个电子 ...
【元胞自动机】元胞自动机地铁火灾疏散模型【含Matlab源码 246期】
⛄一.元胞自动机简介 1 元胞自动机发展历程最初的元胞自动机是由冯 · 诺依曼在 1950 年代为模拟生物细胞的自我复制而提出的. 但是并未受到学术界重视. 1970 年, 剑桥大学的约翰 · 何 ...

【图像分割】基于snake模型的图像分割matlab 源码

一、简介

二、源代码

三、运行结果

【图像分割】基于snake模型的图像分割matlab 源码相关推荐

最新文章

热门文章