for i = 1:row_max
for j = 1:col_max
A = Input_Im(i:i+L-1, j:j+L-1);
A = A+1-A(C,C);
A(A>0) = 1;
LBP_Im(i,j) = A(C,L) + A(L,L)*2 + A(L,C)*4 + A(L,1)*8 + A(C,1)*16 + A(1,1)*32 + A(1,C)*64 + A(1,L)*128;
block_index=fix((i-1)/20)*6+fix((j-1)/20)+1;
bin_index=fix(LBP_Im(i,j)/16)+1;

hist(block_index,bin_index)=hist(block_index,bin_index)+1;

    end;
end;

https://zhidao.baidu.com/question/499581164397800164.html