基于MATLAB的游程平滑算法(RLSA)

function smoothedImage = rlsa(image, minLength, maxLength)
    % image: 输入的二值图像
    % minLength: 合并游程的最小长度
    % maxLength: 分割游程的最大长度
    
    % 获取图像尺寸
    [height, width] = size(image);
    
    % 初始化平滑后的图像
    smoothedImage = false(height, width);
    
    % 遍历图像的每一行
    for i = 1:height
        runStart = -1;
        runLength = 0;
        
        for j = 1:width
            if image(i, j) ~= runStart
                % 开始一个新的游程
                runStart = image(i, j);
                runLength = 1;
            else
                % 继续当前游程
                runLength = runLength + 1;
            end
            
            % 根据游程长度阈值进行处理
            if (runStart == 1 && runLength >= minLength) || (runStart == 0 && runLength <= maxLength)
                smoothedImage(i, j - runLength + 1:j) = runStart;
            end
        end
    end
end