{
"id": 333,
"title": "Benchmark",
"url": "/competitions/WIC2011/discussion/333",
"author_name": "",
"post_date": "2011-03-13T17:45:09Z",
"votes": null,
"comment_count": 2,
"views": 2936,
"content": "Dear participants,
Below is the c++ code used to generate the benchmark.
Comments are welcome.
#include \"stdafx.h\"
#include <fstream>
#include <conio.h>
#include <list>
#include <iostream>
#include <sstream>
#include <vector>
#include <string>
#include <math.h>
using namespace std;
//class handling the provided features
class document_features
{
public:
document_features()
{
ThicknessLengthsCircleHist30=new double[30];
XProjectionHist10=new double[10];
YProjectionHist10=new double[10];
Distribution10x10_100=new double[100];
tortuosityHist10=new double[10];
tortuosityDirectionHist10=new double[10];
tortuosityDerivateHist10=new double[10];
tortuosityDerivateDirectionHist10=new double[10];
nbbranches_1=new double[1];
LengthsOfBranchesHist_10=new double[10];
numberOfPixels1=new double[1];
XFilledProjectionHist10=new double[10];
YFilledProjectionHist10=new double[10];
Barycenter2=new double[2];
DirectionPerpendicular5Hist10=new double[10];
CurvaturePerpendicular5Hist100=new double[100];
luminanceHist256=new double[256];
CurvatureAli5Hist100=new double[100];
CurvaturesDerivateAli5Hist100=new double[100];
CurvatureAli10Hist100=new double[100];
CurvaturesDerivateAli10Hist100=new double[100];
CurvatureAli15Hist100=new double[100];
CurvaturesDerivateAli15Hist100=new double[100];
CurvatureAli20Hist100=new double[100];
CurvaturesDerivateAli20Hist100=new double[100];
chaincodeHist_4=new double[4];
chaincodeHist_8=new double[8];
chaincode8order2_64=new double[64];
chaincode4order2_16=new double[16];
chaincode4order3_64=new double[64];
chaincode8order3_512=new double[512];
chaincode4order4_256=new double[256];
chaincode8order4_4096=new double[4096];
NumberOfConnectedComponents_1=new double[1];
NumberOfHoles_1=new double[1];
Fourier_1=new double[1];
Fourier_5=new double[5];
Fourier_9=new double[9];
Fourier_15=new double[15];
}
double*ThicknessLengthsCircleHist30,*XProjectionHist10,*YProjectionHist10,*Distribution10x10_100,*tortuosityHist10,*tortuosityDirectionHist10,*tortuosityDerivateHist10,*tortuosityDerivateDirectionHist10,*nbbranches_1,*LengthsOfBranchesHist_10,*numberOfPixels1,*XFilledProjectionHist10,*YFilledProjectionHist10,*Barycenter2,*DirectionPerpendicular5Hist10,*CurvaturePerpendicular5Hist100,*luminanceHist256,*CurvatureAli5Hist100,*CurvaturesDerivateAli5Hist100,*CurvatureAli10Hist100,*CurvaturesDerivateAli10Hist100,*CurvatureAli15Hist100,*CurvaturesDerivateAli15Hist100,*CurvatureAli20Hist100,*CurvaturesDerivateAli20Hist100,*chaincodeHist_4,*chaincodeHist_8,*chaincode8order2_64,*chaincode4order2_16,*chaincode4order3_64,*chaincode8order3_512,*chaincode4order4_256,*chaincode8order4_4096,*NumberOfConnectedComponents_1,*NumberOfHoles_1,*Fourier_1,*Fourier_5,*Fourier_9,*Fourier_15;
double SpatialMoment00,SpatialMoment10,SpatialMoment20,SpatialMoment30,SpatialMoment01,SpatialMoment11,SpatialMoment21,SpatialMoment02,SpatialMoment12,SpatialMoment03,CentralMoment00,CentralMoment10,CentralMoment20,CentralMoment30,CentralMoment01,CentralMoment11,CentralMoment21,CentralMoment02,CentralMoment12,CentralMoment03,NormalizedCentralMoment00,NormalizedCentralMoment10,NormalizedCentralMoment20,NormalizedCentralMoment30,NormalizedCentralMoment01,NormalizedCentralMoment11,NormalizedCentralMoment21,NormalizedCentralMoment02,NormalizedCentralMoment12,NormalizedCentralMoment03,HuMoment1,HuMoment2,HuMoment3,HuMoment4,HuMoment5,HuMoment6,HuMoment7;
int writer_id;
};
void main()
{
//training phase
std::ifstream data(\"train.csv\");
int linenum=0;
string line;
document_features**my_document_features=new document_features*[108];
string item;//=new char[4];
//reading the training documents
while (getline (data, line))
{
cout << \"\\nLine #\" << linenum << \":\" << endl;
istringstream linestream(line);
if(linenum>0)
{
int row_index=linenum-1;
my_document_features[row_index]=new document_features();
getline (linestream, item, ',');
my_document_features[row_index]->writer_id=atoi(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NumberOfConnectedComponents_1[0]=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NumberOfHoles_1[0]=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment00=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment10=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment20=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment30=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment01=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment11=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment21=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment02=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment12=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment03=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment00=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment10=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment20=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment30=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment01=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment11=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment21=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment02=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment12=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment03=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment00=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment10=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment20=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment30=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment01=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment11=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment21=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment02=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment12=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment03=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment1=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment2=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment3=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment4=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment5=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment6=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment7=atof(item.c_str());
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->XProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->YProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->XFilledProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->YFilledProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->Distribution10x10_100[x]=atof(item.c_str());
}
getline (linestream, item, ',');
my_document_features[row_index]->Barycenter2[0]=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->Barycenter2[1]=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->numberOfPixels1[0]=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->Fourier_1[0]=atof(item.c_str());
for(int x=0;x<5;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->Fourier_5[x]=atof(item.c_str());
}
for(int x=0;x<9;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->Fourier_9[x]=atof(item.c_str());
}
for(int x=0;x<15;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->Fourier_15[x]=atof(item.c_str());
}
getline (linestream, item, ',');
my_document_features[row_index]->nbbranches_1[0]=atof(item.c_str());
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->LengthsOfBranchesHist_10[x]=atof(item.c_str());
}
for(int x=0;x<30;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->ThicknessLengthsCircleHist30[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->tortuosityHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->tortuosityDirectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->tortuosityDerivateHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->tortuosityDerivateDirectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->DirectionPerpendicular5Hist10[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvaturePerpendicular5Hist100[x]=atof(item.c_str());
}
for(int x=0;x<256;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->luminanceHist256[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvatureAli5Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvaturesDerivateAli5Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvatureAli10Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvaturesDerivateAli10Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvatureAli15Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvaturesDerivateAli15Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvatureAli20Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvaturesDerivateAli20Hist100[x]=atof(item.c_str());
}
for(int x=0;x<4;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincodeHist_4[x]=atof(item.c_str());
}
for(int x=0;x<8;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincodeHist_8[x]=atof(item.c_str());
}
for(int x=0;x<64;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincode8order2_64[x]=atof(item.c_str());
}
for(int x=0;x<16;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincode4order2_16[x]=atof(item.c_str());
}
for(int x=0;x<64;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincode4order3_64[x]=atof(item.c_str());
}
for(int x=0;x<512;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincode8order3_512[x]=atof(item.c_str());
}
for(int x=0;x<256;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincode4order4_256[x]=atof(item.c_str());
}
for(int x=0;x<4096;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincode8order4_4096[x]=atof(item.c_str());
}
}
linenum++;
}
//difference between features of all documents in the training set
std::ofstream differences(\"train_differences.csv\");
differences<<\"SameWriter,NumberOfConnectedComponents_1,NumberOfHoles_1,SpatialMoment00,SpatialMoment10,SpatialMoment20,SpatialMoment30,SpatialMoment01,SpatialMoment11,SpatialMoment21,SpatialMoment02,SpatialMoment12,SpatialMoment03,CentralMoment00,CentralMoment10,CentralMoment20,CentralMoment30,CentralMoment01,CentralMoment11,CentralMoment21,CentralMoment02,CentralMoment12,CentralMoment03,NormalizedCentralMoment00,NormalizedCentralMoment10,NormalizedCentralMoment20,NormalizedCentralMoment30,NormalizedCentralMoment01,NormalizedCentralMoment11,NormalizedCentralMoment21,NormalizedCentralMoment02,NormalizedCentralMoment12,NormalizedCentralMoment03,HuMoment1,HuMoment2,HuMoment3,HuMoment4,HuMoment5,HuMoment6,HuMoment7,XProjectionHist10,YProjectionHist10,XFilledProjectionHist10,YFilledProjectionHist10,Distribution10x10_100,Barycenter2,numberOfPixels1,Fourier_1,Fourier_5,Fourier_9,Fourier_15,nbbranches_1,LengthsOfBranchesHist_10,ThicknessLengthsCircleHist30,tortuosityHist10,tortuosityDirectionHist10,tortuosityDerivateHist10,tortuosityDerivateDirectionHist10,DirectionPerpendicular5Hist10,CurvaturePerpendicular5Hist100,luminanceHist256,CurvatureAli5Hist100,CurvaturesDerivateAli5Hist100,CurvatureAli10Hist100,CurvaturesDerivateAli10Hist100,CurvatureAli15Hist100,CurvaturesDerivateAli15Hist100,CurvatureAli20Hist100,CurvaturesDerivateAli20Hist100,chaincodeHist_4,chaincodeHist_8,chaincode8order2_64,chaincode4order2_16,chaincode4order3_64,chaincode8order3_512,chaincode4order4_256,chaincode8order4_4096\\n\";
//comparing each pair of documents
for(int document_index1=0;document_index1<107;document_index1++)
for(int document_index2=document_index1+1;document_index2<108;document_index2++)
{
double NumberOfConnectedComponents_1=0,NumberOfHoles_1=0,SpatialMoment00=0,SpatialMoment10=0,SpatialMoment20=0,SpatialMoment30=0,SpatialMoment01=0,SpatialMoment11=0,SpatialMoment21=0,SpatialMoment02=0,SpatialMoment12=0,SpatialMoment03=0,CentralMoment00=0,CentralMoment10=0,CentralMoment20=0,CentralMoment30=0,CentralMoment01=0,CentralMoment11=0,CentralMoment21=0,CentralMoment02=0,CentralMoment12=0,CentralMoment03=0,NormalizedCentralMoment00=0,NormalizedCentralMoment10=0,NormalizedCentralMoment20=0,NormalizedCentralMoment30=0,NormalizedCentralMoment01=0,NormalizedCentralMoment11=0,NormalizedCentralMoment21=0,NormalizedCentralMoment02=0,NormalizedCentralMoment12=0,NormalizedCentralMoment03=0,HuMoment1=0,HuMoment2=0,HuMoment3=0,HuMoment4=0,HuMoment5=0,HuMoment6=0,HuMoment7=0,ThicknessLengthsCircleHist30=0,XProjectionHist10=0,YProjectionHist10=0,Distribution10x10_100=0,tortuosityHist10=0,tortuosityDirectionHist10=0,tortuosityDerivateHist10=0,tortuosityDerivateDirectionHist10=0,nbbranches_1=0,LengthsOfBranchesHist_10=0,numberOfPixels1=0,XFilledProjectionHist10=0,YFilledProjectionHist10=0,Barycenter2=0,DirectionPerpendicular5Hist10=0,CurvaturePerpendicular5Hist100=0,luminanceHist256=0,CurvatureAli5Hist100=0,CurvaturesDerivateAli5Hist100=0,CurvatureAli10Hist100=0,CurvaturesDerivateAli10Hist100=0,CurvatureAli15Hist100=0,CurvaturesDerivateAli15Hist100=0,CurvatureAli20Hist100=0,CurvaturesDerivateAli20Hist100=0,chaincodeHist_4=0,chaincodeHist_8=0,chaincode8order2_64=0,chaincode4order2_16=0,chaincode4order3_64=0,chaincode8order3_512=0,chaincode4order4_256=0,chaincode8order4_4096=0,Fourier_1=0,Fourier_5=0,Fourier_9=0,Fourier_15=0;
NumberOfConnectedComponents_1+=fabs(my_document_features[document_index1]->NumberOfConnectedComponents_1[0]-my_document_features[document_index2]->NumberOfConnectedComponents_1[0]);
NumberOfHoles_1+=fabs(my_document_features[document_index1]->NumberOfHoles_1[0]-my_document_features[document_index2]->NumberOfHoles_1[0]);
SpatialMoment00+=fabs(my_document_features[document_index1]->SpatialMoment00-my_document_features[document_index2]->SpatialMoment00);
SpatialMoment10+=fabs(my_document_features[document_index1]->SpatialMoment10-my_document_features[document_index2]->SpatialMoment10);
SpatialMoment20+=fabs(my_document_features[document_index1]->SpatialMoment20-my_document_features[document_index2]->SpatialMoment20);
SpatialMoment30+=fabs(my_document_features[document_index1]->SpatialMoment30-my_document_features[document_index2]->SpatialMoment30);
SpatialMoment01+=fabs(my_document_features[document_index1]->SpatialMoment01-my_document_features[document_index2]->SpatialMoment01);
SpatialMoment11+=fabs(my_document_features[document_index1]->SpatialMoment11-my_document_features[document_index2]->SpatialMoment11);
SpatialMoment21+=fabs(my_document_features[document_index1]->SpatialMoment21-my_document_features[document_index2]->SpatialMoment21);
SpatialMoment02+=fabs(my_document_features[document_index1]->SpatialMoment02-my_document_features[document_index2]->SpatialMoment02);
SpatialMoment12+=fabs(my_document_features[document_index1]->SpatialMoment12-my_document_features[document_index2]->SpatialMoment12);
SpatialMoment03+=fabs(my_document_features[document_index1]->SpatialMoment03-my_document_features[document_index2]->SpatialMoment03);
CentralMoment00+=fabs(my_document_features[document_index1]->CentralMoment00-my_document_features[document_index2]->CentralMoment00);
CentralMoment10+=fabs(my_document_features[document_index1]->CentralMoment10-my_document_features[document_index2]->CentralMoment10);
CentralMoment20+=fabs(my_document_features[document_index1]->CentralMoment20-my_document_features[document_index2]->CentralMoment20);
CentralMoment30+=fabs(my_document_features[document_index1]->CentralMoment30-my_document_features[document_index2]->CentralMoment30);
CentralMoment01+=fabs(my_document_features[document_index1]->CentralMoment01-my_document_features[document_index2]->CentralMoment01);
CentralMoment11+=fabs(my_document_features[document_index1]->CentralMoment11-my_document_features[document_index2]->CentralMoment11);
CentralMoment21+=fabs(my_document_features[document_index1]->CentralMoment21-my_document_features[document_index2]->CentralMoment21);
CentralMoment02+=fabs(my_document_features[document_index1]->CentralMoment02-my_document_features[document_index2]->CentralMoment02);
CentralMoment12+=fabs(my_document_features[document_index1]->CentralMoment12-my_document_features[document_index2]->CentralMoment12);
CentralMoment03+=fabs(my_document_features[document_index1]->CentralMoment03-my_document_features[document_index2]->CentralMoment03);
NormalizedCentralMoment00+=fabs(my_document_features[document_index1]->NormalizedCentralMoment00-my_document_features[document_index2]->NormalizedCentralMoment00);
NormalizedCentralMoment10+=fabs(my_document_features[document_index1]->NormalizedCentralMoment10-my_document_features[document_index2]->NormalizedCentralMoment10);
NormalizedCentralMoment20+=fabs(my_document_features[document_index1]->NormalizedCentralMoment20-my_document_features[document_index2]->NormalizedCentralMoment20);
NormalizedCentralMoment30+=fabs(my_document_features[document_index1]->NormalizedCentralMoment30-my_document_features[document_index2]->NormalizedCentralMoment30);
NormalizedCentralMoment01+=fabs(my_document_features[document_index1]->NormalizedCentralMoment01-my_document_features[document_index2]->NormalizedCentralMoment01);
NormalizedCentralMoment11+=fabs(my_document_features[document_index1]->NormalizedCentralMoment11-my_document_features[document_index2]->NormalizedCentralMoment11);
NormalizedCentralMoment21+=fabs(my_document_features[document_index1]->NormalizedCentralMoment21-my_document_features[document_index2]->NormalizedCentralMoment21);
NormalizedCentralMoment02+=fabs(my_document_features[document_index1]->NormalizedCentralMoment02-my_document_features[document_index2]->NormalizedCentralMoment02);
NormalizedCentralMoment12+=fabs(my_document_features[document_index1]->NormalizedCentralMoment12-my_document_features[document_index2]->NormalizedCentralMoment12);
NormalizedCentralMoment03+=fabs(my_document_features[document_index1]->NormalizedCentralMoment03-my_document_features[document_index2]->NormalizedCentralMoment03);
HuMoment1+=fabs(my_document_features[document_index1]->HuMoment1-my_document_features[document_index2]->HuMoment1);
HuMoment2+=fabs(my_document_features[document_index1]->HuMoment2-my_document_features[document_index2]->HuMoment2);
HuMoment3+=fabs(my_document_features[document_index1]->HuMoment3-my_document_features[document_index2]->HuMoment3);
HuMoment4+=fabs(my_document_features[document_index1]->HuMoment4-my_document_features[document_index2]->HuMoment4);
HuMoment5+=fabs(my_document_features[document_index1]->HuMoment5-my_document_features[document_index2]->HuMoment5);
HuMoment6+=fabs(my_document_features[document_index1]->HuMoment6-my_document_features[document_index2]->HuMoment6);
HuMoment7+=fabs(my_document_features[document_index1]->HuMoment7-my_document_features[document_index2]->HuMoment7);
for(int x=0;x<10;x++)
XProjectionHist10+=fabs(my_document_features[document_index1]->XProjectionHist10[x]-my_document_features[document_index2]->XProjectionHist10[x]);
for(int x=0;x<10;x++)
YProjectionHist10+=fabs(my_document_features[document_index1]->YProjectionHist10[x]-my_document_features[document_index2]->YProjectionHist10[x]);
for(int x=0;x<10;x++)
XFilledProjectionHist10+=fabs(my_document_features[document_index1]->XFilledProjectionHist10[x]-my_document_features[document_index2]->XFilledProjectionHist10[x]);
for(int x=0;x<10;x++)
YFilledProjectionHist10+=fabs(my_document_features[document_index1]->YFilledProjectionHist10[x]-my_document_features[document_index2]->YFilledProjectionHist10[x]);
for(int x=0;x<10;x++)
Distribution10x10_100+=fabs(my_document_features[document_index1]->Distribution10x10_100[x]-my_document_features[document_index2]->Distribution10x10_100[x]);
for(int x=0;x<2;x++)
Barycenter2+=fabs(my_document_features[document_index1]->Barycenter2[x]-my_document_features[document_index2]->Barycenter2[x]);
numberOfPixels1+=fabs(my_document_features[document_index1]->numberOfPixels1[0]-my_document_features[document_index2]->numberOfPixels1[0]);
for(int x=0;x<1;x++)
Fourier_1+=fabs(my_document_features[document_index1]->Fourier_1[x]-my_document_features[document_index2]->Fourier_1[x]);
for(int x=0;x<5;x++)
Fourier_5+=fabs(my_document_features[document_index1]->Fourier_5[x]-my_document_features[document_index2]->Fourier_5[x]);
for(int x=0;x<9;x++)
Fourier_9+=fabs(my_document_features[document_index1]->Fourier_9[x]-my_document_features[document_index2]->Fourier_9[x]);
for(int x=0;x<15;x++)
Fourier_15+=fabs(my_document_features[document_index1]->Fourier_15[x]-my_document_features[document_index2]->Fourier_15[x]);
nbbranches_1+=fabs(my_document_features[document_index1]->nbbranches_1[0]-my_document_features[document_index2]->nbbranches_1[0]);
for(int x=0;x<10;x++)
LengthsOfBranchesHist_10+=fabs(my_document_features[document_index1]->LengthsOfBranchesHist_10[x]-my_document_features[document_index2]->LengthsOfBranchesHist_10[x]);
for(int x=0;x<30;x++)
ThicknessLengthsCircleHist30+=fabs(my_document_features[document_index1]->ThicknessLengthsCircleHist30[x]-my_document_features[document_index2]->ThicknessLengthsCircleHist30[x]);
for(int x=0;x<10;x++)
tortuosityHist10+=fabs(my_document_features[document_index1]->tortuosityHist10[x]-my_document_features[document_index2]->tortuosityHist10[x]);
for(int x=0;x<10;x++)
tortuosityDirectionHist10+=fabs(my_document_features[document_index1]->tortuosityDirectionHist10[x]-my_document_features[document_index2]->tortuosityDirectionHist10[x]);
for(int x=0;x<10;x++)
tortuosityDerivateHist10+=fabs(my_document_features[document_index1]->tortuosityDerivateHist10[x]-my_document_features[document_index2]->tortuosityDerivateHist10[x]);
for(int x=0;x<10;x++)
tortuosityDerivateDirectionHist10+=fabs(my_document_features[document_index1]->tortuosityDerivateDirectionHist10[x]-my_document_features[document_index2]->tortuosityDerivateDirectionHist10[x]);
for(int x=0;x<10;x++)
DirectionPerpendicular5Hist10+=fabs(my_document_features[document_index1]->DirectionPerpendicular5Hist10[x]-my_document_features[document_index2]->DirectionPerpendicular5Hist10[x]);
for(int x=0;x<100;x++)
CurvaturePerpendicular5Hist100+=fabs(my_document_features[document_index1]->CurvaturePerpendicular5Hist100[x]-my_document_features[document_index2]->CurvaturePerpendicular5Hist100[x]);
for(int x=0;x<256;x++)
luminanceHist256+=fabs(my_document_features[document_index1]->luminanceHist256[x]-my_document_features[document_index2]->luminanceHist256[x]);
for(int x=0;x<100;x++)
CurvatureAli5Hist100+=fabs(my_document_features[document_index1]->CurvatureAli5Hist100[x]-my_document_features[document_index2]->CurvatureAli5Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli5Hist100+=fabs(my_document_features[document_index1]->CurvaturesDerivateAli5Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli5Hist100[x]);
for(int x=0;x<100;x++)
CurvatureAli10Hist100+=fabs(my_document_features[document_index1]->CurvatureAli10Hist100[x]-my_document_features[document_index2]->CurvatureAli10Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli10Hist100+=fabs(my_document_features[document_index1]->CurvaturesDerivateAli10Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli10Hist100[x]);
for(int x=0;x<100;x++)
CurvatureAli15Hist100+=fabs(my_document_features[document_index1]->CurvatureAli15Hist100[x]-my_document_features[document_index2]->CurvatureAli15Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli15Hist100+=fabs(my_document_features[document_index1]->CurvaturesDerivateAli15Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli15Hist100[x]);
for(int x=0;x<100;x++)
CurvatureAli20Hist100+=fabs(my_document_features[document_index1]->CurvatureAli20Hist100[x]-my_document_features[document_index2]->CurvatureAli20Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli20Hist100+=fabs(my_document_features[document_index1]->CurvaturesDerivateAli20Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli20Hist100[x]);
for(int x=0;x<4;x++)
chaincodeHist_4+=fabs(my_document_features[document_index1]->chaincodeHist_4[x]-my_document_features[document_index2]->chaincodeHist_4[x]);
for(int x=0;x<8;x++)
chaincodeHist_8+=fabs(my_document_features[document_index1]->chaincodeHist_8[x]-my_document_features[document_index2]->chaincodeHist_8[x]);
for(int x=0;x<64;x++)
chaincode8order2_64+=fabs(my_document_features[document_index1]->chaincode8order2_64[x]-my_document_features[document_index2]->chaincode8order2_64[x]);
for(int x=0;x<16;x++)
chaincode4order2_16+=fabs(my_document_features[document_index1]->chaincode4order2_16[x]-my_document_features[document_index2]->chaincode4order2_16[x]);
for(int x=0;x<64;x++)
chaincode4order3_64+=fabs(my_document_features[document_index1]->chaincode4order3_64[x]-my_document_features[document_index2]->chaincode4order3_64[x]);
for(int x=0;x<512;x++)
chaincode8order3_512+=fabs(my_document_features[document_index1]->chaincode8order3_512[x]-my_document_features[document_index2]->chaincode8order3_512[x]);
for(int x=0;x<256;x++)
chaincode4order4_256+=fabs(my_document_features[document_index1]->chaincode4order4_256[x]-my_document_features[document_index2]->chaincode4order4_256[x]);
for(int x=0;x<4096;x++)
chaincode8order4_4096+=fabs(my_document_features[document_index1]->chaincode8order4_4096[x]-my_document_features[document_index2]->chaincode8order4_4096[x]);
if(my_document_features[document_index1]->writer_id==my_document_features[document_index2]->writer_id)differences<<\"1,\";
else differences<<\"0,\";
differences<<NumberOfConnectedComponents_1<<\",\"<<NumberOfHoles_1<<\",\"<<SpatialMoment00<<\",\"<<SpatialMoment10<<\",\"<<SpatialMoment20<<\",\"<<SpatialMoment30<<\",\"<<SpatialMoment01<<\",\"<<SpatialMoment11<<\",\"<<SpatialMoment21<<\",\"<<SpatialMoment02<<\",\"<<SpatialMoment12<<\",\"<<SpatialMoment03<<\",\"<<CentralMoment00<<\",\"<<CentralMoment10<<\",\"<<CentralMoment20<<\",\"<<CentralMoment30<<\",\"<<CentralMoment01<<\",\"<<CentralMoment11<<\",\"<<CentralMoment21<<\",\"<<CentralMoment02<<\",\"<<CentralMoment12<<\",\"<<CentralMoment03<<\",\"<<NormalizedCentralMoment00<<\",\"<<NormalizedCentralMoment10<<\",\"<<NormalizedCentralMoment20<<\",\"<<NormalizedCentralMoment30<<\",\"<<NormalizedCentralMoment01<<\",\"<<NormalizedCentralMoment11<<\",\"<<NormalizedCentralMoment21<<\",\"<<NormalizedCentralMoment02<<\",\"<<NormalizedCentralMoment12<<\",\"<<NormalizedCentralMoment03<<\",\"<<HuMoment1<<\",\"<<HuMoment2<<\",\"<<HuMoment3<<\",\"<<HuMoment4<<\",\"<<HuMoment5<<\",\"<<HuMoment6<<\",\"<<HuMoment7<<\",\"<<XProjectionHist10<<\",\"<<YProjectionHist10<<\",\"<<XFilledProjectionHist10<<\",\"<<YFilledProjectionHist10<<\",\"<<Distribution10x10_100<<\",\"<<Barycenter2<<\",\"<<numberOfPixels1<<\",\"<<Fourier_1<<\",\"<<Fourier_5<<\",\"<<Fourier_9<<\",\"<<Fourier_15<<\",\"<<nbbranches_1<<\",\"<<LengthsOfBranchesHist_10<<\",\"<<ThicknessLengthsCircleHist30<<\",\"<<tortuosityHist10<<\",\"<<tortuosityDirectionHist10<<\",\"<<tortuosityDerivateHist10<<\",\"<<tortuosityDerivateDirectionHist10<<\",\"<<DirectionPerpendicular5Hist10<<\",\"<<CurvaturePerpendicular5Hist100<<\",\"<<luminanceHist256<<\",\"<<CurvatureAli5Hist100<<\",\"<<CurvaturesDerivateAli5Hist100<<\",\"<<CurvatureAli10Hist100<<\",\"<<CurvaturesDerivateAli10Hist100<<\",\"<<CurvatureAli15Hist100<<\",\"<<CurvaturesDerivateAli15Hist100<<\",\"<<CurvatureAli20Hist100<<\",\"<<CurvaturesDerivateAli20Hist100<<\",\"<<chaincodeHist_4<<\",\"<<chaincodeHist_8<<\",\"<<chaincode8order2_64<<\",\"<<chaincode4order2_16<<\",\"<<chaincode4order3_64<<\",\"<<chaincode8order3_512<<\",\"<<chaincode4order4_256<<\",\"<<chaincode8order4_4096<<\"\\n\";
}
differences.close();
//You might use the obtained \"train_differences.csv\" file to build your model
//I used here the model trained using a logistic regression classifier after keeping the 5 most discriminant features :
//double linear=(-10.667528249636) + -0.0667090106640684 * NumberOfConnectedComponents_1 +9.34910163750235 * ThicknessLengthsCircleHist30 + 30.4471654988069 *CurvaturePerpendicular5Hist100 + 30.5600351874986 * luminanceHist256 +33.6853635484793 * chaincode8order3_512;
//double prob_same_writer=1. / (1. + exp(linear));
//reading the testing documents
std::ifstream test_data(\"test.csv\");
linenum=0;
document_features**my_questioned_document_features=new document_features*[53];
while (getline (test_data, line))
{
cout << \"\\nLine #\" << linenum << \":\" << endl;
istringstream linestream(line);
if(linenum>0)
{
int row_index=linenum-1;
my_questioned_document_features[row_index]=new document_features();
getline (linestream, item, ',');
//my_questioned_document_features[row_index]->writer_id=atoi(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NumberOfConnectedComponents_1[0]=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NumberOfHoles_1[0]=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment00=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment10=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment20=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment30=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment01=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment11=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment21=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment02=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment12=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment03=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment00=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment10=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment20=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment30=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment01=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment11=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment21=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment02=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment12=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment03=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment00=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment10=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment20=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment30=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment01=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment11=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment21=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment02=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment12=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment03=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment1=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment2=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment3=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment4=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment5=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment6=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment7=atof(item.c_str());
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->XProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->YProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->XFilledProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->YFilledProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Distribution10x10_100[x]=atof(item.c_str());
}
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Barycenter2[0]=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Barycenter2[1]=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->numberOfPixels1[0]=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Fourier_1[0]=atof(item.c_str());
for(int x=0;x<5;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Fourier_5[x]=atof(item.c_str());
}
for(int x=0;x<9;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Fourier_9[x]=atof(item.c_str());
}
for(int x=0;x<15;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Fourier_15[x]=atof(item.c_str());
}
getline (linestream, item, ',');
my_questioned_document_features[row_index]->nbbranches_1[0]=atof(item.c_str());
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->LengthsOfBranchesHist_10[x]=atof(item.c_str());
}
for(int x=0;x<30;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->ThicknessLengthsCircleHist30[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->tortuosityHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->tortuosityDirectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->tortuosityDerivateHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->tortuosityDerivateDirectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->DirectionPerpendicular5Hist10[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvaturePerpendicular5Hist100[x]=atof(item.c_str());
}
for(int x=0;x<256;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->luminanceHist256[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvatureAli5Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvaturesDerivateAli5Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvatureAli10Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvaturesDerivateAli10Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvatureAli15Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvaturesDerivateAli15Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvatureAli20Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvaturesDerivateAli20Hist100[x]=atof(item.c_str());
}
for(int x=0;x<4;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincodeHist_4[x]=atof(item.c_str());
}
for(int x=0;x<8;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincodeHist_8[x]=atof(item.c_str());
}
for(int x=0;x<64;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincode8order2_64[x]=atof(item.c_str());
}
for(int x=0;x<16;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincode4order2_16[x]=atof(item.c_str());
}
for(int x=0;x<64;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincode4order3_64[x]=atof(item.c_str());
}
for(int x=0;x<512;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincode8order3_512[x]=atof(item.c_str());
}
for(int x=0;x<256;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincode4order4_256[x]=atof(item.c_str());
}
for(int x=0;x<4096;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincode8order4_4096[x]=atof(item.c_str());
}
}
linenum++;
}
//comparing all the documents of the test set with all the documents of the training set
for(int document_index1=0;document_index1<53;document_index1++)
{
double maxprob=0;
for(int document_index2=0;document_index2<108;document_index2++)
{
double NumberOfConnectedComponents_1=0,NumberOfHoles_1=0,SpatialMoment00=0,SpatialMoment10=0,SpatialMoment20=0,SpatialMoment30=0,SpatialMoment01=0,SpatialMoment11=0,SpatialMoment21=0,SpatialMoment02=0,SpatialMoment12=0,SpatialMoment03=0,CentralMoment00=0,CentralMoment10=0,CentralMoment20=0,CentralMoment30=0,CentralMoment01=0,CentralMoment11=0,CentralMoment21=0,CentralMoment02=0,CentralMoment12=0,CentralMoment03=0,NormalizedCentralMoment00=0,NormalizedCentralMoment10=0,NormalizedCentralMoment20=0,NormalizedCentralMoment30=0,NormalizedCentralMoment01=0,NormalizedCentralMoment11=0,NormalizedCentralMoment21=0,NormalizedCentralMoment02=0,NormalizedCentralMoment12=0,NormalizedCentralMoment03=0,HuMoment1=0,HuMoment2=0,HuMoment3=0,HuMoment4=0,HuMoment5=0,HuMoment6=0,HuMoment7=0,ThicknessLengthsCircleHist30=0,XProjectionHist10=0,YProjectionHist10=0,Distribution10x10_100=0,tortuosityHist10=0,tortuosityDirectionHist10=0,tortuosityDerivateHist10=0,tortuosityDerivateDirectionHist10=0,nbbranches_1=0,LengthsOfBranchesHist_10=0,numberOfPixels1=0,XFilledProjectionHist10=0,YFilledProjectionHist10=0,Barycenter2=0,DirectionPerpendicular5Hist10=0,CurvaturePerpendicular5Hist100=0,luminanceHist256=0,CurvatureAli5Hist100=0,CurvaturesDerivateAli5Hist100=0,CurvatureAli10Hist100=0,CurvaturesDerivateAli10Hist100=0,CurvatureAli15Hist100=0,CurvaturesDerivateAli15Hist100=0,CurvatureAli20Hist100=0,CurvaturesDerivateAli20Hist100=0,chaincodeHist_4=0,chaincodeHist_8=0,chaincode8order2_64=0,chaincode4order2_16=0,chaincode4order3_64=0,chaincode8order3_512=0,chaincode4order4_256=0,chaincode8order4_4096=0,Fourier_1=0,Fourier_5=0,Fourier_9=0,Fourier_15=0;
NumberOfConnectedComponents_1+=fabs(my_questioned_document_features[document_index1]->NumberOfConnectedComponents_1[0]-my_document_features[document_index2]->NumberOfConnectedComponents_1[0]);
NumberOfHoles_1+=fabs(my_questioned_document_features[document_index1]->NumberOfHoles_1[0]-my_document_features[document_index2]->NumberOfHoles_1[0]);
SpatialMoment00+=fabs(my_questioned_document_features[document_index1]->SpatialMoment00-my_document_features[document_index2]->SpatialMoment00);
SpatialMoment10+=fabs(my_questioned_document_features[document_index1]->SpatialMoment10-my_document_features[document_index2]->SpatialMoment10);
SpatialMoment20+=fabs(my_questioned_document_features[document_index1]->SpatialMoment20-my_document_features[document_index2]->SpatialMoment20);
SpatialMoment30+=fabs(my_questioned_document_features[document_index1]->SpatialMoment30-my_document_features[document_index2]->SpatialMoment30);
SpatialMoment01+=fabs(my_questioned_document_features[document_index1]->SpatialMoment01-my_document_features[document_index2]->SpatialMoment01);
SpatialMoment11+=fabs(my_questioned_document_features[document_index1]->SpatialMoment11-my_document_features[document_index2]->SpatialMoment11);
SpatialMoment21+=fabs(my_questioned_document_features[document_index1]->SpatialMoment21-my_document_features[document_index2]->SpatialMoment21);
SpatialMoment02+=fabs(my_questioned_document_features[document_index1]->SpatialMoment02-my_document_features[document_index2]->SpatialMoment02);
SpatialMoment12+=fabs(my_questioned_document_features[document_index1]->SpatialMoment12-my_document_features[document_index2]->SpatialMoment12);
SpatialMoment03+=fabs(my_questioned_document_features[document_index1]->SpatialMoment03-my_document_features[document_index2]->SpatialMoment03);
CentralMoment00+=fabs(my_questioned_document_features[document_index1]->CentralMoment00-my_document_features[document_index2]->CentralMoment00);
CentralMoment10+=fabs(my_questioned_document_features[document_index1]->CentralMoment10-my_document_features[document_index2]->CentralMoment10);
CentralMoment20+=fabs(my_questioned_document_features[document_index1]->CentralMoment20-my_document_features[document_index2]->CentralMoment20);
CentralMoment30+=fabs(my_questioned_document_features[document_index1]->CentralMoment30-my_document_features[document_index2]->CentralMoment30);
CentralMoment01+=fabs(my_questioned_document_features[document_index1]->CentralMoment01-my_document_features[document_index2]->CentralMoment01);
CentralMoment11+=fabs(my_questioned_document_features[document_index1]->CentralMoment11-my_document_features[document_index2]->CentralMoment11);
CentralMoment21+=fabs(my_questioned_document_features[document_index1]->CentralMoment21-my_document_features[document_index2]->CentralMoment21);
CentralMoment02+=fabs(my_questioned_document_features[document_index1]->CentralMoment02-my_document_features[document_index2]->CentralMoment02);
CentralMoment12+=fabs(my_questioned_document_features[document_index1]->CentralMoment12-my_document_features[document_index2]->CentralMoment12);
CentralMoment03+=fabs(my_questioned_document_features[document_index1]->CentralMoment03-my_document_features[document_index2]->CentralMoment03);
NormalizedCentralMoment00+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment00-my_document_features[document_index2]->NormalizedCentralMoment00);
NormalizedCentralMoment10+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment10-my_document_features[document_index2]->NormalizedCentralMoment10);
NormalizedCentralMoment20+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment20-my_document_features[document_index2]->NormalizedCentralMoment20);
NormalizedCentralMoment30+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment30-my_document_features[document_index2]->NormalizedCentralMoment30);
NormalizedCentralMoment01+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment01-my_document_features[document_index2]->NormalizedCentralMoment01);
NormalizedCentralMoment11+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment11-my_document_features[document_index2]->NormalizedCentralMoment11);
NormalizedCentralMoment21+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment21-my_document_features[document_index2]->NormalizedCentralMoment21);
NormalizedCentralMoment02+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment02-my_document_features[document_index2]->NormalizedCentralMoment02);
NormalizedCentralMoment12+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment12-my_document_features[document_index2]->NormalizedCentralMoment12);
NormalizedCentralMoment03+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment03-my_document_features[document_index2]->NormalizedCentralMoment03);
HuMoment1+=fabs(my_questioned_document_features[document_index1]->HuMoment1-my_document_features[document_index2]->HuMoment1);
HuMoment2+=fabs(my_questioned_document_features[document_index1]->HuMoment2-my_document_features[document_index2]->HuMoment2);
HuMoment3+=fabs(my_questioned_document_features[document_index1]->HuMoment3-my_document_features[document_index2]->HuMoment3);
HuMoment4+=fabs(my_questioned_document_features[document_index1]->HuMoment4-my_document_features[document_index2]->HuMoment4);
HuMoment5+=fabs(my_questioned_document_features[document_index1]->HuMoment5-my_document_features[document_index2]->HuMoment5);
HuMoment6+=fabs(my_questioned_document_features[document_index1]->HuMoment6-my_document_features[document_index2]->HuMoment6);
HuMoment7+=fabs(my_questioned_document_features[document_index1]->HuMoment7-my_document_features[document_index2]->HuMoment7);
for(int x=0;x<10;x++)
XProjectionHist10+=fabs(my_questioned_document_features[document_index1]->XProjectionHist10[x]-my_document_features[document_index2]->XProjectionHist10[x]);
for(int x=0;x<10;x++)
YProjectionHist10+=fabs(my_questioned_document_features[document_index1]->YProjectionHist10[x]-my_document_features[document_index2]->YProjectionHist10[x]);
for(int x=0;x<10;x++)
XFilledProjectionHist10+=fabs(my_questioned_document_features[document_index1]->XFilledProjectionHist10[x]-my_document_features[document_index2]->XFilledProjectionHist10[x]);
for(int x=0;x<10;x++)
YFilledProjectionHist10+=fabs(my_questioned_document_features[document_index1]->YFilledProjectionHist10[x]-my_document_features[document_index2]->YFilledProjectionHist10[x]);
for(int x=0;x<10;x++)
Distribution10x10_100+=fabs(my_questioned_document_features[document_index1]->Distribution10x10_100[x]-my_document_features[document_index2]->Distribution10x10_100[x]);
for(int x=0;x<2;x++)
Barycenter2+=fabs(my_questioned_document_features[document_index1]->Barycenter2[x]-my_document_features[document_index2]->Barycenter2[x]);
numberOfPixels1+=fabs(my_questioned_document_features[document_index1]->numberOfPixels1[0]-my_document_features[document_index2]->numberOfPixels1[0]);
for(int x=0;x<1;x++)
Fourier_1+=fabs(my_questioned_document_features[document_index1]->Fourier_1[x]-my_document_features[document_index2]->Fourier_1[x]);
for(int x=0;x<5;x++)
Fourier_5+=fabs(my_questioned_document_features[document_index1]->Fourier_5[x]-my_document_features[document_index2]->Fourier_5[x]);
for(int x=0;x<9;x++)
Fourier_9+=fabs(my_questioned_document_features[document_index1]->Fourier_9[x]-my_document_features[document_index2]->Fourier_9[x]);
for(int x=0;x<15;x++)
Fourier_15+=fabs(my_questioned_document_features[document_index1]->Fourier_15[x]-my_document_features[document_index2]->Fourier_15[x]);
nbbranches_1+=fabs(my_questioned_document_features[document_index1]->nbbranches_1[0]-my_document_features[document_index2]->nbbranches_1[0]);
for(int x=0;x<10;x++)
LengthsOfBranchesHist_10+=fabs(my_questioned_document_features[document_index1]->LengthsOfBranchesHist_10[x]-my_document_features[document_index2]->LengthsOfBranchesHist_10[x]);
for(int x=0;x<30;x++)
ThicknessLengthsCircleHist30+=fabs(my_questioned_document_features[document_index1]->ThicknessLengthsCircleHist30[x]-my_document_features[document_index2]->ThicknessLengthsCircleHist30[x]);
for(int x=0;x<10;x++)
tortuosityHist10+=fabs(my_questioned_document_features[document_index1]->tortuosityHist10[x]-my_document_features[document_index2]->tortuosityHist10[x]);
for(int x=0;x<10;x++)
tortuosityDirectionHist10+=fabs(my_questioned_document_features[document_index1]->tortuosityDirectionHist10[x]-my_document_features[document_index2]->tortuosityDirectionHist10[x]);
for(int x=0;x<10;x++)
tortuosityDerivateHist10+=fabs(my_questioned_document_features[document_index1]->tortuosityDerivateHist10[x]-my_document_features[document_index2]->tortuosityDerivateHist10[x]);
for(int x=0;x<10;x++)
tortuosityDerivateDirectionHist10+=fabs(my_questioned_document_features[document_index1]->tortuosityDerivateDirectionHist10[x]-my_document_features[document_index2]->tortuosityDerivateDirectionHist10[x]);
for(int x=0;x<10;x++)
DirectionPerpendicular5Hist10+=fabs(my_questioned_document_features[document_index1]->DirectionPerpendicular5Hist10[x]-my_document_features[document_index2]->DirectionPerpendicular5Hist10[x]);
for(int x=0;x<100;x++)
CurvaturePerpendicular5Hist100+=fabs(my_questioned_document_features[document_index1]->CurvaturePerpendicular5Hist100[x]-my_document_features[document_index2]->CurvaturePerpendicular5Hist100[x]);
for(int x=0;x<256;x++)
luminanceHist256+=fabs(my_questioned_document_features[document_index1]->luminanceHist256[x]-my_document_features[document_index2]->luminanceHist256[x]);
for(int x=0;x<100;x++)
CurvatureAli5Hist100+=fabs(my_questioned_document_features[document_index1]->CurvatureAli5Hist100[x]-my_document_features[document_index2]->CurvatureAli5Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli5Hist100+=fabs(my_questioned_document_features[document_index1]->CurvaturesDerivateAli5Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli5Hist100[x]);
for(int x=0;x<100;x++)
CurvatureAli10Hist100+=fabs(my_questioned_document_features[document_index1]->CurvatureAli10Hist100[x]-my_document_features[document_index2]->CurvatureAli10Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli10Hist100+=fabs(my_questioned_document_features[document_index1]->CurvaturesDerivateAli10Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli10Hist100[x]);
for(int x=0;x<100;x++)
CurvatureAli15Hist100+=fabs(my_questioned_document_features[document_index1]->CurvatureAli15Hist100[x]-my_document_features[document_index2]->CurvatureAli15Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli15Hist100+=fabs(my_questioned_document_features[document_index1]->CurvaturesDerivateAli15Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli15Hist100[x]);
for(int x=0;x<100;x++)
CurvatureAli20Hist100+=fabs(my_questioned_document_features[document_index1]->CurvatureAli20Hist100[x]-my_document_features[document_index2]->CurvatureAli20Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli20Hist100+=fabs(my_questioned_document_features[document_index1]->CurvaturesDerivateAli20Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli20Hist100[x]);
for(int x=0;x<4;x++)
chaincodeHist_4+=fabs(my_questioned_document_features[document_index1]->chaincodeHist_4[x]-my_document_features[document_index2]->chaincodeHist_4[x]);
for(int x=0;x<8;x++)
chaincodeHist_8+=fabs(my_questioned_document_features[document_index1]->chaincodeHist_8[x]-my_document_features[document_index2]->chaincodeHist_8[x]);
for(int x=0;x<64;x++)
chaincode8order2_64+=fabs(my_questioned_document_features[document_index1]->chaincode8order2_64[x]-my_document_features[document_index2]->chaincode8order2_64[x]);
for(int x=0;x<16;x++)
chaincode4order2_16+=fabs(my_questioned_document_features[document_index1]->chaincode4order2_16[x]-my_document_features[document_index2]->chaincode4order2_16[x]);
for(int x=0;x<64;x++)
chaincode4order3_64+=fabs(my_questioned_document_features[document_index1]->chaincode4order3_64[x]-my_document_features[document_index2]->chaincode4order3_64[x]);
for(int x=0;x<512;x++)
chaincode8order3_512+=fabs(my_questioned_document_features[document_index1]->chaincode8order3_512[x]-my_document_features[document_index2]->chaincode8order3_512[x]);
for(int x=0;x<256;x++)
chaincode4order4_256+=fabs(my_questioned_document_features[document_index1]->chaincode4order4_256[x]-my_document_features[document_index2]->chaincode4order4_256[x]);
for(int x=0;x<4096;x++)
chaincode8order4_4096+=fabs(my_questioned_document_features[document_index1]->chaincode8order4_4096[x]-my_document_features[document_index2]->chaincode8order4_4096[x]);
//using the model previously built
double linear=(-10.667528249636) + -0.0667090106640684 * NumberOfConnectedComponents_1 +9.34910163750235 * ThicknessLengthsCircleHist30 + 30.4471654988069 *CurvaturePerpendicular5Hist100 + 30.5600351874986 * luminanceHist256 +33.6853635484793 * chaincode8order3_512;
double prob_same_writer=1. / (1. + exp(linear));
//the most probable writer
if(prob_same_writer>maxprob)
{
maxprob=prob_same_writer;
my_questioned_document_features[document_index1]->writer_id=my_document_features[document_index2]->writer_id;
}
}
}
//creating the file to submit
//this benchmark does not provide a way of detecting unknown writers
std::ofstream submission_file(\"new_benchmark.csv\");
submission_file<<\",031,034,035,036,037,038,039,041,042,046,047,048,049,050,051,052,053,054,055,056,057,058,059,060,061,062,063,064,065,066,067,068,069,070,071,072,073,074,077,078,080,081,082,083,084,087,088,090,091,092,093,095,096,unknown\\n\";
int writers[]={31,34,35,36,37,38,39,41,42,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,77,78,80,81,82,83,84,87,88,90,91,92,93,95,96};
string questioned_writers[]={\"AA\",\"AB\",\"AC\",\"AD\",\"AE\",\"AF\",\"AG\",\"AH\",\"AI\",\"AJ\",\"AK\",\"AM\",\"AN\",\"AO\",\"AP\",\"AQ\",\"AR\",\"AS\",\"AT\",\"AU\",\"AV\",\"AW\",\"AX\",\"AY\",\"AZ\",\"BA\",\"BB\",\"BC\",\"BD\",\"BE\",\"BF\",\"BG\",\"BH\",\"BI\",\"BK\",\"BL\",\"BM\",\"BN\",\"BO\",\"BP\",\"BQ\",\"BR\",\"BS\",\"BT\",\"BU\",\"BV\",\"BW\",\"BX\",\"BY\",\"BZ\",\"CA\",\"CB\",\"CC\"};
for(int document_index1=0;document_index1<53;document_index1++)
{
submission_file<<questioned_writers[document_index1];
bool at_least_one=false;
for(int writer_index=0;writer_index<53;writer_index++)
{
if(my_questioned_document_features[document_index1]->writer_id==writers[writer_index])
{
at_least_one=true;
submission_file<<\",1\";
}
else submission_file<<\",0\";
}
//unknown writer
if(at_least_one)submission_file<<\",0\\n\";
else submission_file<<\",1\\n\";
}
submission_file.close();
std::cout<<\"Press a key to exit...\";
_getch();
}
",
"messages": [
{
"id": "1993",
"postDate": "03/13/2011 17:45:09",
"content": "Dear participants,
Below is the c++ code used to generate the benchmark.
Comments are welcome.
#include \"stdafx.h\"
#include <fstream>
#include <conio.h>
#include <list>
#include <iostream>
#include <sstream>
#include <vector>
#include <string>
#include <math.h>
using namespace std;
//class handling the provided features
class document_features
{
public:
document_features()
{
ThicknessLengthsCircleHist30=new double[30];
XProjectionHist10=new double[10];
YProjectionHist10=new double[10];
Distribution10x10_100=new double[100];
tortuosityHist10=new double[10];
tortuosityDirectionHist10=new double[10];
tortuosityDerivateHist10=new double[10];
tortuosityDerivateDirectionHist10=new double[10];
nbbranches_1=new double[1];
LengthsOfBranchesHist_10=new double[10];
numberOfPixels1=new double[1];
XFilledProjectionHist10=new double[10];
YFilledProjectionHist10=new double[10];
Barycenter2=new double[2];
DirectionPerpendicular5Hist10=new double[10];
CurvaturePerpendicular5Hist100=new double[100];
luminanceHist256=new double[256];
CurvatureAli5Hist100=new double[100];
CurvaturesDerivateAli5Hist100=new double[100];
CurvatureAli10Hist100=new double[100];
CurvaturesDerivateAli10Hist100=new double[100];
CurvatureAli15Hist100=new double[100];
CurvaturesDerivateAli15Hist100=new double[100];
CurvatureAli20Hist100=new double[100];
CurvaturesDerivateAli20Hist100=new double[100];
chaincodeHist_4=new double[4];
chaincodeHist_8=new double[8];
chaincode8order2_64=new double[64];
chaincode4order2_16=new double[16];
chaincode4order3_64=new double[64];
chaincode8order3_512=new double[512];
chaincode4order4_256=new double[256];
chaincode8order4_4096=new double[4096];
NumberOfConnectedComponents_1=new double[1];
NumberOfHoles_1=new double[1];
Fourier_1=new double[1];
Fourier_5=new double[5];
Fourier_9=new double[9];
Fourier_15=new double[15];
}
double*ThicknessLengthsCircleHist30,*XProjectionHist10,*YProjectionHist10,*Distribution10x10_100,*tortuosityHist10,*tortuosityDirectionHist10,*tortuosityDerivateHist10,*tortuosityDerivateDirectionHist10,*nbbranches_1,*LengthsOfBranchesHist_10,*numberOfPixels1,*XFilledProjectionHist10,*YFilledProjectionHist10,*Barycenter2,*DirectionPerpendicular5Hist10,*CurvaturePerpendicular5Hist100,*luminanceHist256,*CurvatureAli5Hist100,*CurvaturesDerivateAli5Hist100,*CurvatureAli10Hist100,*CurvaturesDerivateAli10Hist100,*CurvatureAli15Hist100,*CurvaturesDerivateAli15Hist100,*CurvatureAli20Hist100,*CurvaturesDerivateAli20Hist100,*chaincodeHist_4,*chaincodeHist_8,*chaincode8order2_64,*chaincode4order2_16,*chaincode4order3_64,*chaincode8order3_512,*chaincode4order4_256,*chaincode8order4_4096,*NumberOfConnectedComponents_1,*NumberOfHoles_1,*Fourier_1,*Fourier_5,*Fourier_9,*Fourier_15;
double SpatialMoment00,SpatialMoment10,SpatialMoment20,SpatialMoment30,SpatialMoment01,SpatialMoment11,SpatialMoment21,SpatialMoment02,SpatialMoment12,SpatialMoment03,CentralMoment00,CentralMoment10,CentralMoment20,CentralMoment30,CentralMoment01,CentralMoment11,CentralMoment21,CentralMoment02,CentralMoment12,CentralMoment03,NormalizedCentralMoment00,NormalizedCentralMoment10,NormalizedCentralMoment20,NormalizedCentralMoment30,NormalizedCentralMoment01,NormalizedCentralMoment11,NormalizedCentralMoment21,NormalizedCentralMoment02,NormalizedCentralMoment12,NormalizedCentralMoment03,HuMoment1,HuMoment2,HuMoment3,HuMoment4,HuMoment5,HuMoment6,HuMoment7;
int writer_id;
};
void main()
{
//training phase
std::ifstream data(\"train.csv\");
int linenum=0;
string line;
document_features**my_document_features=new document_features*[108];
string item;//=new char[4];
//reading the training documents
while (getline (data, line))
{
cout << \"\\nLine #\" << linenum << \":\" << endl;
istringstream linestream(line);
if(linenum>0)
{
int row_index=linenum-1;
my_document_features[row_index]=new document_features();
getline (linestream, item, ',');
my_document_features[row_index]->writer_id=atoi(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NumberOfConnectedComponents_1[0]=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NumberOfHoles_1[0]=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment00=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment10=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment20=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment30=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment01=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment11=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment21=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment02=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment12=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->SpatialMoment03=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment00=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment10=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment20=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment30=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment01=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment11=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment21=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment02=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment12=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->CentralMoment03=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment00=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment10=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment20=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment30=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment01=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment11=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment21=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment02=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment12=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->NormalizedCentralMoment03=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment1=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment2=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment3=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment4=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment5=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment6=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->HuMoment7=atof(item.c_str());
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->XProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->YProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->XFilledProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->YFilledProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->Distribution10x10_100[x]=atof(item.c_str());
}
getline (linestream, item, ',');
my_document_features[row_index]->Barycenter2[0]=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->Barycenter2[1]=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->numberOfPixels1[0]=atof(item.c_str());
getline (linestream, item, ',');
my_document_features[row_index]->Fourier_1[0]=atof(item.c_str());
for(int x=0;x<5;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->Fourier_5[x]=atof(item.c_str());
}
for(int x=0;x<9;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->Fourier_9[x]=atof(item.c_str());
}
for(int x=0;x<15;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->Fourier_15[x]=atof(item.c_str());
}
getline (linestream, item, ',');
my_document_features[row_index]->nbbranches_1[0]=atof(item.c_str());
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->LengthsOfBranchesHist_10[x]=atof(item.c_str());
}
for(int x=0;x<30;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->ThicknessLengthsCircleHist30[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->tortuosityHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->tortuosityDirectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->tortuosityDerivateHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->tortuosityDerivateDirectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->DirectionPerpendicular5Hist10[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvaturePerpendicular5Hist100[x]=atof(item.c_str());
}
for(int x=0;x<256;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->luminanceHist256[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvatureAli5Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvaturesDerivateAli5Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvatureAli10Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvaturesDerivateAli10Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvatureAli15Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvaturesDerivateAli15Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvatureAli20Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->CurvaturesDerivateAli20Hist100[x]=atof(item.c_str());
}
for(int x=0;x<4;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincodeHist_4[x]=atof(item.c_str());
}
for(int x=0;x<8;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincodeHist_8[x]=atof(item.c_str());
}
for(int x=0;x<64;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincode8order2_64[x]=atof(item.c_str());
}
for(int x=0;x<16;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincode4order2_16[x]=atof(item.c_str());
}
for(int x=0;x<64;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincode4order3_64[x]=atof(item.c_str());
}
for(int x=0;x<512;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincode8order3_512[x]=atof(item.c_str());
}
for(int x=0;x<256;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincode4order4_256[x]=atof(item.c_str());
}
for(int x=0;x<4096;x++)
{
getline (linestream, item, ',');
my_document_features[row_index]->chaincode8order4_4096[x]=atof(item.c_str());
}
}
linenum++;
}
//difference between features of all documents in the training set
std::ofstream differences(\"train_differences.csv\");
differences<<\"SameWriter,NumberOfConnectedComponents_1,NumberOfHoles_1,SpatialMoment00,SpatialMoment10,SpatialMoment20,SpatialMoment30,SpatialMoment01,SpatialMoment11,SpatialMoment21,SpatialMoment02,SpatialMoment12,SpatialMoment03,CentralMoment00,CentralMoment10,CentralMoment20,CentralMoment30,CentralMoment01,CentralMoment11,CentralMoment21,CentralMoment02,CentralMoment12,CentralMoment03,NormalizedCentralMoment00,NormalizedCentralMoment10,NormalizedCentralMoment20,NormalizedCentralMoment30,NormalizedCentralMoment01,NormalizedCentralMoment11,NormalizedCentralMoment21,NormalizedCentralMoment02,NormalizedCentralMoment12,NormalizedCentralMoment03,HuMoment1,HuMoment2,HuMoment3,HuMoment4,HuMoment5,HuMoment6,HuMoment7,XProjectionHist10,YProjectionHist10,XFilledProjectionHist10,YFilledProjectionHist10,Distribution10x10_100,Barycenter2,numberOfPixels1,Fourier_1,Fourier_5,Fourier_9,Fourier_15,nbbranches_1,LengthsOfBranchesHist_10,ThicknessLengthsCircleHist30,tortuosityHist10,tortuosityDirectionHist10,tortuosityDerivateHist10,tortuosityDerivateDirectionHist10,DirectionPerpendicular5Hist10,CurvaturePerpendicular5Hist100,luminanceHist256,CurvatureAli5Hist100,CurvaturesDerivateAli5Hist100,CurvatureAli10Hist100,CurvaturesDerivateAli10Hist100,CurvatureAli15Hist100,CurvaturesDerivateAli15Hist100,CurvatureAli20Hist100,CurvaturesDerivateAli20Hist100,chaincodeHist_4,chaincodeHist_8,chaincode8order2_64,chaincode4order2_16,chaincode4order3_64,chaincode8order3_512,chaincode4order4_256,chaincode8order4_4096\\n\";
//comparing each pair of documents
for(int document_index1=0;document_index1<107;document_index1++)
for(int document_index2=document_index1+1;document_index2<108;document_index2++)
{
double NumberOfConnectedComponents_1=0,NumberOfHoles_1=0,SpatialMoment00=0,SpatialMoment10=0,SpatialMoment20=0,SpatialMoment30=0,SpatialMoment01=0,SpatialMoment11=0,SpatialMoment21=0,SpatialMoment02=0,SpatialMoment12=0,SpatialMoment03=0,CentralMoment00=0,CentralMoment10=0,CentralMoment20=0,CentralMoment30=0,CentralMoment01=0,CentralMoment11=0,CentralMoment21=0,CentralMoment02=0,CentralMoment12=0,CentralMoment03=0,NormalizedCentralMoment00=0,NormalizedCentralMoment10=0,NormalizedCentralMoment20=0,NormalizedCentralMoment30=0,NormalizedCentralMoment01=0,NormalizedCentralMoment11=0,NormalizedCentralMoment21=0,NormalizedCentralMoment02=0,NormalizedCentralMoment12=0,NormalizedCentralMoment03=0,HuMoment1=0,HuMoment2=0,HuMoment3=0,HuMoment4=0,HuMoment5=0,HuMoment6=0,HuMoment7=0,ThicknessLengthsCircleHist30=0,XProjectionHist10=0,YProjectionHist10=0,Distribution10x10_100=0,tortuosityHist10=0,tortuosityDirectionHist10=0,tortuosityDerivateHist10=0,tortuosityDerivateDirectionHist10=0,nbbranches_1=0,LengthsOfBranchesHist_10=0,numberOfPixels1=0,XFilledProjectionHist10=0,YFilledProjectionHist10=0,Barycenter2=0,DirectionPerpendicular5Hist10=0,CurvaturePerpendicular5Hist100=0,luminanceHist256=0,CurvatureAli5Hist100=0,CurvaturesDerivateAli5Hist100=0,CurvatureAli10Hist100=0,CurvaturesDerivateAli10Hist100=0,CurvatureAli15Hist100=0,CurvaturesDerivateAli15Hist100=0,CurvatureAli20Hist100=0,CurvaturesDerivateAli20Hist100=0,chaincodeHist_4=0,chaincodeHist_8=0,chaincode8order2_64=0,chaincode4order2_16=0,chaincode4order3_64=0,chaincode8order3_512=0,chaincode4order4_256=0,chaincode8order4_4096=0,Fourier_1=0,Fourier_5=0,Fourier_9=0,Fourier_15=0;
NumberOfConnectedComponents_1+=fabs(my_document_features[document_index1]->NumberOfConnectedComponents_1[0]-my_document_features[document_index2]->NumberOfConnectedComponents_1[0]);
NumberOfHoles_1+=fabs(my_document_features[document_index1]->NumberOfHoles_1[0]-my_document_features[document_index2]->NumberOfHoles_1[0]);
SpatialMoment00+=fabs(my_document_features[document_index1]->SpatialMoment00-my_document_features[document_index2]->SpatialMoment00);
SpatialMoment10+=fabs(my_document_features[document_index1]->SpatialMoment10-my_document_features[document_index2]->SpatialMoment10);
SpatialMoment20+=fabs(my_document_features[document_index1]->SpatialMoment20-my_document_features[document_index2]->SpatialMoment20);
SpatialMoment30+=fabs(my_document_features[document_index1]->SpatialMoment30-my_document_features[document_index2]->SpatialMoment30);
SpatialMoment01+=fabs(my_document_features[document_index1]->SpatialMoment01-my_document_features[document_index2]->SpatialMoment01);
SpatialMoment11+=fabs(my_document_features[document_index1]->SpatialMoment11-my_document_features[document_index2]->SpatialMoment11);
SpatialMoment21+=fabs(my_document_features[document_index1]->SpatialMoment21-my_document_features[document_index2]->SpatialMoment21);
SpatialMoment02+=fabs(my_document_features[document_index1]->SpatialMoment02-my_document_features[document_index2]->SpatialMoment02);
SpatialMoment12+=fabs(my_document_features[document_index1]->SpatialMoment12-my_document_features[document_index2]->SpatialMoment12);
SpatialMoment03+=fabs(my_document_features[document_index1]->SpatialMoment03-my_document_features[document_index2]->SpatialMoment03);
CentralMoment00+=fabs(my_document_features[document_index1]->CentralMoment00-my_document_features[document_index2]->CentralMoment00);
CentralMoment10+=fabs(my_document_features[document_index1]->CentralMoment10-my_document_features[document_index2]->CentralMoment10);
CentralMoment20+=fabs(my_document_features[document_index1]->CentralMoment20-my_document_features[document_index2]->CentralMoment20);
CentralMoment30+=fabs(my_document_features[document_index1]->CentralMoment30-my_document_features[document_index2]->CentralMoment30);
CentralMoment01+=fabs(my_document_features[document_index1]->CentralMoment01-my_document_features[document_index2]->CentralMoment01);
CentralMoment11+=fabs(my_document_features[document_index1]->CentralMoment11-my_document_features[document_index2]->CentralMoment11);
CentralMoment21+=fabs(my_document_features[document_index1]->CentralMoment21-my_document_features[document_index2]->CentralMoment21);
CentralMoment02+=fabs(my_document_features[document_index1]->CentralMoment02-my_document_features[document_index2]->CentralMoment02);
CentralMoment12+=fabs(my_document_features[document_index1]->CentralMoment12-my_document_features[document_index2]->CentralMoment12);
CentralMoment03+=fabs(my_document_features[document_index1]->CentralMoment03-my_document_features[document_index2]->CentralMoment03);
NormalizedCentralMoment00+=fabs(my_document_features[document_index1]->NormalizedCentralMoment00-my_document_features[document_index2]->NormalizedCentralMoment00);
NormalizedCentralMoment10+=fabs(my_document_features[document_index1]->NormalizedCentralMoment10-my_document_features[document_index2]->NormalizedCentralMoment10);
NormalizedCentralMoment20+=fabs(my_document_features[document_index1]->NormalizedCentralMoment20-my_document_features[document_index2]->NormalizedCentralMoment20);
NormalizedCentralMoment30+=fabs(my_document_features[document_index1]->NormalizedCentralMoment30-my_document_features[document_index2]->NormalizedCentralMoment30);
NormalizedCentralMoment01+=fabs(my_document_features[document_index1]->NormalizedCentralMoment01-my_document_features[document_index2]->NormalizedCentralMoment01);
NormalizedCentralMoment11+=fabs(my_document_features[document_index1]->NormalizedCentralMoment11-my_document_features[document_index2]->NormalizedCentralMoment11);
NormalizedCentralMoment21+=fabs(my_document_features[document_index1]->NormalizedCentralMoment21-my_document_features[document_index2]->NormalizedCentralMoment21);
NormalizedCentralMoment02+=fabs(my_document_features[document_index1]->NormalizedCentralMoment02-my_document_features[document_index2]->NormalizedCentralMoment02);
NormalizedCentralMoment12+=fabs(my_document_features[document_index1]->NormalizedCentralMoment12-my_document_features[document_index2]->NormalizedCentralMoment12);
NormalizedCentralMoment03+=fabs(my_document_features[document_index1]->NormalizedCentralMoment03-my_document_features[document_index2]->NormalizedCentralMoment03);
HuMoment1+=fabs(my_document_features[document_index1]->HuMoment1-my_document_features[document_index2]->HuMoment1);
HuMoment2+=fabs(my_document_features[document_index1]->HuMoment2-my_document_features[document_index2]->HuMoment2);
HuMoment3+=fabs(my_document_features[document_index1]->HuMoment3-my_document_features[document_index2]->HuMoment3);
HuMoment4+=fabs(my_document_features[document_index1]->HuMoment4-my_document_features[document_index2]->HuMoment4);
HuMoment5+=fabs(my_document_features[document_index1]->HuMoment5-my_document_features[document_index2]->HuMoment5);
HuMoment6+=fabs(my_document_features[document_index1]->HuMoment6-my_document_features[document_index2]->HuMoment6);
HuMoment7+=fabs(my_document_features[document_index1]->HuMoment7-my_document_features[document_index2]->HuMoment7);
for(int x=0;x<10;x++)
XProjectionHist10+=fabs(my_document_features[document_index1]->XProjectionHist10[x]-my_document_features[document_index2]->XProjectionHist10[x]);
for(int x=0;x<10;x++)
YProjectionHist10+=fabs(my_document_features[document_index1]->YProjectionHist10[x]-my_document_features[document_index2]->YProjectionHist10[x]);
for(int x=0;x<10;x++)
XFilledProjectionHist10+=fabs(my_document_features[document_index1]->XFilledProjectionHist10[x]-my_document_features[document_index2]->XFilledProjectionHist10[x]);
for(int x=0;x<10;x++)
YFilledProjectionHist10+=fabs(my_document_features[document_index1]->YFilledProjectionHist10[x]-my_document_features[document_index2]->YFilledProjectionHist10[x]);
for(int x=0;x<10;x++)
Distribution10x10_100+=fabs(my_document_features[document_index1]->Distribution10x10_100[x]-my_document_features[document_index2]->Distribution10x10_100[x]);
for(int x=0;x<2;x++)
Barycenter2+=fabs(my_document_features[document_index1]->Barycenter2[x]-my_document_features[document_index2]->Barycenter2[x]);
numberOfPixels1+=fabs(my_document_features[document_index1]->numberOfPixels1[0]-my_document_features[document_index2]->numberOfPixels1[0]);
for(int x=0;x<1;x++)
Fourier_1+=fabs(my_document_features[document_index1]->Fourier_1[x]-my_document_features[document_index2]->Fourier_1[x]);
for(int x=0;x<5;x++)
Fourier_5+=fabs(my_document_features[document_index1]->Fourier_5[x]-my_document_features[document_index2]->Fourier_5[x]);
for(int x=0;x<9;x++)
Fourier_9+=fabs(my_document_features[document_index1]->Fourier_9[x]-my_document_features[document_index2]->Fourier_9[x]);
for(int x=0;x<15;x++)
Fourier_15+=fabs(my_document_features[document_index1]->Fourier_15[x]-my_document_features[document_index2]->Fourier_15[x]);
nbbranches_1+=fabs(my_document_features[document_index1]->nbbranches_1[0]-my_document_features[document_index2]->nbbranches_1[0]);
for(int x=0;x<10;x++)
LengthsOfBranchesHist_10+=fabs(my_document_features[document_index1]->LengthsOfBranchesHist_10[x]-my_document_features[document_index2]->LengthsOfBranchesHist_10[x]);
for(int x=0;x<30;x++)
ThicknessLengthsCircleHist30+=fabs(my_document_features[document_index1]->ThicknessLengthsCircleHist30[x]-my_document_features[document_index2]->ThicknessLengthsCircleHist30[x]);
for(int x=0;x<10;x++)
tortuosityHist10+=fabs(my_document_features[document_index1]->tortuosityHist10[x]-my_document_features[document_index2]->tortuosityHist10[x]);
for(int x=0;x<10;x++)
tortuosityDirectionHist10+=fabs(my_document_features[document_index1]->tortuosityDirectionHist10[x]-my_document_features[document_index2]->tortuosityDirectionHist10[x]);
for(int x=0;x<10;x++)
tortuosityDerivateHist10+=fabs(my_document_features[document_index1]->tortuosityDerivateHist10[x]-my_document_features[document_index2]->tortuosityDerivateHist10[x]);
for(int x=0;x<10;x++)
tortuosityDerivateDirectionHist10+=fabs(my_document_features[document_index1]->tortuosityDerivateDirectionHist10[x]-my_document_features[document_index2]->tortuosityDerivateDirectionHist10[x]);
for(int x=0;x<10;x++)
DirectionPerpendicular5Hist10+=fabs(my_document_features[document_index1]->DirectionPerpendicular5Hist10[x]-my_document_features[document_index2]->DirectionPerpendicular5Hist10[x]);
for(int x=0;x<100;x++)
CurvaturePerpendicular5Hist100+=fabs(my_document_features[document_index1]->CurvaturePerpendicular5Hist100[x]-my_document_features[document_index2]->CurvaturePerpendicular5Hist100[x]);
for(int x=0;x<256;x++)
luminanceHist256+=fabs(my_document_features[document_index1]->luminanceHist256[x]-my_document_features[document_index2]->luminanceHist256[x]);
for(int x=0;x<100;x++)
CurvatureAli5Hist100+=fabs(my_document_features[document_index1]->CurvatureAli5Hist100[x]-my_document_features[document_index2]->CurvatureAli5Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli5Hist100+=fabs(my_document_features[document_index1]->CurvaturesDerivateAli5Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli5Hist100[x]);
for(int x=0;x<100;x++)
CurvatureAli10Hist100+=fabs(my_document_features[document_index1]->CurvatureAli10Hist100[x]-my_document_features[document_index2]->CurvatureAli10Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli10Hist100+=fabs(my_document_features[document_index1]->CurvaturesDerivateAli10Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli10Hist100[x]);
for(int x=0;x<100;x++)
CurvatureAli15Hist100+=fabs(my_document_features[document_index1]->CurvatureAli15Hist100[x]-my_document_features[document_index2]->CurvatureAli15Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli15Hist100+=fabs(my_document_features[document_index1]->CurvaturesDerivateAli15Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli15Hist100[x]);
for(int x=0;x<100;x++)
CurvatureAli20Hist100+=fabs(my_document_features[document_index1]->CurvatureAli20Hist100[x]-my_document_features[document_index2]->CurvatureAli20Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli20Hist100+=fabs(my_document_features[document_index1]->CurvaturesDerivateAli20Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli20Hist100[x]);
for(int x=0;x<4;x++)
chaincodeHist_4+=fabs(my_document_features[document_index1]->chaincodeHist_4[x]-my_document_features[document_index2]->chaincodeHist_4[x]);
for(int x=0;x<8;x++)
chaincodeHist_8+=fabs(my_document_features[document_index1]->chaincodeHist_8[x]-my_document_features[document_index2]->chaincodeHist_8[x]);
for(int x=0;x<64;x++)
chaincode8order2_64+=fabs(my_document_features[document_index1]->chaincode8order2_64[x]-my_document_features[document_index2]->chaincode8order2_64[x]);
for(int x=0;x<16;x++)
chaincode4order2_16+=fabs(my_document_features[document_index1]->chaincode4order2_16[x]-my_document_features[document_index2]->chaincode4order2_16[x]);
for(int x=0;x<64;x++)
chaincode4order3_64+=fabs(my_document_features[document_index1]->chaincode4order3_64[x]-my_document_features[document_index2]->chaincode4order3_64[x]);
for(int x=0;x<512;x++)
chaincode8order3_512+=fabs(my_document_features[document_index1]->chaincode8order3_512[x]-my_document_features[document_index2]->chaincode8order3_512[x]);
for(int x=0;x<256;x++)
chaincode4order4_256+=fabs(my_document_features[document_index1]->chaincode4order4_256[x]-my_document_features[document_index2]->chaincode4order4_256[x]);
for(int x=0;x<4096;x++)
chaincode8order4_4096+=fabs(my_document_features[document_index1]->chaincode8order4_4096[x]-my_document_features[document_index2]->chaincode8order4_4096[x]);
if(my_document_features[document_index1]->writer_id==my_document_features[document_index2]->writer_id)differences<<\"1,\";
else differences<<\"0,\";
differences<<NumberOfConnectedComponents_1<<\",\"<<NumberOfHoles_1<<\",\"<<SpatialMoment00<<\",\"<<SpatialMoment10<<\",\"<<SpatialMoment20<<\",\"<<SpatialMoment30<<\",\"<<SpatialMoment01<<\",\"<<SpatialMoment11<<\",\"<<SpatialMoment21<<\",\"<<SpatialMoment02<<\",\"<<SpatialMoment12<<\",\"<<SpatialMoment03<<\",\"<<CentralMoment00<<\",\"<<CentralMoment10<<\",\"<<CentralMoment20<<\",\"<<CentralMoment30<<\",\"<<CentralMoment01<<\",\"<<CentralMoment11<<\",\"<<CentralMoment21<<\",\"<<CentralMoment02<<\",\"<<CentralMoment12<<\",\"<<CentralMoment03<<\",\"<<NormalizedCentralMoment00<<\",\"<<NormalizedCentralMoment10<<\",\"<<NormalizedCentralMoment20<<\",\"<<NormalizedCentralMoment30<<\",\"<<NormalizedCentralMoment01<<\",\"<<NormalizedCentralMoment11<<\",\"<<NormalizedCentralMoment21<<\",\"<<NormalizedCentralMoment02<<\",\"<<NormalizedCentralMoment12<<\",\"<<NormalizedCentralMoment03<<\",\"<<HuMoment1<<\",\"<<HuMoment2<<\",\"<<HuMoment3<<\",\"<<HuMoment4<<\",\"<<HuMoment5<<\",\"<<HuMoment6<<\",\"<<HuMoment7<<\",\"<<XProjectionHist10<<\",\"<<YProjectionHist10<<\",\"<<XFilledProjectionHist10<<\",\"<<YFilledProjectionHist10<<\",\"<<Distribution10x10_100<<\",\"<<Barycenter2<<\",\"<<numberOfPixels1<<\",\"<<Fourier_1<<\",\"<<Fourier_5<<\",\"<<Fourier_9<<\",\"<<Fourier_15<<\",\"<<nbbranches_1<<\",\"<<LengthsOfBranchesHist_10<<\",\"<<ThicknessLengthsCircleHist30<<\",\"<<tortuosityHist10<<\",\"<<tortuosityDirectionHist10<<\",\"<<tortuosityDerivateHist10<<\",\"<<tortuosityDerivateDirectionHist10<<\",\"<<DirectionPerpendicular5Hist10<<\",\"<<CurvaturePerpendicular5Hist100<<\",\"<<luminanceHist256<<\",\"<<CurvatureAli5Hist100<<\",\"<<CurvaturesDerivateAli5Hist100<<\",\"<<CurvatureAli10Hist100<<\",\"<<CurvaturesDerivateAli10Hist100<<\",\"<<CurvatureAli15Hist100<<\",\"<<CurvaturesDerivateAli15Hist100<<\",\"<<CurvatureAli20Hist100<<\",\"<<CurvaturesDerivateAli20Hist100<<\",\"<<chaincodeHist_4<<\",\"<<chaincodeHist_8<<\",\"<<chaincode8order2_64<<\",\"<<chaincode4order2_16<<\",\"<<chaincode4order3_64<<\",\"<<chaincode8order3_512<<\",\"<<chaincode4order4_256<<\",\"<<chaincode8order4_4096<<\"\\n\";
}
differences.close();
//You might use the obtained \"train_differences.csv\" file to build your model
//I used here the model trained using a logistic regression classifier after keeping the 5 most discriminant features :
//double linear=(-10.667528249636) + -0.0667090106640684 * NumberOfConnectedComponents_1 +9.34910163750235 * ThicknessLengthsCircleHist30 + 30.4471654988069 *CurvaturePerpendicular5Hist100 + 30.5600351874986 * luminanceHist256 +33.6853635484793 * chaincode8order3_512;
//double prob_same_writer=1. / (1. + exp(linear));
//reading the testing documents
std::ifstream test_data(\"test.csv\");
linenum=0;
document_features**my_questioned_document_features=new document_features*[53];
while (getline (test_data, line))
{
cout << \"\\nLine #\" << linenum << \":\" << endl;
istringstream linestream(line);
if(linenum>0)
{
int row_index=linenum-1;
my_questioned_document_features[row_index]=new document_features();
getline (linestream, item, ',');
//my_questioned_document_features[row_index]->writer_id=atoi(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NumberOfConnectedComponents_1[0]=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NumberOfHoles_1[0]=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment00=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment10=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment20=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment30=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment01=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment11=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment21=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment02=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment12=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->SpatialMoment03=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment00=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment10=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment20=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment30=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment01=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment11=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment21=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment02=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment12=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CentralMoment03=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment00=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment10=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment20=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment30=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment01=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment11=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment21=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment02=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment12=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->NormalizedCentralMoment03=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment1=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment2=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment3=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment4=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment5=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment6=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->HuMoment7=atof(item.c_str());
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->XProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->YProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->XFilledProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->YFilledProjectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Distribution10x10_100[x]=atof(item.c_str());
}
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Barycenter2[0]=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Barycenter2[1]=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->numberOfPixels1[0]=atof(item.c_str());
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Fourier_1[0]=atof(item.c_str());
for(int x=0;x<5;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Fourier_5[x]=atof(item.c_str());
}
for(int x=0;x<9;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Fourier_9[x]=atof(item.c_str());
}
for(int x=0;x<15;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->Fourier_15[x]=atof(item.c_str());
}
getline (linestream, item, ',');
my_questioned_document_features[row_index]->nbbranches_1[0]=atof(item.c_str());
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->LengthsOfBranchesHist_10[x]=atof(item.c_str());
}
for(int x=0;x<30;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->ThicknessLengthsCircleHist30[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->tortuosityHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->tortuosityDirectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->tortuosityDerivateHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->tortuosityDerivateDirectionHist10[x]=atof(item.c_str());
}
for(int x=0;x<10;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->DirectionPerpendicular5Hist10[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvaturePerpendicular5Hist100[x]=atof(item.c_str());
}
for(int x=0;x<256;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->luminanceHist256[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvatureAli5Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvaturesDerivateAli5Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvatureAli10Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvaturesDerivateAli10Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvatureAli15Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvaturesDerivateAli15Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvatureAli20Hist100[x]=atof(item.c_str());
}
for(int x=0;x<100;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->CurvaturesDerivateAli20Hist100[x]=atof(item.c_str());
}
for(int x=0;x<4;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincodeHist_4[x]=atof(item.c_str());
}
for(int x=0;x<8;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincodeHist_8[x]=atof(item.c_str());
}
for(int x=0;x<64;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincode8order2_64[x]=atof(item.c_str());
}
for(int x=0;x<16;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincode4order2_16[x]=atof(item.c_str());
}
for(int x=0;x<64;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincode4order3_64[x]=atof(item.c_str());
}
for(int x=0;x<512;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincode8order3_512[x]=atof(item.c_str());
}
for(int x=0;x<256;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincode4order4_256[x]=atof(item.c_str());
}
for(int x=0;x<4096;x++)
{
getline (linestream, item, ',');
my_questioned_document_features[row_index]->chaincode8order4_4096[x]=atof(item.c_str());
}
}
linenum++;
}
//comparing all the documents of the test set with all the documents of the training set
for(int document_index1=0;document_index1<53;document_index1++)
{
double maxprob=0;
for(int document_index2=0;document_index2<108;document_index2++)
{
double NumberOfConnectedComponents_1=0,NumberOfHoles_1=0,SpatialMoment00=0,SpatialMoment10=0,SpatialMoment20=0,SpatialMoment30=0,SpatialMoment01=0,SpatialMoment11=0,SpatialMoment21=0,SpatialMoment02=0,SpatialMoment12=0,SpatialMoment03=0,CentralMoment00=0,CentralMoment10=0,CentralMoment20=0,CentralMoment30=0,CentralMoment01=0,CentralMoment11=0,CentralMoment21=0,CentralMoment02=0,CentralMoment12=0,CentralMoment03=0,NormalizedCentralMoment00=0,NormalizedCentralMoment10=0,NormalizedCentralMoment20=0,NormalizedCentralMoment30=0,NormalizedCentralMoment01=0,NormalizedCentralMoment11=0,NormalizedCentralMoment21=0,NormalizedCentralMoment02=0,NormalizedCentralMoment12=0,NormalizedCentralMoment03=0,HuMoment1=0,HuMoment2=0,HuMoment3=0,HuMoment4=0,HuMoment5=0,HuMoment6=0,HuMoment7=0,ThicknessLengthsCircleHist30=0,XProjectionHist10=0,YProjectionHist10=0,Distribution10x10_100=0,tortuosityHist10=0,tortuosityDirectionHist10=0,tortuosityDerivateHist10=0,tortuosityDerivateDirectionHist10=0,nbbranches_1=0,LengthsOfBranchesHist_10=0,numberOfPixels1=0,XFilledProjectionHist10=0,YFilledProjectionHist10=0,Barycenter2=0,DirectionPerpendicular5Hist10=0,CurvaturePerpendicular5Hist100=0,luminanceHist256=0,CurvatureAli5Hist100=0,CurvaturesDerivateAli5Hist100=0,CurvatureAli10Hist100=0,CurvaturesDerivateAli10Hist100=0,CurvatureAli15Hist100=0,CurvaturesDerivateAli15Hist100=0,CurvatureAli20Hist100=0,CurvaturesDerivateAli20Hist100=0,chaincodeHist_4=0,chaincodeHist_8=0,chaincode8order2_64=0,chaincode4order2_16=0,chaincode4order3_64=0,chaincode8order3_512=0,chaincode4order4_256=0,chaincode8order4_4096=0,Fourier_1=0,Fourier_5=0,Fourier_9=0,Fourier_15=0;
NumberOfConnectedComponents_1+=fabs(my_questioned_document_features[document_index1]->NumberOfConnectedComponents_1[0]-my_document_features[document_index2]->NumberOfConnectedComponents_1[0]);
NumberOfHoles_1+=fabs(my_questioned_document_features[document_index1]->NumberOfHoles_1[0]-my_document_features[document_index2]->NumberOfHoles_1[0]);
SpatialMoment00+=fabs(my_questioned_document_features[document_index1]->SpatialMoment00-my_document_features[document_index2]->SpatialMoment00);
SpatialMoment10+=fabs(my_questioned_document_features[document_index1]->SpatialMoment10-my_document_features[document_index2]->SpatialMoment10);
SpatialMoment20+=fabs(my_questioned_document_features[document_index1]->SpatialMoment20-my_document_features[document_index2]->SpatialMoment20);
SpatialMoment30+=fabs(my_questioned_document_features[document_index1]->SpatialMoment30-my_document_features[document_index2]->SpatialMoment30);
SpatialMoment01+=fabs(my_questioned_document_features[document_index1]->SpatialMoment01-my_document_features[document_index2]->SpatialMoment01);
SpatialMoment11+=fabs(my_questioned_document_features[document_index1]->SpatialMoment11-my_document_features[document_index2]->SpatialMoment11);
SpatialMoment21+=fabs(my_questioned_document_features[document_index1]->SpatialMoment21-my_document_features[document_index2]->SpatialMoment21);
SpatialMoment02+=fabs(my_questioned_document_features[document_index1]->SpatialMoment02-my_document_features[document_index2]->SpatialMoment02);
SpatialMoment12+=fabs(my_questioned_document_features[document_index1]->SpatialMoment12-my_document_features[document_index2]->SpatialMoment12);
SpatialMoment03+=fabs(my_questioned_document_features[document_index1]->SpatialMoment03-my_document_features[document_index2]->SpatialMoment03);
CentralMoment00+=fabs(my_questioned_document_features[document_index1]->CentralMoment00-my_document_features[document_index2]->CentralMoment00);
CentralMoment10+=fabs(my_questioned_document_features[document_index1]->CentralMoment10-my_document_features[document_index2]->CentralMoment10);
CentralMoment20+=fabs(my_questioned_document_features[document_index1]->CentralMoment20-my_document_features[document_index2]->CentralMoment20);
CentralMoment30+=fabs(my_questioned_document_features[document_index1]->CentralMoment30-my_document_features[document_index2]->CentralMoment30);
CentralMoment01+=fabs(my_questioned_document_features[document_index1]->CentralMoment01-my_document_features[document_index2]->CentralMoment01);
CentralMoment11+=fabs(my_questioned_document_features[document_index1]->CentralMoment11-my_document_features[document_index2]->CentralMoment11);
CentralMoment21+=fabs(my_questioned_document_features[document_index1]->CentralMoment21-my_document_features[document_index2]->CentralMoment21);
CentralMoment02+=fabs(my_questioned_document_features[document_index1]->CentralMoment02-my_document_features[document_index2]->CentralMoment02);
CentralMoment12+=fabs(my_questioned_document_features[document_index1]->CentralMoment12-my_document_features[document_index2]->CentralMoment12);
CentralMoment03+=fabs(my_questioned_document_features[document_index1]->CentralMoment03-my_document_features[document_index2]->CentralMoment03);
NormalizedCentralMoment00+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment00-my_document_features[document_index2]->NormalizedCentralMoment00);
NormalizedCentralMoment10+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment10-my_document_features[document_index2]->NormalizedCentralMoment10);
NormalizedCentralMoment20+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment20-my_document_features[document_index2]->NormalizedCentralMoment20);
NormalizedCentralMoment30+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment30-my_document_features[document_index2]->NormalizedCentralMoment30);
NormalizedCentralMoment01+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment01-my_document_features[document_index2]->NormalizedCentralMoment01);
NormalizedCentralMoment11+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment11-my_document_features[document_index2]->NormalizedCentralMoment11);
NormalizedCentralMoment21+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment21-my_document_features[document_index2]->NormalizedCentralMoment21);
NormalizedCentralMoment02+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment02-my_document_features[document_index2]->NormalizedCentralMoment02);
NormalizedCentralMoment12+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment12-my_document_features[document_index2]->NormalizedCentralMoment12);
NormalizedCentralMoment03+=fabs(my_questioned_document_features[document_index1]->NormalizedCentralMoment03-my_document_features[document_index2]->NormalizedCentralMoment03);
HuMoment1+=fabs(my_questioned_document_features[document_index1]->HuMoment1-my_document_features[document_index2]->HuMoment1);
HuMoment2+=fabs(my_questioned_document_features[document_index1]->HuMoment2-my_document_features[document_index2]->HuMoment2);
HuMoment3+=fabs(my_questioned_document_features[document_index1]->HuMoment3-my_document_features[document_index2]->HuMoment3);
HuMoment4+=fabs(my_questioned_document_features[document_index1]->HuMoment4-my_document_features[document_index2]->HuMoment4);
HuMoment5+=fabs(my_questioned_document_features[document_index1]->HuMoment5-my_document_features[document_index2]->HuMoment5);
HuMoment6+=fabs(my_questioned_document_features[document_index1]->HuMoment6-my_document_features[document_index2]->HuMoment6);
HuMoment7+=fabs(my_questioned_document_features[document_index1]->HuMoment7-my_document_features[document_index2]->HuMoment7);
for(int x=0;x<10;x++)
XProjectionHist10+=fabs(my_questioned_document_features[document_index1]->XProjectionHist10[x]-my_document_features[document_index2]->XProjectionHist10[x]);
for(int x=0;x<10;x++)
YProjectionHist10+=fabs(my_questioned_document_features[document_index1]->YProjectionHist10[x]-my_document_features[document_index2]->YProjectionHist10[x]);
for(int x=0;x<10;x++)
XFilledProjectionHist10+=fabs(my_questioned_document_features[document_index1]->XFilledProjectionHist10[x]-my_document_features[document_index2]->XFilledProjectionHist10[x]);
for(int x=0;x<10;x++)
YFilledProjectionHist10+=fabs(my_questioned_document_features[document_index1]->YFilledProjectionHist10[x]-my_document_features[document_index2]->YFilledProjectionHist10[x]);
for(int x=0;x<10;x++)
Distribution10x10_100+=fabs(my_questioned_document_features[document_index1]->Distribution10x10_100[x]-my_document_features[document_index2]->Distribution10x10_100[x]);
for(int x=0;x<2;x++)
Barycenter2+=fabs(my_questioned_document_features[document_index1]->Barycenter2[x]-my_document_features[document_index2]->Barycenter2[x]);
numberOfPixels1+=fabs(my_questioned_document_features[document_index1]->numberOfPixels1[0]-my_document_features[document_index2]->numberOfPixels1[0]);
for(int x=0;x<1;x++)
Fourier_1+=fabs(my_questioned_document_features[document_index1]->Fourier_1[x]-my_document_features[document_index2]->Fourier_1[x]);
for(int x=0;x<5;x++)
Fourier_5+=fabs(my_questioned_document_features[document_index1]->Fourier_5[x]-my_document_features[document_index2]->Fourier_5[x]);
for(int x=0;x<9;x++)
Fourier_9+=fabs(my_questioned_document_features[document_index1]->Fourier_9[x]-my_document_features[document_index2]->Fourier_9[x]);
for(int x=0;x<15;x++)
Fourier_15+=fabs(my_questioned_document_features[document_index1]->Fourier_15[x]-my_document_features[document_index2]->Fourier_15[x]);
nbbranches_1+=fabs(my_questioned_document_features[document_index1]->nbbranches_1[0]-my_document_features[document_index2]->nbbranches_1[0]);
for(int x=0;x<10;x++)
LengthsOfBranchesHist_10+=fabs(my_questioned_document_features[document_index1]->LengthsOfBranchesHist_10[x]-my_document_features[document_index2]->LengthsOfBranchesHist_10[x]);
for(int x=0;x<30;x++)
ThicknessLengthsCircleHist30+=fabs(my_questioned_document_features[document_index1]->ThicknessLengthsCircleHist30[x]-my_document_features[document_index2]->ThicknessLengthsCircleHist30[x]);
for(int x=0;x<10;x++)
tortuosityHist10+=fabs(my_questioned_document_features[document_index1]->tortuosityHist10[x]-my_document_features[document_index2]->tortuosityHist10[x]);
for(int x=0;x<10;x++)
tortuosityDirectionHist10+=fabs(my_questioned_document_features[document_index1]->tortuosityDirectionHist10[x]-my_document_features[document_index2]->tortuosityDirectionHist10[x]);
for(int x=0;x<10;x++)
tortuosityDerivateHist10+=fabs(my_questioned_document_features[document_index1]->tortuosityDerivateHist10[x]-my_document_features[document_index2]->tortuosityDerivateHist10[x]);
for(int x=0;x<10;x++)
tortuosityDerivateDirectionHist10+=fabs(my_questioned_document_features[document_index1]->tortuosityDerivateDirectionHist10[x]-my_document_features[document_index2]->tortuosityDerivateDirectionHist10[x]);
for(int x=0;x<10;x++)
DirectionPerpendicular5Hist10+=fabs(my_questioned_document_features[document_index1]->DirectionPerpendicular5Hist10[x]-my_document_features[document_index2]->DirectionPerpendicular5Hist10[x]);
for(int x=0;x<100;x++)
CurvaturePerpendicular5Hist100+=fabs(my_questioned_document_features[document_index1]->CurvaturePerpendicular5Hist100[x]-my_document_features[document_index2]->CurvaturePerpendicular5Hist100[x]);
for(int x=0;x<256;x++)
luminanceHist256+=fabs(my_questioned_document_features[document_index1]->luminanceHist256[x]-my_document_features[document_index2]->luminanceHist256[x]);
for(int x=0;x<100;x++)
CurvatureAli5Hist100+=fabs(my_questioned_document_features[document_index1]->CurvatureAli5Hist100[x]-my_document_features[document_index2]->CurvatureAli5Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli5Hist100+=fabs(my_questioned_document_features[document_index1]->CurvaturesDerivateAli5Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli5Hist100[x]);
for(int x=0;x<100;x++)
CurvatureAli10Hist100+=fabs(my_questioned_document_features[document_index1]->CurvatureAli10Hist100[x]-my_document_features[document_index2]->CurvatureAli10Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli10Hist100+=fabs(my_questioned_document_features[document_index1]->CurvaturesDerivateAli10Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli10Hist100[x]);
for(int x=0;x<100;x++)
CurvatureAli15Hist100+=fabs(my_questioned_document_features[document_index1]->CurvatureAli15Hist100[x]-my_document_features[document_index2]->CurvatureAli15Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli15Hist100+=fabs(my_questioned_document_features[document_index1]->CurvaturesDerivateAli15Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli15Hist100[x]);
for(int x=0;x<100;x++)
CurvatureAli20Hist100+=fabs(my_questioned_document_features[document_index1]->CurvatureAli20Hist100[x]-my_document_features[document_index2]->CurvatureAli20Hist100[x]);
for(int x=0;x<100;x++)
CurvaturesDerivateAli20Hist100+=fabs(my_questioned_document_features[document_index1]->CurvaturesDerivateAli20Hist100[x]-my_document_features[document_index2]->CurvaturesDerivateAli20Hist100[x]);
for(int x=0;x<4;x++)
chaincodeHist_4+=fabs(my_questioned_document_features[document_index1]->chaincodeHist_4[x]-my_document_features[document_index2]->chaincodeHist_4[x]);
for(int x=0;x<8;x++)
chaincodeHist_8+=fabs(my_questioned_document_features[document_index1]->chaincodeHist_8[x]-my_document_features[document_index2]->chaincodeHist_8[x]);
for(int x=0;x<64;x++)
chaincode8order2_64+=fabs(my_questioned_document_features[document_index1]->chaincode8order2_64[x]-my_document_features[document_index2]->chaincode8order2_64[x]);
for(int x=0;x<16;x++)
chaincode4order2_16+=fabs(my_questioned_document_features[document_index1]->chaincode4order2_16[x]-my_document_features[document_index2]->chaincode4order2_16[x]);
for(int x=0;x<64;x++)
chaincode4order3_64+=fabs(my_questioned_document_features[document_index1]->chaincode4order3_64[x]-my_document_features[document_index2]->chaincode4order3_64[x]);
for(int x=0;x<512;x++)
chaincode8order3_512+=fabs(my_questioned_document_features[document_index1]->chaincode8order3_512[x]-my_document_features[document_index2]->chaincode8order3_512[x]);
for(int x=0;x<256;x++)
chaincode4order4_256+=fabs(my_questioned_document_features[document_index1]->chaincode4order4_256[x]-my_document_features[document_index2]->chaincode4order4_256[x]);
for(int x=0;x<4096;x++)
chaincode8order4_4096+=fabs(my_questioned_document_features[document_index1]->chaincode8order4_4096[x]-my_document_features[document_index2]->chaincode8order4_4096[x]);
//using the model previously built
double linear=(-10.667528249636) + -0.0667090106640684 * NumberOfConnectedComponents_1 +9.34910163750235 * ThicknessLengthsCircleHist30 + 30.4471654988069 *CurvaturePerpendicular5Hist100 + 30.5600351874986 * luminanceHist256 +33.6853635484793 * chaincode8order3_512;
double prob_same_writer=1. / (1. + exp(linear));
//the most probable writer
if(prob_same_writer>maxprob)
{
maxprob=prob_same_writer;
my_questioned_document_features[document_index1]->writer_id=my_document_features[document_index2]->writer_id;
}
}
}
//creating the file to submit
//this benchmark does not provide a way of detecting unknown writers
std::ofstream submission_file(\"new_benchmark.csv\");
submission_file<<\",031,034,035,036,037,038,039,041,042,046,047,048,049,050,051,052,053,054,055,056,057,058,059,060,061,062,063,064,065,066,067,068,069,070,071,072,073,074,077,078,080,081,082,083,084,087,088,090,091,092,093,095,096,unknown\\n\";
int writers[]={31,34,35,36,37,38,39,41,42,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,77,78,80,81,82,83,84,87,88,90,91,92,93,95,96};
string questioned_writers[]={\"AA\",\"AB\",\"AC\",\"AD\",\"AE\",\"AF\",\"AG\",\"AH\",\"AI\",\"AJ\",\"AK\",\"AM\",\"AN\",\"AO\",\"AP\",\"AQ\",\"AR\",\"AS\",\"AT\",\"AU\",\"AV\",\"AW\",\"AX\",\"AY\",\"AZ\",\"BA\",\"BB\",\"BC\",\"BD\",\"BE\",\"BF\",\"BG\",\"BH\",\"BI\",\"BK\",\"BL\",\"BM\",\"BN\",\"BO\",\"BP\",\"BQ\",\"BR\",\"BS\",\"BT\",\"BU\",\"BV\",\"BW\",\"BX\",\"BY\",\"BZ\",\"CA\",\"CB\",\"CC\"};
for(int document_index1=0;document_index1<53;document_index1++)
{
submission_file<<questioned_writers[document_index1];
bool at_least_one=false;
for(int writer_index=0;writer_index<53;writer_index++)
{
if(my_questioned_document_features[document_index1]->writer_id==writers[writer_index])
{
at_least_one=true;
submission_file<<\",1\";
}
else submission_file<<\",0\";
}
//unknown writer
if(at_least_one)submission_file<<\",0\\n\";
else submission_file<<\",1\\n\";
}
submission_file.close();
std::cout<<\"Press a key to exit...\";
_getch();
}
",
"rawMarkdown": "",
"votes": null
},
{
"id": "2024",
"postDate": "03/13/2011 17:45:09",
"content": "Hello,
When I submit the output of the provided benchmark program, I get an MAE of 0.00931315 whereas the Qatar University Benchmark is 0.00814901. Any reason why they are different? Also the benchmark source code provided below has a lot of errors due to spaces inserted in the text like this : Distributi on10x10_100, so for convenience I attached the corrected file.
Regards,
Sameh
",
"rawMarkdown": "",
"votes": null
},
{
"id": "2025",
"postDate": "03/13/2011 17:45:09",
"content": "Thanks for your effort.
I have posted two benchmarks and I have chosen the entry which has the lowest MAE.
In the other benchmark, there is a slight difference in the way of computing the sums of absolute differences.
PS: I had to create another Kaggle account because I couldn't access to this topic with my other account!
",
"rawMarkdown": "",
"votes": null
}
],
"comments": [
{
"id": 2024,
"author_name": "drsameh",
"author_url": "",
"post_date": "03/13/2011 17:45:09",
"content": "Hello,
When I submit the output of the provided benchmark program, I get an MAE of 0.00931315 whereas the Qatar University Benchmark is 0.00814901. Any reason why they are different? Also the benchmark source code provided below has a lot of errors due to spaces inserted in the text like this : Distributi on10x10_100, so for convenience I attached the corrected file.
Regards,
Sameh
",
"votes": null,
"replies": []
},
{
"id": 2025,
"author_name": "hassaine",
"author_url": "",
"post_date": "03/13/2011 17:45:09",
"content": "Thanks for your effort.
I have posted two benchmarks and I have chosen the entry which has the lowest MAE.
In the other benchmark, there is a slight difference in the way of computing the sums of absolute differences.
PS: I had to create another Kaggle account because I couldn't access to this topic with my other account!
",
"votes": null,
"replies": []
}
],
"raw_markdown_by_id": {
"1993": "",
"2024": "",
"2025": ""
},
"source": "meta"
}