LBP算法（局部二值模式）

#include<cv.h>
#include<highgui.h>
#include<iostream>
#include<fstream>
using namespace std;
int main()
{
IplImage *iplImage1 = NULL;
IplImage *iplImage2 = NULL;
IplImage *iplImage3 = NULL;
CvSize cvSize;
cvSize.width = 9;
cvSize.height = 9;
iplImage1 = cvLoadImage("test15.jpg");
/*iplImage3 = cvCreateImage(cvSize,IPL_DEPTH_8U, 1);
if(!iplImage3) cout<<"cannot create"<<endl;*/
iplImage2 = cvCreateImage(cvSize,iplImage1->depth,iplImage1->nChannels);
/*cvCvtColor(iplImage1,iplImage3,CV_BGR2GRAY);
cvNamedWindow("test",1);
cvShowImage("test",iplImage3);
cvWaitKey(0);*/
cvResize(iplImage1,iplImage2,CV_INTER_LINEAR);
int n = 5;
for(int i = 4;i < n;i++)
{

double **mat;
mat = new double *[2*i+1];
for (int j=0; j<2*i+1; j++)
{