帮我写一段区域生长的matlab代码

好的，下面是一段MATLAB代码实现区域生长：

function [segImg] = regionGrowing(img, seed, threshold)
% img: 待分割的图像
% seed: 种子点坐标[x,y]
% threshold: 生长的阈值

% 获取图像尺寸
[row, col] = size(img);

% 初始化分割图像
segImg = zeros(row, col);
segImg(seed(1), seed(2)) = 1;

% 初始化生长队列
queue = [seed(1), seed(2)];

% 逐步生长
while ~isempty(queue)
    % 取出队首像素
    curPix = queue(1,:);
    queue(1,:) = [];
    
    % 获取当前像素邻域内所有未被标记的像素
    if curPix(1) > 1 && segImg(curPix(1)-1,curPix(2)) == 0 % 上
        if abs(img(curPix(1),curPix(2)) - img(curPix(1)-1,curPix(2))) <= threshold
            segImg(curPix(1)-1,curPix(2)) = 1;
            queue = [queue;curPix(1)-1,curPix(2)];
        end
    end
    if curPix(1) < row && segImg(curPix(1)+1,curPix(2)) == 0 % 下
        if abs(img(curPix(1),curPix(2)) - img(curPix(1)+1,curPix(2))) <= threshold
            segImg(curPix(1)+1,curPix(2)) = 1;
            queue = [queue;curPix(1)+1,curPix(2)];
        end
    end
    if curPix(2) > 1 && segImg(curPix(1),curPix(2)-1) == 0 % 左
        if abs(img(curPix(1),curPix(2)) - img(curPix(1),curPix(2)-1)) <= threshold
            segImg(curPix(1),curPix(2)-1) = 1;
            queue = [queue;curPix(1),curPix(2)-1];
        end
    end
    if curPix(2) < col && segImg(curPix(1),curPix(2)+1) == 0 % 右
        if abs(img(curPix(1),curPix(2)) - img(curPix(1),curPix(2)+1)) <= threshold
            segImg(curPix(1),curPix(2)+1) = 1;
            queue = [queue;curPix(1),curPix(2)+1];
        end
    end
end

该代码实现了基于种子点的区域生长算法，输入待分割的图像、种子点坐标和生长阈值，输出分割结果。具体实现过程为：首先初始化分割图像，将种子点标记为1，初始化生长队列，将种子点加入队列。然后逐步生长，从队首取出一个像素，依次检查其上、下、左、右邻域内所有未被标记的像素，如果像素值与当前像素值的差小于生长阈值，则将其标记为1并加入生长队列。最终输出分割图像。