1 #include 
     
       2 #include 
      
        3 #include "math.h" 4  5 using namespace cv; 6 using namespace std; 7  8 Mat src, gray_src; 9 10 const char* output_tt = "LBP Result";11 12 int main(int argc, char** argv) {13     src = imread("test.jpg");14     if (src.empty()) {15         printf("could not load image...\n");16         return -1;17     }18     19     namedWindow("input image", CV_WINDOW_AUTOSIZE);20     namedWindow(output_tt, CV_WINDOW_AUTOSIZE);21     imshow("input image", src);22 23     // convert to gray24     cvtColor(src, gray_src, COLOR_BGR2GRAY);25     int width = gray_src.cols;26     int height = gray_src.rows;27 28     // 基本LBP演示29     Mat lbpImage = Mat::zeros(gray_src.rows - 2, gray_src.cols - 2, CV_8UC1);//3*3窗口，边界有左右各1个像素不处理30     for (int row = 1; row < height - 1; row++) {31         for (int col = 1; col < width - 1; col++) {32             uchar c = gray_src.at
       
        (row, col);//获取中心像素值33             uchar code = 0;//特征码34             code |= (gray_src.at
        
         (row - 1, col - 1) > c) << 7;35             code |= (gray_src.at
         
          (row - 1, col) > c) << 6;36             code |= (gray_src.at
          
           (row - 1, col + 1) > c) << 5;37 code |= (gray_src.at
           
            (row, col + 1) > c) << 4;38 code |= (gray_src.at
            
             (row + 1, col + 1) > c) << 3;39 code |= (gray_src.at
             
              (row + 1, col) > c) << 2;40 code |= (gray_src.at
              
               (row + 1, col - 1) > c) << 1;41 code |= (gray_src.at
               
                (row, col - 1) > c) << 0;42 lbpImage.at
                
                 (row - 1, col - 1) = code;//原图1，1=效果图0，043 }44 }45 imshow(output_tt, lbpImage);46 47 waitKey(0);48 return 0;49 }
                
               
              
             
            
           
          
         
        
       
      
     

LBP的扩展（多尺度表达）：

1 #include 
     
       2 #include 
      
        3 #include "math.h" 4  5 using namespace cv; 6 using namespace std; 7  8 Mat src, gray_src; 9 int current_radius = 3;10 int max_count = 10;11 12 const char* output_tt = "LBP Result";13 14 void ELBP_Demo(int, void*);15 int main(int argc, char** argv) {16     src = imread("test.jpg");17     if (src.empty()) {18         printf("could not load image...\n");19         return -1;20     }21     22     namedWindow("input image", CV_WINDOW_AUTOSIZE);23     namedWindow(output_tt, CV_WINDOW_AUTOSIZE);24     imshow("input image", src);25 26     // convert to gray27     cvtColor(src, gray_src, COLOR_BGR2GRAY);28     int width = gray_src.cols;29     int height = gray_src.rows;30 31     // 基本LBP演示32     Mat lbpImage = Mat::zeros(gray_src.rows - 2, gray_src.cols - 2, CV_8UC1);//3*3窗口，边界有左右各1个像素不处理33     for (int row = 1; row < height - 1; row++) {34         for (int col = 1; col < width - 1; col++) {35             uchar c = gray_src.at
       
        (row, col);//获取中心像素值36             uchar code = 0;//特征码37             code |= (gray_src.at
        
         (row - 1, col - 1) > c) << 7;38             code |= (gray_src.at
         
          (row - 1, col) > c) << 6;39             code |= (gray_src.at
          
           (row - 1, col + 1) > c) << 5;40 code |= (gray_src.at
           
            (row, col + 1) > c) << 4;41 code |= (gray_src.at
            
             (row + 1, col + 1) > c) << 3;42 code |= (gray_src.at
             
              (row + 1, col) > c) << 2;43 code |= (gray_src.at
              
               (row + 1, col - 1) > c) << 1;44 code |= (gray_src.at
               
                (row, col - 1) > c) << 0;45 lbpImage.at
                
                 (row - 1, col - 1) = code;46 }47 }48 imshow(output_tt, lbpImage);49 50 // ELBP 演示51 52 namedWindow("ELBP Result", CV_WINDOW_AUTOSIZE);53 createTrackbar("ELBP Radius:", "ELBP Result", &current_radius, max_count, ELBP_Demo);54 ELBP_Demo(0, 0);55 56 waitKey(0);57 return 0;58 }59 60 void ELBP_Demo(int, void*) {61 int offset = current_radius * 2;//输出图无法计算到的边界宽度62 Mat elbpImage = Mat::zeros(gray_src.rows - offset, gray_src.cols - offset, CV_8UC1);63 int width = gray_src.cols;64 int height = gray_src.rows;65 66 int numNeighbors = 8;//应根据半径变化，为了简便这里固定为867 for (int n = 0; n < numNeighbors; n++) {68 float x = static_cast
                 
                  (current_radius) * cos(2.0 * CV_PI*n / static_cast
                  
                   (numNeighbors));//x=r*cos(2*PI*n/8)69 float y = static_cast
                   
                    (current_radius) * -sin(2.0 * CV_PI*n / static_cast
                    
                     (numNeighbors));70 71 int fx = static_cast
                     
                      (floor(x));//floor“向下取整”72 int fy = static_cast
                      
                       (floor(y));73 int cx = static_cast
                       
                        (ceil(x));//ceil“向上取整”74 int cy = static_cast
                        
                         (ceil(y));75 76 //双线性插值77 float ty = y - fy;78 float tx = x - fx;79 80 //双线性插值，周边4个点的权重81 float w1 = (1 - tx)*(1 - ty);82 float w2 = tx * (1 - ty);83 float w3 = (1 - tx)* ty;84 float w4 = tx * ty;85 86 for (int row = current_radius; row < (height - current_radius); row++) {87 for (int col = current_radius; col < (width - current_radius); col++) {88 float t = w1 * gray_src.at
                         
                          (row + fy, col + fx) + w2 * gray_src.at
                          
                           (row + fy, col + cx) +89 w3 * gray_src.at
                           
                            (row + cy, col + fx) + w4 * gray_src.at
                            
                             (row + cy, col + cx);90 elbpImage.at
                             
                              (row - current_radius, col - current_radius) +=91 ((t > gray_src.at
                              
                               (row, col)) && (abs(t - gray_src.at
                               
                                (row, col)) > std::numeric_limits
                                
                                 ::epsilon())) << n;92 }93 }94 }95 imshow("ELBP Result", elbpImage);96 return;97 }