一种基于非线性优化的相机标定方法技术

本发明专利技术涉及一种基于非线性优化的相机标定方法，属于图像处理技术领域。该方法包括以下步骤：S1：获取数据：制作棋盘格标定板，并测量各方格的尺寸，从不同角度拍摄两张以上照片；S2：计算单应矩阵；S3：计算内参数、外参数和畸变因数的初始值；S4：用非线性优化方法精化计算结果。本发明专利技术在张正友标定算法的基础上，引入知识和约束条件，使得搜索范围更小，从而减小陷入局部最优的可能性，使算法更稳定；并且只需要打印一个棋盘格图形，测量其尺寸就可以标定相机，其精度明显优于比经典的张正友算法。

【技术实现步骤摘要】
一种基于非线性优化的相机标定方法
本专利技术属于图像处理
，涉及一种基于非线性优化的少图高精度相机标定方法。
技术介绍
目前，相机标定方法有线性变换法、非线性优化法、两步标定法等方法，(1)线性变换法：计算复杂度小，有唯一的解析解，但是由于相机标定对数据噪声相当敏感，该方法的结果精度较低。(2)非线性优化法：通过迭代优化得到一个最优解，可以得到高精度的结果，但相机标定是一个非凸优化问题，其精度对初始值高度依赖。(3)两步标定法(如张正友标定算法)：通过线性变换求得参数初始值，然后对其进行非线性优化，本质上是由粗到精的方法，以其精度高和标定板制作简单而得到广泛应用，但其在只有一到两张标定图时效果不佳。相机模型就是用表达式模拟相机成像过程，通常采用小孔成像模型，如图1所示，相机模型表达式为：其中，表示图像上的齐次坐标，表示标定点(标定板上的角点)的齐次坐标，I是3×3的单位矩阵，0是3个元素的0向量(列向量)，其中，r11-r33的3×3矩阵为旋转矩阵，太复杂，但可以用ψ、θ和参数表示，两者关系式为当不考虑相机畸变时，相机标定就是求K中的5个参数和R中的6个参数。但实际上在精度要求高的情况下必须考虑相机畸变，也就是说通过计算的结果与实际成像是有差异的，需要对结果进行修正，这个修正的表达式一般采用这里又有2个参数，所以考虑相机畸变时，相机标定就是求13个参数。直接计算非常困难，一般先算出P＝K[I|O]R,然后再分解出11个参数。但由于标定数据存在噪声，用解析方法计算的结果往往不精确，所以，一般只是将其作为非线性优化的初始值。因此，现有相机标定方法中，对少图相机标定方法，要么标定物特别复杂，要么精度不高。故，有必要研究一种少图高精度的相机标定方法。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种基于非线性优化的相机标定方法，在张正友标定算法的基础上，引入知识和约束条件，使得搜索范围更小，从而减小陷入局部最优的可能性，使算法更稳定。该方法只要有两张图甚至一张图就可以算出所有参数，然后通过图中的已经尺寸物体精确确定未知尺寸的物体等数据，适用于刑侦等特殊场合。为达到上述目的，本专利技术提供如下技术方案：一种基于非线性优化的相机标定方法，该方法具体包括以下步骤：S1：获取数据：制作棋盘格标定板，并测量各方格的尺寸，从不同角度拍摄两张以上照片；S2：计算单应矩阵；S3：计算内参数、外参数和畸变因数的初始值；S4：用非线性优化方法精化计算结果。进一步，所述步骤S2具体包括：令齐次坐标为的空间点在图像平面上的齐次坐标为有上式结果为2×12矩阵；将标定点(角点)坐标和及其图像上的像点坐标代入上式解方程组，得到单应矩阵P(是一个3×4的矩阵)；进一步，所述步骤S3具体包括：采用无畸变情况下的相机模型，如下式：其中，f表示焦距；uo、vo表示主点在图像中的坐标；r11～33表示旋转矩阵的各元素；t1～3表示平移向量的各元素、x、y表示三维空间中目标点的坐标，u、v表示目标点的像的图像坐标；zc表示相机坐标系中目标点的Z坐标。已知P，计算下式等号右侧中的参数：其中，P是3X4的矩阵，只取其中的第1、2、4列，即删除第3列后，得到另外一个3X3的矩阵H，其中的第i个列向量表示为：hi＝[hi1hi2hi3]T；(1)计算内参数：f和uo、vo求解如下方程组：即Vb＝0其中，vij＝[hi3hj3-(hi3hj1+hi1hj3)-(hi2hj3+hi3hj2)hi1hj1+hi2hj2]，其解为VTV的最小特征值对应的特征向量，即b向量的值；通过上式求解出第一个参数矩阵，称为内参数矩阵；(2)计算外参数：这个矩阵中的12个参数；令则通过下式求出第二个参数矩阵，称为外参数矩阵；其中，τ＝1/‖A-1h1‖，r1、r2、r3和t分别为外参数矩阵R中的各个列向量；(3)计算畸变因数：k＝(DTD)-1DTd其中，k＝[k1k2]T，u,v是用无畸变相机模型算出来的图像点，是实际的像点位置，u0,v0是内参数矩阵中的两个元素，k1,k2分别表示一阶和二阶畸变因数。进一步，所述步骤S4具体包括：采用Levenberg-Marquardt方法进一步优化内外参数和畸变因子，优化的目标函数为：其中，是通过相机模型算出的点，mij是实际的图像坐标，d(k,mij)是角点离其所在的行和列方向上的线的距离，k是控制因子，m,n分别表示每行和列的标定点数量。本专利技术的有益效果在于：本专利技术所述方法只需要打印一个棋盘格图形，测量其尺寸就可以标定相机，其精度明显优于经典的张正友标定算法。附图说明为了使本专利技术的目的、技术方案和效果说明更加清楚，本专利技术提供如下附图：图1为小孔成像模型图；图2为棋盘格标定板；图3为本专利技术所述相机标定方法流程图；图4为单图无噪声的成功率对比图；图5为本专利技术所述方法与张氏算法误差的均值对比图(单图无噪声)；图6为本专利技术所述方法与张氏算法误差的标准差对比图(单图无噪声)；图7为双图无噪声的成功率对比图；图8为本专利技术所述方法与张氏算法误差的均值对比图(双图无噪声)；图9为本专利技术所述方法与张氏算法误差的标准差对比图(双图无噪声)；图10为双图量化噪声的成功率对比图；图11为本专利技术所述方法与张氏算法误差的均值对比图(双图量化噪声)；图12为本专利技术所述方法与张氏算法误差的标准差对比图(双图量化噪声)。具体实施方式下面将结合附图，对本专利技术的优选实施例进行详细的描述。一种基于非线性优化的相机标定方法，如图3所示，具体包括以下步骤：S1：获取数据：如图2所示，制作棋盘格标定板，并测量各方格的尺寸，从不同角度拍摄两张以上照片；S2：计算单应矩阵；S3：计算内参数、外参数和畸变因数的初始值；S4：用非线性优化方法精化计算结果。仿真实验：分别对单图、双图以及实图(两张图)使用本专利技术所述相机标定方法和张正友标定方法进行相机标定处理，分别从成功率、准确度和精确度三方面进行对比分析两种方法的效果。成功率用标定结果为具有物理意义的实数的次数与总实验次数的比值表示。准确度计算结果与参考值的差的均值来度量。精确度用标定结果的标准差来度量。1、单图无噪声效果(1)成功率：如图4所示，图4a中灰色表示成功，黑色表示不成功，上图为本专利技术提出的算法，下图为经典的张正友算法；图4b为不同角点数下的结果；图4c为标定板与相机不同角度下的结果；图4b和4c中，实线为本专利技术提出的算法；虚线为张正友算法。(2)准确度：如图5所示，图5表示本专利技术所述方法与张氏算法的误差均值对比图，反映准确度。(3)精确度：如图6所示，图6表示本专利技术所述方法与张氏算法的误差标准差对比图，反映精确度。2、双图无噪声效果(1)成功率：如图7所示，图7a中灰色表示成功，黑色表示不成功，上图为本专利技术提出的算法，下图为经典的张正友算法；图7b为不同角点数下的结果；图7c为标定板与相机不同角度下的结果；图7b和7c中，实线为本专利技术提出的算法；虚线为张正友算法。(2)准确度：如图8所示，图8表示本专利技术所述方法与张氏算法误差的均值对比图，反映准确度。(3)精确度：如图9所示，图9表示本专利技术所述方法与张氏算法误差的标准差对比图，反映精确度。3、双图量化噪声效果(1)成功率：如图10所示，图10a中灰色表示成功，黑色表示本文档来自技高网...

【技术保护点】
1.一种基于非线性优化的相机标定方法，其特征在于，该方法具体包括以下步骤：S1：获取数据：制作棋盘格标定板，并测量各方格的尺寸，从不同角度拍摄两张以上照片；S2：计算单应矩阵；S3：计算内参数、外参数和畸变因数的初始值；S4：用非线性优化方法精化计算结果。

【技术特征摘要】
1.一种基于非线性优化的相机标定方法，其特征在于，该方法具体包括以下步骤：S1：获取数据：制作棋盘格标定板，并测量各方格的尺寸，从不同角度拍摄两张以上照片；S2：计算单应矩阵；S3：计算内参数、外参数和畸变因数的初始值；S4：用非线性优化方法精化计算结果。2.根据权利要求1所述的一种基于非线性优化的相机标定方法，其特征在于，所述步骤S3具体包括：采用无畸变情况下的相机模型，如下式：其中，f表示焦距；uo、vo表示主点在图像中的坐标；r11～33表示旋转矩阵的各元素；t1～3表示平移向量的各元素；x、y表示三维空间中目标点的坐标；u、v表示目标点的像的图像坐标；zc表示相机坐标系中目标点的Z坐标；已知P，计算下式等号右侧中的参数：其中，P是3X4的矩阵，只取其中的第1、2、4列，即删除第3列后，得到另外一个3X3的矩阵H，其中的第i个列向量表示为：hi＝[hi1hi2hi3]T；(1)计算内参数：f和uo、vo求解如下方程组：即Vb＝0其中，vij＝[ji3hj3-(ji3hj1+h...

【专利技术属性】
技术研发人员：朱红军，钱鹰，刘歆，聂顺超，夏璨，杨鹏程，张俊，
申请(专利权)人：重庆邮电大学，
类型：发明
国别省市：重庆,50