function [ output ] = LBP( I )
%UNTITLED2 Summary of this function goes here
% Detailed explanation goes here
I = rgb2gray(I);
[m,n] = size(I);
count = (m-2)*(n-2); %count为需要计算LBP的像素个数。
P = 8;
%--------------------------------------------------------------------------
%取3*3的周长为8，半径为1的区域。左下角为第一个，顺时针方向依次排序。
%--------------------------------------------------------------------------
label = zeros(count,8);
k = 1;
for i=2:m-1
for j=2:n-1
label(k,1) = S(I(i+1,j-1),I(i,j));
label(k,2) = S(I(i,j-1),I(i,j));
label(k,3) = S(I(i-1,j-1),I(i,j));
label(k,4) = S(I(i-1,j),I(i,j));
label(k,5) = S(I(i-1,j+1),I(i,j));
label(k,6) = S(I(i,j+1),I(i,j));
label(k,7) = S(I(i+1,j+1),I(i,j));
label(k,8) = S(I(i+1,j),I(i,j));
k = k + 1;
end
end
%--------------------------------------------------------------------------
%计算Ulbp和countS
%--------------------------------------------------------------------------
Ulbp = zeros(count,1);
countS = zeros(count,1);
for i=1:count
Ulbp(i) = computUlbp(label(i,:));
countS(i) = sumS(label(i,:));
end
%--------------------------------------------------------------------------
%计算LBPpr
%--------------------------------------------------------------------------
LBPpr = zeros(count,1);
for i=1:count
if Ulbp(i)<=2 && countS(i)<P
LBPpr(i) = countS(i);
elseif Ulbp(i)<=2 && countS(i)==P
LBPpr(i) = 0;
else
LBPpr(i) = P;
end
end
%--------------------------------------------------------------------------
%计算H直方图
%--------------------------------------------------------------------------
H = zeros(count,P+1);
for i=1:count
H(i,:) = sigma(LBPpr(i),P);
end
temp = zeros(1,P+1);
for i=1:count
for j=1:P+1
temp(j) = temp(j) + H(i,j);
end
end
%--------------------------------------------------------------------------
%计算LBP纹理特征向量
%--------------------------------------------------------------------------
output = zeros(1,P+1);
for i=1:P+1
output(i) = temp(i)/count;
end
end
