www.gusucode.com > 《matlab图像处理与界面编程宝典》秦襄培 编著，每章的MATLAB源代码程序 > 第23章/代码23-6.txt

    
f1=imread('peppers.png');                 % 读入图像
f1=rgb2gray(f1);                          % 图像灰度转换
f=double(f1);                             % 数据类型转换
s=255;
t=65;  
if numel(s)==1
    si=f==s;
    s1=s;
else
    si=bwmorph(s,'shrink',Inf);
    j=find(si);
    s1=f(j);
end
ti=false(size(f));
for k=1:length(s1)
    sv=s1(k);
    s=abs(f-sv)<=t;
    ti=ti|s;
end
[g nr]=bwlabel(imreconstruct(si,ti));      % 图像标识
figure,imshow(f1);                         % 显示原始图像 
title('Original Image');                   % 设置图像标题
figure,imshow(g);                          % 显示处理后的图像
title('Segmented Image - Region Growing'); % 设置图像标题
display('No. of regions');                 % 显示文字
nr