在OpenCV里使用特征匹配和单映射变换的代码详解


Posted in Python onOctober 23, 2019

前面已经学习特征查找和对应匹配,接着下来在特征匹配之后,再使用findHomography函数来找出对应图像的投影矩阵。首先使用一个查询图片,然后在另外一张图片里找到目标对象,其实就是想在图片里查找所需要目标的一部分区域。为了实现这样的功能,需要使用calib3d库里的一个函数cv.findHomography(),把从两张图片里找到的特征点当作参数,传送给这个函数,然后这个函数返回一个投影变换矩阵,我们就可以使用 cv.perspectiveTransform()函数来对查找的目标进行投影,这样就可以在复杂图片里标记出相应的目标位置。

我们已经看到,在匹配时可能会有一些错误,这可能会影响结果。为了解决这个问题,需要使用RANSAC 或 LEAST_MEDIAN算法。所以提供正确估计的良好匹配称为内聚,其余的称为外联。cv.findHomography()函数返回一个值表示内聚还是外联的点。

在例子里,先使用ORB来寻找两个图片的特征点,接着根据设置条件为10个匹配特征,如果满足就会计算投影变换矩阵,一旦获得3x3的矩阵,就可以把寻找的目标对象在图片里标记出来。最后在复杂的图片里用白色线条标记出来。

在OpenCV里使用特征匹配和单映射变换的代码详解

参数详解:

srcPoints 源平面中点的坐标矩阵,可以是CV_32FC2类型,也可以是vector<Point2f>类型

dstPoints 目标平面中点的坐标矩阵,可以是CV_32FC2类型,也可以是vector<Point2f>类型

method 计算单应矩阵所使用的方法。不同的方法对应不同的参数,具体如下:

0 - 利用所有点的常规方法

RANSAC - RANSAC-基于RANSAC的鲁棒算法

LMEDS - 最小中值鲁棒算法

RHO - PROSAC-基于PROSAC的鲁棒算法

ransacReprojThreshold

将点对视为内点的最大允许重投影错误阈值(仅用于RANSAC和RHO方法)。如果

则点被认为是个外点(即错误匹配点对)。若srcPoints和dstPoints是以像素为单位的,则该参数通常设置在1到10的范围内。

mask

可选输出掩码矩阵,通常由鲁棒算法(RANSAC或LMEDS)设置。 请注意,输入掩码矩阵是不需要设置的。

maxIters RANSAC算法的最大迭代次数,默认值为2000。

confidence 可信度值,取值范围为0到1.

该函数能够找到并返回源平面和目标平面之间的转换矩阵H,以便于反向投影错误率达到最小。

在OpenCV里使用特征匹配和单映射变换的代码详解

演示使用的例子如下:

#python 3.7.4,opencv4.1
#蔡军生 https://blog.csdn.net/caimouse/article/details/51749579
#
import numpy as np
import cv2
from matplotlib import pyplot as plt
 
MIN_MATCH_COUNT = 10
 
#读取文件
img1 = cv2.imread('rmb3.png')
img2 = cv2.imread('rmb4.png')
 
#初始化ORB检测器
orb = cv2.ORB_create()
 
#用ORB查找关键点
kp1, des1 = orb.detectAndCompute(img1,None)
kp2, des2 = orb.detectAndCompute(img2,None)
 
# FLANN参数
FLANN_INDEX_LSH = 6
index_params = dict(algorithm = FLANN_INDEX_LSH,
          table_number = 6, 
          key_size = 12,   
          multi_probe_level = 1)
search_params = dict(checks=50)  #或者使用一个空的字典
flann = cv2.FlannBasedMatcher(index_params,search_params)
 
matches = flann.knnMatch(des1,des2,k=2)
 
# 比率
good = []
for m,n in matches:
  if m.distance < 0.7*n.distance:
    good.append(m)
if len(good)>MIN_MATCH_COUNT:
  src_pts = np.float32([ kp1[m.queryIdx].pt for m in good ]).reshape(-1,1,2)
  dst_pts = np.float32([ kp2[m.trainIdx].pt for m in good ]).reshape(-1,1,2)
  #找到投影变换矩阵
  M, mask = cv2.findHomography(src_pts, dst_pts, cv2.RANSAC,5.0)
  matchesMask = mask.ravel().tolist()
  #进行投影变换
  h,w,d = img1.shape
  pts = np.float32([ [0,0],[0,h-1],[w-1,h-1],[w-1,0] ]).reshape(-1,1,2)
  dst = cv2.perspectiveTransform(pts,M)
  #画变换后的外形
  img2 = cv2.polylines(img2,[np.int32(dst)],True,(255,255,255),3, cv2.LINE_AA)
 
else:
  print( "Not enough matches are found - {}/{}".format(len(good), MIN_MATCH_COUNT) )
  matchesMask = Non
 
draw_params = dict(matchColor = (0,255,0), # draw matches in green color
          singlePointColor = None,
          matchesMask = matchesMask, # draw only inliers
          flags = 2)
 
img3 = cv2.drawMatches(img1,kp1,img2,kp2,good,None,**draw_params)
 
#显示图片
cv2.imshow('img3',img3)
      
cv2.waitKey(0)
cv2.destroyAllWindows()

结果输出如下:

在OpenCV里使用特征匹配和单映射变换的代码详解

总结

以上所述是小编给大家介绍的在OpenCV里使用特征匹配和单映射变换的代码详解,希望对大家有所帮助,如果大家有任何疑问请给我留言,小编会及时回复大家的。在此也非常感谢大家对三水点靠木网站的支持!
如果你觉得本文对你有帮助,欢迎转载,烦请注明出处,谢谢!

Python 相关文章推荐
python使用xlrd模块读写Excel文件的方法
May 06 Python
python查看zip包中文件及大小的方法
Jul 09 Python
Python实现配置文件备份的方法
Jul 30 Python
Python脚本实现Web漏洞扫描工具
Oct 25 Python
Python中int()函数的用法浅析
Oct 17 Python
Numpy掩码式数组详解
Apr 17 Python
python 限制函数调用次数的实例讲解
Apr 21 Python
python3.6中@property装饰器的使用方法示例
Aug 17 Python
浅谈python print(xx, flush = True) 全网最清晰的解释
Feb 21 Python
python实现mean-shift聚类算法
Jun 10 Python
pycharm-professional-2020.1下载与激活的教程
Sep 21 Python
详解pycharm的python包opencv(cv2)无代码提示问题的解决
Jan 29 Python
手把手教你Python yLab的绘制折线图的画法
Oct 23 #Python
Python之Numpy的超实用基础详细教程
Oct 23 #Python
Python从列表推导到zip()函数的5种技巧总结
Oct 23 #Python
Python箱型图绘制与特征值获取过程解析
Oct 22 #Python
Python使用贪婪算法解决问题
Oct 22 #Python
python元组和字典的内建函数实例详解
Oct 22 #Python
Python List列表对象内置方法实例详解
Oct 22 #Python
You might like
PHP的加密方式及原理
2012/06/14 PHP
php实现俄罗斯乘法实例
2015/03/07 PHP
php编程每天必学之表单验证
2016/03/01 PHP
Yii2实现ajax上传图片插件用法
2016/04/28 PHP
Laravel框架实现定时发布任务的方法
2018/08/16 PHP
js的event详解。
2006/09/06 Javascript
任意位置显示html菜单
2007/02/01 Javascript
最新28个很棒的jQuery 教程
2011/05/28 Javascript
ajax不执行success回调而是执行了error回调
2012/12/10 Javascript
分享Javascript中最常用的55个经典小技巧
2013/11/29 Javascript
js获取和设置属性的方法
2014/02/20 Javascript
Node.js 学习笔记之简介、安装及配置
2015/03/03 Javascript
第三章之Bootstrap 表格与按钮功能
2016/04/25 Javascript
AngularJS基础 ng-srcset 指令简单示例
2016/08/03 Javascript
Dropzone.js实现文件拖拽上传功能(附源码下载)
2016/11/22 Javascript
D3.js中强制异步文件读取同步的几种方法
2017/02/06 Javascript
详解用vue-cli来搭建vue项目和webpack
2017/04/20 Javascript
详解webpack 打包文件体积过大解决方案(code splitting)
2018/04/10 Javascript
highCharts提示框中显示当前时间的方法
2019/01/18 Javascript
python+opencv实现动态物体识别
2018/01/09 Python
python数据封装json格式数据
2018/03/04 Python
ubuntu系统下使用pm2设置nodejs开机自启动的方法
2018/05/12 NodeJs
Python requests设置代理的方法步骤
2020/02/23 Python
python梯度下降算法的实现
2020/02/24 Python
tensorflow dataset.shuffle、dataset.batch、dataset.repeat顺序区别详解
2020/06/03 Python
Python切片列表字符串如何实现切换
2020/08/06 Python
python利用线程实现多任务
2020/09/18 Python
HTTP状态码详解
2021/03/18 杂记
css实例教程 一款纯css3实现的超炫动画背画特效
2014/11/05 HTML / CSS
凯普林包包西班牙官网:Kipling西班牙
2019/04/12 全球购物
学雷锋标语
2014/06/25 职场文书
2014小学二年级班主任工作总结
2014/12/05 职场文书
小学教师师德师风自我评价
2015/03/04 职场文书
2016反腐倡廉警示教育心得体会
2016/01/13 职场文书
2016年党员干部廉政承诺书
2016/03/24 职场文书
Django使用redis配置缓存的方法
2021/06/01 Redis