zhang-suen细化算法matlab代码

以下是张氏细化算法的Matlab代码示例：

function img_thin = zhangSuenThinning(img)
    img = imbinarize(img); % 将图像二值化
    img_thin = img;

    [rows, cols] = size(img);
    has_change = true;

    while has_change
        has_change = false;
        
        % 第一次迭代
        for i = 2:rows-1
            for j = 2:cols-1
                if img_thin(i, j) == 0 % 当前像素为黑色
                    continue;
                end
                
                neighbors = [img_thin(i-1, j-1), img_thin(i-1, j), img_thin(i-1, j+1), ...
                             img_thin(i, j+1), img_thin(i+1, j+1), img_thin(i+1, j), ...
                             img_thin(i+1, j-1), img_thin(i, j-1)];

                transitions = 0;
                for k = 2:numel(neighbors)
                    if neighbors(k-1) == 0 && neighbors(k) == 1
                        transitions = transitions + 1;
                    end
                end

                if transitions == 1 && sum(neighbors) >= 2 && sum(neighbors) <= 6
                    if neighbors(2) * neighbors(4) * neighbors(6) == 0 && ...
                       neighbors(4) * neighbors(6) * neighbors(8) == 0
                        img_thin(i, j) = 0;
                        has_change = true;
                    end
                end
            end
        end
        
        % 第二次迭代
        for i = 2:rows-1
            for j = 2:cols-1
                if img_thin(i, j) == 0 % 当前像素为黑色
                    continue;
                end
                
                neighbors = [img_thin(i-1, j-1), img_thin(i-1, j), img_thin(i-1, j+1), ...
                             img_thin(i, j+1), img_thin(i+1, j+1), img_thin(i+1, j), ...
                             img_thin(i+1, j-1), img_thin(i, j-1)];

                transitions = 0;
                for k = 2:numel(neighbors)
                    if neighbors(k-1) == 0 && neighbors(k) == 1
                        transitions = transitions + 1;
                    end
                end

                if transitions == 1 && sum(neighbors) >= 2 && sum(neighbors) <= 6
                    if neighbors(2) * neighbors(4) * neighbors(8) == 0 && ...
                       neighbors(2) * neighbors(6) * neighbors(8) == 0
                        img_thin(i, j) = 0;
                        has_change = true;
                    end
                end
            end
        end
    end
end

请注意，此代码只是一个简单的示例，可能无法处理所有情况。实际使用时，请根据您的需求进行适当修改和调整。