PCA提取特征脸（vs2013+opencv249）

最新推荐文章于 2022-05-11 10:06:50 发布

zijingping

最新推荐文章于 2022-05-11 10:06:50 发布

阅读量1.1k

点赞数

分类专栏： opencv基础知识文章标签： PCA eigenface

本文链接：https://blog.csdn.net/bettyshasha/article/details/48957791

版权

opencv基础知识专栏收录该内容

4 篇文章 0 订阅

订阅专栏

#include "opencv\cv.h"
#include "opencv\highgui.h"

using namespace std;
using namespace cv;

Mat normalize(const Mat& src) {
 Mat srcnorm;
 normalize(src, srcnorm, 0, 255, NORM_MINMAX, CV_8UC1);
 return srcnorm;
}

int main(int argc, char *argv[]) {
 vector<Mat> db;


 string prefix = "E:/CSolution/2013test/test/att_faces/";db.push_back(imread(prefix + "s1/1.pgm", 0));
 //db.push_back(imread(prefix + "s1\\1.pgm", 0));
 db.push_back(imread(prefix + "s1/2.pgm", 0));
 db.push_back(imread(prefix + "s1/3.pgm", 0));

 db.push_back(imread(prefix + "s2/1.pgm", 0));
 db.push_back(imread(prefix + "s2/2.pgm", 0));
 db.push_back(imread(prefix + "s2/3.pgm", 0));

 db.push_back(imread(prefix + "s3/1.pgm", 0));
 db.push_back(imread(prefix + "s3/2.pgm", 0));
 db.push_back(imread(prefix + "s3/3.pgm", 0));

 db.push_back(imread(prefix + "s4/1.pgm", 0));
 db.push_back(imread(prefix + "s4/2.pgm", 0));
 db.push_back(imread(prefix + "s4/3.pgm", 0));
 double alpha = 1;
 int nprincipal = 10;
 double beta = 0;
 size_t n = db.size();
 size_t d = db[0].total(); 

Mat src(n, d, CV_32FC1);
 for (int i = 0; i < n; i++)
 {
  Mat xi = src.row(i);
  //db[i].reshape(1, 1).row(0).convertTo(xi, CV_32FC1, alpha, beta);
  db[i].reshape(1, 1).row(0).convertTo(xi, CV_32FC1, alpha, beta);
 }
 PCA pca(src, Mat(), CV_PCA_DATA_AS_ROW, nprincipal);
 imshow("pca", normalize(pca.eigenvectors.row(0).reshape(1, db[0].rows)));
 waitKey();

 return 0;
}