CN114359405A - 一种离轴沙姆3d线激光相机的标定方法 - Google Patents

一种离轴沙姆3d线激光相机的标定方法 Download PDF

Info

Publication number
CN114359405A
CN114359405A CN202111600729.9A CN202111600729A CN114359405A CN 114359405 A CN114359405 A CN 114359405A CN 202111600729 A CN202111600729 A CN 202111600729A CN 114359405 A CN114359405 A CN 114359405A
Authority
CN
China
Prior art keywords
coordinate system
camera
laser
calibration
coordinates
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111600729.9A
Other languages
English (en)
Inventor
翟晓彤
彭思龙
汪雪林
顾庆毅
周大勇
王克岩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhongke Xingzhi Jinan Intelligent Technology Co ltd
Suzhou Zhongke Xingzhi Intelligent Technology Co ltd
Original Assignee
Zhongke Xingzhi Jinan Intelligent Technology Co ltd
Suzhou Zhongke Xingzhi Intelligent Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhongke Xingzhi Jinan Intelligent Technology Co ltd, Suzhou Zhongke Xingzhi Intelligent Technology Co ltd filed Critical Zhongke Xingzhi Jinan Intelligent Technology Co ltd
Priority to CN202111600729.9A priority Critical patent/CN114359405A/zh
Publication of CN114359405A publication Critical patent/CN114359405A/zh
Pending legal-status Critical Current

Links

Images

Abstract

本发明公开了一种离轴沙姆3D线激光相机的标定方法,构建离轴沙姆3D线激光相机的成像模型:构建与透镜面平行且成像中心过光轴的虚拟像面坐标系,实现图像像素坐标与世界坐标系下坐标之间的转换;先计算相机内参和外参初值,后利用非线性优化算法计算最优解,获得相机标定;采集不少于2组有无激光用于激光平面标定的标定板图片,通过组内标定板图片做差的方法,得到激光线数据,计算每组激光线中心像素坐标在相机坐标系下的三维坐标,拟合得到激光平面方程,得到激光平面的标定。解决了离轴沙姆3D线激光相机难以标定的问题。

Description

一种离轴沙姆3D线激光相机的标定方法
技术领域
本发明涉及3D线激光测量技术领域,尤其涉及一种离轴沙姆3D线激光相机的标定方法。
背景技术
3D线激光相机主要由相机和激光器组成,是一种快速发展的非接触式测量装置,具有灵活性好、速度快、精度高以及智能化等优点,在逆向工程和复杂曲面零部件等测量领域有着广泛应用。但在实际应用中,对于倾斜目标,传统镜头存在由于景深限制,会出现视野两端离焦模糊的问题。为解决此问题,具有增大景深,使得倾斜目标全视野清晰成像优势的沙姆相机,逐渐得到了重视与青睐。如图1所示,沙姆相机是指拍摄时满足沙姆定律的相机,即目标聚焦面、透镜面、成像面三者延长线相交于一线,且相交线唯一,此相交线即称为沙姆线。离轴是指主光轴不通过成像面中心,具有降低主光线角度CRA,增加相机接收强度的作用。
相机标定是求解相机参数的过程,即通过建立相机模型,构建相机图像像素位置与场景点位置之间的关系,由已知特征点的图像坐标求解相机的模型参数。相机标定的精度及算法的稳定性直接影响相机工作产生结果的准确性,做好相机标定是做好后续工作的前提,提高标定精度是实现相机测量的重点所在。
目前,对于普通相机标定问题已提出了很多成熟的方法,如传统相机标定法、主动视觉相机标定方法、相机自标定法等。但对于离轴沙姆3D线激光相机,其与普通相机的成像模型不同,因此传统成熟的相机标定方法对离轴沙姆3D线激光相机标定无法适用。
发明内容
本发明实施例提供了一种离轴沙姆3D线激光相机的标定方法,用于实现离轴沙姆3D线激光相机的标定。
一种离轴沙姆3D线激光相机的标定方法,包括:
步骤1:构建离轴沙姆3D线激光相机的成像模型:将图像坐标系转换到实际像面坐标系,构建与透镜面平行且成像中心过光轴的虚拟像面坐标系,将实际像面坐标系转换到虚拟像面坐标系,在虚拟像面坐标系进行畸变校正,得到虚拟像面坐标系下无失真的虚拟像面坐标,再将无失真的虚拟像面坐标转换到实际像面坐标系下无失真的实际像面坐标,再将实际像面坐标转换到图像坐标系对应的无失真的图像像素坐标,再近似按照小孔成像模型,将无失真的图像像素坐标转换到相机坐标系下的坐标,再将相机坐标系下坐标转换到世界坐标系下的坐标。从而得到图像像素坐标与世界坐标系下坐标的转换关系,实现离轴沙姆3D线激光相机成像模型的构建;
步骤2:基于成像模型进行相机标定:先计算相机内参和外参初值,后利用非线性优化算法计算最优解,进行相机标定获得相机标定参数:相机内参、畸变系数、沙姆角、偏移量;
步骤3:对激光平面进行标定:采集不少于2组用于激光平面标定的标定板图片,每组中标定板位于同一位置,包含有激光线图片和无激光线图片;通过组内标定板图片做差的方法,得到激光线数据,通过中心线提取算法提取激光线中心像素坐标;对每组无激光线标定板图片,提取特征点坐标,以距离激光线最邻近的特征点计算标定板图片对应的外参,根据标定板图片对应的外参,计算每组激光线中心像素坐标在相机坐标系下的三维坐标,通过拟合得到激光平面方程,实现对激光平面的标定。
一种较佳的实施方式,步骤1具体为:
(1)图像坐标系下的一点qi=(ui,vi)T转换到实际像面坐标系下的点qs=(xs,ys)T,实际像面坐标系为倾斜的像面;
Figure BDA0003431654870000031
其中,sx、sy分别表示像素在x方向和y方向的实际物理尺寸,(u0,v0)表示图像坐标系原点在图像坐标系中的坐标;
(2)将实际像面坐标系下的点qs转换到虚拟像面坐标系下的点qp=(xp,yp)T,虚拟像面垂直于光轴且成像中心在光轴上;
Figure BDA0003431654870000032
其中,ox、oy分别表示实际像面坐标系中心在x方向和y方向的偏移量,Sx、Sy分别表示倾斜的实际像面中心与虚拟像面的x方向和y方向的沙姆角;
(3)在虚拟像面坐标系上引入畸变模型,将虚拟像面坐标系下的点qp变换到无失真的点
Figure BDA0003431654870000033
Figure BDA0003431654870000034
其中,k1、k2、k3表示径向畸变参数,p1、p2表示切向畸变参数,r2=xp 2+yp 2
(4)将虚拟像面坐标系上无失真的点
Figure BDA0003431654870000035
转换到实际像面坐标系上无失真的点
Figure BDA0003431654870000041
Figure BDA0003431654870000042
(5)将实际像面坐标系上无失真的点
Figure BDA0003431654870000043
转换到图像坐标系下无失真的像素点
Figure BDA0003431654870000044
Figure BDA0003431654870000045
(6)将实际像面坐标系上无失真的像素点
Figure BDA0003431654870000046
转换到相机坐标系下的点qc=(xc,yc,zc)T
Figure BDA0003431654870000047
其中,αx、αy表示x和y方向的有效焦距;
(7)将相机坐标系下的点qc转换为世界坐标系下的点qw=(xw,yw,zw)T
Figure BDA0003431654870000048
其中,R表示由世界坐标系到相机坐标系的旋转矩阵,T表示由世界坐标系到相机坐标系的平移向量。
一种较佳的实施方式,步骤2中先计算相机内参和外参初值,具体为:通过张氏解析算法计算相机的内参和外参初值。
一种较佳的实施方式,步骤2中利用非线性优化算法计算最优解,进行相机标定,具体为:建立内参、畸变参数、外参以及沙姆角和偏移量的参数的非线性优化方程,通过迭代循环求解进行优化,最终的最优解即为相机标定参数。
一种较佳的实施方式,步骤2具体为:
(1)对不少于N1张的标定板图片,提取每张标定板图片上的N2个特征点;
(2)采用张氏解析算法,计算相机内参αx,αy,u0,v0,以及外参Ri,Ti
(3)利用非线性优化算法对内参、畸变系数、外参以及沙姆角和偏移量进行优化,对应的非线性优化模型如下:
Figure BDA0003431654870000051
L为以重投影误差最小为准则的目标方程,i表示标定板第i个位姿,j表示第j个特征点,N1为标定板位姿个数,N1≥3,N2为特征点数目,N2≥2,A为相机内参,Ri、Ti为外参,K=(k1,k2,k3,p1,p2)为畸变系数,(Sx,Sy)表示x和y方向的沙姆角,(ox,oy)表示成像中心在成像面x和y方向的偏移量,mij表示第i个位姿下第j个特征点的实际像面的像素坐标,Mij表示第i个位姿下第j个特征点的重投影的实际成像面的像素坐标。
一种较佳的实施方式,步骤3具体为:
(1)将标定板图片放置于测量的景深范围内,与激光平面呈一定倾角,控制激光器在标定板图片表面形成激光线,采集有激光线的标定板图像,同时保持标定板位置不变,采集同一位置无激光线的标定板图像;重复采集不少于2组不同位置的有激光线和无激光线的标定板图像;
(2)对标定板图像进行特征检测,得到标定板图像特征点的图像坐标;
(3)对每组有激光线和无激光线的标定板图像进行做差,对做差后的图像通过中心线提取算法,获取每组图片对应的激光线的图像坐标(ui,vi),i=1,2,…,N为每组标定板图片的编号,N为采集的标定板图片的组数;
(4)对得到的激光线的图像坐标做线性拟合得到对应直线,找出每组中特征点距离直线最近的M个特征点;
(5)结合相机标定得到相机内参,利用每组图片中的M个特征点
Figure BDA0003431654870000061
计算与相机坐标系之间的旋转矩阵Ri和平移向量Ti
Figure BDA0003431654870000062
其中,
Figure BDA0003431654870000063
为第i组图片的M个特征点的像素坐标,
Figure BDA0003431654870000064
为该点在第i组图片的M个特征点在世界坐标系下的三维坐标,Ri和Ti为计算得到的第i组图片世界坐标系与相机坐标系之间的旋转矩阵和平移向量,s为缩放变量,A为相机内参;
(6)根据得到的标定板图片对应的外参,计算每组M个特征点在相机坐标系下的三维坐标
Figure BDA0003431654870000065
并拟合得到N个平面方程Aix+Biy+CiZ+Di=0;
Figure BDA0003431654870000066
(7)计算每条激光线在相机坐标系下的三维坐标
Figure BDA0003431654870000067
Figure BDA0003431654870000071
(8)对所有激光线在相机坐标系下的三维坐标
Figure BDA0003431654870000072
进行最小二乘平面拟合,得到激光平面方程PA*x+PB*y+PC*z+PD=0,其中方程中的各个系数通过拟合获得。
本发明实施例的有益效果包括:构建了一种离轴沙姆相机的成像模型,降低了离轴沙姆3D线激光相机的标定难度;通过获得的相机内参以及畸变系数,并使用非线性优化算法对沙姆角、偏移量进行迭代优化,提高了标定结果的精度。
附图说明
图1为离轴沙姆3D线激光相机结构图;
图2为本发明实施例的离轴沙姆相机的标定模型示意图;
图3为本发明实施例的离轴沙姆3D线激光相机的标定流程图;
图4为本发明实施例标定相机的某一位置的标定板图片;
图5为本发明实施例的标定激光平面的一组标定板图片;
图6为本发明实施例的获取到距离激光线最近的20个特征点。
具体实施方式
下面结合附图2、附图3、附图4、附图5和附图6,对本发明的具体实施方式进行详细说明。下面所描述的实施例仅仅是本发明的一部分实施例,而不是全部实施例。基于本发明的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图2为本发明实施例的离轴沙姆相机的标定模型示意图,构建图像像素坐标位置与世界坐标位置之间的关系,由已知特征点的图像坐标求解相机的模型参数。
如图3所示,一种离轴沙姆3D线激光相机的标定方法,包括:
步骤1:构建离轴沙姆3D线激光相机的成像模型:将图像坐标系转换到实际像面坐标系,构建与透镜面平行且成像中心过光轴的虚拟像面坐标系,将实际像面坐标系转换到虚拟像面坐标系,在虚拟像面坐标系进行畸变校正,得到虚拟像面坐标系下无失真的虚拟像面坐标,再将无失真的虚拟像面坐标转换到实际像面坐标系下无失真的实际像面坐标,再将实际像面坐标转换到图像坐标系对应的无失真的图像像素坐标,再近似按照小孔成像模型,将无失真的图像像素坐标转换到相机坐标系下的坐标,再将相机坐标系下坐标转换到世界坐标系下的坐标。从而得到图像像素坐标与世界坐标系下坐标的转换关系,实现离轴沙姆3D线激光相机成像模型的构建;
步骤2:基于成像模型进行相机标定:先计算相机内参和外参初值,后利用非线性优化算法计算最优解,进行相机标定获得相机标定参数:相机内参、畸变系数、沙姆角、偏移量;
步骤3:对激光平面进行标定:采集不少于2组用于激光平面标定的标定板图片,每组中标定板位于同一位置,包含有激光线图片和无激光线图片;通过组内标定板图片做差的方法,得到激光线数据,通过中心线提取算法提取激光线中心像素坐标;对每组无激光线标定板图片,提取特征点坐标,以距离激光线最邻近的特征点计算标定板图片对应的外参,根据标定板图片对应的外参,计算每组激光线中心像素坐标在相机坐标系下的三维坐标,通过拟合得到激光平面方程,实现对激光平面的标定。
其中,步骤1具体为:
(1)图像坐标系下的一点qi=(ui,vi)T转换到实际像面坐标系下的点qs=(xs,ys)T,实际像面坐标系为倾斜的像面;
Figure BDA0003431654870000091
其中,sx、sy分别表示像素在x方向和y方向的实际物理尺寸,(u0,v0)表示图像坐标系原点在图像坐标系中的坐标;
(2)将实际像面坐标系下的点qs转换到虚拟像面坐标系下的点qp=(xp,yp)T,虚拟像面垂直于光轴且成像中心在光轴上;
Figure BDA0003431654870000092
其中,ox、oy分别表示实际像面坐标系中心在x方向和y方向的偏移量,Sx、Sy分别表示倾斜的实际像面中心与虚拟像面的x方向和y方向的沙姆角;
(3)在虚拟像面坐标系上引入畸变模型,将虚拟像面坐标系下的点qp变换到无失真的点
Figure BDA0003431654870000093
Figure BDA0003431654870000094
其中,k1、k2、k3表示径向畸变参数,p1、p2表示切向畸变参数,r2=xp 2+yp 2
(4)将虚拟像面坐标系上无失真的点
Figure BDA0003431654870000095
转换到实际像面坐标系上无失真的点
Figure BDA0003431654870000101
Figure BDA0003431654870000102
(5)将实际像面坐标系上无失真的点
Figure BDA0003431654870000103
转换到图像坐标系下无失真的像素点
Figure BDA0003431654870000104
Figure BDA0003431654870000105
(6)将实际像面坐标系上无失真的像素点
Figure BDA0003431654870000106
转换到相机坐标系下的点qc=(xc,yc,zc)T
Figure BDA0003431654870000107
其中,αx、αy表示x和y方向的有效焦距;
(7)将相机坐标系下的点qc转换为世界坐标系下的点qw=(xw,yw,zw)T
Figure BDA0003431654870000108
其中,R表示由世界坐标系到相机坐标系的旋转矩阵,T表示由世界坐标系到相机坐标系的平移向量。
步骤2中先计算相机内参和外参初值,具体为:通过张氏解析算法计算相机的内参和外参初值。
步骤2中利用非线性优化算法计算最优解,进行相机标定,具体为:建立内参、畸变参数、外参以及沙姆角和偏移量的参数的非线性优化方程,通过迭代循环求解进行优化,最终的最优解即为相机标定参数。
步骤2具体转化过程为:
(1)对不少于N1张的标定板图片,提取每张标定板图片上的N2个特征点;
(2)采用张氏解析算法,计算相机内参αx,αy,u0,v0,以及外参Ri,Ti
(3)利用非线性优化算法对内参、畸变系数、外参以及沙姆角和偏移量进行优化,对应的非线性优化模型如下:
Figure BDA0003431654870000111
L为以重投影误差最小为准则的目标方程,i表示标定板第i个位姿,j表示第j个特征点,N1为标定板位姿个数,N1≥3,N2为特征点数目,N2≥2,A为相机内参,Ri、Ti为外参,K=(k1,k2,k3,p1,p2)为畸变系数,(Sx,Sy)表示x和y方向的沙姆角,(ox,oy)表示成像中心在成像面x和y方向的偏移量,mij表示第i个位姿下第j个特征点的实际像面的像素坐标,Mij表示第i个位姿下第j个特征点的重投影的实际成像面的像素坐标。
具体对相机标定时的一个具体实施方式,使用9*11的圆形图案标定板,将标定板放置激光平面上,在相机的视野范围内,上下以及左右移动,采集不同位置的30张相机标定图片,采集的某一位置的图片如图4所示。
步骤3具体标定过程为:
(1)将标定板图片放置于测量的景深范围内,与激光平面呈一定倾角,控制激光器在标定板图片表面形成激光线,采集有激光线的标定板图像,同时保持标定板位置不变,采集同一位置无激光线的标定板图像;重复采集不少于2组不同位置的有激光线和无激光线的标定板图像。具体一个实施方式是重复直至采集4组激光平面标定图片,采集的某一位置的两种图片如图5所示。
(2)对标定板图像进行特征检测,得到标定板图像特征点的图像坐标;
(3)对每组有激光线和无激光线的标定板图像进行做差,对做差后的图像通过中心线提取算法,获取每组图片对应的激光线的图像坐标(ui,vi),i=1,2,…,N为每组标定板图片的编号,N为采集的标定板图片的组数;采用灰度重心法对激光线进行检测,可以清晰得到激光线的图像坐标。
(4)对得到的激光线的图像坐标做线性拟合得到对应直线,找出每组中特征点距离直线最近的M个特征点;具体一个实施方式,如图6所示,可以找到与激光线相邻最近的20个圆点中标记圆圈的特征点。
(5)结合相机标定得到相机内参,利用每组图片中的M个特征点
Figure BDA0003431654870000121
计算与相机坐标系之间的旋转矩阵Ri和平移向量Ti
Figure BDA0003431654870000122
其中,
Figure BDA0003431654870000123
为第i组图片的M个特征点的像素坐标,
Figure BDA0003431654870000124
为该点在第i组图片的M个特征点在世界坐标系下的三维坐标,Ri和Ti为计算得到的第i组图片世界坐标系与相机坐标系之间的旋转矩阵和平移向量,s为缩放变量,A为相机内参;
(6)根据得到的标定板图片对应的外参,计算每组M个特征点在相机坐标系下的三维坐标
Figure BDA0003431654870000131
并拟合得到N个平面方程Aix+Biy+CiZ+Di=0;
Figure BDA0003431654870000132
(7)计算每条激光线在相机坐标系下的三维坐标
Figure BDA0003431654870000133
Figure BDA0003431654870000134
(8)对所有激光线在相机坐标系下的三维坐标
Figure BDA0003431654870000135
进行最小二乘平面拟合,得到激光平面方程PA*x+PB*y+PC*z+PD=0,其中方程中的各个系数通过拟合获得。
具体一个实施方式中,使用的相机分辨率为1920*1080,x方向和y方向的沙姆角分别为[0°,5.71°],x方向和y方向的偏移量为[0,4.15]单位为毫米。
按照本发明实施例中的方法,采用张氏解析算法得到的相机内参为:
Figure BDA0003431654870000136
经非线性优化后所得到的的沙姆角和偏移量的优化结果为:
[0.57804,6.42792]
[-0.01706,4.16547]
重投影误差为:0.17647
激光平面的方程系数为:0.0042*x-2.3758y+z-320.1950=0
从而,可以获得离轴沙姆3D线激光相机的标定:相机的内参、畸变系数、沙姆角、偏移量以及激光平面在相机坐标系下的平面方程系数。
本发明实施例的离轴沙姆3D线激光相机的标定方法,构建了离轴沙姆相机的成像模型,实现了离轴沙姆3D线激光相机的标定;并且,通过获得的相机内参以及畸变系数,并使用非线性优化算法对沙姆角、偏移量进行迭代优化,提高了标定结果的精度。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (6)

1.一种离轴沙姆3D线激光相机的标定方法,其特征在于,包括:
步骤1:构建所述离轴沙姆3D线激光相机的成像模型:将图像坐标系转换到实际像面坐标系,构建与透镜面平行且成像中心过光轴的虚拟像面坐标系,将实际像面坐标系转换到虚拟像面坐标系,在虚拟像面坐标系进行畸变校正,得到虚拟像面坐标系下无失真的虚拟像面坐标,再将无失真的虚拟像面坐标转换到实际像面坐标系下无失真的实际像面坐标,再将实际像面坐标转换到图像坐标系对应的无失真的图像像素坐标,再近似按照小孔成像模型,将无失真的图像像素坐标转换到相机坐标系下的坐标,再将相机坐标系下坐标转换到世界坐标系下的坐标。从而得到图像像素坐标与世界坐标系下坐标的转换关系,实现离轴沙姆3D线激光相机成像模型的构建;
步骤2:基于所述成像模型进行相机标定:先计算相机内参和外参初值,后利用非线性优化算法计算最优解,进行相机标定获得所述相机标定参数:相机内参、畸变系数、沙姆角、偏移量;
步骤3:对激光平面进行标定:采集不少于2组用于激光平面标定的标定板图片,每组中所述标定板位于同一位置,包含有激光线图片和无激光线图片;通过组内标定板图片做差的方法,得到激光线数据,通过中心线提取算法提取激光线中心像素坐标;对每组所述无激光线标定板图片,提取特征点坐标,以距离激光线最邻近的特征点计算所述标定板图片对应的外参,根据所述标定板图片对应的外参,计算每组所述激光线中心像素坐标在相机坐标系下的三维坐标,通过拟合得到激光平面方程,实现对所述激光平面的标定。
2.根据权利要求1所述的方法,其特征在于,所述步骤1具体为:
(1)所述图像坐标系下的一点qi=(ui,vi)T转换到所述实际像面坐标系下的点qs=(xs,ys)T,所述实际像面坐标系为倾斜的像面;
Figure FDA0003431654860000021
其中,sx、sy分别表示像素在x方向和y方向的实际物理尺寸,(u0,v0)表示所述图像坐标系原点在所述图像坐标系中的坐标;
(2)将所述实际像面坐标系下的点qs转换到所述虚拟像面坐标系下的点qp=(xp,yp)T,所述虚拟像面垂直于所述光轴且成像中心在所述光轴上;
Figure FDA0003431654860000022
其中,ox、oy分别表示所述实际像面坐标系中心在x方向和y方向的偏移量,Sx、Sy分别表示倾斜的所述实际像面中心与所述虚拟像面的x方向和y方向的沙姆角;
(3)在所述虚拟像面坐标系上引入畸变模型,将所述虚拟像面坐标系下的点qp变换到无失真的点
Figure FDA0003431654860000023
Figure FDA0003431654860000024
其中,k1、k2、k3表示径向畸变参数,p1、p2表示切向畸变参数,r2=xp 2+yp 2
(4)将所述虚拟像面坐标系上无失真的点
Figure FDA0003431654860000025
转换到所述实际像面坐标系上无失真的点
Figure FDA0003431654860000026
Figure FDA0003431654860000027
(5)将所述实际像面坐标系上无失真的点
Figure FDA0003431654860000031
转换到图像坐标系下无失真的像素点
Figure FDA0003431654860000032
Figure FDA0003431654860000033
(6)将所述实际像面坐标系上无失真的像素点qi,转换到相机坐标系下的点qc=(xc,yc,zc)T
Figure FDA0003431654860000034
其中,αx、αy表示x和y方向的有效焦距;
(7)将相机坐标系下的点qc转换为世界坐标系下的点qw=(xw,yw,zw)T
Figure FDA0003431654860000035
其中,R表示由世界坐标系到相机坐标系的旋转矩阵,T表示由世界坐标系到相机坐标系的平移向量。
3.根据权利要求1所述的方法,其特征在于,所述步骤2中所述先计算相机内参和外参初值,具体为:通过张氏解析算法计算相机的内参和外参初值。
4.根据权利要求1所述的方法,其特征在于,所述步骤2中所述利用非线性优化算法计算最优解,进行相机标定,具体为:建立内参、畸变参数、外参以及沙姆角和偏移量的参数的非线性优化方程,通过迭代循环求解进行优化,最终的最优解即为相机标定参数。
5.根据权利要求1、3或4所述的方法,其特征在于,所述步骤2具体为:
(1)对不少于N1张的标定板图片,提取每张所述标定板图片上的N2个特征点;
(2)采用张氏解析算法,计算相机内参αx,αy,u0,v0,以及外参Ri,Ti
(3)利用非线性优化算法对内参、畸变系数、外参以及沙姆角和偏移量进行优化,对应的非线性优化模型如下:
Figure FDA0003431654860000041
L为以重投影误差最小为准则的目标方程,i表示标定板第i个位姿,j表示第j个特征点,N1为标定板位姿个数,N1≥3,N2为特征点数目,N2≥2,A为相机内参,Ri、Ti为外参,K=(k1,k2,k3,p1,p2)为畸变系数,(Sx,Sy)表示x和y方向的沙姆角,(ox,oy)表示成像中心在成像面x和y方向的偏移量,mij表示第i个位姿下第j个特征点的实际像面的像素坐标,Mij表示第i个位姿下第j个特征点的重投影的实际成像面的像素坐标。
6.根据权利要求1所述的方法,其特征在于,所述步骤3具体为:
(1)将所述标定板图片放置于测量的景深范围内,与激光平面呈一定倾角,控制激光器在所述标定板图片表面形成激光线,采集有激光线的标定板图像,同时保持标定板位置不变,采集同一位置无激光线的标定板图像;重复采集不少于2组不同位置的有激光线和无激光线的标定板图像;
(2)对标定板图像进行特征检测,得到标定板图像特征点的图像坐标;
(3)对每组有激光线和无激光线的标定板图像进行做差,对做差后的图像通过中心线提取算法,获取每组图片对应的激光线的图像坐标(ui,vi),i=1,2,…,N为每组标定板图片的编号,N为采集的标定板图片的组数;
(4)对得到的所述激光线的图像坐标做线性拟合得到对应直线,找出每组中特征点距离直线最近的M个特征点;
(5)结合相机标定得到相机内参,利用每组图片中的M个特征点
Figure FDA0003431654860000051
计算与相机坐标系之间的旋转矩阵Ri和平移向量Ti
Figure FDA0003431654860000052
其中,
Figure FDA0003431654860000053
为第i组图片的M个特征点的像素坐标,
Figure FDA0003431654860000054
为该点在第i组图片的M个特征点在世界坐标系下的三维坐标,Ri和Ti为计算得到的第i组图片世界坐标系与相机坐标系之间的旋转矩阵和平移向量,s为缩放变量,A为相机内参;
(6)根据得到的所述标定板图片对应的外参,计算每组M个特征点在相机坐标系下的三维坐标
Figure FDA0003431654860000055
并拟合得到N个平面方程Aix+Biy+CiZ+Di=0;
Figure FDA0003431654860000056
(7)计算每条激光线在相机坐标系下的三维坐标
Figure FDA0003431654860000057
Figure FDA0003431654860000061
(8)对所有激光线在相机坐标系下的三维坐标
Figure FDA0003431654860000062
进行最小二乘平面拟合,得到激光平面方程PA*x+PB*y+PC*z+PD=0,其中方程中的各个系数通过拟合获得。
CN202111600729.9A 2021-12-24 2021-12-24 一种离轴沙姆3d线激光相机的标定方法 Pending CN114359405A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111600729.9A CN114359405A (zh) 2021-12-24 2021-12-24 一种离轴沙姆3d线激光相机的标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111600729.9A CN114359405A (zh) 2021-12-24 2021-12-24 一种离轴沙姆3d线激光相机的标定方法

Publications (1)

Publication Number Publication Date
CN114359405A true CN114359405A (zh) 2022-04-15

Family

ID=81101235

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111600729.9A Pending CN114359405A (zh) 2021-12-24 2021-12-24 一种离轴沙姆3d线激光相机的标定方法

Country Status (1)

Country Link
CN (1) CN114359405A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114993207A (zh) * 2022-08-03 2022-09-02 广东省智能机器人研究院 基于双目测量系统的三维重建方法
CN116071240A (zh) * 2023-03-07 2023-05-05 广东利元亨智能装备股份有限公司 图像拼接方法、装置、计算机设备和存储介质
CN116188594A (zh) * 2022-12-31 2023-05-30 梅卡曼德(北京)机器人科技有限公司 沙姆相机的标定方法、标定系统、装置和电子设备
CN116342718A (zh) * 2023-05-26 2023-06-27 合肥埃科光电科技股份有限公司 一种线激光3d相机的标定方法、装置、存储介质及设备

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114993207A (zh) * 2022-08-03 2022-09-02 广东省智能机器人研究院 基于双目测量系统的三维重建方法
CN116188594A (zh) * 2022-12-31 2023-05-30 梅卡曼德(北京)机器人科技有限公司 沙姆相机的标定方法、标定系统、装置和电子设备
CN116188594B (zh) * 2022-12-31 2023-11-03 梅卡曼德(北京)机器人科技有限公司 沙姆相机的标定方法、标定系统、装置和电子设备
CN116071240A (zh) * 2023-03-07 2023-05-05 广东利元亨智能装备股份有限公司 图像拼接方法、装置、计算机设备和存储介质
CN116071240B (zh) * 2023-03-07 2023-06-20 广东利元亨智能装备股份有限公司 图像拼接方法、装置、计算机设备和存储介质
CN116342718A (zh) * 2023-05-26 2023-06-27 合肥埃科光电科技股份有限公司 一种线激光3d相机的标定方法、装置、存储介质及设备
CN116342718B (zh) * 2023-05-26 2023-08-01 合肥埃科光电科技股份有限公司 一种线激光3d相机的标定方法、装置、存储介质及设备

Similar Documents

Publication Publication Date Title
CN114359405A (zh) 一种离轴沙姆3d线激光相机的标定方法
CN107121109B (zh) 一种基于前镀膜平面镜的结构光参数标定装置及方法
CN111243033B (zh) 一种优化双目相机外参数的方法
CN105716542B (zh) 一种基于柔性特征点的三维数据拼接方法
CN111486802A (zh) 基于自适应距离加权的旋转轴标定方法
CN107025670A (zh) 一种远心相机标定方法
CN112465912B (zh) 一种立体相机标定方法及装置
CN106971408B (zh) 一种基于时空转换思想的摄像机标定方法
CN109325981B (zh) 基于聚焦像点的微透镜阵列型光场相机几何参数标定方法
CN113160339A (zh) 一种基于沙姆定律的投影仪标定方法
CN109544643A (zh) 一种摄像机图像校正方法及装置
CN110889829A (zh) 一种基于鱼眼镜头的单目测距方法
CN110738608B (zh) 一种平面图像校正方法及系统
CN113012234B (zh) 基于平面变换的高精度相机标定方法
CN113724337B (zh) 一种无需依赖云台角度的相机动态外参标定方法及装置
CN108154536A (zh) 二维平面迭代的相机标定法
CN106600549A (zh) 鱼眼图像的矫正方法和装置
CN112598747A (zh) 一种单目相机与投影仪的联合标定方法
CN110108230B (zh) 基于图像差分与lm迭代的二值光栅投影离焦程度评估方法
CN116740187A (zh) 一种无重叠视场多相机联合标定方法
CN115546289A (zh) 一种基于机器人的复杂结构件三维形貌测量方法
CN114998448A (zh) 一种多约束双目鱼眼相机标定与空间点定位的方法
CN117115272A (zh) 用于降水粒子多角度成像的远心相机标定和三维重构方法
CN112489141B (zh) 车载摄像头单板单图带中继镜的产线标定方法及装置
CN115439541A (zh) 用于折射成像系统的玻璃方位标定系统及方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination