Compute gradient vector flow.

snake模型的改进算法，本程序相对于opencv自带的cvSnakeImage能够得到更好的获取轮廓的效果，但是经测试本程序存在着Bug，在稳定性方面不如cvSnakeImage，仅仅用于测试效果还是可行的，如果谁能够去掉这个程序中的Bug，欢迎联系我。

GVF-Snake代码,里面一份是C++写的，需要opencv库支持才能运行，另一份是matlab代码可以通过matlab直接运行

GVF snake模型的matlab代码

这是snake模型中GVF模型的matlab源代码,希望对大家有帮助! 这是snake模型中GVF模型的matlab源代码,希望对大家有帮助!

实现snake算法的源程序，对源程序的英文注释做了翻译。 程序同时实现了GVFsnake算法，同时加入了气球力对snake算法进行改进。

能够实现，希望谁能在此基础上添加一个有效外力。是的算法的速度更快。

基于GVF的snake算法模型，用matlab实现。有例子和源代码及demo

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

使用GVF域和VFC域进行图片分割 % Vector field convolution (VFC) external force field example. % % See also AMT, EXAMPLE_PIG, AM_VFC, AM_VFK, AC_DISPLAY. % % Reference % [1] Bing Li and Scott T. Acton, "Active contour external force using % vector field convolution for image segmentation," Image Processing, % IEEE Trans. on, vol. 16, pp. 2096-2106, 2007. % [2] Bing Li and Scott T. Acton, "Automatic Active Model % Initialization via Poisson Inverse Gradient," Image Processing, % IEEE Trans. on, vol. 17, pp. 1406-1420, 2008. % % (c) Copyright Bing Li 2005 - 2009. clear all disp('======================================') disp('Vector field convolution (VFC) example') %% parameter settings disp('Initializing parameters ...') SAVE_AVI = 0; % set it to 1 if you want to save the process as .avi movie DISPLAY_STREAMLINE = 0; % set it to 1 if you want to plot streamlines, note that it takes a while mu = .2; GVF_ITER = 100; normalize = 1; alpha = .5; beta = 0; tau = .5; SNAKE_ITER = 5; SNAKE_ITER1 = 60; RES = .5; clr = {'b' 'b' 'r'}; %% Read images disp('Reading images ...') U = imread('im_U.bmp'); noisyU=imread('im_Unoisy.bmp'); figure(1) %% compare 3 different cases for cs = 1:3, %% compute external force fields switch cs, case 1, % traditional GVF with Gaussian filter disp('--------------------------------------------------') disp('Case 1: GVF snake with initial circle close to FOI') disp('Computing the external force field ...') h = fspecial('gaussian',[5 5],5); f = imfilter(double(noisyU),h); titl = 'GVF'; Fext = AM_GVF(f, mu, GVF_ITER, normalize); R = 20; case 2, % traditional GVF with Gaussian filter disp('--------------------------------------------------') disp('Case 2: GVF snake with initial circle far away from FOI') disp('Computing the external force field ...