python实现识别相似图片小结


Posted in Python onFebruary 22, 2016

文章简介

在网上看到python做图像识别的相关文章后,真心感觉python的功能实在太强大,因此将这些文章总结一下,建立一下自己的知识体系。
当然了,图像识别这个话题作为计算机科学的一个分支,不可能就在本文简单几句就说清,所以本文只作基本算法的科普向。

如有错误,请多包涵和多多指教。

参考的文章和图片来源会在底部一一列出。

以及本篇文章所用的代码都会在底下给出github地址。

安装相关库

python用作图像处理的相关库主要有openCV(C++编写,提供了python语言的接口),PIL,但由于PIL很早就停了,所以不支持python3.x,所以建议使用基于PIL的pillow,本文也是在python3.4和pillow的环境下进行实验。

pillow下载地址
PIL的下载地址
openCV的官网

至于opencv,在做人脸识别的时候会用到,但本文不会涉及到,在本专栏的后续中会谈及openCV的人脸识别和基于此的python图片爬虫,有兴趣的朋友可以关注本专栏。

相关背景

要识别两张相似图像,我们从感性上来谈是怎么样的一个过程?首先我们会区分这两张相片的类型,例如是风景照,还是人物照。风景照中,是沙漠还是海洋,人物照中,两个人是不是都是国字脸,还是瓜子脸(还是倒瓜子脸……哈哈……)。

那么从机器的角度来说也是这样的,先识别图像的特征,然后再相比。

很显然,在没有经过训练的计算机(即建立模型),那么计算机很难区分什么是海洋,什么是沙漠。但是计算机很容易识别到图像的像素值。

因此,在图像识别中,颜色特征是最为常用的。(其余常用的特征还有纹理特征、形状特征和空间关系特征等)

其中又分为

  1. 直方图
  2. 颜色集
  3. 颜色矩
  4. 聚合向量
  5. 相关图

直方图计算法

这里先用直方图进行简单讲述。

先借用一下恋花蝶的图片,

python实现识别相似图片小结
python实现识别相似图片小结

从肉眼来看,这两张图片大概也有八成是相似的了。

在python中可以依靠Image对象的histogram()方法获取其直方图数据,但这个方法返回的结果是一个列表,如果想得到下图可视化数据,需要另外使用 matplotlib,这里因为主要介绍算法思路,matplotlib的使用这里不做介绍。

python实现识别相似图片小结

是的,我们可以明显的发现,两张图片的直方图是近似重合的。所以利用直方图判断两张图片的是否相似的方法就是,计算其直方图的重合程度即可。

计算方法如下:

python实现识别相似图片小结

其中gi和si是分别指两条曲线的第i个点。

最后计算得出的结果就是就是其相似程度。

不过,这种方法有一个明显的弱点,就是他是按照颜色的全局分布来看的,无法描述颜色的局部分布和色彩所处的位置。

也就是假如一张图片以蓝色为主,内容是一片蓝天,而另外一张图片也是蓝色为主,但是内容却是妹子穿了蓝色裙子,那么这个算法也很可能认为这两张图片的相似的。

缓解这个弱点有一个方法就是利用Image的crop方法把图片等分,然后再分别计算其相似度,最后综合考虑。

图像指纹与汉明距离

在介绍下面其他判别相似度的方法前,先补充一些概念。第一个就是图像指纹

图像指纹和人的指纹一样,是身份的象征,而图像指纹简单点来讲,就是将图像按照一定的哈希算法,经过运算后得出的一组二进制数字。

说到这里,就可以顺带引出汉明距离的概念了。

假如一组二进制数据为101,另外一组为111,那么显然把第一组的第二位数据0改成1就可以变成第二组数据111,所以两组数据的汉明距离就为1

简单点说,汉明距离就是一组二进制数据变成另一组数据所需的步骤数,显然,这个数值可以衡量两张图片的差异,汉明距离越小,则代表相似度越高。汉明距离为0,即代表两张图片完全一样。

如何计算得到汉明距离,情况下面三种哈希算法

平均哈希法(aHash)

此算法是基于比较灰度图每个像素与平均值来实现的

一般步骤

1.缩放图片,可利用Image对象的resize(size)改变,一般大小为8*8,64个像素值。
2.转化为灰度图
转灰度图的算法。
1.浮点算法:Gray=Rx0.3+Gx0.59+Bx0.11
2.整数方法:Gray=(Rx30+Gx59+Bx11)/100
3.移位方法:Gray =(Rx76+Gx151+Bx28)>>8;
4.平均值法:Gray=(R+G+B)/3;
5.仅取绿色:Gray=G;

在python中,可用Image的对象的方法convert('L')直接转换为灰度图

3.计算平均值:计算进行灰度处理后图片的所有像素点的平均值。
4.比较像素灰度值:遍历灰度图片每一个像素,如果大于平均值记录为1,否则为0.
5.得到信息指纹:组合64个bit位,顺序随意保持一致性。
最后比对两张图片的指纹,获得汉明距离即可。

感知哈希算法(pHash)

平均哈希算法过于严格,不够精确,更适合搜索缩略图,为了获得更精确的结果可以选择感知哈希算法,它采用的是DCT(离散余弦变换)来降低频率的方法

一般步骤:

  1. 缩小图片:32 * 32是一个较好的大小,这样方便DCT计算
  2. 转化为灰度图:把缩放后的图片转化为256阶的灰度图。(具体算法见平均哈希算法步骤)
  3. 计算DCT:DCT把图片分离成分率的集合
  4. 缩小DCT:DCT是32 * 32,保留左上角的8 * 8,这些代表的图片的最低频率
  5. 计算平均值:计算缩小DCT后的所有像素点的平均值。
  6. 进一步减小DCT:大于平均值记录为1,反之记录为0.
  7. 得到信息指纹:组合64个信息位,顺序随意保持一致性。

最后比对两张图片的指纹,获得汉明距离即可。

这里给出别人的DCT的介绍和计算方法(离散余弦变换的方法)

相比pHash,dHash的速度要快的多,相比aHash,dHash在效率几乎相同的情况下的效果要更好,它是基于渐变实现的。

步骤:

  1. 缩小图片:收缩到9*8的大小,一遍它有72的像素点
  2. 转化为灰度图:把缩放后的图片转化为256阶的灰度图。(具体算法见平均哈希算法步骤)
  3. 计算差异值:dHash算法工作在相邻像素之间,这样每行9个像素之间产生了8个不同的差异,一共8行,则产生了64个差异值
  4. 获得指纹:如果左边的像素比右边的更亮,则记录为1,否则为0.

最后比对两张图片的指纹,获得汉明距离即可。

总结

这几种算法是识别相似图像的基础,显然,有时两图中的人相似比整体的颜色相似更重要,所以我们有时需要进行人脸识别,
然后在脸部区进行局部哈希,或者进行其他的预处理再进行哈希,这里涉及其他知识本文不作介绍。

下一次将讲述利用opencv和以训练好的模型来进行人脸识别。

本文算法的实现在下面,点一下下面的连接就好

github仓库

Python 相关文章推荐
Python的字典和列表的使用中一些需要注意的地方
Apr 24 Python
KMP算法精解及其Python版的代码示例
Jun 01 Python
python中map()函数的使用方法示例
Sep 29 Python
在Python程序员面试中被问的最多的10道题
Dec 05 Python
Python使用zip合并相邻列表项的方法示例
Mar 17 Python
python3爬虫怎样构建请求header
Dec 23 Python
对python周期性定时器的示例详解
Feb 19 Python
Python实现钉钉发送报警消息的方法
Feb 20 Python
python sklearn库实现简单逻辑回归的实例代码
Jul 01 Python
python、PyTorch图像读取与numpy转换实例
Jan 13 Python
python GUI库图形界面开发之PyQt5选项卡控件QTabWidget详细使用方法与实例
Mar 01 Python
Python装饰器详细介绍
Mar 25 Python
python脚本设置系统时间的两种方法
Feb 21 #Python
python脚本设置超时机制系统时间的方法
Feb 21 #Python
Go语言基于Socket编写服务器端与客户端通信的实例
Feb 19 #Python
使用C#配合ArcGIS Engine进行地理信息系统开发
Feb 19 #Python
Python中使用OpenCV库来进行简单的气象学遥感影像计算
Feb 19 #Python
Python实现以时间换空间的缓存替换算法
Feb 19 #Python
Python使用爬虫猜密码
Feb 19 #Python
You might like
第六节--访问属性和方法
2006/11/16 PHP
实用函数7
2007/11/08 PHP
php常用ODBC函数集(详细)
2013/06/24 PHP
php基于CodeIgniter实现图片上传、剪切功能
2016/05/14 PHP
深入浅出讲解:php的socket通信原理
2016/12/03 PHP
Mootools 1.2教程 排序类和方法简介
2009/09/15 Javascript
Extjs Ext.MessageBox.confirm 确认对话框详解
2010/04/02 Javascript
使用JavaScript动态设置样式实现代码及演示动画
2013/01/25 Javascript
JavaScript怎么判断图片是否加载完成以便获取其尺寸
2014/05/08 Javascript
node.js中的fs.statSync方法使用说明
2014/12/16 Javascript
javascript实现类似java中getClass()得到对象类名的方法
2015/07/27 Javascript
jQuery实现无限往下滚动效果代码
2016/04/16 Javascript
js实现页面a向页面b传参的方法
2016/05/29 Javascript
jquery使用on绑定a标签无效 只能用live解决
2016/06/02 Javascript
移动端H5开发 Turn.js实现很棒的翻书效果
2016/06/20 Javascript
element UI upload组件上传附件格式限制方法
2018/09/04 Javascript
vue webpack重写cookie路径的方法
2019/07/10 Javascript
JavaScript 函数用法详解【函数定义、参数、绑定、作用域、闭包等】
2020/05/12 Javascript
[15:23]教你分分钟做大人:虚空假面
2014/10/30 DOTA
Python中List.index()方法的使用教程
2015/05/20 Python
详谈套接字中SO_REUSEPORT和SO_REUSEADDR的区别
2018/04/28 Python
Python对象属性自动更新操作示例
2018/06/15 Python
Python和Sublime整合过程图示
2019/12/25 Python
Pycharm自动添加文件头注释和函数注释参数的方法
2020/10/23 Python
CSS实现进度条和订单进度条的示例
2020/11/05 HTML / CSS
新西兰演唱会和体育门票网站:Ticketmaster新西兰
2017/10/07 全球购物
机关门卫岗位职责
2013/12/30 职场文书
拖鞋店创业计划书
2014/01/15 职场文书
人力资源主管岗位职责
2014/01/29 职场文书
市场营销工作计划书
2014/05/06 职场文书
人身损害赔偿协议书范本
2014/09/27 职场文书
大学生党员自我剖析材料
2014/10/06 职场文书
陈斌强事迹观后感
2015/06/17 职场文书
疾病证明书
2015/06/19 职场文书
运动会致辞稿
2015/07/29 职场文书
baselines示例程序train_cartpole.py的ImportError
2022/05/20 Python