Python opencv相机标定实现原理及步骤详解


Posted in Python onApril 09, 2020

相机标定相机标定的目的

获取摄像机的内参和外参矩阵(同时也会得到每一幅标定图像的选择和平移矩阵),内参和外参系数可以对之后相机拍摄的图像就进行矫正,得到畸变相对很小的图像。

相机标定的输入

标定图像上所有内角点的图像坐标,标定板图像上所有内角点的空间三维坐标(一般情况下假定图像位于Z=0平面上)。

相机标定的输出

摄像机的内参、外参系数。

拍摄的物体都处于三维世界坐标系中,而相机拍摄时镜头看到的是三维相机坐标系,成像时三维相机坐标系向二维图像坐标系转换。不同的镜头成像时的转换矩阵不同,同时可能引入失真,标定的作用是近似地估算出转换矩阵和失真系数。为了估算,需要知道若干点的三维世界坐标系中的坐标和二维图像坐标系中的坐标,也就是拍摄棋盘的意义。

相机成像

相机的成像原理:小孔成像

Python opencv相机标定实现原理及步骤详解

相机的内参相机的外参

在实际由于设计工艺问题、相机安装环境或物体摆放位置等影响,会照成成像与实际图像不一样的现象。

由于设计工艺照成的影响是无法改变的事实,所以这将是相机的内参;

由环境或安装方式照成的影响是可以改变的,这就是相机的外参。

张正友标定相机原理

1.求得相机内参数:

用于标定的棋盘格是特制的,其角点坐标已知。标定棋盘格是三维场景中的一个平面∏,棋盘格在成像平面为π(知道了∏与π的对应点坐标之后,可求解两个平面1对应的单应矩阵H)。

根据相机成像模型,P为标定的棋盘坐标,p为其像素点坐标。则Python opencv相机标定实现原理及步骤详解,通过对应的点坐标求解H后,可用于求K,R,T。

2.设棋盘格所在平面为世界坐标系上XOY平面,则棋盘格上任一角点P世界坐标系为(X,Y,0)。

Python opencv相机标定实现原理及步骤详解

3、内参约束条件

Python opencv相机标定实现原理及步骤详解

Python opencv相机标定实现原理及步骤详解

Python opencv相机标定实现原理及步骤详解

Python opencv相机标定实现原理及步骤详解

实验步骤打印棋盘图片(网上找一张)

Python opencv相机标定实现原理及步骤详解

将打印出的纸固定放到一个平面上,使用同一相机从不同的位置,不同的角度,拍摄标定板的多张照片(我拍了15张)手机型号是华为mate9

Python opencv相机标定实现原理及步骤详解

提取标定板的世界坐标

标定板的大小是标定板在水平和竖直方向上内角点的个数。内角点指的是,标定板上不挨着边界的角点。

我打印的是6x9的标定板。

Python opencv相机标定实现原理及步骤详解

标定相机

Python opencv相机标定实现原理及步骤详解

Python opencv相机标定实现原理及步骤详解

Python opencv相机标定实现原理及步骤详解

mtx -->内参数矩阵

dist --> 畸变系数

rvecs --> 旋转向量

tvecs --> 平移向量

我们可以通过反投影误差来评估结果的好坏,越接近0,说明结果越理想。

通过之前计算的内参数矩阵、畸变系数、旋转矩阵和平移向量,使用cv2.projectPoints()计算三维点到二维图像的投影,然后计算反投影得到的点与图像上检测到的点的误差,最后计算一个对于所有标定图像的平均误差即反投影误差

我的棋盘打印出来有些不平整,可能是打印的纸张没有放正,导致有些地方翘着,效果不是很好,误差值有些大了,把纸张贴平整应该会好很多。而且我可能拍照的角度变化不是太大,可以试着把拍照的角度更加差异些,结果会更明显。

import cv2
import numpy as np
import glob
# 设置寻找亚像素角点的参数,采用的停止准则是最大循环次数30和最大误差容限0.001
criteria = (cv2.TERM_CRITERIA_MAX_ITER | cv2.TERM_CRITERIA_EPS, 30, 0.001)
# 获取标定板角点的位置
objp = np.zeros((6 * 9, 3), np.float32)
objp[:, :2] = np.mgrid[0:9, 0:6].T.reshape(-1, 2) # 将世界坐标系建在标定板上,所有点的Z坐标全部为0,所以只需要赋值x和y
obj_points = [] # 存储3D点
img_points = [] # 存储2D点
images = glob.glob("E:/test_pic/qipan/*.jpg")
for fname in images:
  img = cv2.imread(fname)
  cv2.imshow('img',img)
  gray = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY)
  size = gray.shape[::-1]
  ret, corners = cv2.findChessboardCorners(gray, (6, 9), None)
  print(ret)
  if ret:
    obj_points.append(objp)
    corners2 = cv2.cornerSubPix(gray, corners, (5, 5), (-1, -1), criteria) # 在原角点的基础上寻找亚像素角点
    #print(corners2)
    if [corners2]:
      img_points.append(corners2)
    else:
      img_points.append(corners)

    cv2.drawChessboardCorners(img, (8, 6), corners, ret) # 记住,OpenCV的绘制函数一般无返回值
    cv2.imshow('img', img)
    cv2.waitKey(2000)
print(len(img_points))
cv2.destroyAllWindows()
# 标定
ret, mtx, dist, rvecs, tvecs = cv2.calibrateCamera(obj_points, img_points, size, None, None)
print("ret:", ret)
print("mtx:\n", mtx) # 内参数矩阵
print("dist:\n", dist) # 畸变系数  distortion cofficients = (k_1,k_2,p_1,p_2,k_3)
print("rvecs:\n", rvecs) # 旋转向量 # 外参数
print("tvecs:\n", tvecs ) # 平移向量 # 外参数
print("-----------------------------------------------------")

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

Python 相关文章推荐
介绍Python的Urllib库的一些高级用法
Apr 30 Python
Python实现Logger打印功能的方法详解
Sep 01 Python
详解python函数传参是传值还是传引用
Jan 16 Python
对python中GUI,Label和Button的实例详解
Jun 27 Python
运用PyTorch动手搭建一个共享单车预测器
Aug 06 Python
Python判断字符串是否xx开始或结尾的示例
Aug 08 Python
python实现智能语音天气预报
Dec 02 Python
Python异常继承关系和自定义异常实现代码实例
Feb 20 Python
Python拼接字符串的7种方式详解
Mar 19 Python
Visual Studio code 配置Python开发环境
Sep 11 Python
selenium.webdriver中add_argument方法常用参数表
Apr 08 Python
Python中OpenCV实现查找轮廓的实例
Jun 08 Python
浅谈python 中的 type(), dtype(), astype()的区别
Apr 09 #Python
利用python绘制数据曲线图的实现
Apr 09 #Python
利用python生成照片墙的示例代码
Apr 09 #Python
Python 改变数组类型为uint8的实现
Apr 09 #Python
pandas 强制类型转换 df.astype实例
Apr 09 #Python
python 给图像添加透明度(alpha通道)
Apr 09 #Python
Python日志logging模块功能与用法详解
Apr 09 #Python
You might like
《PHP编程最快明白》第七讲:php图片验证码与缩略图
2010/11/01 PHP
php中inlcude()性能对比详解
2012/09/16 PHP
ThinkPHP3.2.2的插件控制器功能简述
2014/07/09 PHP
各种常用浏览器getBoundingClientRect的解析
2009/05/21 Javascript
checkbox设置复选框的只读效果不让用户勾选
2013/08/12 Javascript
jQuery 三击事件实现代码
2013/09/11 Javascript
JS按字节截取字符长度实例
2013/11/20 Javascript
JavaScript继承基础讲解(原型链、借用构造函数、混合模式、原型式继承、寄生式继承、寄生组合式继承)
2014/08/16 Javascript
AngularJS基础 ng-keypress 指令简单示例
2016/08/02 Javascript
详解Node.js:events事件模块
2016/11/24 Javascript
AngularJS 教程及实例代码
2017/10/23 Javascript
vue2.0使用swiper组件实现轮播的示例代码
2018/03/03 Javascript
vue2.0 实现页面导航提示引导的方法
2018/03/13 Javascript
基于vue cli 通过命令行传参实现多环境配置
2018/07/12 Javascript
Quasar Input:type="number" 去掉上下小箭头 实现加减按钮样式功能
2020/04/09 Javascript
Vue 组件注册全解析
2020/12/17 Vue.js
[02:40]DOTA2英雄基础教程 炼金术士
2013/12/23 DOTA
[43:35]TI4 循环赛第二日Liquid vs Fnatic
2014/07/11 DOTA
python基础入门学习笔记(Python环境搭建)
2016/01/13 Python
代码分析Python地图坐标转换
2018/02/08 Python
运用TensorFlow进行简单实现线性回归、梯度下降示例
2018/03/05 Python
pandas 将索引值相加的方法
2018/11/15 Python
python如何制作英文字典
2019/06/25 Python
Flask教程之重定向与错误处理实例分析
2019/08/01 Python
Django文件存储 自己定制存储系统解析
2019/08/02 Python
Pytorch反向求导更新网络参数的方法
2019/08/17 Python
关于阿里云oss获取sts凭证 app直传 python的实例
2019/08/20 Python
使用Python进行防病毒免杀解析
2019/12/13 Python
python3跳出一个循环的实例操作
2020/08/18 Python
python多线程爬取西刺代理的示例代码
2021/01/30 Python
python中time.ctime()实例用法
2021/02/03 Python
工业自动化专业毕业生推荐信
2013/11/18 职场文书
幼儿园家长评语大全
2014/04/16 职场文书
2015年建筑工程工作总结
2015/05/13 职场文书
高中生社会实践心得体会
2016/01/14 职场文书
简单介绍 http请求响应参数、无连接无状态、MIME、状态码、端口、telnet、curl
2021/03/31 HTML / CSS