在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 相关文章推荐
Java多线程编程中ThreadLocal类的用法及深入
Jun 21 Python
Python使用QQ邮箱发送Email的方法实例
Feb 09 Python
Python简单网络编程示例【客户端与服务端】
May 26 Python
python装饰器深入学习
Apr 06 Python
Python装饰器的执行过程实例分析
Jun 04 Python
Python中new方法的详解
Jan 15 Python
python绘制评估优化算法性能的测试函数
Jun 25 Python
python+jinja2实现接口数据批量生成工具
Aug 28 Python
Python中的延迟绑定原理详解
Oct 11 Python
Python CSV文件模块的使用案例分析
Dec 21 Python
Django生成数据库及添加用户报错解决方案
Oct 09 Python
使用Python通过oBIX协议访问Niagara数据的示例
Dec 04 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
德生S2000收音机更换“钕铁硼”全频扬声器
2021/03/02 无线电
php下保存远程图片到本地的办法
2010/08/08 PHP
用PHP书写安全的脚本代码
2012/02/05 PHP
PHP简单实现数字分页功能示例
2016/08/24 PHP
php Session无效分析资料整理
2016/11/29 PHP
Laravel框架源码解析之模型Model原理与用法解析
2020/05/14 PHP
转一个日期输入控件,支持FF
2007/04/27 Javascript
浏览器的JavaScript引擎的识别方法
2013/10/20 Javascript
教你如何使用firebug调试功能了解javascript闭包和this
2015/03/04 Javascript
用JavaScript实现页面重定向功能的教程
2015/06/04 Javascript
基于Jquery制作图片文字排版预览效果附源码下载
2015/11/18 Javascript
值得分享的轻量级Bootstrap Table表格插件
2016/05/30 Javascript
原生JS查找元素的方法(推荐)
2016/11/22 Javascript
react-native-tab-navigator组件的基本使用示例代码
2017/09/07 Javascript
使用jQuery实现两个div中按钮互换位置的实例代码
2017/09/21 jQuery
vue实现导航栏效果(选中状态刷新不消失)
2017/12/13 Javascript
vue登录路由验证的实现
2017/12/13 Javascript
vue.js数据绑定操作详解
2018/04/23 Javascript
Vue实现微信支付功能遇到的坑
2019/06/05 Javascript
angularjs请求数据的方法示例
2019/08/06 Javascript
Nodejs环境实现socket通信过程解析
2020/07/03 NodeJs
Vue移动端项目实现使用手机预览调试操作
2020/07/18 Javascript
python3实现短网址和数字相互转换的方法
2015/04/28 Python
Python functools模块学习总结
2015/05/09 Python
Tensorflow简单验证码识别应用
2017/05/25 Python
Pycharm保存不能自动同步到远程服务器的解决方法
2019/06/27 Python
Python中and和or如何使用
2020/05/28 Python
Python Charles抓包配置实现流程图解
2020/09/29 Python
倩碧香港官方网站:Clinique香港
2017/11/13 全球购物
世界排名第一的万圣节服装店:Spirit Halloween
2018/10/16 全球购物
Vans(范斯)新西兰官方网站:美国原创极限运动品牌
2020/09/19 全球购物
行政人员岗位职责
2013/12/08 职场文书
主持人演讲稿范文
2013/12/28 职场文书
初三化学教学反思
2014/01/23 职场文书
新年联欢会主持词
2014/03/27 职场文书
成都人事代理协议书
2014/10/25 职场文书