python opencv 直方图反向投影的方法


Posted in Python onFebruary 24, 2018

本文介绍了python opencv 直方图反向投影的方法,分享给大家,具体如下:

目标:

直方图反向投影

原理:

反向投影可以用来做图像分割,寻找感兴趣区间。它会输出与输入图像大小相同的图像,每一个像素值代表了输入图像上对应点属于目标对象的概率,简言之,输出图像中像素值越高的点越可能代表想要查找的目标。直方图投影经常与camshift(追踪算法)算法一起使用。

算法实现的方法,首先要为包含我们感兴趣区域的图像建立直方图(样例要找一片草坪,其他的不要)。被查找的对象最好是占据整个图像(图像里全是草坪)。最好使用颜色直方图,物体的颜色信息比灰度图像更容易被分割和识别。再将颜色直方图投影到输入图像查找目标,也就是找到输入图像中每一个像素点的像素值在直方图中对应的概率,这样就得到一个概率图像,最后设置适当的阈值对概率图像进行二值化。

numpy算法:

建立两幅颜色直方图,目标图像直方图(M),输入图像直方图(I)

import cv2
import numpy as np
from matplotlib import pyplot as plt
#roi图片,就想要找的的图片
roi = cv2.imread('3.jpg')
hsv = cv2.cvtColor(roi,cv2.COLOR_BGR2HSV)
#目标搜索图片
target = cv2.imread('33.jpg')
hsvt = cv2.cvtColor(target,cv2.COLOR_BGR2HSV)
#创建直方图
M = cv2.calcHist([hsv],[0,1],None,[180,256],[0,180,0,256])
I = cv2.calcHist([hsvt],[0,1],None,[180,256],[0,180,0,256])

计算比值:R=MI。反向投影R,根据R这个调色板创建新图像,每一个像素代表这个点事目标的概率。例如,B(x,y)=R[h(x,y),s(x,y),其中H为点(x,y)的色调(hue)值,s为点(x,y)的饱和度(saturation)。最后加入条件B(x,y)=min([B(x,y),1]

h,s,v = cv2.split(hsvt)
B = R[h.ravel(),s.ravel()]
B = np.minimum(B,1)
B = B.reshape(hsvt.shape[:2])

使用圆盘算子做卷积,B=D×B,其中D为卷积核

disc = cv2.getStructuringElement(cv2.MORPH_ELLIPSE,(5,5))#定义结构形状,5×5的椭圆
B=cv2.filter2D(B,-1,disc)#对图像进行卷积运算
B = np.uint8(B)
cv2.normalize(B,B,0,255,cv2.NORM_MINMAX)

其中函数cv2.getStructuringElement是定义结构元素,例如element = cv2.getStructuringElement(cv2.MORPH_CROSS,(5,5)) 定义了一个十字形,5×5的结构

python opencv 直方图反向投影的方法 

输出图像中灰度最大的地方就是目标位置。如果要找的是一个区域,可以使用一个阈值对图像二值化,这样能得到不错的结果。

ret,thresh = cv2.threshold(B,50,255,0)

opencv反向投影

函数 cv2.calcBackProject()直接实现反向投影,参数与cv2.calcHist基本一致。其中一个参数是要查找的目标的直方图。在使用目标直方图反向投赢钱应该进行归一化处理。返回结果是一个概率图像,然后进行圆盘形状卷积操作,再二值化。

roi区域图片

python opencv 直方图反向投影的方法 

待搜索图片

python opencv 直方图反向投影的方法 

结果

python opencv 直方图反向投影的方法

import cv2
import numpy as np
from matplotlib import pyplot as plt
#roi图片,就想要找的的图片
roi = cv2.imread('33.jpg')
hsv = cv2.cvtColor(roi,cv2.COLOR_BGR2HSV)

#目标搜索图片
target = cv2.imread('3.jpg')
hsvt = cv2.cvtColor(target,cv2.COLOR_BGR2HSV)

#计算目标直方图
roihist = cv2.calcHist([hsv],[0,1],None,[180,256],[0,180,0,256])
#归一化,参数为原图像和输出图像,归一化后值全部在2到255范围
cv2.normalize(roihist,roihist,0,255,cv2.NORM_MINMAX)
dst = cv2.calcBackProject([hsvt],[0,1],roihist,[0,180,0,256],1)

#卷积连接分散的点
disc = cv2.getStructuringElement(cv2.MORPH_ELLIPSE,(5,5))
dst = cv2.filter2D(dst,-1,disc)

ret,thresh = cv2.threshold(dst,50,255,0)
#使用merge变成通道图像
thresh = cv2.merge((thresh,thresh,thresh))

#蒙板
res = cv2.bitwise_and(target,thresh)
#矩阵按列合并,就是把target,thresh和res三个图片横着拼在一起
res = np.hstack((target,thresh,res))

cv2.imwrite('res.jpg',res)
#显示图像
cv2.imshow('1',res)
cv2.waitKey(0)

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

Python 相关文章推荐
python使用urllib2模块获取gravatar头像实例
Dec 18 Python
Python多进程同步简单实现代码
Apr 27 Python
pycharm恢复默认设置或者是替换pycharm的解释器实例
Oct 29 Python
python模块之subprocess模块级方法的使用
Mar 26 Python
python取余运算符知识点详解
Jun 27 Python
python3.x+pyqt5实现主窗口状态栏里(嵌入)显示进度条功能
Jul 04 Python
python文字转语音实现过程解析
Nov 12 Python
Python简单实现区域生长方式
Jan 16 Python
Python random模块制作简易的四位数验证码
Feb 01 Python
Python递归实现打印多重列表代码
Feb 27 Python
python+OpenCV实现图像拼接
Mar 05 Python
pd.DataFrame中的几种索引变换的实现
Jun 16 Python
python爬虫爬取淘宝商品信息
Feb 23 #Python
python爬取淘宝商品详情页数据
Feb 23 #Python
Python如何抓取天猫商品详细信息及交易记录
Feb 23 #Python
python列表生成式与列表生成器的使用
Feb 23 #Python
1分钟快速生成用于网页内容提取的xslt
Feb 23 #Python
python使用xslt提取网页数据的方法
Feb 23 #Python
Python爬虫使用Selenium+PhantomJS抓取Ajax和动态HTML内容
Feb 23 #Python
You might like
PHP 字符截取 解决中文的截取问题,不用mb系列
2009/09/29 PHP
php实现简单的语法高亮函数实例分析
2015/04/27 PHP
PHP伪造来源HTTP_REFERER的方法实例详解
2015/07/06 PHP
javascript检查日期格式的函数[比较全]
2008/10/17 Javascript
jquery 弹出登录窗口实现代码
2009/12/24 Javascript
分享几个超级震憾的图片特效
2012/01/08 Javascript
js实现通用的微信分享组件示例
2014/03/10 Javascript
果断收藏9个Javascript代码高亮脚本
2016/01/06 Javascript
indexedDB bootstrap angularjs之 MVC DOMO (应用示例)
2016/06/20 Javascript
JS类的定义与使用方法深入探索
2016/11/26 Javascript
AngularJS入门教程之Helloworld示例
2016/12/25 Javascript
微信小程序实现倒计时60s获取验证码
2020/04/17 Javascript
JavaScript使用FileReader实现图片上传预览效果
2020/03/27 Javascript
使用 webpack 插件自动生成 vue 路由文件的方法
2019/08/20 Javascript
浅析vue-cli3配置webpack-bundle-analyzer插件【推荐】
2019/10/23 Javascript
jQuery实现王者荣耀手风琴效果
2020/01/17 jQuery
vue中watch和computed的区别与使用方法
2020/08/23 Javascript
vue+elementUI中表格高亮或字体颜色改变操作
2020/11/02 Javascript
js制作提示框插件
2020/12/24 Javascript
[10:54]Team Spirit vs Navi
2018/06/07 DOTA
Python3实现从指定路径查找文件的方法
2015/05/22 Python
基于python的七种经典排序算法(推荐)
2016/12/08 Python
python实现生命游戏的示例代码(Game of Life)
2018/01/24 Python
Python向Excel中插入图片的简单实现方法
2018/04/24 Python
python 实现批量xls文件转csv文件的方法
2018/10/23 Python
在Python 中同一个类两个函数间变量的调用方法
2019/01/31 Python
Python 经典算法100及解析(小结)
2019/09/13 Python
django ajax发送post请求的两种方法
2020/01/05 Python
解决python3.x安装numpy成功但import出错的问题
2020/11/17 Python
小狗电器官方商城:中国高端吸尘器品牌
2017/03/29 全球购物
Boden美国官网:英伦原创时装品牌
2017/07/03 全球购物
运动会加油口号
2014/06/07 职场文书
陈胜吴广起义口号
2014/06/20 职场文书
推广普通话标语
2014/06/27 职场文书
党员带头倡议书
2015/04/29 职场文书
如何起草一份正确的合伙创业协议书?
2019/07/04 职场文书