LBP C++实现代码

// FaceRecognize.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"

#include "opencv2/opencv.hpp"
#include "CossValidtion.h"
using namespace std;
using namespace cv;
using namespace kl1p;
//#define LEN 30

int PixelToBLP(Mat image,arma::Col<klab::DoubleReal> &lbp){

int i,j,k,g=0;
double lbp_value,a[58][2],sum;
a[0][0]=0;a[1][0]=1;a[2][0]=2;a[3][0]=3;a[4][0]=4;a[5][0]=6;a[6][0]=7;a[7][0]=8;a[8][0]=12;a[9][0]=14;
a[10][0]=15;a[11][0]=16;a[12][0]=24;a[13][0]=28;a[14][0]=30;a[15][0]=31;a[16][0]=32;a[17][0]=48;a[18][0]=56;a[19][0]=60;
a[20][0]=62;a[21][0]=63;a[22][0]=64;a[23][0]=96;a[24][0]=112;a[25][0]=120;a[26][0]=124;a[27][0]=126;a[28][0]=127;a[29][0]=128;
a[30][0]=129;a[31][0]=131;a[32][0]=135;a[33][0]=143;a[34][0]=159;a[35][0]=191;a[36][0]=192;a[37][0]=193;a[38][0]=195;a[39][0]=199;
a[40][0]=207;a[41][0]=223;a[42][0]=224;a[43][0]=225;a[44][0]=227;a[45][0]=231;a[46][0]=239;a[47][0]=240;a[48][0]=241;a[49][0]=243;
a[50][0]=247;a[51][0]=248;a[52][0]=249;a[53][0]=251;a[54][0]=252;a[55][0]=253;a[56][0]=254;a[57][0]=255;

sum=0;
for(i=0;i<58;i++)//cell1
a[i][1]=0;
for(i=1;i<10;i++)
for(j=1;j<10;j++)
{
lbp_value=0;
if(image.at<uchar>(i-1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+128;
if(image.at<uchar>(i-1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+64;
if(image.at<uchar>(i-1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+32;
if(image.at<uchar>(i,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+16;
if(image.at<uchar>(i+1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+8;
if(image.at<uchar>(i+1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+4;
if(image.at<uchar>(i+1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+2;
if(image.at<uchar>(i,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+1;
for(k=0;k<58;k++)
if(a[k][0]==lbp_value)//若是uniform
{
a[k][1]++;
sum++;
}

}
for(i=0;i<58;i++)
lbp(g++)=a[i][1]/sum;//一范式归一化

sum=0;
for(i=0;i<58;i++)//cell2
a[i][1]=0;
for(i=1;i<10;i++)
for(j=10;j<20;j++)
{
lbp_value=0;
if(image.at<uchar>(i-1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+128;
if(image.at<uchar>(i-1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+64;
if(image.at<uchar>(i-1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+32;
if(image.at<uchar>(i,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+16;
if(image.at<uchar>(i+1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+8;
if(image.at<uchar>(i+1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+4;
if(image.at<uchar>(i+1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+2;
if(image.at<uchar>(i,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+1;
for(k=0;k<58;k++)
if(a[k][0]==lbp_value)//若是uniform
{
a[k][1]++;
sum++;
}
}
for(i=0;i<58;i++)
lbp(g++)=a[i][1]/sum;//一范式归一化

sum=0;
for(i=0;i<58;i++)//cell3
a[i][1]=0;
for(i=1;i<10;i++)
for(j=20;j<29;j++)
{
lbp_value=0;
if(image.at<uchar>(i-1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+128;
if(image.at<uchar>(i-1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+64;
if(image.at<uchar>(i-1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+32;
if(image.at<uchar>(i,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+16;
if(image.at<uchar>(i+1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+8;
if(image.at<uchar>(i+1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+4;
if(image.at<uchar>(i+1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+2;
if(image.at<uchar>(i,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+1;
for(k=0;k<58;k++)
if(a[k][0]==lbp_value)//若是uniform
{
a[k][1]++;
sum++;
}
}
for(i=0;i<58;i++)
lbp(g++)=a[i][1]/sum;//一范式归一化

sum=0;
for(i=0;i<58;i++)//cell4
a[i][1]=0;
for(i=10;i<20;i++)
for(j=1;j<10;j++)
{
lbp_value=0;
if(image.at<uchar>(i-1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+128;
if(image.at<uchar>(i-1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+64;
if(image.at<uchar>(i-1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+32;
if(image.at<uchar>(i,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+16;
if(image.at<uchar>(i+1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+8;
if(image.at<uchar>(i+1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+4;
if(image.at<uchar>(i+1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+2;
if(image.at<uchar>(i,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+1;
for(k=0;k<58;k++)
if(a[k][0]==lbp_value)//若是uniform
{
a[k][1]++;
sum++;
}
}
for(i=0;i<58;i++)
lbp(g++)=a[i][1]/sum;//一范式归一化

sum=0;
for(i=0;i<58;i++)//cell5
a[i][1]=0;
for(i=10;i<20;i++)
for(j=10;j<20;j++)
{
lbp_value=0;
if(image.at<uchar>(i-1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+128;
if(image.at<uchar>(i-1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+64;
if(image.at<uchar>(i-1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+32;
if(image.at<uchar>(i,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+16;
if(image.at<uchar>(i+1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+8;
if(image.at<uchar>(i+1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+4;
if(image.at<uchar>(i+1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+2;
if(image.at<uchar>(i,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+1;
for(k=0;k<58;k++)
if(a[k][0]==lbp_value)//若是uniform
{
a[k][1]++;
sum++;
}
}
for(i=0;i<58;i++)
lbp(g++)=a[i][1]/sum;//一范式归一化

sum=0;
for(i=0;i<58;i++)//cell6
a[i][1]=0;
for(i=10;i<20;i++)
for(j=20;j<29;j++)
{
lbp_value=0;
if(image.at<uchar>(i-1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+128;
if(image.at<uchar>(i-1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+64;
if(image.at<uchar>(i-1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+32;
if(image.at<uchar>(i,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+16;
if(image.at<uchar>(i+1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+8;
if(image.at<uchar>(i+1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+4;
if(image.at<uchar>(i+1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+2;
if(image.at<uchar>(i,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+1;
for(k=0;k<58;k++)
if(a[k][0]==lbp_value)//若是uniform
{
a[k][1]++;
sum++;
}
}
for(i=0;i<58;i++)
lbp(g++)=a[i][1]/sum;//一范式归一化

sum=0;
for(i=0;i<58;i++)//cell7
a[i][1]=0;
for(i=20;i<29;i++)
for(j=1;j<10;j++)
{
lbp_value=0;
if(image.at<uchar>(i-1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+128;
if(image.at<uchar>(i-1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+64;
if(image.at<uchar>(i-1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+32;
if(image.at<uchar>(i,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+16;
if(image.at<uchar>(i+1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+8;
if(image.at<uchar>(i+1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+4;
if(image.at<uchar>(i+1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+2;
if(image.at<uchar>(i,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+1;
for(k=0;k<58;k++)
if(a[k][0]==lbp_value)//若是uniform
{
a[k][1]++;
sum++;
}
}
for(i=0;i<58;i++)
lbp(g++)=a[i][1]/sum;//一范式归一化

sum=0;
for(i=0;i<58;i++)//cell8
a[i][1]=0;
for(i=20;i<29;i++)
for(j=10;j<20;j++)
{
lbp_value=0;
if(image.at<uchar>(i-1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+128;
if(image.at<uchar>(i-1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+64;
if(image.at<uchar>(i-1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+32;
if(image.at<uchar>(i,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+16;
if(image.at<uchar>(i+1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+8;
if(image.at<uchar>(i+1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+4;
if(image.at<uchar>(i+1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+2;
if(image.at<uchar>(i,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+1;
for(k=0;k<58;k++)
if(a[k][0]==lbp_value)//若是uniform
{
a[k][1]++;
sum++;
}
}
for(i=0;i<58;i++)
lbp(g++)=a[i][1]/sum;//一范式归一化

sum=0;
for(i=0;i<58;i++)//cell9
a[i][1]=0;
for(i=20;i<29;i++)
for(j=20;j<29;j++)
{
lbp_value=0;
if(image.at<uchar>(i-1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+128;
if(image.at<uchar>(i-1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+64;
if(image.at<uchar>(i-1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+32;
if(image.at<uchar>(i,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+16;
if(image.at<uchar>(i+1,j+1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+8;
if(image.at<uchar>(i+1,j) >= image.at<uchar>(i,j))
lbp_value=lbp_value+4;
if(image.at<uchar>(i+1,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+2;
if(image.at<uchar>(i,j-1) >= image.at<uchar>(i,j))
lbp_value=lbp_value+1;
for(k=0;k<58;k++)
if(a[k][0]==lbp_value)//若是uniform
{
a[k][1]++;
sum++;
}
}
for(i=0;i<58;i++)
lbp(g++)=a[i][1]/sum;//一范式归一化

//cout<<g<<endl;
return 0;
}

int _tmain(int argc, _TCHAR* argv[])
{
arma::Col<klab::DoubleReal> *lbp = new arma::Col<klab::DoubleReal>(522);//58*9=522
Mat image = imread("face_database\\001\\meitu_00001.jpg"); //图片为30*30=900的图片，我要把它划成10*10的cell，划9个cell
PixelToBLP(image,*lbp);

cout<<"over"<<endl;
getchar();
return 0;
}