一种基于二维标定板的三维标定方法技术

技术编号:35453506 阅读:21 留言:0更新日期:2022-11-03 12:09
本发明专利技术公开了一种基于二维标定板的三维标定方法,涉及计算机视觉技术领域,包括:保持标定板位置不变,通过相机获取m幅不同位置、不同角度的标定板图像;对获取的各个标定板图像的所有角点进行检测;获取各个标定板图像的所有角点的世界坐标;对各个标定板图像的所有角点的像素坐标以及获取的各个标定板图像的所有角点的世界坐标进行相关运算获取相机的内部参数矩阵;根据获取的相机内部参数矩阵和求得的单应矩阵分离出相机的外部参数;将各个标定板图像的所有角点的像素坐标和获取的相机的内部参数待入至畸变方程中,获取畸变系数;通过构建似然函数对畸变系数进行优化,本发明专利技术采用了二维标定板进行三维标定,在不失标定板便携易维护的情况下,保留了标定板的维度信息,从而提高了标定算法的精度。从而提高了标定算法的精度。从而提高了标定算法的精度。

【技术实现步骤摘要】
一种基于二维标定板的三维标定方法


[0001]本专利技术涉及计算机视觉
,尤其涉及一种基于二维标定板的三维标定方法。

技术介绍

[0002]随着计算机图像处理和计算机视觉的发展,为适应不同的场景与产品的需求,工业界与科研工作者提出了许多相机标定技术,相机标定是机器视觉领域的重要基础,它在空间物体表面某点的三维几何坐标与其在二维图像的像素点之间建立关系,根据相机成像模型求解相机模型的内部参数和外部参数,从而为视觉计算、三维重建、无人驾驶等其它计算机视觉领域提供基础运算,在某些领域,相机标定的精度对其后续工作结果的精度有着直接影响,因此完成相机标定是后续工作的前提,提高相机标定精度是一项重要的科研议题。
[0003]常见的相机标定算法有传统三维标定法、张氏标定法、相机自标定法等,其中相机自标定法无需使用尺寸已知的标定物,而是利用相机运动的约束以求解相机内外参数,但在实际应用中由于相机的运动约束条件太强,算法并不具有强鲁棒性;传统三维标定板标定法和张氏标定法都需要特定的标定板,传统三维标定法需要一块精度极高的立体标定板,对制作工艺有较高的要求,且维护较为困难;张氏标定法克服了传统标定法需要高精度标定物的缺点,仅需一个打印的平面棋盘格作为标定物,操作简单易于维护,但该方法中采用的标定角点均位于Z=0的平面,没有充分利用标定板的三维信息,进而有时达不到某些应用需求的精度。
[0004]因此,提供一种新的技术方案改善上述问题,是本领域技术人员急需解决的问题。

技术实现思路

[0005]有鉴于此,本专利技术提供一种基于二维标定板的三维标定方法,以解决上述技术问题。
[0006]为实现上述目的,本专利技术提供如下技术方案:
[0007]一种基于二维标定板的三维标定方法,其特征在于,包括:
[0008]保持标定板位置不变,通过相机获取m幅不同位置、不同角度的标定板图像;
[0009]对获取的各个标定板图像的所有角点进行检测;
[0010]获取各个标定板图像的所有角点的世界坐标;
[0011]对各个标定板图像的所有角点的像素坐标以及获取的各个标定板图像的所有角点的世界坐标进行相关运算获取相机的内部参数矩阵;
[0012]根据获取的相机内部参数矩阵和求得的单应矩阵分离出相机的外部参数;
[0013]将各个标定板图像的所有角点的像素坐标和获取的相机的内部参数待入至畸变方程中,获取畸变系数;
[0014]通过构建似然函数对畸变系数进行优化。
[0015]在上述的方案中,所述获取各个标定板图像的所有角点的世界坐标包括:
[0016]定义标定板图像中心角点为世界坐标系原点;
[0017]获取标定板图像在世界坐标系上的倾斜角度;
[0018]根据世界坐标系原点以及标定板图像在世界坐标系上的倾斜角度获取各个标定板图像的所有角点的世界坐标。
[0019]在上述的方案中,所述对各个标定板图像的所有角点的像素坐标以及获取的各个标定板图像的所有角点的世界坐标进行相关运算获取相机的内部参数矩阵包括:
[0020]根据映射公式建立相机空间点到像素点的映射关系;
[0021]定义并求得单应矩阵;
[0022]根据求得的单应矩阵获取相机的内部参数矩阵。
[0023]在上述的方案中,所述映射公式为:其中,s为固定实数,为单个标定板图像的所有角点的像素坐标矩阵,为单个标定板图像的所有角点的世界坐标矩阵,
[0024]为相机的内参矩阵,为相机的外参矩阵,R由三个互相垂直的单位方向矢量构成。
[0025]在上述的方案中,所述定义并求得单应矩阵包括:
[0026]定义单应矩阵
[0027]将单个标定板图像的所有角点的像素坐标矩阵和单个标定板图像的所有角点的世界坐标矩阵待入至映射公式中获取单应矩阵中各项元素的数值,并令相机的内参矩阵中γ=0,求得单应矩阵
[0028]取u0=单个标定板图像高度的一半,v0=单个标定板图像高度的一半,对求得的单应矩阵矩阵进行多次初等行变换,得到:
[0029][0030]在上述的方案中,所述根据求得的单应矩阵获取相机的内部参数矩阵包括:
[0031]重新定义单应矩阵:
[0032]其中,
h0=αr
11
,v0=αr
21
,k0=αr
31
,j0=αt
13
,h1=βr
12
,v1=βr
22
,k1=βr
32
,j1=βt2,h2=r
13
,v2=r
23
,k2=r
33
,j2=t3;
[0033]对重新定义的单应矩阵中的部分列向量进行相关运算:
[0034]令令
[0035]对重新定义的单应矩阵中的部分元素以及对重新定义的单应矩阵中的部分列向量进行相关运算获取的各个向量中的部分元素进行运算以及组合,获取新的矩阵,并将获取的新矩阵与对相机内参矩阵中的α和β分别进行相关运算后组合成的矩阵进行乘积运算:
[0036]令新矩阵为
[0037]对相机内参矩阵中的α和β分别进行相关运算后组合成的矩阵为
[0038]将获取的新矩阵与对相机内参矩阵中的α和β分别进行相关运算后组合成的矩阵进行乘积运算得到
[0039][0040]在上述的方案中,所述根据求得的单应矩阵获取相机的内部参数矩阵还包括:
[0041]根据求得的单应矩中的元素对新矩阵中各个元素和乘积运算得到的矩阵中各个元素进行赋值:
[0042][0043]根据R由三个互相垂直的单位方向矢量构成,得到矩阵R中元素的相关关系式:
[0044]根据矩阵R中元素的相关关系式对赋值处理的乘积运算得到的矩阵进行再次赋值:
[0045][0046]根据赋值结果求解相机内参矩阵中的α和β。
[0047]在上述的方案中,所述根据获取的相机内部参数矩阵和求得的单应矩阵分离出相机的外部参数包括:
[0048]根据单应矩阵H=A[R t],得到[R t]=A
‑1H;
[0049]将得到的相机内参矩阵A进行矩阵求逆运算,得到A
‑1;
[0050]将A
‑1以及求得的单应矩阵获取相机的外部参数矩阵[R t]。
[0051]在上述的方案中,所述畸变方程为:
[0052][0053]其中,k1与k2为畸变系数,(x,y)为理想无畸变点,为实际畸变点。
[0054]在上述的方案中,所述通过构建似然函数对畸变系数进行优化包括:
[0055]构建似然函数其中,m
ij
为单个标定板图像中其中一个角点的高斯噪声,为角点预设高斯噪声;
[0056]求得似然函数的最小值,获取畸变参数优化值。
[0057]综上所述,本专利技术的有益效果是:采用了二维标定板进行本文档来自技高网
...

【技术保护点】

【技术特征摘要】
1.一种基于二维标定板的三维标定方法,其特征在于,包括:保持标定板位置不变,通过相机获取m幅不同位置、不同角度的标定板图像;对获取的各个标定板图像的所有角点进行检测;获取各个标定板图像的所有角点的世界坐标;对各个标定板图像的所有角点的像素坐标以及获取的各个标定板图像的所有角点的世界坐标进行相关运算获取相机的内部参数矩阵;根据获取的相机内部参数矩阵和求得的单应矩阵分离出相机的外部参数;将各个标定板图像的所有角点的像素坐标和获取的相机的内部参数待入至畸变方程中,获取畸变系数;通过构建似然函数对畸变系数进行优化。2.根据权利要求1所述的基于二维标定板的三维标定方法,其特征在于,所述获取各个标定板图像的所有角点的世界坐标包括:定义标定板图像中心角点为世界坐标系原点;获取标定板图像在世界坐标系上的倾斜角度;根据世界坐标系原点以及标定板图像在世界坐标系上的倾斜角度获取各个标定板图像的所有角点的世界坐标。3.根据权利要求1所述的基于二维标定板的三维标定方法,其特征在于,所述对各个标定板图像的所有角点的像素坐标以及获取的各个标定板图像的所有角点的世界坐标进行相关运算获取相机的内部参数矩阵包括:根据映射公式建立相机空间点到像素点的映射关系;定义并求得单应矩阵;根据求得的单应矩阵获取相机的内部参数矩阵。4.根据权利要求3所述的基于二维标定板的三维标定方法,其特征在于,所述映射公式为:其中,s为固定实数,为单个标定板图像的所有角点的像素坐标矩阵,为单个标定板图像的所有角点的世界坐标矩阵,为相机的内参矩阵,为相机的外参矩阵,R由三个互相垂直的单位方向矢量构成。5.根据权利要求3所述的基于二维标定板的三维标定方法,其特征在于,所述定义并求得单应矩阵包括:定义单应矩阵将单个标定板图像的所有角点的像素坐标矩阵和单个标定板图像的所有角点的世界坐标矩阵待入至映射公式中获取单应矩阵中各项元素的数值,并令相机的内参矩阵中γ=0,求得单应矩阵
取u0=单个标定板图像高度的一半,v0=单个标定板图像高度的一半,对求得的单应矩阵矩阵进行多次初等行变换,得到:6.根据权利要求3所述的基于二维标定板的三维标定方法,其特征在于,所述根据求得的单应矩阵获取相机的内部参数矩阵包括:重新定义单应矩阵:其中,其中,h0=αr
11
,v0=αr
21<...

【专利技术属性】
技术研发人员:李毅林绿开黄辉钱乐旦张笑钦
申请(专利权)人:温州大学大数据与信息技术研究院
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1