snake活动轮廓边界提取源码

matlab轮廓提取代码，主要用于医学图像轮廓的提取。

function B = BoundMirrorEnsure（A）
% Ensure确保 mirror镜像 boundary边界 condition条件
%
% The number of rows and columns of A must be greater than 2
%
% for example （X means value that is not of interest）
% 
% A = [
%     X  X  X  X  X   X
%     X  1  2  3  11  X
%     X  4  5  6  12  X 
%     X  7  8  9  13  X 
%     X  X  X  X  X   X
%     ]
%
% B = BoundMirrorEnsure（A） will yield
%
%     5  4  5  6  12  6
%     2  1  2  3  11  3
%     5  4  5  6  12  6 
%     8  7  8  9  13  9 
%     5  4  5  6  12  6
%

% Chenyang Xu and Jerry L. Prince 9/9/1999
% http://iacl.ece.jhu.edu/projects/gvf

[mn] = size（A）;

if （m<3 | n<3） 
    error（‘either the number of rows or columns is smaller than 3‘）;
end

yi = 2:m-1;
xi = 2:n-1;
B = A;
B（[1 m][1 n]） = B（[3 m-2][3 n-2]）;  % mirror corners
B（[1 m]xi） = B（[3 m-2]xi）;          % mirror left and right boundary
B（yi[1 n]） = B（yi[3 n-2]）;          % mirror top and bottom boundary

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

     文件        932  2007-09-18 22:31  snake\BoundMirrorEnsure.m

     文件        703  2007-09-18 22:33  snake\BoundMirrorExpand.m

     文件        492  2007-09-18 22:33  snake\BoundMirrorShrink.m

     文件       1252  2007-09-18 22:42  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

     文件        268  2007-09-18 22:44  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

     文件        551  1997-08-28 12:32  snake\snakedisp.m

     文件        194  2007-09-16 14:53  snake\snakeindex.m

     文件        894  1997-08-28 12:32  snake\snakeinit.m

     文件       1593  1997-08-28 12:32  snake\snakeinterp.m

     文件       1310  1997-08-28 12:32  snake\snakeinterp1.m

     文件        913  1997-08-28 12:32  snake\xconv2.m

    ..AD...         0  2007-09-24 10:32  snake

----------- ---------  ---------- -----  ----

                21152                    21