相机标定解析

点击上方“小白学视觉”，选择加"星标"或“置顶”

重磅干货，第一时间送达

转载自：3d tof
编辑：东岸因为@一点人工一点智能
01  为什么要进行相机标定随着机器视觉的迅猛发展，我们已经不满足于使用摄像机进行监控、抓拍这种较为简单的功能。更多的用户青睐于它在非接触三维尺寸测量上的应用。我们所谓的三维测量是广义的三维测量，它不仅包括三维物体的重构与测量，还包括在三维空间中识别任意二维平面上的尺寸以及位置。这种技术目前已被应用在高精度的工业模具以及装配测量中，其中任意二维平面上的尺寸检测技术应用得更为广泛。

图一如图1当被测平面和像平面平行且成像模型为理想的小孔成像模型，我们设焦距为、工作距离为，则被测物和它的像关系可简单的表示为：但是在实际应用中并非如此，我们无法严格控制像平面和被测平面的位置，所用的镜头也不是严格的小孔模型。如果直接使用【1】式计算将会产生极大的误差。因此，为了获取更高的测量精度，我们需要通过标定来实现坐标平面的转换以及图像的校正。02  什么是相机标定在实际应用中，被测平面的不确定性以及镜头的畸变使我们已经无法简单的使用【1】式计算出实际距离，但是我们可以将目前能够获得的数据进行转换，使这些数据符合【1】式的使用条件。也就是将任意坐标平面通过旋转和平移映射到理想坐标平面上，对有畸变的图像进行校正，让它成为符合小孔成像模型的像平面。有了这种方法，我们只要确定转换算法、校正算法以及【1】式中的参数就可以实现三维空间中任意平面上尺寸与位置的测量。我们将这种确定参数的过程称之为标定。03  相机单目标定相机标定的方法根据摄像机的数目可分为单目标定、双目标定以及多目标定。其中单目相机标定是双目标定的基础，而多目相机的标定则是双目相机的扩展。因此，我们今天首先来为大家介绍单目标定。在平面测量中影响我们拍摄图像形变的因素有两个：镜头和相机姿态。根据这两个因素我们将摄像机的参数分为两组，相机内参和相机外参。
3.1 相机内参内参一般包括镜头的焦距、镜头畸变参数、光轴中心坐标以及像元尺寸，当摄像机和镜头确定时，这些参数唯一确定。下面我们来详细介绍一下各参数的数学模型。1) 焦距根据镜头类型不同焦距的计算可分为针孔模型和远心模型。如图3我们假设世界坐标系有任一点，在摄像机靶面所成的像为，根据不同的光路模型它们之间有如下的对应关系a) 针孔模型b) 远心模型由于远心镜头特殊的光路设计使得像的大小与拍摄距离无关，因此表达式比针孔模型更为简单。2) 镜头畸变受到镜头的制作和安装精度的影响，我们所获得图像会产生非线性失真。我们称这种失真为镜头畸变。镜头畸变产生的误差使得理想针孔模型已不再适用。因此我们需要先将所得图像进行校正，再应用理想的针孔模型。假设我们所获取的原始图像坐标、校正后的结果，其畸变模型坐标关系为：a) 径向畸变
径向畸变主要由透镜制作过程中表面曲率引起，它会使图像发生桶形畸变和枕形畸变（如图4）。其数学模型如下：其中，如果对精度要求不是很高我们可以令将上式简化为如下表达式：b) 离心畸变离心变量又称偏心变量，它的误差来源于透镜的安装精度，这主要是因为所有镜片的光学中心并不能严格的保证在同一条直线上。这种误差除了在引入径向畸变同时还会引入切向畸变。由于之前我们已经进行了径向畸变的校正，因此我们在此基础上只需加入切向畸变校正即可。其数学模型如下：c) 薄棱镜畸变
影响薄棱镜畸变的主要因素的是透镜以及摄像机靶面的平行度，镜片与摄像机靶面夹角越大畸变就越严重。其数学表达式如下：d) 畸变校正在实际的应用中，大多数工业摄像机的厂商可以通过摄像机接口螺纹的机械精度来保证镜头透镜与靶面的平行性，而且这种畸变产生的误差较小，因此在一般的图像标定中不作考虑。至此，我们已经基本掌握了大多数情况下畸变产生的原因以及数学模型。结合【5】、【7】两式我们可以推导出镜头畸变校正模型：
3.2 相机外参相机外参是指摄像机坐标系与世界坐标系的转换参数它主要由旋转矩阵R和平移矩阵T组成。对于任意三维坐标系，我们都可以通过这两个矩阵将其转换到摄像机坐标系中。其数学模型为 【10】
1) 旋转矩阵如图5所示，我们通过沿坐标轴分别旋转来实现坐标系的转换。因此，旋转矩阵R可分解为三个矩阵相乘的形式。我们以轴为例，假设与轴夹角为，且到原点距离为，通过旋转矩阵坐标系沿轴旋转后得到点，我们可得方程组:由三角函数展开得:由与轴夹角为得：将【15】式带入【13】式、【16】式带入【14】式得：因此
以此类推求得和将它们与相乘得：2) 平移矩阵通过旋转矩阵运算后，世界坐标系的三个坐标轴会与摄像机坐标系对应的坐标轴相平行。此时我们已经离我们的目标又近了一步。如图7所示我们现在只要沿各坐标轴做平移运算即可，由此得：3) 参数求解根据【2】、【9】、【20】、【21】这几个数学模型，我们可以得知，若想确定一个摄像机与被测平面的相对位置，则需要确定包括内参、外参在内的14个参数 ，其中是已知的。因此，我们至少需要9个坐标点，构成9个方程才可以解出剩余的9个未知数。在通常情况下，点的分布以覆盖大部分视场为准，获得的数据点越多，统计的参数就越准确。我们一般采用最小二乘法或者线性规划等统计算法来求解相应参数。下面我们以HALCON为例演示一个标定的全过程。HALCON是德国MVTEC Software GmbH公司开发的一套完善的机器视觉算法软件包。它除了拥有亚像素精度的算法以及高效的处理性能外，在三维重构方面它也有卓越的表现。它的开发环境中自带相机标定工具，可以轻松的完成相机标定工作。同时，您还可以使用HALCON生成可打印的标定板文件。下面我们就来介绍一下相机标定的整个流程。a) 生成标定板i. 创建标定板使用HALCON开发环境HDevelop创建标定板在选择尺寸选择时推荐大家使用边长为视野1/3左右的标定板                ii. 打印标定板
通过GSView等高精度打印软件打印标定板。b) 相机参数设置输入所使用的摄像机以及镜头的相关参数 c) 拍摄标定板图像我们采用平移和倾斜的方式使得拍摄图像中的标定板尽量覆盖整个视场。一般情况我们需要保存15幅不同位置的图像，具体位置如下：d) 标定图像载入通过HALCON我们可以实时拍摄图像也可以读取我们事先拍摄好的图像进行标定。如果标定板识别成功，图像上将绘制出标定板坐标系。e) 标定标定后我们会获得相机内参和外参。我们还可以将它们保存起来用于坐标转换或图像校正。  f) 验证标定结果在HALCON中不仅有相机标定工具还有测量工具，下面我们使用HDevelop自带的一维测量助手来验证一下我们的标定结果。首先我们要加载我们刚才标定的数据。成功加载标定数据后我们就可以使用我们随机拍摄的一张图像进行测量。通过以上几步操作，现在我们就已经得到了像素点的实际距离。复杂的标定工作在HALCON的帮助下就轻松的完成了。不仅如此，以上的所有步骤都可以导出为代码，我们可以将这些代码集成到我们自己的程序中。
备注：若采用高精度的标定板将会得到更高精度的标定结果。04  应用单目相机标定技术适用于被测表面曲率较小且需要获取实际数据的应用场合，比如食品、机械制造以及半导体等。理论上，在检测过程中所有的被测物尺寸的判断都可以使用像素尺寸来完成。虽然像素尺寸和实际尺寸对于计算机来说几乎是一样的，一般只用于比较大小，但是对于我们来说实际尺寸更加直观。
图15（左图）；图16（右图）图15为速冻鱼柳检测系统的图像处理示例。由于实际应用中不同批次的鱼柳宽度不同，因此，最好的办法是使用标定后获取的实际数据来进行计算。这种判定方法使得系统参数更加直观，设置更加方便。图16为工件尺寸检测系统，通过标定，我们不仅可以获取实际数据，而且可以将它们直接和CAD中数据进行比较，提高检测效率。05  结语在机器视觉产业的蓬勃发展的今天，相机标定已经逐步应用到医疗、食品、磨具生产、半导体生产等诸多检测系统中。随着像HALCON这样的算法库用户群的不断扩大，相机标定在作为我们研究课题的同时，也会不断的走进工业应用中。它将为我们提供更完善，更精准的二维以及三维空间的解决方案，成为提高工业生产效率和产品质量的推动力。

声明：部分内容来源于网络，仅供读者学习、交流之目的。文章版权归原作者所有。如有不妥，请联系删除。

下载1：OpenCV-Contrib扩展模块中文版教程

在「小白学视觉」公众号后台回复：扩展模块中文教程，即可下载全网第一份OpenCV扩展模块教程中文版，涵盖扩展模块安装、SFM算法、立体视觉、目标跟踪、生物视觉、超分辨率处理等二十多章内容。


下载2：Python视觉实战项目52讲
在「小白学视觉」公众号后台回复：Python视觉实战项目，即可下载包括图像分割、口罩检测、车道线检测、车辆计数、添加眼线、车牌识别、字符识别、情绪检测、文本内容提取、面部识别等31个视觉实战项目，助力快速学校计算机视觉。


下载3：OpenCV实战项目20讲
在「小白学视觉」公众号后台回复：OpenCV实战项目20讲，即可下载含有20个基于OpenCV实现20个实战项目，实现OpenCV学习进阶。


交流群

欢迎加入公众号读者群一起和同行交流，目前有SLAM、三维视觉、传感器、自动驾驶、计算摄影、检测、分割、识别、医学影像、GAN、算法竞赛等微信群（以后会逐渐细分），请扫描下面微信号加群，备注：”昵称+学校/公司+研究方向“，例如：”张三 + 上海交大 + 视觉SLAM“。请按照格式备注，否则不予通过。添加成功后会根据研究方向邀请进入相关微信群。请勿在群内发送广告，否则会请出群，谢谢理解~