基于分块的摄像机标定方法
【专利摘要】本发明基于分块的摄像机标定方法属于图像处理和计算机视觉检测领域，特别涉及大型锻件尺寸测量系统中摄像机内、外参数的现场标定方法。标定方法首先利用变焦的方法计算出摄像机主点坐标，在不同焦距下拍摄标定模板的7幅图像，通过缩比不变性求得摄像机主点坐标，而后将大视场划分成n×n个子空间，分别对每个区域进行标定；在测量所需的焦距下拍摄标定模不同方位的图像，然后对模板图像进行直线提取；利用所得直线求取对应正交消隐点，即可求得摄像机内参数；最后，再标定双目视觉测量系统的外参数。本发明解决了摄像机图像不同区域标定与校正参数不一致的难题，有效提高了标定的精度，特别适用于锻造现场的准确标定。
【专利说明】基于分块的摄像机标定方法
【技术领域】
[0001] 本发明属于图像处理和计算机视觉检测领域，特别涉及大型锻件尺寸测量系统中 摄像机内、外参数的现场标定方法。
【背景技术】
[0002] 计算机视觉处理的基本任务之一是根据二维图像信息恢复物体的三维几何信息。 要实现利用图像点求取相应的空间物体表面点的任务，需要确定摄像机参数，即对摄像机 进行标定。摄像机标定过程是建立摄像机所拍摄场景中物体上特征点的图像坐标与物体上 特征点世界坐标的关系。测量系统中需标定的摄像机参数分为摄像机的畸变等内参数和被 测平面相对于摄像机坐标系的外参数，其中内参数包括摄像机的主点坐标以及尺度因子， 外参数包括两摄像机间的旋转矩阵和平移矢量。标定时，由场景中已知特征点的图像坐标 和世界坐标即可求解摄像机的内外参数。摄像机标定是构建测量系统地一项关键技术，标 定结果对系统的测量准确度及效果有至关重要的作用。
[0003] 摄像机的标定方法较多，Tsai提出的基于径向约束的两步法及张正友提出的基于 棋盘格靶标标定方法均可以获得较高的精度。但是对于大视场下的摄像机标定而言，由于 传统标定方法所占有视场较小，导致特征点提取误差大，通过标定获得的摄像机内参数存 在误差，并对外参数的标定精度产生影响。在实际的成像系统中，由于系统复杂非线性的影 响，在视场的不同区域，畸变参数并不一致，用单一的一组内外参数表示被测平面相对于摄 像机的姿态进行平面图像测量，获得的测量结果精度较差。
【发明内容】
[0004] 本发明所要解决的技术问题是克服现有技术的不足，针对在锻造现场传动标定方 法精度较差的问题，发明一种基于分块的摄像机标定方法。根据实际成像系统非线性畸变 的特点，将整个视场划分为ηΧη个子空间，该方法首先求取主点坐标，利用分块的思想将 视场划分成ηΧη个子空间，而后分别对ηΧη个子空间进行标定：利用标定模板特征点计 算图像的正交消隐点，并计算相机的尺度因子，最后，以基本矩阵和本质矩阵分别得出ηΧη 个子空间的参数。在每个子空间上利用基于主动视觉系统的摄像机标定方法分别进行其余 内参数和外参数的标定，能较好地补偿实际成像系统中存在的误差，并取得较高的图像测 量精度。所以，这种标定方法有效提1? 了标定的精度。
[0005] 本发明采用的技术方案是一种基于分块的摄像机标定方法，其特征是，标定方法 首先利用变焦的方法计算出摄像机主点坐标，在不同焦距下拍摄标定模板的7幅图像，通 过缩比不变性求得摄像机主点坐标，而后将大视场划分成ηΧη个子空间，分别对每个区域 进行标定，具体标定方法如下：在测量所需的焦距下拍摄标定模不同方位的图像，然后对模 板图像进行直线提取；利用所得直线求取对应正交消隐点，即可求得摄像机内参数；最后， 同时标定双目视觉测量系统的外参数，具体步骤如下：
[0006] 步骤1 :求取摄像机主点坐标
[0007] 摄像机标定实际上是要求出内参数矩阵K，两摄像机间的旋转矩阵R和平移矢量 t。通用的摄像机模型可表示为：
[0008]
【权利要求】
1. 一种基于分块的摄像机标定方法，其特征是，标定方法首先利用变焦的方法计算出 摄像机主点坐标，在不同焦距下拍摄标定模板的7幅图像，通过缩比不变性求得摄像机主 点坐标，而后将大视场划分成nXn个子空间，分别对每个区域进行标定；在测量所需的焦 距下拍摄标定模不同方位的图像，然后对模板图像进行直线提取；利用所得直线求取对应 正交消隐点，即可求得摄像机内参数；最后，再标定双目视觉测量系统的外参数，标定方法 的具体步骤如下： 步骤1 :求取摄像机主点坐标 摄像机标定实际上是要求出内参数矩阵K，两摄像机间的旋转矩阵R和平移矢量t;通 用的摄像机模型可表示为：
其中，（Xw，Yw，Zw，1)T为空间点在世界坐标系下的齐次坐标，（x，y，1)T为空间点在图像坐 标系下的对应像点坐标，fx，fy为摄像机尺度因子，为像机的主点坐标，z为比例因 子； (1) 采集变焦距图像 在待测目标的同一侧布置两台摄像机构成双目视觉测量系统，按照摄像的位置关系， 分别为左摄像机（2a)和右摄像机（2b);标定左摄像机（2a)时，调整标定模板的位置，保证 标定模板在不同拍摄焦距时，均在摄像机视场内，在最大焦距下，对摄像机合焦；标定模板 静止，逐渐减小摄像机焦距，在不同焦距下拍摄7幅模板图像；使用同样的方法采集右摄像 机（2b)的7张图像； (2) 分别计算左、右摄像机主点坐标 改变焦距时，摄像机视场呈现比例扩缩现象，因此在不同焦距拍摄的图像中，特征点连 线的交点即为主点；计算摄像机主点坐标，利用拍摄的七张图片，将不同图像上标定模板相 同特征点对应的7个像点采用线性回归的方法拟合成n条直线，第i条直线的拟合方程为： Ii:ajU+biV+Cj= 0,(i=I, 2. . .n) (2) 以到像点连线距离平方和最小为目标函数S，其中，
计算主点坐标（uQ,V。），分别求出左、右摄像机的主点坐标（u1Q,v1(l)和（U1^vki); 步骤2 :大视场子空间划分 将双目测量系统中的左、右摄像机合理放置，调整摄像机的视场，使左、右摄像机视场 重叠，并根据测量距离选择合适的焦距，进行合焦操作；按照分块的思想，将大视场划分 成nXn个等分，得到nXn个子空间，按照行列子空间分别编号为\_，（i=l，2...n，j= I, 2.. .n), 其中，i为子空间所在的行，j为子空间所在的列； 步骤3:单个子空间标定 根据变焦距计算的主点坐标，分别在各个子空间标定出尺度因子和摄像机的外部参 数； (1) 求取摄像机尺度因子 利用已经得出的主点坐标，基于正交消隐点计算摄像机的尺度因子匕和&，以左像机 为例，说明标定的具体过程： (a) 图像采集，分别在每个子空间Sij,(i= 1，2. . .n,j= 1，2. . .n)区域内，拍摄标定 模板在不同位姿下的图像10张，作为求取摄像机尺度因子的图像； (b) 根据拍摄得到标定模板的不同方位的图像，分别求出左摄像机的10组正交消隐点 对Pi (uPi，vPi)，Qi (uQi，vQi)，i= 1，2. ? ? 10 ; (c) 将10组正交消隐点对分别代入非线性方程：
其中，（uPi，vPi)、（uQi，vQi)分别为第i幅图像的正交消隐点对，联立方程组，利用非线性 最小二乘方法迭代求解非线性方程组获得尺度因子，即fy，fx ;其中，（uw，Vltl)为左摄像机的 主点坐标，fy，fx为摄像机等效焦距，由此可得左摄像机内参数矩阵
采用相同的步骤（a)、（b)和（c)可得到右摄像机的尺度因子，由步骤1中的主点坐标 构成右摄像机的内参数矩阵Kk ; (2) 求取摄像机的外参数 利用基本矩阵求取摄像机的外参数，首先在子空间区域摆放两次标定模板，并拍摄图 像，摄像机间外部参数的标定的具体步骤如下： (a) 求取基本矩阵F，根据拍摄的两张图像，提取对应像点坐标，利用对应像点的齐次 坐标获得线性方程： [uu' ，vu' ,U1 ，uv' ,W1 ,V1 ,u,v, 1]f= 0 (6) 其中，（u，v)，（u'，v')分别为左右摄像机图像中对应的像点坐标，为9维列向量， 为f= [fn，f12,f13,f21，f22,f23,f31，f32,f33]T基本矩阵F的第i行、第j列上的元素；提取所 采集的两幅图像中n(n彡8)组对应点构成对应点矩阵A:
其中，（Ui，Vi)，（u/，v/ )为两幅图像中对应点坐标；利用Af= 0,即可求得基本矩 阵F; (b) 根据左、右摄像机的内参数矩阵KuKk，由公式⑶即可得到本质矩阵E; E=KetFKl (8) (c) 本质矩阵由双目立体视觉的内外参数决定，对公式⑶中得到的本质矩阵E进行奇 异值分解：
其中，OO2,O3为特征值且Oi彡O2彡O3 ;从本质矩阵算出摄像机间的外部参数OIG，其中，f=M,Z为相差的比例因子，t为平移矢量； 外部参数具有下述4种组合： [UWVtIu3] ； [UWVtI-U3]； [UWtVt |u3]； [UWtVt |-u3] (10) 其中，
》U和V为分解出的矩阵，U3为U的第3列向量，|f| =l,对以 上4种参数组合进行判断，方法如下： 由于重建点应在两台摄像机前方，因此将世界坐标系分别与左、右摄像机坐标系重合， 重建出的点应取三维坐标中Z向分量均应为正值的组合；由上述参数组合的形式判断,计 算出旋转矩阵R和单位平移矢量？； (d)计算摄像机的比例因子，采用长度为L的高精度标尺作为长度基准，使用上述双目 测量系统对其进行拍摄并重建出三维长度1，则摄像机的比例因子
；根据f=zr确定 平移矢量t，就得到左右摄像机的外参数； 通过上述步骤计算出双目视觉测量系统的内外参数，完成对单个子空间的标定过程； 步骤4 :大视场整体标定 按照步骤3中的标定方法，分别对nXn个子空间进行标定，确定每个分块子空间的内 外参数，将所得的内外参数对应其所在子空间进行记录；在双目视觉测量的后续工作中，根 据拍摄特征所在大视场中子空间的位置对内外参数进行调用。
【文档编号】G06T7/00GK104268876SQ201410505306
【公开日】2015年1月7日 申请日期:2014年9月26日 优先权日:2014年9月26日
【发明者】贾振元, 刘巍, 杨景豪, 樊超楠, 刘阳, 王灵丽 申请人:大连理工大学