%LBP?returns?the?local?binary?pattern?image?or?LBP?histogram?of?an?image.
%??J?=?LBP（IRNMAPPINGMODE）?returns?either?a?local?binary?pattern
%??coded?image?or?the?local?binary?pattern?histogram?of?an?intensity
%??image?I.?The?LBP?codes?are?computed?using?N?sampling?points?on?a?
%??circle?of?radius?R?and?using?mapping?table?defined?by?MAPPING.?
%??See?the?getmapping?function?for?different?mappings?and?use?0?for
%??no?mapping.?Possible?values?for?MODE?are
%???????‘h‘?or?‘hist‘??to?get?a?histogram?of?LBP?codes
%???????‘nh‘???????????to?get?a?normalized?histogram
%??Otherwise?an?LBP?code?image?is?returned.
%
%??J?=?LBP（I）?returns?the?original?（basic）?LBP?histogram?of?image?I
%
%??J?=?LBP（ISPMAPPINGMODE）?computes?the?LBP?codes?using?n?sampling
%??points?defined?in?（n?*?2）?matrix?SP.?The?sampling?points?should?be
%??defined?around?the?origin?（coordinates?（00））.
%
%??Examples
%??--------
%???????I=imread（‘rice.png‘）;
%???????mapping=getmapping（8‘u2‘）;?
%???????H1=lbp（I18mapping‘h‘）;?%LBP?histogram?in?（81）?neighborhood
%??????????????????????????????????using?uniform?patterns
%???????subplot（211）stem（H1）;
%?
%???????H2=lbp（I）;
%???????subplot（212）stem（H2）;
%?
%???????SP=[-1?-1;?-1?0;?-1?1;?0?-1;?-0?1;?1?-1;?1?0;?1?1];
%???????I2=LBP（ISP0‘i‘）;?%LBP?code?image?using?sampling?points?in?SP
%???????????????????????????and?no?mapping.?Now?H2?is?equal?to?histogram
%???????????????????????????of?I2.

function?result?=?lbp（varargin）?%?imageradiusneighborsmappingmode）
%?Version?0.3.3
%?Authors:?Marko?Heikkil?and?Timo?Ahonen

%?Changelog
%?Version?0.3.2:?A?bug?fix?to?enable?using?mappings?together?with?a
%?predefined?spoints?array
%?Version?0.3.1:?Changed?MAPPING?input?to?be?a?struct?containing?the?mapping
%?table?and?the?number?of?bins?to?make?the?function?run?faster?with?high?number
%?of?sampling?points.?Lauge?Sorensen?is?acknowledged?for?spotting?this?problem.


%?Check?number?of?input?arguments.
error（nargchk（15nargin））;

image=varargin{1};
d_image=double（image）;

if?nargin==1
????spoints=[-1?-1;?-1?0;?-1?1;?0?-1;?-0?1;?1?-1;?1?0;?1?1];
????neighbors=8;
????mapping=0;
????mode=‘h‘;
end

if?（nargin?==?2）?&&?（length（varargin{2}）?==?1）
????error（‘Input?arguments‘）;
end

if?（nargin?>?2）?&&?（length（varargin{2}）?==?1）
????radius=varargin{2};
????neighbors=varargin{3};
????
????spoints=zeros（neighbors2）;

????%?Angle?step.
????a?=?2*pi/neighbors;
????
????for?i?=?1:neighbors
????????spoints（i1）?=?-radius*sin（（i-1）*a）;
????????spoints（i2）?=?radius*cos（（i-1）*a）;
????end
????
????if（nargin?>=?4）
????????mapping=varargin{4};
????????if（isstruct（mapping）?&&?mapping.samples?~=?neighbors）
????????????error（‘Incompatible?mapping‘）;
????????end
????else
????????mapping=0;
????end
????
????if（nargin?>=?5）
????????mode=varargin{5};
????else
????????mode=‘h‘;
????end
end

if?（nargin?>?1）?&&?（length（varargin{2}）?>?1）
????spoints=varargin{2};
????neighbors=siz