最简单的针孔相机模型标定方法,使用Matlab进行快速可视化相机标定,小白专用,超详细

本文主要是介绍最简单的针孔相机模型标定方法,使用Matlab进行快速可视化相机标定,小白专用,超详细,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

相机标定是机器视觉、SLAM等绕不开的环节,今天介绍一种最为简单的可视化标定方法。
  
  这篇文章只是讲解如果标定,关于相机的标定原理,可以参考别的文章,这里推荐一篇《【机器视觉】张氏法相机标定》

一. 准备工作

  1. 方格标定板一块(实际上没有也行);
  2. 已经安装Matlab的电脑一台(Matlab最低版本要求是R2013b,推荐使用最新版本R2019b);

注意:
1.如果你有正版Matlab账号,或者你是学生,你可以使用网页版的Matlab;
2.如果你电脑上没有Matlab,而且你也没有正版账号,那么你可以选择这个网站提供的云端Matlab,你只需要免费注册个账号,然后将图片上传云端就可以按照如下方法进行标定。

二. 开始标定

1.第一步 打印标定板

如果你还没有标定,那么你需要打印一块标定板出来,但是你放心,准备标定板的过程是非常简单的。你先将我提供的标定板图片下载——下载标定板(提取码:5yun),然后打印出来,只要你保持图片的横纵比,打印成任何尺寸都可以,我推荐打印在A4纸上(如果你不方便打印,其实你也可以直接下载下来显示在电脑屏幕上也是一样的,也不需要固定大小,只要横纵比不变就行)。打印之后的样式,如下图:
在这里插入图片描述

2. 拍摄包含标定板的图片

使用将要被标定的相机,拍摄一些包含标定板的照片,根据我的经验20到30张图片已经够了。当然图片越多结果也会越准确,但是会多用一些时间。拍摄图片的技巧是,要多个角度大范围的去拍摄标定板,让它成像在你相机屏幕的各个地方,大家可以参考我下面的拍摄方式。
注意: 拍摄照片时有几点需要注意:标定板的摆放没有要求,横竖都行。但是拍照片时对相机的姿态有要求,你相机如果是最终横着使用,那么拍照片时相机也要横着拍。至于这里面的原因,可以参考我的另外一篇博客:《相机内参fx,fy,cx,cy的方向问题!》
在这里插入图片描述

3. 使用Matlab进行相机内参数标定

① 打开已经安装好的Matlab软件,菜单栏中依次点击 App -> Camera Calibrator,如下图:在这里插入图片描述
②在打开的Camera Calibrator工具中,依次点击Add Image -> From File,如下图:
在这里插入图片描述
③然后选择你刚刚已经拍摄好的照片,点击Open按钮,如下图:
在这里插入图片描述
④弹出的Checkboard Square Size窗口中,直接点击OK按钮,如下图:
在这里插入图片描述
注意:蓝色框中的数字是标定板一个小方格的边长,这里不用管它,数值对最终的标定结果影响不大。
⑤ 选择相机模型,如果你使用的是普通的针孔相机模型,请选择Standard,如果是鱼眼相机,就选择Fisheye,如下图:
在这里插入图片描述
⑥ 然后选择需要标定的参数。点击Options选择想要的畸变参数形式,然后点击,Calibrate按钮开始标定,如下图:
注意: 在畸变参数中,径向畸变是默认勾选2 Coefficients,也就是说径向畸变只使用两个系数,如果你不是很明白其中的道理,那你就用默认状态;切向畸变是可选择的,如果你需要切向畸变,请勾选Tangential Distortiorskew很少用到,它表示的是是单个像素x,y轴方向的夹角系数。在这里插入图片描述
⑦ 稍等一会,就能看到标定工具界面发生了变化,右边多出来两个窗口,红框中表示的是每一张图片投影产生的误差是几个像素,横坐标是图片的序号,纵坐标是误差。绿色框中表示的是通过照片还原得到的拍摄角度。
在这里插入图片描述
⑧ 鼠标点击红框中图标上的红线,将它往下拉,拉到0.8以内就可以了,那些误差大的照片就会被选中,如下图:
在这里插入图片描述
⑨ 上一步拉动红线被选中的图片,会在标定工具左侧窗口中选中,蓝色状态的图片就是误差较大的图片,如下图:
在这里插入图片描述
⑩ 在任意一张被选中的图片上,单击右键,选择Remove and Recalibrate删掉这些图片,然后就会自动重新标定,如下图:
在这里插入图片描述
⑪等待标定完成之后,依次点击Export Camera Parameters -> OK,将标定结果输出到Matlab命令行中,如下图:
在这里插入图片描述
⑫关闭相机标定工具页面,在Matlab主界面命令行中,可以看到最终的标定结果,如下图:
在这里插入图片描述

三. 参数介绍:

相机的内参矩阵:
K = [ f x 0 c x 0 f y c y 0 0 1 ] K = \left[ \begin{matrix}f_x & 0 & c_x \\ 0 & f_y & c_y \\ 0 & 0 & 1 \end{matrix} \right] K=fx000fy0cxcy1其中 FocalLength分别对应 f x fx fx f y f_y fy
   PrincipalPoint:分别对应 c x c_x cx c y c_y cy
  
  相机的径向畸变公式 x d i s t o r t e d = x ( 1 + k 1 r 2 + k 2 r 4 + k 3 r 6 ) x_{distorted}=x(1+k_1r^2+k_2r^4+k_3r^6) xdistorted=x(1+k1r2+k2r4+k3r6) y d i s t o r t e d = y ( 1 + k 1 r 2 + k 2 r 4 + k 3 r 6 ) y_{distorted}=y(1+k_1r^2+k_2r^4+k_3r^6) ydistorted=y(1+k1r2+k2r4+k3r6)其中RadialDistortion分别对应 k 1 k_1 k1 k 2 k_2 k2
  
  相机的切向畸变公式 x d i s t o r t e d = x + 2 p 1 x y + p 2 ( r 2 + 2 x 2 ) x_{distorted} = x+2p_1xy+p_2(r^2+2x^2) xdistorted=x+2p1xy+p2(r2+2x2) y d i s t o r t e d = y + p 1 ( r 2 + 2 y 2 ) + 2 p 2 x y y_{distorted} = y+p_1(r^2+2y^2)+2p_2xy ydistorted=y+p1(r2+2y2)+2p2xy其中TangentialDistortion分别对应 p 1 p_1 p1 p 2 p_2 p2

这篇关于最简单的针孔相机模型标定方法,使用Matlab进行快速可视化相机标定,小白专用,超详细的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

C++使用栈实现括号匹配的代码详解

《C++使用栈实现括号匹配的代码详解》在编程中,括号匹配是一个常见问题,尤其是在处理数学表达式、编译器解析等任务时,栈是一种非常适合处理此类问题的数据结构,能够精确地管理括号的匹配问题,本文将通过C+... 目录引言问题描述代码讲解代码解析栈的状态表示测试总结引言在编程中,括号匹配是一个常见问题,尤其是在

Python调用Orator ORM进行数据库操作

《Python调用OratorORM进行数据库操作》OratorORM是一个功能丰富且灵活的PythonORM库,旨在简化数据库操作,它支持多种数据库并提供了简洁且直观的API,下面我们就... 目录Orator ORM 主要特点安装使用示例总结Orator ORM 是一个功能丰富且灵活的 python O

Nginx设置连接超时并进行测试的方法步骤

《Nginx设置连接超时并进行测试的方法步骤》在高并发场景下,如果客户端与服务器的连接长时间未响应,会占用大量的系统资源,影响其他正常请求的处理效率,为了解决这个问题,可以通过设置Nginx的连接... 目录设置连接超时目的操作步骤测试连接超时测试方法:总结:设置连接超时目的设置客户端与服务器之间的连接

Java中String字符串使用避坑指南

《Java中String字符串使用避坑指南》Java中的String字符串是我们日常编程中用得最多的类之一,看似简单的String使用,却隐藏着不少“坑”,如果不注意,可能会导致性能问题、意外的错误容... 目录8个避坑点如下:1. 字符串的不可变性:每次修改都创建新对象2. 使用 == 比较字符串,陷阱满

Java判断多个时间段是否重合的方法小结

《Java判断多个时间段是否重合的方法小结》这篇文章主要为大家详细介绍了Java中判断多个时间段是否重合的方法,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录判断多个时间段是否有间隔判断时间段集合是否与某时间段重合判断多个时间段是否有间隔实体类内容public class D

Python使用国内镜像加速pip安装的方法讲解

《Python使用国内镜像加速pip安装的方法讲解》在Python开发中,pip是一个非常重要的工具,用于安装和管理Python的第三方库,然而,在国内使用pip安装依赖时,往往会因为网络问题而导致速... 目录一、pip 工具简介1. 什么是 pip?2. 什么是 -i 参数?二、国内镜像源的选择三、如何

使用C++实现链表元素的反转

《使用C++实现链表元素的反转》反转链表是链表操作中一个经典的问题,也是面试中常见的考题,本文将从思路到实现一步步地讲解如何实现链表的反转,帮助初学者理解这一操作,我们将使用C++代码演示具体实现,同... 目录问题定义思路分析代码实现带头节点的链表代码讲解其他实现方式时间和空间复杂度分析总结问题定义给定

IDEA编译报错“java: 常量字符串过长”的原因及解决方法

《IDEA编译报错“java:常量字符串过长”的原因及解决方法》今天在开发过程中,由于尝试将一个文件的Base64字符串设置为常量,结果导致IDEA编译的时候出现了如下报错java:常量字符串过长,... 目录一、问题描述二、问题原因2.1 理论角度2.2 源码角度三、解决方案解决方案①:StringBui

Linux使用nload监控网络流量的方法

《Linux使用nload监控网络流量的方法》Linux中的nload命令是一个用于实时监控网络流量的工具,它提供了传入和传出流量的可视化表示,帮助用户一目了然地了解网络活动,本文给大家介绍了Linu... 目录简介安装示例用法基础用法指定网络接口限制显示特定流量类型指定刷新率设置流量速率的显示单位监控多个

Java覆盖第三方jar包中的某一个类的实现方法

《Java覆盖第三方jar包中的某一个类的实现方法》在我们日常的开发中,经常需要使用第三方的jar包,有时候我们会发现第三方的jar包中的某一个类有问题,或者我们需要定制化修改其中的逻辑,那么应该如何... 目录一、需求描述二、示例描述三、操作步骤四、验证结果五、实现原理一、需求描述需求描述如下:需要在