相机内参的标定方法

2023-05-16

简介

摄像机标定(Camera calibration)简单来说是从世界坐标系换到图像坐标系的过程,也就是求最终的投影矩阵 PP 的过程,下面相关的部分主要参考UIUC的计算机视觉的课件(网址Spring 2016 CS543 / ECE549 Computer vision)。

基本的坐标系:

  • 世界坐标系(world coordinate system);
  • 相机坐标系(camera coordinate system);
  • 图像坐标系(image coordinate system);

一般来说,标定的过程分为两个部分:

  • 第一步是从世界坐标系转换为相机坐标系,这一步是三维点到三维点的转换,包括 RR,tt (相机外参)等参数;
  • 第二部是从相机坐标系转为图像坐标系,这一步是三维点到二维点的转换,包括 KK(相机内参)等参数;

相机坐标系 转换到 图像坐标系

坐标系介绍

相机到图片

如上图所示(图片来自UIUC计算机视觉课件),是一个小孔成像的模型,其中:

  • CC 点表示camera centre,即相机的中心点,也是相机坐标系的中心点;
  • ZZ 轴表示principal axis,即相机的主轴;
  • pp 点所在的平面表示image plane,即相机的像平面,也就是图片坐标系所在的二维平面;
  • pp 点表示principal point,即主点,主轴与像平面相交的点;
  • CC 点到 pp 点的距离,也就是右边图中的 ff 表示focal length,即相机的焦距;
  • 像平面上的 xx 和 yy 坐标轴是与相机坐标系上的 XX 和 YY 坐标轴互相平行的;
  • 相机坐标系是以 XX, YY, ZZ(大写)三个轴组成的且原点在 CC 点,度量值为米(m);
  • 像平面坐标系是以 xx,yy(小写)两个轴组成的且原点在 pp 点,度量值为米(m);
  • 图像坐标系一般指图片相对坐标系,在这里可以认为和像平面坐标系在一个平面上,不过原点是在图片的角上,而且度量值为像素的个数(pixel);

相机 转换到 像平面

知道上面的简单知识后,如果知道相机坐标系中的一个点 XX(现实三维世界中的点),在像平面坐标系对应的点是 xx,要求求从相机坐标系转为像平面坐标系的转换,也就是从 XX 点的(X,Y,Z)(X,Y,Z)通过一定的转换变为 xx 点的(x,y)(x,y)。注意:(X,Y,Z)(X,Y,Z)(大写)是在相机坐标系,而(x,y)(x,y)(小写)是在像平面坐标系(还不是图像坐标系,原点不同。)观察第二个图,很简单的可以得到这个转换: 

x=fX/Z

y=fY/Zy

(X,Y,Z)↦(fX/Z,fY/Z)

加入偏移量

通过上面,可以把相机坐标系转换到像平面坐标系,但是像平面坐标系和图像坐标系虽然在同一个平面上,但是原点并不是同一个,而目标是要转换到图像坐标系下,所以还需要一步操作,如下图:

这里写图片描述

如上图所示(图片来自UIUC计算机视觉课件),其中主点 pp 是像平面坐标系的原点,但在图像坐标系中的位置为(px,py)(px,py),在这里,图形坐标系的原点是图片的左下角,所以可以得到: 

(X,Y,Z)↦(fX/Z+px,fY/Z+py)

整理后可以得到 K,也就是平时所说的相机内参(Intrinsic parameters): 

K=\begin{bmatrix} f & &p_x\\ & f & p_y\\ & & 1 \end{bmatrix}

f为焦距,即相机镜头光心到成像平面的距离。

像素坐标

前面也提到了在图像坐标系中用的不是现实生活中的m来度量,而是用的 pixel 的个数,所以在上面转换到图像坐标系中还有个问题,就是坐标的表示还是m,并没有转换到像素坐标系统;在这里需要引入一个新概念就是:

  • mx表示在水平方向1m的长度包含的像素的个数;
  • my 表示在竖直方向1m的长度包含的像素的个数;

可能有人奇怪为啥不是一个值,还需要分别指定 mxmx 和 mymy 呀,这是因为通过上面可以得到一个像素点的大小(m度量)为: 

\frac{1}{m_x{}}\times \frac{1}{m_y{}}

但是需要说明的是像素并不一定是一个正方形,有时候可能也是一个矩形,所以要分别指定。 

一般来说,在使用相机内参K计算坐标系转换时,提供的都是已经变换后的值;例如会提供 fx,fy,cx,cy 四个值代表相机内参K,其实 fx就是这里的 αx,同理 fy 是 αy,cx 是 βx,cy 是 βy。

世界坐标系 转换到 图像坐标系

坐标系介绍

如上图所示(图片来自UIUC计算机视觉课件),从世界坐标系转换到相机坐标系是三维空间到三维空间的变换,一般来说需要一个平移操作和一个旋转操作就可以完成这个转换,用公式表示如下(可以理解为世界坐标系原点先平移到相机坐标系的位置然后在做一次坐标系旋转,使坐标轴对齐。): 

X~cam=R(X~−C~)

  • RR 表示旋转矩阵;
  • X˜X~ 表示 XX 点在世界坐标系中的位置;
  • C˜C~ 表示相机原点 CC 在世界坐标系中的位置;
  • X˜camX~cam 表示 XX 点在相机坐标系中的位置;

世界 转换到 相机

根据上面的公式可以得到从一个三维点从世界坐标系转换到相机坐标的变换公式如下(也是用的齐次坐标的表示方式): 

Xcam=(X˜cam1)=[R0−RC˜1](X˜1)=[R0−RC˜1]X

世界 转换到 图像

根据上面的讨论知道了怎样从世界坐标系转换到相机坐标系(平移和旋转)以及从相机坐标系转换到图像坐标系(相机内参变换),所以带入上面的矩阵计算,可以得到: 

x=K[I0]Xcam=K[R−RC˜]X

这样就得到了最终的投影矩阵 PP : 

P=K[Rt]

其中: 

t=−RC~

 

在这里,KK 一般称为相机内参(intrinsic parameters),描述了相机的内部参数,包括焦距 ff、主点 pp 的位置、以及像素与真实环境的大小比例等,这个是固有属性,是提供好的;RR 和 tt 称为相机外参(extrinsic parameters),RR 在这里是旋转矩阵,可以转换为三维的旋转向量,分别表示绕xx,yy,zz 三个轴的旋转角度,tt 目前就是一个平移向量,分别表示在xx,yy,zz 三个方向上的平移量。

畸变参数(distortion parameters)

在几何光学和阴极射线管(CRT)显示中,畸变(distortion) 是对直线投影(rectilinear projection)的一种偏移。简单来说直线投影是场景内的一条直线投影到图片上也保持为一条直线。那畸变简单来说就是一条直线投影到图片上不能保持为一条直线了,这是一种光学畸变(optical aberration)。可能由于摄像机镜头的原因,这里不讨论,有兴趣的可以查阅光学畸变的相关的资料。 
畸变一般可以分为两大类,包括径向畸变和切向畸变。主要的一般径向畸变有时也会有轻微的切向畸变。

径向畸变(Radial distortion

径向畸变的效应有三种,一种是桶形畸变(barrel distortion),另一种是枕形畸变(pincushion distortion),还有一种是两种的结合叫做胡子畸变(mustache distortion),从图片中可以很容易看出区别,具体见下图(图片来自wikipedia):

径向畸变可以用如下公式修正: 

xcorr=xdis(1+k1r2+k2r4+k3r6)

ycorr=ydis(1+k1r2+k2r4+k3r6)

切向畸变(tangential distortion

切向畸变是由于透镜与成像平面不严格的平行,其可以用如下公式修正: 

xcorr=xdis+[2p1xy+p2(r2+2x2)]

ycorr=ydis+[p1(r2+2y2)+2p2xy]

其中:

  • xdisxdis 和 ydisydis 表示有畸变的坐标;
  • xcorrxcorr 和 ycorrycorr 表示修复后的坐标;
  • k1k1,k2k2,k3k3 表示径向畸变参数;
  • p1p1,p2p2 表示切向畸变参数;

所以最终得到5个畸变参数: 

D=(k1,k2,p1,p2,k3)

 

相机标定

那么可以利用这些来进行最终的任务相机标定,简单的过程可以描述为通过标定板,如下图,可以得到n个对应的世界坐标三维点 XiXi 和对应的图像坐标二维点 xixi,这些三维点到二维点的转换都可以通过上面提到的相机内参 KK,相机外参 RR和 tt,以及畸变参数 DD 经过一系列的矩阵变换得到。现在就用这些对应关系来求解这些相机参数。最后就是用线性方法求解方程式,这里就不做讨论了。

那为什么要做相机标定呢? 
每个镜头的畸变程度各不相同,通过相机标定可以校正这种镜头畸变。其实可以认为用这种标定的方式来求解相机内参和畸变参数,相当于一种相机校准,然后这些参数就可以用于后面的求解。例如求解新拍的两幅图片相对的 RR 和 tt,求解这个外参用到就是标定得到的相机内参和畸变参数。

齐次坐标 
就是将一个原本是n维的向量用一个n+1维向量来表示。 
许多图形应用涉及到几何变换,主要包括平移、旋转、缩放。以矩阵表达式来计算这些变换时,平移是矩阵相加,旋转和缩放则是矩阵相乘,综合起来可以表示为 x=R∗X+tx=R∗X+t(注:因为习惯的原因,实际使用时一般使用变化矩阵左乘向量)(RR 旋转缩放矩阵,tt 为平移矩阵,XX 为原向量,xx 为变换后的向量)。引入齐次坐标的目的主要是合并矩阵运算中的乘法和加法,表示为 x=P∗Xx=P∗X 的形式。即它提供了用矩阵运算把二维、三维甚至高维空间中的一个点集从一个坐标系变换到另一个坐标系的有效方法。和上面的计算过程是对应的。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

相机内参的标定方法 的相关文章

随机推荐

  • 【EHub_tx1_tx2_E100】Ubuntu18.04 + ROS_ Melodic + WHEELTEC 9轴IMU测试(如何在该环境下配置驱动/RVIZ查看)

    简介 xff1a 介绍WHEELTEC 9轴IMU 在EHub tx1 tx2 E100载板 xff0c TX1核心模块环境 xff08 Ubuntu18 04 xff09 下测试ROS驱动 xff0c 配置驱动环境 xff0c RVIZ查
  • #kali网卡配置以及乱码问题解决 #title1:kali-linux-2019.4-i386 #title2:network

    kali网卡配置以及乱码问题解决 title1 kali linux 2019 4 i386 title2 network 0x01 改网卡配置 xff0c 重启网卡服务 保证虚拟机系统桥接 vi etc network interface
  • git错误记录:git-lfs filter-process: git-lfs: command not found

    参考文章 The Git init amp fetch approach breaks Git LFS 问题描述 在pull远程主分支代码时 xff0c 因为合作开发的伙伴上传了个静态库 xff0c 大文件拉取不下来 xff0c 出现了 g
  • 智能车调试

    首先 xff1a 摄像头采集图像 xff0c 采集是一张一张的 xff0c 在每一张图片然后根据二维数组识别黑白边界 xff0c 进而通过每一行得到黑白边界求出中值点 xff0c 车随中值点构成的中值线驱动电机的行驶 xff0c 电机有加速
  • STM32串口不定长度传输策略

    目录 一 接收中断 43 空闲中断 二 接收中断 43 外加定时器 三 接收中断 43 自定义结束符 四 串口DMA接收不定长数据 一 接收中断 43 空闲中断 RXNE中断和IDLE中断的区别 xff1f 当接收到1个字节 xff0c 就
  • 【AUTOSAR】【信息安全】CSM

    目录 一 概述 二 依赖模块 三 功能描述 3 1 基本体系结构 3 2 通用行为 3 2 1 正常操作 3 2 2 设计说明 3 3 错误分类 3 3 1 开发错误 3 3 2 运行时错误 四 API接口 4 1 通用接口 4 2 加密接
  • ::在c++中什么意思

    34 34 在C 43 43 中表示作用域 xff0c 和所属关系 34 34 是运算符中等级最高的 xff0c 它分为三种 xff0c 分别如下 xff1a 一 作用域符号 xff1a 作用域符号 的前面一般是类名称 xff0c 后面一般
  • 在eclipse上配置使用tomcat

    在eclipse上配置使用tomcat 1 确认自己已经下载好tomcat后 xff0c 打开eclipse后选择菜单栏中的windows列表中的首选项配置 2 找到Server选项中的Runtime Envirnm选项 xff0c 如图所
  • 本科学完C语言、C++、python(学透点),还有必要学别的语言吗?

    原作者是一名高校的信息技术类的教师 xff0c 主教程序设计类课程 这样的问题 xff0c 作者的学生也会经常问他 本篇文章意于为各位大学生提供一些在编程上的疑惑 xff0c 希望能够对大家有帮助 作者 xff1a 悟空问答丨EXCEL进阶
  • C/C++编程笔记:C/C++中的strrchr()函数,到底该怎么用?

    在C 43 43 中 xff0c strrchr xff08 xff09 是用于字符串处理的预定义函数 cstring是字符串函数所需的头文件 此函数返回一个指针 xff0c 该指针指向字符串中最后一次出现的字符 我们想要找到的最后一个出现
  • C++编程书籍推荐:零基础入门书籍,学C++看它们就够了!

    如果你是一个没有编程经验的C 43 43 零基础小白 xff0c 或者有其它语言经验的C 43 43 初学者 xff0c 那么强烈推荐下面的十本零基础小白入门C 43 43 书籍 1 C 43 43 Primer 作者 xff1a Stan
  • 【ROS2 入门】虚拟机环境 ubuntu 18.04 ROS2 安装

    大家好 xff0c 我是虎哥 xff0c 从今天开始 xff0c 我将花一段时间 xff0c 开始将自己从ROS1切换到ROS2 xff0c 做为有别于ROS1的版本 xff0c 做了很多更新和改变 xff0c 我还是很期待自己逐步去探索R
  • 如何解压.gz的压缩文件

    如何解压 gz的压缩文件 gzip d xxx gz tar命令 root 64 linux tar cxtzjvfpPN 文件与目录 参数 xff1a c xff1a 建立一个压缩文件的参数指令 create 的意思 xff1b x xf
  • GPS经纬度坐标与XY坐标相互转换的python程序

    文章目录 前言一 说明二 函数1 import 和 常数2 GPS经纬度转XY坐标3 XY坐标转GPS经纬度 总结 前言 室外定位常用的是GPS xff0c 故编队队形 设定轨迹都是基于GPS经纬度坐标 而在仿真中我们通常会在XY坐标系下进
  • AD20 原理图设计流程

    Altium Designer 20 的原理图设计大致可以分为 9 个步骤 xff1a xff08 1 xff09 新建原理图 这是原理图设计的第一步 xff08 2 xff09 图纸设置 图纸设置就是要设置图纸的大小 xff0c 方向等信
  • JavaScript基础——DOM节点操作学习笔记

    目录 笔记 方法的使用 案例一 动态生成表格 案例二 下拉菜单 xff0c 鼠标经过和离开实现 案例全部代码 笔记 节点概述 1 网页中的任何内容都是节点 文字 标签 元素 文档等 节点至少有nodeType 节点类型 nodeName 节
  • MAVLINK包的校验方法

    这段时间做一个项目要进行MAVLINK的解包校验 xff0c 但有一个叫做 CRC EXTRA的位导致这个校验码怎么算结果都不对 xff0c 后来找了好久还是在github的论坛上看见别人讨论才找到方法的 1 先上从官网上拿的mavlink
  • 机器人工程专业课程

    1 机器人工程专业的课程主要有 xff1a 高级语言程序设计 电路分析 机械设计基础 模拟电路技术 数字电子技术 自动控制原理 微机原理及接口技术 电机与电气控制技术 单片机原理及其应用 机械制造基础 工业机器人控制系统 运动控制系统 工业
  • python获取当前执行py文件的绝对路径

    python获取当前执行py文件的绝对路径 python3 home appuser test py span class token comment 获取当前执行py文件的绝对路径 span py file path span class
  • 相机内参的标定方法

    简介 摄像机标定 Camera calibration 简单来说是从世界坐标系换到图像坐标系的过程 xff0c 也就是求最终的投影矩阵 PP 的过程 xff0c 下面相关的部分主要参考UIUC的计算机视觉的课件 xff08 网址Spring