lbp特征提取算法

%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