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 相关文章推荐
pygame学习笔记(4):声音控制
Apr 15 Python
python读取TXT到数组及列表去重后按原来顺序排序的方法
Jun 26 Python
Python入门_条件控制(详解)
May 16 Python
python去掉空白行的多种实现代码
Mar 19 Python
selenium+python设置爬虫代理IP的方法
Nov 29 Python
python面向对象实现名片管理系统文件版
Apr 26 Python
梅尔频率倒谱系数(mfcc)及Python实现
Jun 18 Python
画pytorch模型图,以及参数计算的方法
Aug 17 Python
基于Python批量生成指定尺寸缩略图代码实例
Nov 20 Python
python 解决微分方程的操作(数值解法)
May 26 Python
Python下opencv库的安装过程及问题汇总
Jun 11 Python
浅谈Python中的正则表达式
Jun 28 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 数组教程 定义数组
2009/10/23 PHP
浅析Yii2中GridView常见操作
2016/04/22 PHP
PHP中ID设置自增后不连续的原因分析及解决办法
2016/08/21 PHP
由Javascript实现的页面日历
2011/11/04 Javascript
jquery中的事件处理详细介绍
2013/06/24 Javascript
javaScript 逻辑运算符使用技巧整理
2017/05/03 Javascript
Js利用console计算代码运行时间的方法示例
2017/09/24 Javascript
React-Native之定时器Timer的实现代码
2017/10/04 Javascript
微信小程序使用map组件实现检索(定位位置)周边的POI功能示例
2019/01/23 Javascript
[01:15:16]DOTA2-DPC中国联赛 正赛 Elephant vs Aster BO3 第一场 1月26日
2021/03/11 DOTA
Python实现SVN的目录周期性备份实例
2015/07/17 Python
python向已存在的excel中新增表,不覆盖原数据的实例
2018/05/02 Python
Python3 中sorted() 函数的用法
2020/03/24 Python
python mysql中in参数化说明
2020/06/05 Python
利用python3筛选excel中特定的行(行值满足某个条件/行值属于某个集合)
2020/09/04 Python
Python 实现键盘鼠标按键模拟
2020/11/18 Python
CSS3 RGBA色彩模式使用实例讲解
2016/04/26 HTML / CSS
埃弗顿足球俱乐部官方网上商店:Everton Direct
2018/01/13 全球购物
德国Discount-Apotheke中文官网:DC德式康线上药房
2020/02/18 全球购物
屈臣氏俄罗斯在线商店:Watsons俄罗斯
2020/08/03 全球购物
综合素质的自我鉴定
2013/10/07 职场文书
普罗米修斯教学反思
2014/02/06 职场文书
课程改革实施方案
2014/03/16 职场文书
大学生就业策划书范文
2014/04/04 职场文书
情人节寄语大全
2014/04/11 职场文书
责任书格式范文
2014/07/28 职场文书
2014年新教师工作总结
2014/11/08 职场文书
十佳少年事迹材料
2014/12/25 职场文书
生死抉择观后感
2015/06/09 职场文书
亲戚关系证明
2015/06/24 职场文书
搞笑婚礼主持词开场白
2015/11/24 职场文书
2016优秀班主任个人先进事迹材料
2016/02/26 职场文书
2016年社区文体活动总结
2016/04/06 职场文书
详解MongoDB的条件查询和排序
2021/06/23 MongoDB
html form表单基础入门案例讲解
2021/07/21 HTML / CSS
详细聊聊浏览器是如何看闭包的
2021/11/11 Javascript