CN103426168B - 基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法 - Google Patents
基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法 Download PDFInfo
- Publication number
- CN103426168B CN103426168B CN201310303466.4A CN201310303466A CN103426168B CN 103426168 B CN103426168 B CN 103426168B CN 201310303466 A CN201310303466 A CN 201310303466A CN 103426168 B CN103426168 B CN 103426168B
- Authority
- CN
- China
- Prior art keywords
- centerdot
- camera
- prime
- point
- theta
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000005457 optimization Methods 0.000 claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims description 27
- 238000003384 imaging method Methods 0.000 claims description 11
- 238000000354 decomposition reaction Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000013519 translation Methods 0.000 claims description 8
- 230000003287 optical effect Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 3
- 238000002474 experimental method Methods 0.000 abstract description 3
- 241000251468 Actinopterygii Species 0.000 abstract 2
- 241000208340 Araliaceae Species 0.000 abstract 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 abstract 1
- 235000003140 Panax quinquefolius Nutrition 0.000 abstract 1
- 235000008434 ginseng Nutrition 0.000 abstract 1
- 230000000007 visual effect Effects 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 8
- 238000012360 testing method Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Image Analysis (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
一种基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法,该方法有四大步骤:步骤一:摄像机内参的初始化;步骤二:摄像机外参的初始化;步骤三:优化目标函数的建立;步骤四:摄像机内外参的非线性优化。由于现有的标定算法不能够满足鱼眼镜头摄像机立体视觉***的准确性与方便性,因此本发明即在此通用相机模型的基础上设计了对于普通、广角、鱼眼镜头摄像机通用的立体视觉***标定方法。通过Matlab仿真和真实的实验可以验证本发明是可行的并且有效的。它在机器视觉技术领域里有较好的实用价值和应用前景。
Description
技术领域
本发明涉及一种基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法,该发明属于机器视觉领域。
背景技术
在机器视觉领域,立体视觉是重要的组成部分,被广泛应用于导航、三维测量及场景建模等领域。当前大部分的研究是集中于标定普通和广角立体摄像机***,这些方法并不能直接用于标定鱼眼摄像机。另外,普通和广角摄像机***只能提供比较有限的公共视场,而鱼眼立体摄像机***的公共视场却能覆盖摄像机***整个区域。因此,为了更好的发挥鱼眼摄像机的功能,提出一种精确、高效、方便的标定方法是非常重要和有意义的。
发明内容
本发明提供了一种基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法。它解决了鱼眼立体摄像机***乃至多摄像机混合***的标定复杂以及灵活性不高的问题。
本发明中采用的通用摄像机模型如下:
r(θ)=k1θ+k2θ3(1)
这一模型仅有两个参数k1和k2,由于r(θ)在[0,π/2]区间单调递增,因此模型可以反向求解。
如图1所示,一个三维空间点P通过一个鱼眼摄像机成像点为p,若通过***成像点则是p′,p′′为P在以Oc为球心的半球体的投影点。Oc-XcYcZc表示摄像机坐标系,o-xy表示图像坐标系。通过下式计算p的图像像素坐标为:
其中r(θ)可由式(1)获得,是径向方向与o-xy的x正轴方向的夹角,(x,y)为图像坐标(单位:mm),(u0,v0)是主点像素坐标,dx和dy分别是水平和垂直方向的像素大小。p′′点的坐标可表示为对应的准齐次坐标p′′′为:
对于任意空间点及其成像像素坐标可通过下列计算式求得θ和
得到径向距离为:
由式(1)通过卡丹公式求解一元三次方程组选择符合θ∈[0,π/2]的解。
径向夹角可通过反向求解获得
对于立体视觉而言,需要标定的参数包括各摄像机内参和摄像机间外参。摄像机内参即为摄像机模型中{k1,k2,u0,v0,dx,dy}六个参数,摄像机间外参即为摄像机坐标系之间的位置和姿态的关系,一般用3×3旋转矩阵R和3×1平移向量T表示:
本发明提出了一种基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法,其基本思路为:首先利用摄像机产品说明书和摄像机成像特征点对信息计算出立体摄像机内外参数,然后利用计算得到的摄像机内外参数作为优化初始值,对建立的目标函数进行两步优化,从而得到摄像机内外参数的最优值。
本发明一种基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法,该方法具体步骤如下:
步骤一:摄像机内参的初始化
如上所述,每个摄像机的内参包括{k1,k2,u0,v0,dx,dy}六个参数。
1、主点(u0,v0)是光轴通过图像坐标系交点的坐标,一般可直接用图像中心点坐标初始化。
2、像素大小dx和dy可以摄像机生产商提供值作为初始值。
3、对于普通和广角摄像机,k1和k2可以通过摄像机生产商提供的标称焦距f和最大视角θmax利用透视投影模型r=ftanθ来得到。具体实现方法为:将0~θmax等分取m(m≥2)个θi(i=1,…,m),分别代入透视投影模型r=ftanθ和(1)有
此时,即得到关于k1和k2的m个线性方程,可通过线性最小二乘方法求得k1和k2。
其中
而对于鱼眼相机而言,则可令k1=f,k2=0。
步骤二:摄像机外参的初始化
摄像机外参的初始化可通过对摄像机***的本质矩阵E进行奇异值分解得到。
1、对初始特征点进行归一化处理
给定n对两幅图像中的特征点集,根据式(4)~(6)以及各摄像机初始内参可得各特征点成像时对应的θi和通过式(3)计算出特征点在各摄像机模型的半球面投影点的准齐次坐标为其中Xi=(xi,yi,1)T,Xi′=(xi′,yi′,1)T。
令
其中,
2、计算本质矩阵
将归一化后的数据Yi和Yi′(i=1,…,n)代入
对n×9的矩阵Ae进行奇异值分解得
其中Se=diag(σ1,…,σ9),且σ1≥…≥σ9,diag表示对角阵,0(n-9)×9为(n-9)×9的全0矩阵。Ve为9×9矩阵,可由Ve的第9列元素得到:
采用奇异值分解将Enorm分解为
其中,D1≥D2≥D3。令D3=0,则有:
可得本质矩阵为:
3、利用本质矩阵求取摄像机外参的可能解R′和T′
当通过式(17)获得本质矩阵E后,摄像机外参的可能解R′和T′即可通过本质矩阵E的奇异值分解来获得。令
E=USVT(18)
此时R′和T′有四个可能的组合:
这里
为了最终得到正确的组合,需遵循以下准则:1)R′的行列式值为1;2)每一个重建的三维点X必然位于两个摄像机的前方,也就是说,X的Z坐标必须全为正值。如因噪声影响不全为正,应取使得重建的X的Z坐标为正值最多的组合作为最终正确组合。
4、重建三维点得到真实解
假设三维点M在左右两个摄像机O0-X0Y0Z0坐标系和O1-X1Y1Z1坐标系的坐标分别为M0=(X0,Y0,Z0)T和M1=RM0+T,且点在两相机的成像像素坐标点为(um,vm)和(u′m,v′m),由式(4)~(6)及初始内参可求得对应的θ0,θ1,
此时即可利用线性最小二乘方法可直接求解出X0,Y0,Z0:
其中
因此可将(19)中四组可能解分别替换(R,T)代入式(21)中进行三维点重建计算,利用判断准则对可能解进行筛选,直到得到最终的唯一解,此时有R=R′。由于||T′||=1,因此计算得到的平移向量T'与真实的平移向量T相差一个常值因子。令L′表示两个重建三维特征点的距离,L表示真实距离。为了最小化计算误差,常值因子λ需要使用两个相机所有的图像点对计算得到,此时有
其中N表示图像点对的数量。最终求得的转移向量为:
步骤三:目标函数的建立
本发明中一维标定杆上特征点的个数并没有限制,只要能保证每一对图像上同一三维点的不同成像点能够准确匹配就能完成标定。在本发明中在一维标定杆上装载了三个特征点,如图2所示,若因应用需要添加特征点,依次类推计算即可。由于噪声,一维标定杆上的特征点重建以后会产生误差。令d(A,B)=||A-B||(||·||表示欧几里得向量范数),得到距离和共线误差可表示为:
其中 k1,k2,u0,v0,dx,dy和k1′,k2′,u0′,v0′,dx′,dy′分别为左右摄像机的内部参数;可以通过Rodrigues公式与旋转矩阵R互相转化,即 为单位向量;为第k幅图像对中通过(21)计算的A,B,C的重建点;L,L1,L2分别为真实特征点AC、AB及BC之间的实际距离。
假设a0,k,b0,k,c0,k和a1,k,b1,k,c1,k是采集第k幅图像点对时通过式(3)求得的三维特征点A,B,C在左右摄像机半球面模型上投影点对应的准齐次坐标,有:
其中,[T]×为由T定义的反对称矩阵:
结合式(24)~(30),得到的优化目标函数为
其中ρ1∈(800,1200),ρ2∈(0.5,1.5)。
步骤四:内外参数的非线性优化
由步骤三知,共有18个参数需要优化,可以将优化过程分为两步:
1、固定y为y0,即y=y0,y0可由步骤一的初始值给定。设定x的估计为初始值,利用LM(Levenberg-Marquardt)方法来求解如下优化问题
得到最优解x*。
2、设定y0为初始值,利用LM求解如下优化问题
得到最优解y*。
此时,即得到了摄像机内外参数的最优解x*和y*。
其中,在步骤四中所述的“LM方法”是一种使用广泛的非线性最小二乘算法,中文为列文伯格-马夸尔特法,即Levenberg-Marquardt法。
优点与功效:
本发明一种基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法,其优点在于,解决了当前标定算法操作的复杂性以及针对鱼眼相机标定的不精确性等问题。本发明通过采用一维标定杆作为标定特征物以及通用的相机模型,使得在标定中实现了对不同相机进行同时标定,大大减小了实际应用中标定的工作量,同时具有精度高,操作方便,简单易行等特点。
附图说明
图1:鱼眼相机模型
图2:球形立体视觉示意图
图3:不同视觉***标定仿真结果图
图4:实际试验使用的一维标定杆图
图5:为本发明所述方法流程图
图中符号说明如下:
图1中的符号说明:O表示图像坐标系原点,x,y为图像坐标系的坐标轴。Oc表示摄像机坐标系的原点,Xc,Yc,Zc为摄像机坐标系的坐标轴。P为空间一个三维点,p为P在鱼眼摄像机下的成像点,p′为P在***下的成像点,p′′为P在以Oc为球心的半球体上的投影点。f为摄像机的焦距,r为鱼眼摄像机成像点p距离主点的径向距离。θ为入射线与摄像机光轴之间的夹角,是径向距离r与图像坐标系x正轴方向的夹角。
图2中的符号说明:A,B,C(Ak,Bk,Ck分别表示它们在第k幅图像对中的位置)分别表示一维标定杆上的三个共线特征点;O0-X0Y0Z0和O1-X1Y1Z1分别表示左右两个摄像机所处的摄像机坐标系,(R,T)分别为两摄像机之间的旋转矩阵和平移向量;a0,k,b0,k,c0,k和a1,k,b1,k,c1,k分别表示三维空间点Ak,Bk,Ck分别投影到以O0为球心的单位半球体上的和以O1为球心的单位半球体上的投影点。
具体实施方式
本发明一种基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法,见图1~图5所示。
仿真过程是在主频3.07GHz、内存4.00GB的计算机上,WindowsXP环境下的Matlab2009b上进行的。实验对象采用的是BaslerscA640-120gc相机,视场角为185°的FujinonFE185C057HA-1鱼眼镜头,视场角为86.77°的PentaxC60402KP广角镜头。
首先介绍一下仿真中立体视觉摄像机的各个参数,这些参数根据真实摄像机参数假设。采集图像对数量N=300,图像主点u0=u0′=320,v0=v0′=240,像素大小dx=dx′=0.0056mm,dy=dy′=0.0056mm。普通镜头的视场角为45°,标称焦距为5.5mm;广角镜头的视场角为100°,标称焦距为2.2mm;鱼眼镜头的视场角为185°,标称焦距为1.8mm。仿真实验分为四组不同的立体视觉仿真:1)普通镜头立体摄像机组;2)广角镜头立体摄像机组;3)鱼眼镜头立体摄像机组;4)广角镜头与鱼眼镜头混合组成的立体摄像机组;
(1)本发明一种基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法,见图5,该方法具体步骤如下:
步骤一:摄像机内参的初始化
根据仿真参数,所有***中u0=u0′=320,v0=v0′=240,dx=dx′=0.0056,dy=dy′=0.0056;
对于普通相机,根据式(9)得k1=5.3727,k2=2.6344;对于广角相机,根据式(9)得k1=1.4586,k2=2.2714;
对于鱼眼相机,根据仿真设置有k1=f=1.8,k2=0;
步骤二:摄像机外参的初始化
利用随机序列产生N对图像对仿真实际过程中摄像机采集到的一维标定杆特征点,根据该图像对数据和初始化内参利用式(4)~(6)计算得到对应的θ和,其次利用(10)~(17)求得本质矩阵E,通过式(18)~(21)对本质矩阵E进行奇异值分解并利用立体视觉特征筛选出正确的摄像机内外参数(R′,T′),此时有R=R′,通过式(22)~(23)最小化计算误差,得到平移向量T;
步骤三:目标函数的建立
利用式(24)~(30)计算各类误差项,其中L1=200mm,L2=400mm,L=600mm,建立的目标函数为:
其中ρ1=1000,ρ2=1;
步骤四:内外参数的非线性优化
给定上述图像对数据以及内外参数的初始值,对步骤三中的目标函数进行两步优化,直接使用matlab的lsqnonlin函数,得到最终的最优内外参数x*和y*。
(2)仿真结果分析
仿真过程中,一维标定杆在每个立体视觉***中随机挥舞300次,得到了N=300组图像对数据。为了描述偏差,在得到的图像点中添加均值μ=0、标准差σ=0~3个像素的高斯白噪声,每次试验执行10次去均值作为结果。图3显示了不同的立体视觉***的标定结果,这里RMS表示线段AC的真实值与所有图像对重建值的均方根误差。
由图可知,即使当噪声标准差σ=3个像素时,RMS误差依然小于4%,表明本发明提出的算法有高的准确性和稳定性。
(3)实物试验验证
为了进一步验证本发明的可行性,我们在真实的立体视觉***中进行了测试。真实试验主要完成了三个不同的立体视觉***:1)鱼眼镜头立体视觉***;2)广角镜头立体视觉***;3)鱼眼和广角镜头混合立体视觉***。试验采用的一维标定杆如图4所示,它是一根镶有三个共线LED灯的空心棒,这三个LED灯之间的距离分别为600mm,400mm和200mm。同样以实际标志点距离AC与标定后重建距离作为衡量标准,不同立体视觉***的实际实验结果如下表1‐3:
表1.鱼眼镜头立体视觉***中测量值与实际值的比较
表2.广角镜头立体视觉***中测量值与实际值的比较
表3.混合镜头立体视觉***中测量值与实际值的比较
Claims (1)
1.一种基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法,其特征在于:该方法具体作法如下:
本通用标定方法所采用的通用摄像机模型如下:
r(θ)=k1θ+k2θ3(1)
这一模型仅有两个参数k1和k2,由于r(θ)在[0,π/2]区间单调递增,因此模型是反向求解;
一个三维空间点P通过一个鱼眼摄像机成像点为p,若通过***成像点则是p′,p″为P在以Oc为球心的半球体的投影点,Oc-XcYcZc表示摄像机坐标系,o-xy表示图像坐标系,通过下式计算p的图像像素坐标为:
其中r(θ)由式(1)获得,是径向方向与o-xy的x正轴方向的夹角,(x,y)为图像坐标,(u0,v0)是主点像素坐标,dx和dy分别是水平和垂直方向的像素大小;p″点的坐标表示为对应的准齐次坐标p″′为:
对于任意空间点及其成像像素坐标通过下列计算式求得θ和
得到径向距离为:
由式(1)通过卡丹公式求解一元三次方程组选择符合θ∈[0,π/2]的解;
径向夹角通过反向求解获得
对于立体视觉而言,需要标定的参数包括各摄像机内参和摄像机间外参,摄像机内参即为摄像机模型中{k1,k2,u0,v0,dx,dy}六个参数,摄像机间外参即为摄像机坐标系之间的位置和姿态的关系,一般用3×3旋转矩阵R和3×1平移向量T表示:
该方法具体步骤如下:
步骤一:摄像机内参的初始化
每个摄像机的内参包括{k1,k2,u0,v0,dx,dy}六个参数;
1、主点(u0,v0)是光轴通过图像坐标系交点的坐标,一般直接用图像中心点坐标初始化;
2、像素大小dx和dy以摄像机生产商提供值作为初始值;
3、对于普通和广角摄像机,k1和k2通过摄像机生产商提供的标称焦距f和最大视角θmax利用透视投影模型r=ftanθ来得到;具体实现方法为:将0~θmax等分取m(m≥2)个θi(i=1,…,m),分别代入透视投影模型r=ftanθ和式(1),得到
此时,即得到关于k1和k2的m个线性方程,通过线性最小二乘方法求得k1和k2;
其中
而对于鱼眼相机而言,则令k1=f,k2=0;
步骤二:摄像机外参的初始化
摄像机外参的初始化通过对摄像机***的本质矩阵E进行奇异值分解得到;
1、对初始特征点进行归一化处理
给定n对两幅图像中的特征点集,根据式(4)~(6)以及各摄像机初始内参得各特征点成像时对应的θi和通过式(3)计算出特征点在各摄像机模型的半球面投影点的准齐次坐标为其中Xi=(xi,yi,1)T,X′i=(x′i,y′i,1)T;
令
其中,
2、计算本质矩阵
将归一化后的数据Yi和Y′i(i=1,…,n)代入
对n×9的矩阵Ae进行奇异值分解得
其中Se=diag(σ1,…,σ9),且σ1≥…≥σ9,diag表示对角阵,0(n-9)×9为(n-9)×9的全0矩阵;Ve为9×9矩阵,由Ve的第9列元素得到:
采用奇异值分解将Enorm分解为
其中,D1≥D2≥D3,令D3=0,则有:
得到本质矩阵为:
E=(T′norm)TE′normTnorm(17)
3、利用本质矩阵求取摄像机外参的可能解R′和T′
当通过式(17)获得本质矩阵E后,摄像机外参的可能解R′和T′即通过本质矩阵E的奇异值分解来获得;令
E=USVT(18)
此时R′和T′有四个组合:
这里
为了最终得到正确的组合,需遵循以下准则:1)R′的行列式值为1;2)每一个重建的三维点X必然位于两个摄像机的前方,也就是说,X的Z坐标必须全为正值,如因噪声影响不全为正,应取使得重建的X的Z坐标为正值最多的组合作为最终正确组合;
4、重建三维点得到真实解
假设三维点M在左右两个摄像机O0-X0Y0Z0坐标系和O1-X1Y1Z1坐标系的坐标分别为M0=(X0,Y0,Z0)T和M1=RM0+T,且点在两相机的成像像素坐标点为(um,vm)和(u′m,v′m),由式(4)~(6)及初始内参求得对应的θ0,θ1,
此时利用线性最小二乘方法直接求解出X0,Y0,Z0:
其中
因此将式(19)中四组可能解分别替换(R,T)代入式(21)中进行三维点重建计算,利用判断准则对可能解进行筛选,直到得到最终的唯一解,此时有R=R′;由于||T′||=1,因此计算得到的平移向量T'与真实的平移向量T相差一个常值因子;令L′表示两个重建三维特征点的距离,L表示真实距离;为了最小化计算误差,常值因子λ需要使用两个相机所有的图像点对计算得到,此时有
其中N表示图像点对的数量,最终求得的转移向量为:
步骤三:目标函数的建立
本方法中一维标定杆上特征点的个数并没有限制,只要能保证每一对图像上同一三维点的不同成像点能够准确匹配就能完成标定;本方法中在一维标定杆上装载了三个特征点,若因应用需要添加特征点,依次类推计算即行;由于噪声,一维标定杆上的特征点重建以后会产生误差,令d(A,B)=||A-B||,||·||表示欧几里得向量范数,得到距离和共线误差表示为:
其中k1,k2,u0,v0,dx,dy和k′1,k′2,u′0,v′0,dx′,dy′分别为左右摄像机的内部参数;通过Rodrigues公式与旋转矩阵R互相转化,即 为单位向量;为第k幅图像对中通过式(21)计算的A,B,C的重建点;L,L1,L2分别为真实特征点AC、AB及BC之间的实际距离;
假设a0,k,b0,k,c0,k和a1,k,b1,k,c1,k是采集第k幅图像点对时通过式(3)求得的三维特征点A,B,C在左右摄像机半球面模型上投影点对应的准齐次坐标,有:
其中,[T]×为由T定义的反对称矩阵:
结合式(24)~(30),得到的优化目标函数为
其中ρ1∈(800,1200),ρ2∈(0.5,1.5);
步骤四:内外参数的非线性优化
由步骤三知,共有18个参数需要优化,将优化过程分为两步:
1、固定y为y0,即y=y0,y0由步骤一的初始值给定,设定x的估计为初始值,利用LM方法来求解如下优化问题
得到最优解x*;
2、设定y0为初始值,利用LM求解如下优化问题
得到最优解y*;
此时,即得到了摄像机内外参数的最优解x*和y*。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310303466.4A CN103426168B (zh) | 2013-07-18 | 2013-07-18 | 基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310303466.4A CN103426168B (zh) | 2013-07-18 | 2013-07-18 | 基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103426168A CN103426168A (zh) | 2013-12-04 |
CN103426168B true CN103426168B (zh) | 2016-04-27 |
Family
ID=49650858
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310303466.4A Active CN103426168B (zh) | 2013-07-18 | 2013-07-18 | 基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103426168B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104376553B (zh) * | 2014-09-28 | 2017-04-05 | 北京航空航天大学 | 一种基于移动相机和双靶标的单轴旋转角的视觉测量方法 |
CN105139393B (zh) * | 2015-08-19 | 2018-03-02 | 中国航空工业集团公司北京长城航空测控技术研究所 | 线阵相机内参数的标定方法 |
CN105488766B (zh) * | 2015-11-23 | 2018-06-19 | 上海汽车集团股份有限公司 | 鱼眼镜头图像校正方法及装置 |
CN106267734A (zh) * | 2016-08-30 | 2017-01-04 | 石家庄铁路职业技术学院 | 一种数字化健身*** |
CN108288291A (zh) * | 2018-06-07 | 2018-07-17 | 北京轻威科技有限责任公司 | 基于单点标定物的多相机标定 |
CN109448060B (zh) * | 2018-09-04 | 2020-11-06 | 惠州市德赛西威智能交通技术研究院有限公司 | 一种基于蝙蝠算法的相机标定参数优化方法 |
EP3967969B1 (en) | 2019-05-08 | 2024-05-08 | Sichuan Visensing Technology Co., Ltd. | Fisheye camera calibration system, method and apparatus, electronic device, and storage medium |
CN110782498B (zh) * | 2019-09-26 | 2022-03-15 | 北京航空航天大学 | 一种视觉传感网络的快速通用标定方法 |
CN110677599B (zh) * | 2019-09-30 | 2021-11-05 | 西安工程大学 | 一种重建360度全景视频图像的***和方法 |
CN111598954A (zh) * | 2020-04-21 | 2020-08-28 | 哈尔滨拓博科技有限公司 | 一种快速高精度摄像头参数计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101299270A (zh) * | 2008-05-27 | 2008-11-05 | 东南大学 | 三维扫描***中的多个摄像机同步快速标定方法 |
CN102968794A (zh) * | 2012-11-22 | 2013-03-13 | 上海交通大学 | 小型二维靶大视场双目立体视觉***结构参数在位标定法 |
-
2013
- 2013-07-18 CN CN201310303466.4A patent/CN103426168B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101299270A (zh) * | 2008-05-27 | 2008-11-05 | 东南大学 | 三维扫描***中的多个摄像机同步快速标定方法 |
CN102968794A (zh) * | 2012-11-22 | 2013-03-13 | 上海交通大学 | 小型二维靶大视场双目立体视觉***结构参数在位标定法 |
Non-Patent Citations (2)
Title |
---|
A Flexible New Technique for Camera Calibration;Zhengyou Zhang;《Microsoft Research》;20080813;1-22 * |
双目鱼眼摄像机噪声分析与外参标定;李学东等;《Processing of the 30th Chinese Control Conference》;20110724;5017-5022 * |
Also Published As
Publication number | Publication date |
---|---|
CN103426168A (zh) | 2013-12-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103426168B (zh) | 基于一维标定杆的普通、广角、鱼眼立体摄像机的通用标定方法 | |
CN110296691B (zh) | 融合imu标定的双目立体视觉测量方法与*** | |
CN107705252B (zh) | 适用于双目鱼眼图像拼接展开校正的方法及*** | |
CN109272570A (zh) | 一种基于立体视觉数学模型的空间点三维坐标求解方法 | |
CN108288292A (zh) | 一种三维重建方法、装置及设备 | |
CN102750704B (zh) | 一种摄像机分步自标定方法 | |
CN107886546B (zh) | 利用球像及公共自极三角形标定抛物折反射摄像机的方法 | |
CN104376552A (zh) | 一种3d模型与二维图像的虚实配准算法 | |
CN101216296A (zh) | 双目视觉转轴标定方法 | |
CN103983186B (zh) | 双目视觉***校正方法及校正设备 | |
CN110782498B (zh) | 一种视觉传感网络的快速通用标定方法 | |
CN103278138A (zh) | 一种复杂结构薄部件三维位置及姿态的测量方法 | |
CN102855620B (zh) | 基于球形投影模型的纯旋转摄像机自标定方法 | |
CN104463791A (zh) | 一种基于球面模型的鱼眼图像校正法 | |
WO2019056782A1 (zh) | 一种基于球体投影公切线的多相机标定及参数优化方法 | |
CN101887585A (zh) | 基于非共面特征点的摄像机标定方法 | |
CN101354796B (zh) | 基于泰勒级数模型的全向立体视觉三维重建方法 | |
CN106157246A (zh) | 一种全自动的快速柱面全景图像拼接方法 | |
CN110044374A (zh) | 一种基于图像特征的单目视觉测量里程的方法及里程计 | |
CN1897715A (zh) | 立体视觉半实物仿真***及方法 | |
CN104504691B (zh) | 基于低秩纹理的摄像机位置和姿态测量方法 | |
CN105931261A (zh) | 一种双目立体相机外参数修正方法及装置 | |
CN107507232A (zh) | 基于多尺度迭代的立体匹配算法 | |
Lu et al. | Sphere-based calibration method for trinocular vision sensor | |
CN101661625B (zh) | 基于一幅标准图像生成各视角下样本图像的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |