-1

我在将二进制转换为十进制时遇到问题（看起来很长）。

%% Read

clear all;
close all;
clc;
I=imread('test.png');
imshow(I);

%% Crop
I2 = imcrop(I);
figure, imshow(I2)
w=size(I2,1);
h=size(I2,2);
%% LBP
for i=2:w-1
    for j=2:h-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^8+I3(i-1,j)*2^7+I3(i-1,j+1)*2^6+I3(i,j+1)*2^5+I3(i+1,j+1)*2^4+I3(i+1,j)*2^3+I3(i+1,j-1)*2^2+I3(i,j-1)*2^1;
    end
end
figure,imshow(I3)
figure,imhist(LBP)

是否可以更改此行

LBP(i,j)=I3(i-1,j-1)*2^8+I3(i-1,j)*2^7+I3(i-1,j+1)*2^6+I3(i,j+1)*2^5+I3(i+1,j+1)*2^4+I3(i+1,j)*2^3+I3(i+1,j-1)*2^2+I3(i,j-1)*2^1;

更短的东西？

4

3 回答 3

2

选项1：

简化您正在做的事情的一种方法是首先创建一个 3×3 比例因子矩阵（即 2 的幂）并在循环之前对其进行初始化：

scale = 2.^[8 7 6; 1 -inf 5; 2 3 4];

然后，您可以用这些矢量化操作替换循环内的所有内容：

temp = (I2(i-1:i+1,j-1:j+1) > I2(i,j)).*scale;
LBP(i,j) = sum(temp(:));

选项 2：

或者，我相信您可以完全删除两个循环并用对NLFILTER的一次调用替换它们以获取您的矩阵LBP：

LBP = nlfilter(I2,[3 3],@(x) sum((x(:) > x(5)).*scale(:)));
于 2011-09-12T15:53:33.930 回答
1

我不完全确定你在那里bin2dec做什么，但做你想做的事吗？

于 2011-09-12T15:32:26.040 回答
0

我正在使用COLFILT函数添加另一个解决方案。

它涉及将所有滑动窗口放入矩阵的列中，我们使用自定义函数对其进行处理，然后将结果重新排列到原始矩阵中。它在内部使用IM2COL和COL2IM函数。

这是一个例子：

I = imread('coins.png');
fun = @(b) sum( bsxfun(@times, 2.^(8:-1:1)', ...
           bsxfun(@gt, b([1 4 7 8 9 6 3 2],:), b(5,:))) );
II = colfilt(I, [3 3], 'sliding', fun);
imshow(II, [])

截屏

与@gnovice的第二个答案相比，请从NLFILTER文档中阅读此提示：

nlfilter 处理大图像可能需要很长时间。在某些情况下，colfilt 函数可以更快地执行相同的操作。

于 2011-09-13T01:41:29.667 回答