双目立体视觉摄像机标定参数的整体优化方法
【专利摘要】本发明公开了一种双目立体视觉摄像机标定参数的整体优化方法，是一种通过迭代极小化重投影像点与理想像点偏差来优化双目视觉摄像机参数的方法其针对双目立体视觉测量系统，利用两台摄像机同时拍摄二维棋盘格标定板的图像，提出了一种将所有参数进行整体优化的优化方法；视标定板上所有特征点为共面特征点，以位于标定板左上角的特征点为原点建立标定板坐标系，标定板上所有特征点在标定板坐标系下已知，通过坐标变换分别得到标定板坐标系同左右摄像机坐标系的旋转矩阵和平移向量，再通过迭代极小化重投影像点与理想像点偏差来优化双目视觉摄像机参数。
【专利说明】双目立体视觉摄像机标定参数的整体优化方法
【技术领域】
[0001]本发明属于计算机视觉检测领域，特别涉及一种双目立体视觉摄像机标定参数的整体优化方法。
【背景技术】
[0002]双目立体视觉是根据双目摄像机同时拍摄的两张图片的二维图像信息得到物体的几何三维信息。从二维信息恢复三维信息依赖于双目视觉系统的模型，最重要的是需要确定摄像机成像几何模型，该几何模型的参数称为摄像机参数。
[0003]摄像机内参数为摄像机与几何和光学特性有关的参数，内参数只依赖于摄像机自身的制造特点，外参数为摄像机坐标系相对于建立好的世界坐标系的旋转和平移。摄像机标定过程就是确定摄像机内、外参数的过程，标定方法的精度直接影响到计算机视觉测量的精度。因此，对摄像机进行快速、简捷、精准标定的研究无疑具有重大的意义。
[0004]目前摄像机标定方法可以分为三类:传统标定方法、自标定法和基于主动视觉的标定方法。传统标定方法标定精度高，但标定过程繁琐、鲁棒性差。自标定法效率高，但依赖严重特征点提取，精度无法保证。基于主动视觉标定方法精度较高，适应性强，能应用于复杂工况，是进行标定较好选择。
[0005]目前对标定的优化方法大多为复杂非线性算法，引入很多参数作为迭代变量，需要大量计算，若模型设计不当，更容易加大计算量，同时导致结果并不稳定，容易导致结果不收敛，无法得到正确结果。
【发明内容】
[0006]本发明主要解决的是各参数在摄像机标定过程中的误差问题，针对双目立体视觉测量系统，利用两台摄像机同时拍摄二维棋盘格标定板的图像，提出了一种将所有参数进行整体优化的优化方法；视标定板上所有特征点为共面特征点，以位于标定板左上角的特征点为原点建立标定板坐标系，标定板上所有特征点在标定板坐标系下已知，通过坐标变换分别得到标定板坐标系同左右摄像机坐标系的旋转矩阵和平移向量，再通过得到的旋转矩阵及平移向量优化摄像机内、外参数。
[0007]本发明通过迭代极小化重投影像点与理想像点偏差来优化双目视觉摄像机参数，利用双目摄像机同时拍摄二维棋盘格标定板，根据标定板和摄像机的坐标对应关系分别得到左、右摄像机坐标系和标定板坐标系的旋转矩阵及平移向量，再通过迭代极小化重投影像点与理想像点偏差来优化双目视觉摄像机参数；具体步骤如下:
[0008]步骤1:拍摄棋盘格标定板并建立标定板坐标系
[0009]利用两台摄像机同时拍摄棋盘格标定板，得到两张标定板图像；将标定板上所有特征点视为共面特征点，以左上角的特征点为坐标原点，以标定板平面为XOY平面，建立空间立体坐标系，称为标定板坐标系；由于标定板方格实际尺寸已知，所以各个特征点在标定板坐标系下的X向坐标和Y向坐标已知，Z向坐标均为O。
[0010]步骤2:建立两摄像机坐标系与标定板坐标系之间的外参数估计模型
[0011]利用初始摄像机参数以及上一步骤中得到的两张标定板图像可以重建出标定板上特征点在世界坐标系下的三维坐标，一般将世界坐标系与双目摄像机其中一个摄像机坐标系重合，所以重建的三维坐标可以视作该选定摄像机坐标系下的三维坐标，利用初始外参数可以得到重建点在另一摄像机下的三维坐标，此时得到了同一组特征点在左摄像机坐标系、右摄像机坐标系和标定板坐标系下的坐标。
[0012]利用部分特征点来求取标定板坐标系到左摄像机坐标系的旋转矩阵&和平移向量标定板坐标系到右摄像机坐标系的旋转矩阵Rk和平移向量tK。
[0013]选取至少三个不共线的特征点，假设选取的特征点在标定板坐标系和摄像机坐标系中对应的点集为KPpP/ )，(P2，PZ(pn, Pni )}，其中PiS特征点在标定板坐标系坐标,Pi'为特征点在摄像机坐标系坐标,i = I, 2，3…η,求取旋转矩阵Rci和平移向量tQ,点集中对应空间点的关系为Pi' = RciPJtci，建立如下模型估计参数(Rc^tci):
[0014]
【权利要求】
1.一种双目立体视觉摄像机标定参数的整体优化方法，是一种通过迭代极小化重投影像点与理想像点偏差来优化双目视觉摄像机参数的方法，其特征在于如下步骤: 步骤1:拍摄棋盘格标定板并建立标定板坐标系 利用两台摄像机同时拍摄棋盘格标定板，得到两张标定板图像；将标定板上所有特征点视为共面特征点，以左上角的特征点为坐标原点，以标定板平面为XOY平面，建立空间立体坐标系，称为标定板坐标系；由于标定板方格实际尺寸已知，所以各个特征点在标定板坐标系下的X向坐标和Y向坐标已知，Z向坐标均为O ; 步骤2:建立两摄像机坐标系与标定板坐标系之间的外参数估计模型 利用初始摄像机参数以及上一步骤中得到的两张标定板图像重建出标定板上特征点在世界坐标系下的三维坐标，将世界坐标系与双目摄像机其中一个摄像机坐标系重合，重建的三维坐标视作该选定摄像机坐标系下的三维坐标；利用初始外参数得到重建点在另一摄像机下的三维坐标，得到了同一组特征点在左摄像机坐标系、右摄像机坐标系和标定板坐标系下的坐标； 利用部分特征点来求取标定板坐标系到左摄像机坐标系的旋转矩阵&和平移向量标定板坐标系到右摄像机坐标系的旋转矩阵Rk和平移向量tK ； 选取至少三个不共线的特征点，假设选取的特征点在标定板坐标系和摄像机坐标系中对应的点集为{(P1, P/ )，(P2, )，...，(Pn, P1/ )}，其中Pi为特征点在标定板坐标系坐标，Pi'为特征点在摄像机坐标系坐标，1 = 1,2, 3…η,求取旋转矩阵Rci和平移向量点集中对应空间点的关系为Pi' = RciPftci，建立如下模型估计参数(Rc^tci)；
步骤3:求取旋转矩阵Rtl和平移向量h 首先求取点集质心所在坐标，具体计算为:
其中户，尹为各特征点各轴分量的均值，Pi为特征点在标定板坐标系坐标，P/为特征点在摄像机坐标系坐标，1 = 1,2, 3…η ； 将点集质心移动至坐标系原点处，平移后的空间点为:
其中月为在标定板坐标系下质心移至坐标原点后的特征点坐标，P力在摄像机坐标系下质心移至坐标原点后的特征点坐标； 由公式(I) (2) (3)，得到
，目标函数简化为只估计旋转矩阵，具体形式如下:
得到旋转矩阵之后，利用质心坐标的求出平移向量如下:
步骤4:建立优化模型优化标定参数 利用摄像机的内部参数、摄像机坐标系与标定板坐标系的旋转矩阵和平移向量求解标定板上除上一步骤选取的特征点之外的所有特征点重投影坐标ρ.,.,?))，具体算法如下:
其中，rjj为旋转矩阵R的第i行、第j列上的元素,平移向量t = Ct1, t2, t3)T, fx为摄像机横向尺度因子，fy为摄像机纵向尺度因子，U0为主点在像素坐标系下的横坐标，V0为主点在像素坐标系下的纵坐标，(Xff, Yff, Zff)为特征点在世界坐标系下的坐标； 根据已知畸变系数，将实际拍摄获得的像点坐标(Ui' ,V/ )校正为相应的理想像点坐标(Ui, Vi);建立优化模型通过迭代极小化重投影像点坐标和理想像点坐标的偏差，目标优化函数为:
采用LM非线性优化算法，将Hessian阵变为两台的对称正定阵，进行求解，当偏差最小时对应的参数即为优化后的双目立体视觉系统摄像机参数。
【文档编号】G06T7/00GK104182982SQ201410427812
【公开日】2014年12月3日 申请日期:2014年8月27日 优先权日:2014年8月27日
【发明者】贾振元, 刘巍, 王亮, 杨景豪, 刘阳, 王灵丽, 张弛 申请人:大连理工大学