如何实现统一LBP？

Question

可能的重复：
MATLAB 中的本地二进制模式

我想实现统一的 LBP。这是维基百科给出的统一 LBP 的定义。

如果二进制模式是一致的，则局部二进制模式被称为统一的 最多包含两个从 0 到 1 的按位转换，反之亦然 当位模式循环遍历时。例如， 模式 00000000（0 个转换）、01110000（2 个转换）和 11001111（2 个过渡）是均匀的，而图案 11001001（4 转换）和 01010010（6 个转换）不是。在计算中 在 LBP 标签中，使用统一的模式，以便有 每个统一图案和所有非统一图案都有单独的标签 图案用单个标签标记。例如，当使用 (8,R)邻域，共有256种模式，其中58种是 统一，产生 59 种不同的标签。

我已经为 LBP 编写了代码，但不确定如何将其转换为统一的 LBP。下面是 LBP 的代码。

for i=2:m-1
    for j=2:n-1
        J0=I2(i,j);
        I3(i-1,j-1)=I2(i-1,j-1)>J0;
        I3(i-1,j)=I2(i-1,j)>J0;
        I3(i-1,j+1)=I2(i-1,j+1)>J0; 
        I3(i,j+1)=I2(i,j+1)>J0;
        I3(i+1,j+1)=I2(i+1,j+1)>J0; 
        I3(i+1,j)=I2(i+1,j)>J0; 
        I3(i+1,j-1)=I2(i+1,j-1)>J0; 
        I3(i,j-1)=I2(i,j-1)>J0;
        LBP(i,j)=I3(i-1,j-1)*2^7+I3(i-1,j)*2^6+I3(i-1,j+1)*2^5+I3(i,j+1)*2^4+I3(i+1,j+1)*2^3+I3(i+1,j)*2^2+I3(i+1,j-1)*2^1+I3(i,j-1)*2^0;

    end
end
figure,imshow(uint8(LBP))

任何帮助将不胜感激。我正在使用 MATLAB。

Answer 1

步骤

1. 下一步是为存储在 LBP 中的值构建一个查找表。
   • 查找表将 256 种可能的组合映射到 59 个不同的标签。
   • 如果速度不重要，可以使用 for 循环构建表格。
2. 使用表格将 LBP 映射到 59 个标签。
   • labeled = table(LBP) % 这称为表查找或 MATLAB 索引。
3. 使用这 59 个标签执行任何其他工作。

---

建议（虽然对于实现来说不是必需的）

1. I3 不需要使用 2D 矩阵。当前像素的八个邻居是您正在处理的当前像素的本地像素；因此您可以简单地将它们分配给 I3(1), I3(2), ... I3(8) 因为每次您移动到下一个中​​心像素 ( i，j）。

---

function table = BitwiseToLBP
    % we reserve label 0 for non-uniform
    table = zeros(1, 256);
    nextLabel = 1;
    for k = 1:256,
        bits = bitand(k, 2.^(0:7)) > 0;
        if IsUniformLBP(bits),
            table(k) = nextLabel;
            nextLabel = nextLabel + 1;
        else
            table(k) = 0;
        end
    end
end

---

function IsUniformLBP(bits)
    nnz(diff(bits([1:end, 1]))) == 2;
end

Answer 2

我尝试了这段代码：

 sum = abs(I3(k-1,l-1)-I3(k-1,l))+ abs(I3(k-1,l)-I3(k-1,l+1))+ abs(I3(k-1,l+1)-I3(k,l+1))+ abs(I3(k,l+1)-I3(k+1,l+1))+abs(I3(k+1,l+1)-I3(k+1,l))+abs(I3(k+1,l)-I3(k+1,l-1))+abs(I3(k+1,l-1)-I3(k,l-1));
    if(sum <=2)
            UniformHist = [UniformHist Hist];
    end