fisherface 的c++源码

SOURCE CODE FISHERFACE:

void trainFISHER(){
numeigen=numtrain-numclass; //number of eigenfaces to use
for PCA
pfile1=trainPCA(); //perform PCA on the training
images, hte projected images are stored in the array weight[]

CvMat *meanclass[200]; // mean of each class
CvMat *mean; // mean of all classes
CvMat *Sb; // between-class scatter
matrix
CvMat *Sw; // within-class scatter matrix
CvMat *Sw1; // inverse within-class
scatter matrix

//** find mean of each class
for(i=0;i<numclass;i++)
meanclass[i]=cvCloneMat(weightclass[i]);
//weightclass[i] is average of the weights of class i

//** find mean of all classes
mean=cvCreateMat(1,numeigen,CV_32F);
mean=cvCloneMat(weight[0]);
for(i=1;i<numtrain;i++)
cvAdd(mean,weight[i],mean);
cvConvertScale(mean,mean,1.0/numtrain,0);

//** find within class scatter matrix
Sw = cvCreateMat(numeigen,numeigen,CV_32F);