matlab区域生长算法代码

MATLAB区域生长算法是一种基于像素相似性的图像分割方法。其基本思想是以一或多个像素为种子，通过逐步合并相邻像素，将一张图像分成若干个区域。常用的MATLAB区域生长算法代码如下：

function [segmented_img,labels] = region_growing(img,threshold)
% 初始化变量
[M,N] = size(img);
segmented_img = zeros(M,N);
labels = zeros(M,N);
label = 0;
% 对每个像素进行判断
for i = 1:M
    for j = 1:N
        % 如果该像素没有被标记，则表示为新的区域
        if (segmented_img(i,j) == 0)
            label = label + 1;
            segmented_img(i,j) = 1;
            labels(i,j) = label;
            % 递归遍历周围的像素
            segmented_img = region_growing_recursive(img,segmented_img,labels,i,j,threshold,label);
        end
    end
end

function segmented_img = region_growing_recursive(img,segmented_img,labels,i,j,threshold,label)
[M,N] = size(img);
% 遍历周围的像素
for x = max(i-1,1):min(i+1,M)
    for y = max(j-1,1):min(j+1,N)
        % 判断周围像素是否符合条件
        if (segmented_img(x,y) == 0 && abs(img(x,y)-img(i,j))<=threshold)
            segmented_img(x,y) = 1;
            labels(x,y) = label;
            % 递归遍历周围的像素
            segmented_img = region_growing_recursive(img,segmented_img,labels,x,y,threshold,label);
        end
    end
end

该代码实现了一个简单的区域生长算法，其输入参数为图像和阈值，输出参数为分割后的图像和标签。具体实现方式是遍历每个像素，对于未标记的像素，以该像素为种子递归遍历周围的像素，将与其相似的像素标记为同一区域。