OpenCV3.3中主成分分析(Principal Components Analysis, PCA)接口简介及使用
OpenCV3.3中给出了主成分分析(Principal Components Analysis, PCA)的实现，即cv::PCA类，类的声明在include/opencv2/core.hpp文件中，实现在modules/core/src/pca.cpp文件中,其中：
(1)、cv::PCA::PCA：构造函数；
(2)、cv::PCA::operator()：函数调用运算符；
(3)、cv::PCA::project：将输入数据投影到PCA主成分空间；
(4)、cv::PCA::backProject：重建原始数据；
(5)、cv::PCA::write：将特征值、特征向量、均值写入指定的文件；
(6)、cv::PCA::read：从指定文件读入特征值、特征向量、均值；
(7)、cv::PCA::eigenvectors：协方差矩阵的特征向量；
(8)、cv::PCA::eigenvalues：协方差矩阵的特征值；
(9)、cv::PCA::mean：均值。
关于PCA的介绍可以参考： http://blog.csdn.net/fengbingchun/article/details/78977202
以下是使用ORL Faces Database作为测试图像。关于ORL Faces Database的介绍可以参考： http://blog.csdn.net/fengbingchun/article/details/79008891
测试代码如下：
#include "opencv.hpp"
#include <string>
#include <vector>
#include <memory>
#include <algorithm>
#include <opencv2/opencv.hpp>
#include <opencv2/ml.hpp>
#include "common.hpp"// PCA(Principal Component Analysis) ///
int test_opencv_pca()
{// reference: opencv-3.3.0/samples/cpp/pca.cppconst std::string image_path{ "E:/GitCode/NN_Test/data/database/ORL_Faces/" };const std::string image_name{ "1.pgm" };std::vector<cv::Mat> images;for (int i = 1; i <= 15; ++i) {std::string name = image_path + "s" + std::to_string(i) + "/" + image_name;cv::Mat mat = cv::imread(name, 0);if (!mat.data) {fprintf(stderr, "read image fail: %s\n", name.c_str());return -1;}images.emplace_back(mat);}cv::Mat data(images.size(), images[0].rows * images[0].cols, CV_32FC1);for (int i = 0; i < images.size(); ++i) {cv::Mat image_row = images[i].clone().reshape(1, 1);cv::Mat row_i = data.row(i);image_row.convertTo(row_i, CV_32F);}cv::PCA pca(data, cv::Mat(), cv::PCA::DATA_AS_ROW, 0.95f);std::vector<cv::Mat> result(images.size());for (int i = 0; i < images.size(); ++i) {// Demonstration of the effect of retainedVariance on the first imagecv::Mat point = pca.project(data.row(i)); // project into the eigenspace, thus the image becomes a "point"cv::Mat reconstruction = pca.backProject(point); // re-create the image from the "point"reconstruction = reconstruction.reshape(images[i].channels(), images[i].rows); // reshape from a row vector into image shapecv::normalize(reconstruction, reconstruction, 0, 255, cv::NORM_MINMAX, CV_8UC1);reconstruction.copyTo(result[i]);}save_images(result, "E:/GitCode/NN_Test/data/pca_result_.jpg", 5);// save fileconst std::string save_file{ "E:/GitCode/NN_Test/data/pca.xml" }; // .xml, .yaml, .jsonscv::FileStorage fs(save_file, cv::FileStorage::WRITE);pca.write(fs);fs.release();// read fileconst std::string& read_file = save_file;cv::FileStorage fs2(read_file, cv::FileStorage::READ);cv::PCA pca2;pca2.read(fs2.root());fs2.release();return 0;
}结果图像如下：
GitHub： https://github.com/fengbingchun/NN_Test
相关文章：
Spring MVC配置
为什么80%的码农都做不了架构师？>>> 一、传统方式配置Spring MVC （1）导入jar包 需要导入如下的jar包 junit-3.8.1.jar spring-core-3.0.5.RELEASE.jar commons-logging-1.1.1.jar spring-context-3.0.5.REL…
主成分分析(PCA)Python代码实现
主成分分析(Principal Components Analysis, PCA)简介可以参考： http://blog.csdn.net/fengbingchun/article/details/78977202 这里是参照 http://sebastianraschka.com/Articles/2014_pca_step_by_step.html 文章中的code整理的Python代码，实现过程为…
AI发展这一年：不断衍生的技术丑闻与抵制声潮
作者 | AI Now学院译者 | Raku编辑 | Jane出品 | AI科技大本营（ID: rgznai100）【导读】10月2日，纽约大学AI Now学院在纽约大学斯克博剧院（Skirball Theatre）组织召开了第四届年度AI Now研讨会。研讨会邀请了业内组织者…
Distributed Configuration Management Platform(分布式配置管理平台)
2019独角兽企业重金招聘Python工程师标准>>> 专注于各种 分布式系统配置管理 的通用组件/通用平台, 提供统一的配置管理服务。 主要目标： 部署极其简单：同一个上线包，无须改动配置，即可在 多个环境中(RD/QA/PRODUCTION…
如何利用zendstudio新建 或导入php项目
为什么80%的码农都做不了架构师？>>> 一、利用ZendStudio创建 PHP Project 1. 打开ZendStudio, 选择：File à New à PHP Project, 如下图所示: 于是弹出如下界面： 在”Project name”后输入工程名(比如我这里…
一文读懂GoogLeNet神经网络 | CSDN博文精选
作者 | .NY&XX来源 | CSDN博客本文介绍的是著名的网络结构GoogLeNet，目的是试图领会其中结构设计思想。GoogLeNet特点优化网络质量的生物学原理GoogLeNet网络结构的动机GoogLeNet架构细节Inception模块和普通卷积结构的差异辅助分类器GoogLeNet网络架构GoogLeNe…
C++中的函数签名
C中的函数签名(function signature)：包含了一个函数的信息，包括函数名、参数类型、参数个数、顺序以及它所在的类和命名空间。普通函数签名并不包含函数返回值部分，如果两个函数仅仅只有函数返回值不同，那么系统是无法区分这两个函…
MyEclipse断点调试
2019独角兽企业重金招聘Python工程师标准>>> 1、在编辑的程序的左边，你会看到一条浅浅的灰色编带，在这里设置断点。 2、设置断点的方法有很多 方法：1）、双击 ； 2）、右键，选择“Toggl…
C primer plus 练习题 第三章
5. 1 #include <stdio.h>2 3 int main()4 {5 float you_sec;6 printf("请输入你的年龄:");7 scanf("%f", &you_sec);8 printf("年龄合计:%e 秒!\n", you_sec * 3.156e7);9 getchar(); 10 return 0; 11 }
Echache整合Spring缓存实例讲解
2019独角兽企业重金招聘Python工程师标准>>> 摘要：本文主要介绍了EhCache，并通过整合Spring给出了一个使用实例。 一、EhCache 介绍 EhCache 是一个纯Java的进程内缓存框架，具有快速、精干等特点，是Hibernate中默认的C…
灰度图像直方图均衡化公式及实现
图像的直方图：直方图是图像中像素强度分布的图形表达方式。它统计了每一个强度值所具有的像素个数。直方图均衡化：是通过拉伸像素强度分布范围来增强图像对比度的一种方法。是图像处理领域中利用图像直方图对对比度进行调整的方法。均衡化指的是把一个分…
超越最新无监督域自适应方法，研究人员提轻量CNN新架构OSNet
作者 | Kaiyang Zhou, Xiatian Zhu, Yongxin Yang, Andrea Cavallaro, and Tao Xiang 译者 | TroyChang 编辑 | Jane 出品 | AI科技大本营（ID：rgznai100） CNN新架构OSNet 【导读】今天推荐论文《Learning Generalisable Omni-Scale Repre…
在一台机器上搭建多个redis实例
2019独角兽企业重金招聘Python工程师标准>>> 默认Redis程序安装在/usr/local/redis目录下； 配置文件：/usr/local/redis/redis.conf，该配置文件中配置的端口为默认端口：6379； Redis的启动命令路径࿱…
使用kaptcha生成验证码
2019独角兽企业重金招聘Python工程师标准>>> kaptcha是一个简单好用的验证码生成工具，通过配置，可以自己定义验证码大小、颜色、显示的字符等等。下面就来讲一下如何使用kaptcha生成验证码以及在服务器端取出验证码进行校验。 一、搭建测试环…
主成分分析(PCA) C++ 实现
主成分分析(Principal Components Analysis, PCA)简介可以参考： http://blog.csdn.net/fengbingchun/article/details/78977202以下是PCA的C实现，参考OpenCV 3.3中的cv::PCA类。使用ORL Faces Database作为测试图像。关于ORL Faces Database的介绍可以参…
为何Google、微软、华为将亿级源代码放一个仓库？从全球最大代码管理库说起...
作者 | 夕颜编辑 | Just出品 | AI 科技大本营（ID:rgznai100）【导读】2017 年，在当时微软的一篇官方博客中，时任微软云开发服务副总裁的 Brian Harry 表示微软内部代码开始向 Git 迁移，宣布推出针对大规模 repo 的“Git…
jquery mobie导致超链接不可用
在a标签中添加rel"external"即可转载于:https://blog.51cto.com/here2142/1435434
编译器GCC与Clang的异同
GCC：GNU(Gnus Not Unix)编译器套装(GNU Compiler Collection，GCC)，指一套编程语言编译器，以GPL及LGPL许可证所发行的自由软件，也是GNU项目的关键部分，也是GNU工具链的主要组成部分之一。GCC(特别是其中的C语…
如何正确选择聚类算法？ | CSDN博文精选
作者 | Josh Thompson翻译 | 张睿毅校对 | 王雨桐来源 | 数据派THU（ID:DatapiTHU）本文将介绍四种基本的聚类算法—层次聚类、基于质心的聚类、最大期望算法和基于密度的聚类算法，并讨论不同算法的优缺点。聚类算法十分容易上手，但…
Python工具 | 9个用来爬取网络站点的 Python 库
1️⃣Scrapy 一个开源和协作框架，用于从网站中提取所需的数据。 以快速，简单，可扩展的方式。 官网2️⃣cola 一个分布式爬虫框架。 GitHub3️⃣Demiurge 基于 PyQuery 的爬虫微型框架。 官网4️⃣feedparser 通用 feed 解析器。 官网5️⃣Gra…
Python并发编程实例教程
有关Python中的并发编程实例,主要是对Threading模块的应用,文中自定义了一个Threading类库。 一、简介 我们将一个正在运行的程序称为进程。每个进程都有它自己的系统状态，包含内存状态、打开文件列表、追踪指令执行情况的程序指针以及一个保存局部变量的调用栈。…
Nginx负载均衡之TCP端口高可用(二)
在前面我们实现了基本的HTTP反向代理，从互联网过来的请求已经可以分发到后端多台网站服务器上，但不是所有的业务都是网络类型的，此篇文章我们主要讨论的是TCP 端口的负载均衡做法，昨天也有小伙伴提到了，在HTTP反向代理…
语音识别大牛Daniel Povey为何加入小米？“手机+AIoT”强大生态，开源战略是关键...
整理 | 夕颜出品 | AI科技大本营（ID:rgznai100）10 月 17 日，语音识别开源工具 Kaldi 创始人，语音和 AI 领域大牛 Daniel Povey 在10 月 19 日，小米集团副总裁、集团技术委员会主席崔宝秋发布微博，欢迎 Dani…
Python运维项目中用到的redis经验及数据类型
先感叹下，学东西一定要活学活用！ 我用redis也有几年的历史了，今个才想到把集合可以当python list用。 最近做了几个项目都掺杂了redis， 遇到了一些个问题和开发中提高性能的方法，这都分享出来，共同学习。…
图像边缘检测之拉普拉斯(Laplacian)C++实现
拉普拉斯算子(Laplacian)可应用到图像边缘检测中。在OpenCV中当kernel大小为3*3时，支持两种kernel算子，分别为：在OpenCV中默认的计算方式如下，假设有一个5*5的小图像，原始值依次为1,2,…25,如下图红色部分，…
从0到1，Airbnb的深度学习实践经验总结
作者 | Haldar译者 | 陆离出品 | AI科技大本营（ID: rgznai100）此前，AI科技大本营发布了关于希望通过介绍的研究成果为读者提供一些有用的帮助和指引。模型中的生态系统本文要讨论的机器学习现实应用，是关于根据用户预约的可能性来…
高并发大流量专题---8、动态语言的并发处理
高并发大流量专题---8、动态语言的并发处理 一、总结 一句话总结： 和本科毕业论文连起来了：基于消息中间件Rocket MQ的研究；用于并发处理的消息队列 1、什么是进程、线程、协程？ 进程（Process）是计算机中的…
1. 文件系统——磁盘分区、各目录功能、硬盘
一、磁盘分区及文件访问入口在前文中介绍过，Linux的整个文件系统像一棵倒置的数，最顶层的是根文件系统，其下有很多一级子目录，一级子目录下面是二级子目录，依此类推：/：根目录/bin，/s…
吴恩达老师深度学习视频课笔记：逻辑回归公式推导及C++实现
逻辑回归(Logistic Regression)是一个二分分类算法。逻辑回归的目标是最小化其预测与训练数据之间的误差。为了训练逻辑回归模型中的参数w和b，需要定义一个成本函数(cost function)。成本函数(cost function)：它是针对整个训练集的。衡量参数w和b在整个训…
网络运行时间提高100倍，Google使用的AI视频理解架构有多强？
译者 | 刘畅出品 | AI科技大本营（ID:rgznai100）视频理解是一个很有挑战性的问题。由于视频包含时空数据，因此图像的特征表示需要同时提取图像和运动信息。这不仅对自动理解视频语义内容有重要性，还对机器人的感知和学习也至关重要…