function LBPhistogram = FnLBP(img, R, P, numRowPatch, numColPatch)
LBPimg = zeros(size(img, 1)-2*R, size(img,2)-2*R);
% The Angle of Neighbour Sampling Points
theta = ((P:-1:1)*(2*pi/P))+(pi/2);
% The Amount of Displacements Required to Find the Neighbour Sampling Points
dxy = round([-R*sin(theta); R*cos(theta)]);
o = img(R+1:size(img,1)-R, R+1:size(img,2)-R);
for i = 1 : P
%
compMat = img(R+dxy(1,i)+1:size(img,1)-R+dxy(1,i), R+dxy(2,i)+1:size(img,2)-R+dxy(2,i)) - o;
LBPimg = (compMat >= 0) * (2^(P-i)) + LBPimg;
end
%%
rowPatchSize = ceil(size(img,1)/numRowPatch);
colPatchSize = ceil(size(img,2)/numColPatch);
LBPimg = padarray(LBPimg,[R+rowPatchSize*numRowPatch-size(img,1) R+(colPatchSize*numColPatch-size(img,2))/2], 256);
LBPimgCells = mat2cell(LBPimg, rowPatchSize * ones(1, numRowPatch), colPatchSize * ones(1, numColPatch));
bins = [0,1,2,3,4,6,7,8,12,14,15,16,24,28,30,31,32,48,56,60,62,63,64,96,112,120,124,126,127,128,129,131,135,143,159,191,192,193,195,199,207,223,224,225,227,231,239,240,241,243,247,248,249,251,252,253,254,255,256];
LBPhistogram = cellfun(@(x) imhistogram(x, bins), LBPimgCells, 'UniformOutput',0);
LBPfeats = cell2mat(cellfun(@(x) x(1:58), LBPhistogram', 'UniformOutput',0));
LBPhistogram = LBPfeats(:)';
