区域生长法matlab代码

以下是一个简单的区域生长法 Matlab 代码示例：

function [segmImage] = regionGrow(inputImage, row, col, threshold)
% inputImage: 输入图像
% row, col: 种子点的行列坐标
% threshold: 阈值

[rows,cols] = size(inputImage);
visited = zeros(rows,cols);

% 设置种子点为已访问并添加到队列中
visited(row,col) = 1;
queue = [row, col];

% 区域生长
while ~isempty(queue)
    % 取出队首点
    currentPixel = queue(1,:);
    row = currentPixel(1);
    col = currentPixel(2);
    queue(1,:) = [];
    
    % 处理8邻域
    for i = -1:1
        for j = -1:1
            % 跳过中心点和越界的点
            if (i == 0 && j == 0) || row+i < 1 || col+j < 1 || row+i > rows || col+j > cols
                continue;
            end
            
            % 如果点没有被访问过，并且与当前点的差值小于阈值，则将其添加到队列中
            if visited(row+i, col+j) == 0 && abs(inputImage(row+i, col+j) - inputImage(row, col)) < threshold
                visited(row+i, col+j) = 1;
                queue = [queue; row+i, col+j];
            end
        end
    end
end

% 构建分割图像
segmImage = visited;
end

在主函数中，可以通过调用 `regionGrow` 函数来实现区域生长。例如：

inputImage = imread('input.png');
seedRow = 100;
seedCol = 150;
threshold = 10;
segmImage = regionGrow(inputImage, seedRow, seedCol, threshold);
imshow(segmImage);

这个例子将在 `input.png` 中以 `(100,150)` 为种子点，以阈值 `10` 进行区域生长，并将结果显示出来。