Question

我正在尝试实现本文中描述的光速标记技术的代码（我不能使用图像处理工具箱）：https://pdfs.semanticscholar.org/ef31/7c257603004d818ca1e2a2aa67d36d40147e.pdf（参见第7页第2节）。

这是我用于LSL等价构造的Matlab代码（算法14，步骤2）。

function [EQ,ERAi,nea] = LSL_equivalence(EQ,ERim1,RLCi,ERAim1,ERAi,NERi,nea,lImg)
    % LSL_EQUIVALENCE build the associative table between er and ea
    % GOAL: to create a Look Up Table to be applied to ERi to create EAi.

    for er = 1:2:NERi % check segments one by one
        % read the boundaries of each segment to obtain de relative labels of every agacent segment in the prev line
        j0 = RLCi(er); 
        j1 = RLCi(er+1);

        er0 = ERim1(j0+1); % label of the first segment
        er1 = ERim1(j1+1); % the label of the last segment

        % check label parity: segments are odd, background is even
        % bitand([1 2 3 4 5],1) == [1 0 1 0 1]
        if bitand(er0,1) == 0 % if er0 is even
            er0 = er0 + 1;
        end
        if bitand(er1,1) == 0 % if er1 is even
            er1 = er1 -1;
        end
        if er1 >= er0 % if there is an adjacency
            ea = ERAim1(er0+1); % absolute label of the first segment 
            a = EQ(ea+1);     % a is the ancestor (smallest label of the equivalence class) 
                for erk = (er0+2):2:er1
                    eak = ERAim1(erk+1);
                    ak = EQ(eak+1);
                    % min extraction and propagation
                    if a < ak
                        EQ(eak+1) = a;
                    else
                        a = ak;
                        EQ(ea+1) = a;
                        ea = eak;
                    end
                end
            ERAi(er+1) = a; % the global min of all ak ancestors 
        else % if there are no adjacent labels make a new label
            nea = nea + 1;
            ERAi(er+1) = nea;
        end
    end
    end

我在使用索引时遇到了一些问题，因为本文中描述的伪代码的索引从0开始，Matlab与1一起使用。我已经在此Stack Overflow帖子Implementing LSL for Connected Component Labeling/Blob Extraction中找到了一些C ++代码（我申请了）建议的更改）以及此git仓库https://github.com/prittt/YACCLAB/blob/master/include/labeling_lacassagne_2016_code.inc。我没有看到差异。

另外，我在理解等价类是什么（矩阵EQ中的内容）时遇到了一些麻烦。提前谢谢！

Answer 1

我意识到这有点迟了，但是我只写了一段类似于光速标记算法的代码。我将简要说明如何解决索引问题。

LSL的第一步是获取像素的每一列，并找到连续标记的像素的开始和结束位置。您可以在Matlab中这样做：

I = imread('text.png');
[rows,cols] = find(xor(I(2:end,:),I(1:end-1,:)));

这为您提供的是列中每行像素的开始和结束位置的行和列，但不包括索引。非包容索引是指每个像素行（n）的像素索引从I(r(2*n-1),c(2*n-1))到I(r(2*n)-1,c(2*n))。您应注意，本文沿行操作，而上面的代码沿列操作，但原理相同。您还应该注意，以上代码并未涵盖图像边缘上带有标记像素的情况。

如果您想查看完整的实现，我将代码发布在Matlab文件交换上。我并不是说它可以完全复制LSL，但可以在许多相同的原理上工作。

https://www.mathworks.com/matlabcentral/fileexchange/70323-ftllabel?s_tid=prof_contriblnk

Matlab实现光速标签

1 个答案: