%?[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.533544?Feb.?2010
function?[houtldpim]=LDP（im?subsz?order）
numBinsPerBlock?=?32;???????%?number?of?bins?per?direction?per?block

assert（order?>=2?&&?order?<=?10）;
assert（all（floor（size（im）?./?subsz（:）‘）?==?（size（im）?./?subsz（:）‘）））;

im?=?double（im）;
der0?=?im;