相机标定:张正友标定原理

2023-11-11 00:12
文章标签 原理 相机 标定 张正友

本文主要是介绍相机标定:张正友标定原理,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

本文来自公众号“AI大道理”

 ——————

计算机视觉的源头是相机,因此我们有必要对相机有所了解。

原始相机拍摄的图像一般都会有所畸变,导致画面和实际观测的有所排查,为了让相机拍摄的图像和肉眼观察的一致,就需要进行相机标定,获得相机参数,从而进行校正。

1、相机模型

相机是如何成像的?

光束从物体表面反射,经过相机镜头,到达感光原件,这一系列物理过程可以通过数学公式表达,最终变成一个简单的矩阵操作将三维空间中的点对应到图片的一个像素。

2、为什么要相机标定?

我们知道实际使用的相机镜头都是透镜,物理上只有通过光心的光线才是沿直线传播的,而大部分的光线在通过透镜后会发生折射,从而在一定程度上改变传播的角度。越靠近透镜的边缘,改变的角度也就越大,这会造成相机所成的像产生距离上的拉伸以及形状的改变,这个现象称为相机畸变。

而标定操作其实就是通过一系列的计算校准后得到修正参数,通过这些参数修正后就可以得到与我们人眼看到的景象相同的图像,也就是将三维景象转换成去除畸变后的二维图像。

进行摄像机标定的目的:求出相机的内、外参数,以及畸变参数。
标定相机后通常是想做两件事:一个通过相机标定可以校正这种镜头畸变矫正畸变,生成矫正后的图像;另一个是根据获得的图像重构三维场景。  

3、相机的畸变

通常分为三大类,径向畸变、离心畸变以及薄棱镜畸变。

(1)径向畸变

使像点产生径向位置的偏差称为径向畸变。这也是相机中最常见的一种畸变,因为光线通过透镜一定会产生折射。径向畸变又可以具体分为正向枕形畸变、和负向桶形畸变。

(2)切向畸变

相机的光学中心与几何中心不一致(镜头各器件的光学中心)所造成的畸变称为离心畸变。广义上讲,离心畸变既包含径向畸变,又包含镜头主光轴不对称所造成的切向畸变(不过我们通常将径向畸变单独分成一类)。

(3)薄棱镜畸变

镜头设计缺陷与加工安装误差所造成的畸变称为薄棱镜畸变,这种畸变会同时引起径向畸变和切向畸变(想象一下相机的镜头安装的不正,或者镜头不够平整)。高价位的镜头可以忽略薄棱镜畸变。

4、标定方法

  • 相机自标定法

  • 基于主动视觉的标定方法

  • 基于标定物的相机标定方法

5、张正友标定法

张氏标定法属于第三种方法。

张氏标定,简单来说就是个通过单应矩阵的约束求解相机内参的过程。

张氏标定法只关注径向畸变,通过五个左右的畸变系数进行修正,通常为k1,k2,p1,p2,k3。

k的个数通过计算的结果与修正的效果进行调整,可多可少,实际计算过程中,如果考虑太多高阶的畸变参数,会导致标定求解的不稳定。

矩阵A包含5个元素,需要3组H才能解出A的唯一封闭解,因此在标定时需拍摄3组以上的图片。

1)图像点获取

张正友的方法是,拿一个棋盘格图像,通过角点检测,得到角点在图像里的坐标。

2)棋盘的排放

算法是基于2D模型的,如果棋盘摆放的不平整,肯定会造成很大的影像。

3)图像多少

标定原理上至少3张,一般要大于10张。

4)图片的角度

太大的角度对于角点提取的精度影像比较大,所以保持在45度以内比较好。

5)张正友标定的空间坐标系

以标定板左上角第一个角点为原点,水平为x坐标,垂直为y坐标,标定板到相机的距离为z坐标。

6、相机标定参数

标定模型:

世界坐标->相机坐标->图像坐标->像素坐标:

内参:

fx、fy: 单位长度的像素个数

cx、cy: 主点坐标

k1、k2、k3、p1、p2:  畸变系数, 径向三个,切向两个

外参:
R、T: 旋转和平移矩阵

实例:

现以NiKon D700相机为例进行求解其内参数矩阵:   
焦距 f=35mm ,最高分辨率:4256×2832 ,传感器尺寸:36.0×23.9mm。
根据以上定义可以有:
u0=4256/2=2128, v0=2832/2=1416, dx=36.0/4256, dy=23.9/2832, 

fx = f/dx = 4137.8  , fy = f/dy = 4147.3。

其中相机的内参和外参可以通过张正友标定获取。

7、相机标定实验

步骤1:固定相机,不同角度移动标定板获取图片,至少三张

步骤2:python相机标定或者matlab相机标定

matlab:

python:

1)特征点检测

获取像素点imgpoints,对棋盘进行角点检测,映射到图片中获取像素点。共35个,因为5*7。

[1,0,0]就是黑白格的坐标,以1为单位。

2)获取特征点对应的像素点的坐标

[1417,435]是特征点在图像的坐标。

3)空间坐标与像素坐标一一对应的标定

cv2.calibrateCamera(objpoints, imgpoints, gray.shape[::-1], None, None)

步骤3:旋转向量->旋转矩阵

相机标定得到的R和T如下:

可见并非我们所熟悉的3*3,而是每张图片都有一个对应的R和T,这个并不是我们想要的结果,因为我们想要唯一的一个参数使得可以应对所有情况。

外参还需要单独求解。

外参标定算法及流程:

求解外参是一个PNP(Perspective-N-Point)问题。

如果场景的三维结构已知,利用多个控制点在三维场景中的坐标及其在图像中的透视投影坐标即可求解出摄像机坐标系与表示三维场景结构的世界坐标系之间的绝对位姿关系,包括绝对平移向量T以及旋转矩阵R

对于透视投影来说,要使得PNP问题有确定解,需要至少三组控制点。

 

# 求解旋转矩阵

retval, rvec, tvec = cv2.solvePnP(objpoints1, imgpoints1, mtx,

dist)

的确得到了唯一的一个参数,但发现这个参数就是上面的第一个图片的参数。

 

但是还不是3*3的,还需要进行转换。此时旋转参数叫旋转向量,要转为旋转矩阵。

转换方法为罗德里格斯旋转公式。

 

# 旋转向量和旋转矩阵的互相转换

rotationMatrix, _ = cv2.Rodrigues(rvec)

终于得到了一个适应全部图片的3*3的旋转矩阵,和1*3的平移矩阵。

这两个参数可以将世界坐标和相机坐标进行转换。

(灵魂的拷问:至于为什么就是第一张图片的参数呢?莫非第一张图片就能代表全部了?或者任何一种标定的图片的旋转矩阵就可以当做全局的旋转矩阵?还是说实验过程有问题?)

8、图像校正

使用上面获得的摄像头矩阵(内参矩阵)和畸变系数 进行畸变矫正。

校正公式:

径向畸变:

切向畸变:

综合以上两种畸变,得到畸变坐标:

输入图像:

校正方法:

cv2.undistort(gray, mtx, dist, None, mtx)

校正效果:

9、总结

标定相机后通常是想做两件事:

一个通过相机标定可以校正这种镜头畸变矫正畸变,生成矫正后的图像;

另一个是根据获得的图像重构三维场景。  

下来将利用相机内外参进行三维重构,即图像坐标转世界坐标。 

 ——————

浅谈则止,细致入微AI大道理

扫描下方“AI大道理”,选择“关注”公众号

—————————————————————

 

—————————————————————

投稿吧   | 留言吧

这篇关于相机标定:张正友标定原理的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/386236

相关文章

数据库原理与安全复习笔记(未完待续)

1 概念 产生与发展:人工管理阶段 → \to → 文件系统阶段 → \to → 数据库系统阶段。 数据库系统特点:数据的管理者(DBMS);数据结构化;数据共享性高,冗余度低,易于扩充;数据独立性高。DBMS 对数据的控制功能:数据的安全性保护;数据的完整性检查;并发控制;数据库恢复。 数据库技术研究领域:数据库管理系统软件的研发;数据库设计;数据库理论。数据模型要素 数据结构:描述数据库

计算机组成原理——RECORD

第一章 概论 1.固件  将部分操作系统固化——即把软件永恒存于只读存储器中。 2.多级层次结构的计算机系统 3.冯*诺依曼计算机的特点 4.现代计算机的组成:CPU、I/O设备、主存储器(MM) 5.细化的计算机组成框图 6.指令操作的三个阶段:取指、分析、执行 第二章 计算机的发展 1.第一台由电子管组成的电子数字积分和计算机(ENIAC) 第三章 系统总线

GaussDB关键技术原理:高性能(二)

GaussDB关键技术原理:高性能(一)从数据库性能优化系统概述对GaussDB的高性能技术进行了解读,本篇将从查询处理综述方面继续分享GaussDB的高性能技术的精彩内容。 2 查询处理综述 内容概要:本章节介绍查询端到端处理的执行流程,首先让读者对查询在数据库内部如何执行有一个初步的认识,充分理解查询处理各阶段主要瓶颈点以及对应的解决方案,本章以GaussDB为例讲解查询执行的几个主要阶段

【计算机组成原理】部分题目汇总

计算机组成原理 部分题目汇总 一. 简答题 RISC和CICS 简要说明,比较异同 RISC(精简指令集)注重简单快速的指令执行,使用少量通用寄存器,固定长度指令,优化硬件性能,依赖软件(如编译器)来提升效率。 CISC(复杂指令集)包含多样复杂的指令,能一条指令完成多步操作,采用变长指令,减少指令数但可能增加执行时间,倾向于硬件直接支持复杂功能减轻软件负担。 两者均追求高性能,但RISC

MySQL数据库锁的实现原理

MySQL数据库的锁实现原理主要涉及到如何确保在多用户并发访问数据库时,保证数据的完整性和一致性。以下是MySQL数据库锁实现原理的详细解释: 锁的基本概念和目的 锁的概念:在数据库中,锁是用于管理对公共资源的并发控制的机制。当多个用户或事务试图同时访问或修改同一数据时,数据库系统通过加锁来确保数据的一致性和完整性。 锁的目的:解决多用户环境下保证数据库完整性和一致性的问题。在并发的情况下,会

线性回归(Linear Regression)原理详解及Python代码示例

一、线性回归原理详解         线性回归是一种基本的统计方法,用于预测因变量(目标变量)与一个或多个自变量(特征变量)之间的线性关系。线性回归模型通过拟合一条直线(在多变量情况下是一条超平面)来最小化预测值与真实值之间的误差。 1. 线性回归模型         对于单变量线性回归,模型的表达式为:         其中: y是目标变量。x是特征变量。β0是截距项(偏置)。β1

标准分幅下的图幅号转换成经纬度坐标【原理+源代码】

最近要批量的把标准分幅下的图幅号转换成经纬度坐标,所以这两天写了个程序来搞定这件事情。 先举个例子说明一下这个程序的作用。 例如:计算出图幅号I50G021040的经纬度范围,即最大经度、最小经度、最大纬度、最小纬度。 运用我编写的这个程序,可以直接算出来,这个图幅号的经纬度范围,最大经度为115.3125°,最小经度为115.25°,最大纬度为31.167°,最小纬度为31.125°。

SpingBoot原理

配置优先级 SpringBoot配置的优先级从高到低依次为命令行参数、JNDI属性、Java系统属性、操作系统环境变量、外部配置文件、内部配置文件、注解指定的配置文件和编码中直接指定的默认属性。具体如下: 命令行参数:启动应用时,通过命令行指定的参数拥有最高优先级。例如,使用--server.port=8081会直接改变应用程序的端口,无论在什么配置文件中定义过该值。JNDI属性:这些属性由当

HashMap 的工作原理及其在 Java 中的应用?

在Java的数据结构中,HashMap是最常见且最重要的一个数据结构之一。HashMap是Java集合框架中的一部分,它存储的是键值对(Key-value)映射,也就是说,你可以通过键(Key)找到对应的值(Value)。让我们来详细地看一下HashMap的工作原理。 HashMap的工作原理 HashMap内部有一个数组,数组中的每个元素又是一个链表。当我们将一个键值对存入HashM

Ajax及其工作原理

Ajax及其工作原理 AJAX 是一种与服务器交换数据无需刷新网页的技术,最早由Google公司在谷歌地图里使用,并迅速风靡。 AJAX是不能跨域的,如需跨域,可以使用document.domain='a.com';或者使用服务器代理,代理XMLHttpRequest文件 AJAX是基于现有的Internet标准,并且联合使用它们: XMLHttpRequest 对象 (异步的与服