人脸识别经典算法一 特征脸方法(Eigenface)


Posted in Python onMarch 13, 2018

这篇文章是撸主要介绍人脸识别经典方法的第一篇,后续会有其他方法更新。特征脸方法基本是将人脸识别推向真正可用的第一种方法,了解一下还是很有必要的。特征脸用到的理论基础PCA在另一篇博客里:特征脸(Eigenface)理论基础-PCA(主成分分析法) 。本文的参考资料附在最后了^_^

步骤一:获取包含M张人脸图像的集合S。在我们的例子里有25张人脸图像(虽然是25个不同人的人脸的图像,但是看着怎么不像呢,难道我有脸盲症么),如下图所示哦。每张图像可以转换成一个N维的向量(是的,没错,一个像素一个像素的排成一行就好了,至于是横着还是竖着获取原图像的像素,随你自己,只要前后统一就可以),然后把这M个向量放到一个集合S里,如下式所示。

人脸识别经典算法一 特征脸方法(Eigenface)

人脸识别经典算法一 特征脸方法(Eigenface)

步骤二:在获取到人脸向量集合S后,计算得到平均图像Ψ,至于怎么计算平均图像,公式在下面。就是把集合S里面的向量遍历一遍进行累加,然后取平均值。得到的这个Ψ其实还挺有意思的,Ψ其实也是一个N维向量,如果再把它还原回图像的形式的话,可以得到如下的“平均脸”,是的没错,还他妈的挺帅啊。那如果你想看一下某计算机学院男生平均下来都长得什么样子,用上面的方法就可以了。

人脸识别经典算法一 特征脸方法(Eigenface)

人脸识别经典算法一 特征脸方法(Eigenface)

步骤三:计算每张图像和平均图像的差值Φ,就是用S集合里的每个元素减去步骤二中的平均值。

人脸识别经典算法一 特征脸方法(Eigenface)

步骤四:找到M个正交的单位向量un,这些单位向量其实是用来描述Φ(步骤三中的差值)分布的。un里面的第k(k=1,2,3...M)个向量uk是通过下式计算的,

人脸识别经典算法一 特征脸方法(Eigenface)

当这个λk(原文里取了个名字叫特征值)取最小的值时,uk基本就确定了。补充一下,刚才也说了,这M个向量是相互正交而且是单位长度的,所以啦,uk还要满足下式:

人脸识别经典算法一 特征脸方法(Eigenface)

上面的等式使得uk为单位正交向量。计算上面的uk其实就是计算如下协方差矩阵的特征向量:

人脸识别经典算法一 特征脸方法(Eigenface)

其中

人脸识别经典算法一 特征脸方法(Eigenface)

对于一个NxN(比如100x100)维的图像来说,上述直接计算其特征向量计算量实在是太大了(协方差矩阵可以达到10000x10000),所以有了如下的简单计算。

步骤四另解:如果训练图像的数量小于图像的维数比如(M<N^2),那么起作用的特征向量只有M-1个而不是N^2个(因为其他的特征向量对应的特征值为0),所以求解特征向量我们只需要求解一个NxN的矩阵。这个矩阵就是步骤四中的AAT,我们可以设该矩阵为L,那么L的第m行n列的元素可以表示为:

人脸识别经典算法一 特征脸方法(Eigenface)

一旦我们找到了L矩阵的M个特征向量vl,那么协方差矩阵的特征向量ul就可以表示为:

人脸识别经典算法一 特征脸方法(Eigenface)

这些特征向量如果还原成像素排列的话,其实还蛮像人脸的,所以称之为特征脸(如下图)。图里有二十五个特征脸,数量上和训练图像相等只是巧合。有论文表明一般的应用40个特征脸已经足够了。论文Eigenface for recognition里只用了7个特征脸来表明实验。

人脸识别经典算法一 特征脸方法(Eigenface)

步骤五:识别人脸。OK,终于到这步了,别绕晕啦,上面几步是为了对人脸进行降维找到表征人脸的合适向量的。首先考虑一张新的人脸,我们可以用特征脸对其进行标示:

人脸识别经典算法一 特征脸方法(Eigenface)

其中k=1,2...M,对于第k个特征脸uk,上式可以计算其对应的权重,M个权重可以构成一个向量:

人脸识别经典算法一 特征脸方法(Eigenface)

perfect,这就是求得的特征脸对人脸的表示了!

那如何对人脸进行识别呢,看下式:

人脸识别经典算法一 特征脸方法(Eigenface)

其中Ω代表要判别的人脸,Ωk代表训练集内的某个人脸,两者都是通过特征脸的权重来表示的。式子是对两者求欧式距离,当距离小于阈值时说明要判别的脸和训练集内的第k个脸是同一个人的。当遍历所有训练集都大于阈值时,根据距离值的大小又可分为是新的人脸或者不是人脸的两种情况。根据训练集的不同,阈值设定并不是固定的。

后续会有对PCA理论的补充^_^.已补充理论:特征脸(Eigenface)理论基础-PCA(主成分分析法)

参考资料:

1、Eigenface for Recognition

2、特征脸维基百科

3、Eigenface_tutorial

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持三水点靠木。

Python 相关文章推荐
用Python制作在地图上模拟瘟疫扩散的Gif图
Mar 31 Python
Python socket编程实例详解
May 27 Python
详解python之多进程和进程池(Processing库)
Jun 09 Python
python微信跳一跳系列之色块轮廓定位棋盘
Feb 26 Python
利用python为运维人员写一个监控脚本
Mar 25 Python
使用pandas批量处理矢量化字符串的实例讲解
Jul 10 Python
python编辑用户登入界面的实现代码
Jul 16 Python
python 实现语音聊天机器人的示例代码
Dec 02 Python
基于Django的乐观锁与悲观锁解决订单并发问题详解
Jul 31 Python
Python处理session的方法整理
Aug 29 Python
Python类反射机制使用实例解析
Dec 30 Python
Python快速优雅的批量修改Word文档样式
May 20 Python
Python cookbook(数据结构与算法)从序列中移除重复项且保持元素间顺序不变的方法
Mar 13 #Python
Python+Django搭建自己的blog网站
Mar 13 #Python
python实现人脸识别经典算法(一) 特征脸法
Mar 13 #Python
tensorflow实现softma识别MNIST
Mar 12 #Python
wxpython实现图书管理系统
Mar 12 #Python
人生苦短我用python python如何快速入门?
Mar 12 #Python
tensorflow实现KNN识别MNIST
Mar 12 #Python
You might like
PHP中Memcache操作类及用法实例
2014/12/12 PHP
PHP使用xmllint命令处理xml与html的方法
2014/12/15 PHP
PHP获取photoshop写入图片文字信息的方法
2015/03/31 PHP
PHP中的随机性 你觉得自己幸运吗?
2016/01/22 PHP
Yii2框架使用计划任务的方法
2016/05/25 PHP
PHP序列化操作方法分析
2016/09/28 PHP
thinkPHP框架中执行事务的方法示例
2018/05/31 PHP
PHP删除字符串中非字母数字字符方法总结
2019/01/20 PHP
Laravel框架源码解析之入口文件原理分析
2020/05/14 PHP
将jQuery应用于login页面的问题及解决
2009/10/17 Javascript
juqery 学习之三 选择器 层级 基本
2010/11/25 Javascript
一个级联菜单代码学习及removeClass与addClass的应用
2013/01/24 Javascript
JavaScript中字符串分割函数split用法实例
2015/04/07 Javascript
JavaScript实现级联菜单的方法
2015/06/29 Javascript
javascript执行环境及作用域详解
2016/05/05 Javascript
基于nodejs 的多页面爬虫实例代码
2017/05/31 NodeJs
vue生成token并保存到本地存储中
2018/07/17 Javascript
nodejs使用node-xlsx生成excel的方法示例
2019/08/22 NodeJs
Vue的click事件防抖和节流处理详解
2019/11/13 Javascript
vue模块移动组件的实现示例
2020/05/20 Javascript
解决vue组件没显示,没起作用,没报错,但该显示的组件没显示问题
2020/09/02 Javascript
[50:45]2018DOTA2亚洲邀请赛 4.6 淘汰赛 VP vs TNC 第一场
2018/04/10 DOTA
Python使用sftp实现上传和下载功能(实例代码)
2017/03/14 Python
python编程实现希尔排序
2017/04/13 Python
使用Python进行AES加密和解密的示例代码
2018/02/02 Python
python字符串循环左移
2019/03/08 Python
flask框架json数据的拿取和返回操作示例
2019/11/28 Python
Python txt文件常用读写操作代码实例
2020/08/03 Python
土耳其时尚潮流在线购物网站:Trendyol
2017/10/10 全球购物
五年级数学教学反思
2014/02/11 职场文书
农村婚礼主持词
2014/03/13 职场文书
《春到梅花山》教学反思
2014/04/16 职场文书
会计系毕业求职信
2014/08/07 职场文书
详解Vue的sync修饰符
2021/05/15 Vue.js
MySQL 全文索引使用指南
2021/05/25 MySQL
使用nginx配置访问wgcloud的方法
2021/06/26 Servers