<span style="font-size:18px;">I=imread('E:6.jpg');
I=rgb2gray(I);
[a b]=LDP(I,[1 1],2);</span>
<span style="font-size:18px;">% [HIST, LDPIM] = LDP(IMAGE, SUBSZ, ORDER)
% Calculates the LDP images and histogram of IMAGE. The concatenation
% of the LDP histograms of each subblock (of size SUBSZ) in the image
% is returned in HIST, and the LDP code images are returned in LDPIM.
% The LDP order is specified by ORDER.
% This code is based on B Zhang, Y Gao, S Zhao & J Liu, "Local Derivative Pattern Versus Local Binary Pattern: Face Recognition With High-Order Local Pattern Descriptor," Image Processing, IEEE Transactions on, vol.19, no.2, pp.533,544, Feb. 2010
function [hout,ldpim]=LDP(im, subsz, order)
numBinsPerBlock = 32; % number of bins per direction per bl