MATLAB实现区域生长算法详解

>> 
I = imread('E:\图像处理\区域生长\区域生长.png');
if( ~( size(I,3)-3 ))%  ~x,表示对表达式取反
    I = rgb2gray(I);%转化为单通道灰度图
end
I = im2double(I); %图像灰度值归一化到[0,1]之间
Ireshape = imresize(I,[600,800]);%将i图像缩放为600*800大小图像
I = Ireshape(51:475,200:699);
gausFilter = fspecial('gaussian',[5 5],0.5);%建立高斯低通滤波掩模，大小5*5，标准差0.5
I = imfilter(I,gausFilter,'replicate');%线性滤波 使用之前建立的掩模，这样可以降低噪声
 
%种子点的交互式选择
if( exist('x','var') == 0 && exist('y','var') == 0)%exist('x','var')若x变量存在返回true，否则返回false
    subplot(2,2,1),imshow(I,[]);%显示原图
    hold on;
    [y,x] = getpts;%鼠标取点  回车确定
    x = round(x(1));%选择种子点round()为四舍五入取整
    y = round(y(1));
end
if( nargin == 0)
    reg_maxdist = 0.1;
    %nargin是matlab代码编写中常用的一个技巧，主要用于计算当前主函数的输入参数个
    %数，一般可以根据nargin的返回值来确定主函数输入参数的缺省值。在实现中，如果
    %用户输入的参数个数为零，那么默认为0.2
end
J = zeros(size(I)); % 主函数的返回值，记录区域生长所得到的区域
Isizes = size(I);%size(I)返回一个2*1的矩阵存放行数、列数
reg_mean = I(x,y);%表示分割好的区域内的平均值，初始化为种子点的灰度值
reg_size = 1;%分割的到的区域，初始化只有种子点一个
neg_free = 10000; %动态分配内存的时候每次申请的连续空间大小