区域生长算法matlab代码

区域生长算法是一种基于像素相似性的图像分割算法，它通过选择种子点并根据一定的相似性准则将相邻像素逐步合并成为一个区域。以下是一个简单的区域生长算法的MATLAB代码示例：

function segmented_image = region_growing(image, seed_point, threshold)
    [rows, cols] = size(image);
    visited = false(rows, cols);
    segmented_image = zeros(rows, cols);
    
    % 获取种子点的像素值
    seed_value = image(seed_point(1), seed_point(2));
    
    % 定义种子点队列
    queue = [];
    queue = [queue; seed_point];
    
    % 区域生长
    while ~isempty(queue)
        % 取出队列中的第一个点
        current_point = queue(1,:);
        queue(1,:) = [];
        
        % 检查当前点是否已经访问过
        if visited(current_point(1), current_point(2))
            continue;
        end
        
        % 检查当前点与种子点的像素值差异是否小于阈值
        if abs(image(current_point(1), current_point(2)) - seed_value) <= threshold
            % 将当前点标记为已访问，并将其加入到分割图像中
            visited(current_point(1), current_point(2)) = true;
            segmented_image(current_point(1), current_point(2)) = 1;
            
            % 将当前点的相邻点加入到队列中
            neighbors = get_neighbors(current_point, rows, cols);
            queue = [queue; neighbors];
        end
    end
end

function neighbors = get_neighbors(point, rows, cols)
    % 获取当前点的相邻点
    x = point(1);
    y = point(2);
    
    neighbors = [x-1, y; x+1, y; x, y-1; x, y+1];
    
    % 去除超出图像范围的点
    out_of_range = neighbors(:,1) < 1 | neighbors(:,1) > rows | neighbors(:,2) < 1 | neighbors(:,2) > cols;
    neighbors(out_of_range,:) = [];
end

使用该代码，你可以将输入的图像进行区域生长分割。其中，`image`是输入的灰度图像，`seed_point`是种子点的坐标，`threshold`是像素值差异的阈值。函数返回的`segmented_image`是分割后的二值图像，其中区域内的像素值为1，其他区域的像素值为0。