CN108921901B - 一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法 - Google Patents
一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法 Download PDFInfo
- Publication number
- CN108921901B CN108921901B CN201810420919.4A CN201810420919A CN108921901B CN 108921901 B CN108921901 B CN 108921901B CN 201810420919 A CN201810420919 A CN 201810420919A CN 108921901 B CN108921901 B CN 108921901B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- turntable
- camera
- tracker
- reference point
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明涉及一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法,包括:精密二轴转台与激光跟踪仪靶球一同做旋转运动,解算跟踪仪坐标系到转台坐标系的转换矩阵;光学反射球作为相机标定的光学参考点,固定于转台前方,通过光学反射球与跟踪仪靶球精密互换及跟踪仪测量,建立光学反射球中心在转台坐标系下的位置;相机与转台一同做二维转动,拍摄光学参考点,并同步记录各站位的角度值,建立虚拟标定控制场;拟合图像中反射球的椭圆轮廓,计算标定靶标圆心的图像位置;基于相机成像模型建立最小化目标函数,进行相机标定。本发明适合在大视场环境下,完成高精度相机内参标定;尤其适合相机测量深度变化大,传统标定靶标难以设计制作的场合。
Description
技术领域
本发明涉及传感器标定技术,具体涉及一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法。
背景技术
大视场视觉测量的测量空间范围从几米到几十米乃至上百米,是精密航空视觉导航、大型机械装备定位、大型构件制造与装配等先进技术领域中的重要测量手段。大视场视觉测量以其测量范围大,测量过程非接触等优势,在测量领域具有不可替代的地位。
大视场相机标定方法及其应用的研究,旨在解决大视场视觉测量中相机标定困难、标定精度不高的问题。相机标定技术作为视觉测量的核心技术,其精度直接影响测量精度,但常规相机标定方法由于受到视场范围有限、视场空间变化、靶标制作困难等诸多因素的制约,已无法满足大视场高精度的测量需求。
大视场相机标定与普通视场相机标定的理论基础相同,而大视场标定过程中,不可避免的问题是获取大尺寸的精密靶标变得较为困难,即使能够制作这样的靶标或构造虚拟靶标,也面临成本问题,以及难以适应相机测量空间和测量视场的变化问题,难以满足高精度需求的问题。一维靶标以其加工成本低,加工精度高等优点而得到广泛研究,但是一维靶标所含信息量少,精度会受到制约;采用小平面靶标多位置标定的方法,尽管能够灵活布置充满相机测量空间,但是各个位置之间缺乏联系与约束,不能有效描述整个测量空间的相机模型,影响标定精度,同时操作流程也往往较为复杂。在测绘领域,对于中高空工作的航空和卫星测绘相机一般在实验室内采取精密转台,平行光管或星点板的精密测角法来标定相机的内方位元素,这些方法需要精密调校平行光管的光轴和转台、相机的位姿关系,使用不够便捷。在大视场视觉测量和摄影测量领域常采用自标定方法实现摄像标定,该方法虽然布置较为灵活,但运算量较大,方程求解的鲁棒性较差,且标定精度通常无法与光学实验室校准方法相比。此外,基于非参数模型的标定方法不使用传统相机模型,而是用类似穷举的方式对相机靶面上的每个测量角度进行一一映射,虽然精度高,但工作量极为庞大。
申请号为CN201611094763.2的发明专利公开了一种基于单LED发光点和二维转台的复眼系统标定装置及标定方法,通过一维平移调整LED位置的方式实现标定;论文“二轴转台测角法用于线阵相机几何参数标定”提出一种基于二轴转台测角法的线阵相机标定方法,通过简单的设备高精度求取线阵相机的内参数和畸变参数,但是两者并未使用激光跟踪仪,无法建立大视场标定控制场,无法实现大视场下相机高精度的标定;论文“基于虚拟平面靶标的大视场摄像机标定技术研究”和论文“构建虚拟立体靶标的大视场高精度视觉标定”都构建出大尺度的虚拟靶标,但是没有用到转台,需要多次移动标定特征以布满测量空间才能进一步完成标定。
发明内容
本发明技术解决问题:克服现有技术的不足,提供一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法,解决在大视场测量深度变化范围大、标定精度要求高的场合中,难以设计标定靶标与视场相适应的问题。
为达到上述目的,本发明的技术方案是这样实现的:
一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法,该方法包括:
a、精密二轴转台与激光跟踪仪靶球一同做旋转运动,激光跟踪仪测量多个位置的靶球球心坐标,解算激光跟踪仪坐标系到精密二轴转台坐标系的转换矩阵;
b、跟踪仪靶球放置于精密二轴转台和激光跟踪仪前方,激光跟踪仪测量跟踪仪靶球坐标;光学反射球与跟踪仪靶球能实现精度无损互换,作为相机标定的光学参考点,则根据激光跟踪仪坐标系到精密二轴转台坐标系的转换矩阵,建立光学参考点在精密二轴转台坐标系下的位置;
c、相机与精密二轴转台一同做二维转动,拍摄光学参考点,并同步记录各转动位置的角度值;根据光学参考点在精密二轴转台初始位置下的位置以及转台转动角度,建立起虚拟标定控制场;
d、在虚拟标定控制场各转动位置下采集的光学参考点图像中,拟合光学参考点的椭圆轮廓,获取椭圆中心,消除透视投影畸变,计算标定光学参考点圆心的图像位置;
e、基于相机成像理论模型以及各转动位置下光学参考点圆心的图像位置,建立最小化目标函数,通过非线性优化方法得到目标函数的最优解。
步骤a中精密二轴转台与激光跟踪仪靶球一同做旋转运动,激光跟踪仪测量多个位置的靶球球心坐标,解算激光跟踪仪坐标系到精密二轴转台坐标系的转换矩阵;实现步骤如下:
(1)将跟踪仪靶球置于精密二轴转台任意位置,控制转台以设定角度间隔做旋转运动,激光跟踪仪采集并存储各转动位置下靶球球心坐标,得到一系列球心坐标,拟合空间圆,建立转台坐标系;
(2)根据转台坐标系原点在跟踪仪坐标系中的坐标位置,解算跟踪仪坐标系到转台坐标系的转换矩阵。
步骤b中跟踪仪靶球放置于精密二轴转台和激光跟踪仪前方,激光跟踪仪测量跟踪仪靶球坐标;光学反射球与跟踪仪靶球能实现精度无损互换,作为相机标定的光学参考点,则根据激光跟踪仪坐标系到精密二轴转台坐标系的转换矩阵,建立光学参考点在精密二轴转台坐标系下的位置;具体如下:
(1)将跟踪仪靶球置于精密二轴转台与激光跟踪仪前方,距离和位置可以根据应用要求自由调整,激光跟踪仪测量跟踪仪靶球球心的三维坐标;
(2)将跟踪仪靶球替换为直径相同,球心尺寸和位置精度相同,并带有高反光圆形特征平面的光学反射球,实现无精度损失的替换,作为相机标定的光学参考点,通过激光跟踪仪坐标系到精密二轴转台坐标系的转换矩阵,建立光学参考点中心位置在转台坐标系下的三维坐标。
步骤c中相机与精密二轴转台一同做二维转动,拍摄光学参考点,并同步记录各转动位置的角度值;根据光学参考点在精密二轴转台初始位置下的位置以及转台转动角度,建立起虚拟标定控制场,过程为:
(1)按照设定的角度间隔,相机与精密二轴转台一同做二维转动,在各转动位置下拍摄光学参考点,并同步记录各转动位置的角度值;
(2)根据标定光学参考点在精密二轴转台初始位置下的三维坐标,以及记录的转台转动角度,建立起虚拟的大尺度三维标定控制场。
步骤d中在虚拟标定控制场各转动位置下采集的光学参考点图像中,拟合光学参考点的椭圆轮廓,获取椭圆中心,消除透视投影畸变,计算标定光学参考点圆心的图像位置,过程为:
(1)在虚拟标定控制场各转动位置下采集的光学参考点图像中,采用Canny算子和Steger方法相结合的方式来提取光学参考点的边缘点;
(2)剔除无法形成封闭圆环的噪声点之后,使用椭圆拟合算法,拟合光学参考点子像素级的椭圆轮廓,进而计算椭圆中心位置在图像坐标系下的坐标;
(3)根据空间圆中心在相机平面上的畸变误差模型,消除相机透视投影畸变引起的偏心误差,计算标定参考点平面圆圆心精确的图像位置。
步骤e中基于相机成像理论模型以及各转动位置下光学参考点圆心的精确图像位置,通过非线性优化方法得到目标函数最优解,采用LM非线性优化方法求解精密二轴转台到相机的旋转矩阵、平移向量和相机内方位元素在最大似然准则下的最优解。
本发明与现有技术相比的优点在于:本发明提出的一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法,使用精密二轴转台作为角度溯源基准,激光跟踪仪作为长度溯源基准,转台带动相机作二维运动,采用单个固定不动的光学参考点作为参照,由相机在每个转动位置下进行拍摄,最终实现标定。期间相机在转台内框自由安装,并且转台转动控制程序、图像采集程序与标定程序由PC端一键式MFC自动完成,中间过程无需人工干预。转台转动幅度以及参考点与相机距离可根据应用视场需求任意灵活配置。本发明适合在大视场环境下,完成高精度相机内参标定;尤其适合相机测量空间尺度大,传统标定靶标难以设计制作的场合。
附图说明
图1为本发明实施例中基于精密二轴转台和激光跟踪仪的大视场相机标定方法流程图;
图2为本发明实施例中大视场相机标定系统示意图;
图3为本发明实施例中跟踪仪坐标系到转台坐标系的标定示意图;
图4为本发明实施例中作为标定参考可与激光跟踪仪靶球机械互换的光学反射球的机械尺寸图;
图5为本发明实施例中大视场虚拟三维靶标效果示意图。
具体实施方式
本发明的基本思想是:通过转台及跟踪仪虚拟出满足大视场要求的精密大尺寸立体靶标,进而根据相机成像模型建立最小化目标函数,采用LM非线性优化方法求解转台到相机的旋转矩阵、平移向量和相机内参最优解,即实现相机的标定。
下面结合附图及由精密二轴转台、相机、激光跟踪仪和圆形光学参考点组成的具体实施例,对本发明作进一步详细说明。
如图1所示,本发明基于精密二轴转台和激光跟踪仪的大视场相机标定方法主要包括以下步骤:
步骤11:精密二轴转台与激光跟踪仪靶球一同做旋转运动,激光跟踪仪测量多个位置的靶球球心坐标,解算激光跟踪仪坐标系到精密二轴转台坐标系的转换矩阵。
具体包括以下步骤:
步骤111:将跟踪仪靶球置于精密二轴转台任意位置,控制转台以设定角度间隔做旋转运动,激光跟踪仪采集并存储各转动位置下靶球球心坐标,得到一系列球心坐标,拟合空间圆,建立转台坐标系。
标定系统示意图如图2所示,其中精密二轴转台包括三个部分,即固定基座4、外框2和内框3,1为内框转轴,5为外框转轴,标定参考点6放置于固定架7上,8为激光跟踪仪。其中所用的精密二轴转台的步长值应小于相机单像素所占视场角。在图中将相机部分放大,具体建立转台坐标系(ORXRYRZR)在转台内框上,与内框固连,随内框转动而转动,将世界坐标系与转台坐标系重合,相机坐标系(OCXCYCZC),跟踪仪坐标系(OTXTYTZT)。跟踪仪靶球分别固定于精密二轴转台外框的任意位置A处和内框的任意位置B处,控制转台外框做旋转运动,内框不转动,每隔一定角度记录跟踪仪坐标系下靶球球心坐标,转动结束后将转台回归零位。用同样的方式控制内框做旋转运动,并每隔设定角度记录靶球球心坐标。两次单轴转动得到两组点,依次进行空间圆拟合,得到的两个圆分别称为圆OA、圆OB,并可得出跟踪仪坐标系下两圆圆心坐标OA=[xTA,yTA,zTA]T和OB=[xTB,yTB,zTB]T、两圆平面通过各自圆心的单位法向量nx和ny。将两个法向量的交点作为转台坐标系的原点OC,nx作为转台坐标系X轴方向,ny作为转台坐标系Y轴方向,两个法向量的叉乘nx×ny=nz作为转台坐标系Z轴方向,正向如图。
步骤112:根据转台坐标系原点在跟踪仪坐标系中的坐标位置,解算跟踪仪坐标系到转台坐标系的转换矩阵。
跟踪仪坐标系到转台坐标系的标定示意图如图3所示,设初始位置处跟踪仪坐标系到转台坐标系的变换关系由旋转矩阵RTR和平移向量tTR表示,转台坐标系到跟踪仪坐标系的变换关系由旋转矩阵RRT和平移向量tRT表示。根据拟合的精密二轴转台外框圆形转动轨迹10、内框圆形转动轨迹9,计算圆平面单位法向量直线的交点坐标并将其作为转台坐标系原点,另外设该交点在跟踪仪坐标系下的坐标则有:
显然tRT=[xTC,yTC,zTC]T,又转台坐标系到跟踪仪坐标系的旋转矩阵为RRT=[nx,ny,nz],可计算得到RTR=[nx,ny,nz]-1,tTR=-[nx,ny,nz]-1[xTC,yTC,zTC]T。
步骤12:跟踪仪靶球放置于精密二轴转台和激光跟踪仪前方,激光跟踪仪测量跟踪仪靶球坐标;光学反射球与跟踪仪靶球能实现精度无损互换,作为相机标定的光学参考点,则根据激光跟踪仪坐标系到精密二轴转台坐标系的转换矩阵,建立光学参考点在精密二轴转台坐标系下的位置。
具体包括以下步骤:
步骤121:将跟踪仪靶球置于精密二轴转台与激光跟踪仪前方,距离和位置可以根据应用要求自由调整,激光跟踪仪测量跟踪仪靶球球心的三维坐标。
将跟踪仪靶球置于精密二轴转台与激光跟踪仪前方,激光跟踪仪测量固定不动的跟踪仪靶球三维精密坐标,具体测量场朝向前方,测得初始位置处跟踪仪靶球中心在跟踪仪坐标系下的坐标,设为pT0=[xT0,yT0,zT0]T,测量方法在激光跟踪仪的使用说明中有详细描述。
步骤122:将跟踪仪靶球替换为直径相同,球心尺寸和位置精度相同,并带有高反光圆形特征平面的光学反射球,实现无精度损失的替换,作为相机标定的光学参考点,通过激光跟踪仪坐标系到精密二轴转台坐标系的转换矩阵,建立光学参考点中心位置在转台坐标系下的三维坐标。
具体将相机固定安装在精密二轴转台的内框上,其测量视场与激光跟踪仪测量场的方向大致相同。将跟踪仪靶球替换为直径相同,球心尺寸和位置精度相同,并带有高反光圆形特征平面的光学反射球,实现无精度损失的精密替换,作为相机标定光学参考点。作为标定参考点可与激光跟踪仪靶球机械互换的光学反射球的机械尺寸图如图4所示,左图为正视图,反光区域的位置误差为右图为侧视图,反光区域的尺寸误差为跟踪仪靶球和光学反射球都属于精密标准器件,忽略加工误差,理论上光学反射球与跟踪仪靶球球心坐标相同,即零位处光学反射球球心在跟踪仪坐标系下的坐标为pT0=[xT0,yT0,zT0]T,又定义转台坐标系下光学反射球球心的坐标为pR0=[xR0,yR0,zR0]T,两坐标之间的关系为:
步骤13:相机与精密二轴转台一同做二维转动,拍摄光学参考点,并同步记录各转动位置的角度值;根据光学参考点在精密二轴转台初始位置下的位置以及转台转动角度,建立起虚拟标定控制场。
具体包括以下步骤:
步骤131:按照设定的角度间隔,相机与精密二轴转台一同做二维转动,在各转动位置下拍摄光学参考点,并同步记录各转动位置的角度值。
相机在精密二轴转台的带动下实现水平方向和垂直方向的精密旋转,根据镜头焦距f、相机分辨率PIX按照公式计算视场角,设计实验转动站位数量,进而得到转动间隔。大视场虚拟三维靶标效果示意图如图5所示,使标定光学参考点轨迹覆盖相机全幅视场,同时记录转动角度,该过程等价于生成一块大的三维立体靶标,且靶标精度很高,由跟踪仪精度和转台精度共同决定。
步骤132:根据标定光学参考点在精密二轴转台初始位置下的三维坐标,以及记录的转台转动角度,建立起虚拟的大尺度三维标定控制场。
其中,
表示转台坐标系相对零位或固定基座变化的旋转矩阵。带入可得旋转后标定参考点在转台坐标系下的坐标为:
根据标定光学参考点在精密二轴转台初始位置下的三维坐标,以及记录的转台转动角度,建立起虚拟的大尺度三维标定控制场。
步骤14:在虚拟标定控制场各转动位置下采集的光学参考点图像中,拟合光学参考点的椭圆轮廓,获取椭圆中心,消除透视投影畸变,计算标定光学参考点圆心的图像位置。
具体包括以下步骤:
步骤141:在虚拟标定控制场各转动位置下采集的光学参考点图像中,采用Canny算子和Steger方法相结合的方式来提取光学参考点的边缘点。
在虚拟标定控制场各转动位置下采集的光学参考点图像中,标定光学参考点具有高反光效果,圆形特征在图像上显示为白色区域,而周边具有黑色圆圈,需要对黑白区域交界的边缘进行检测,利用灰度的梯度信息检测非常合适,采用Canny算子和Steger方法相结合的方式来提取边缘中心点。首先对光学参考点的局部区域使用Canny算子进行梯度计算,梯度计算后数值范围不定,以梯度的最大值和最小值为界,进行线性变化,将梯度值变化至[0-255]的灰度范围内。使用Steger方法对线性变化后的梯度图像进行子像素级提取。
步骤142:剔除无法形成封闭圆环的噪声点之后,使用椭圆拟合算法,拟合光学参考点子像素级的椭圆轮廓,进而计算椭圆中心位置在图像坐标系下的坐标。
从提取的子像素级图像可以看到,边缘光条中心点包括一些噪声点,宫政的博士论文“现场分布式视觉测量关键技术研究”第三章中设计简单的剔除噪声方法:设置阈值处理下限,当光条中心点所在的像素值小于该阈值时,认为该点是噪声点;阈值处理后仍然存在高梯度值的噪声点,以像素灰度值最高的光条中心为起点,沿着梯度的垂直方向进行搜索,依次连接光条中心点,若无法封闭则认为起点为高梯度值的噪声点。剔除无法形成封闭圆环的噪声点之后,提取到完整的椭圆轮廓,使用椭圆拟合算法,得到椭圆中心位置。
步骤143:根据空间圆中心在相机平面上的畸变误差模型,消除相机透视投影畸变引起的偏心误差,计算标定参考点平面圆圆心精确的图像位置。
相机像面与空间圆所在平面不平行情形下,空间圆在图像平面上的投影不是一个标准的圆,存在透视投影畸变,基于透视投影变换和空间解析几何理论,建立了透视投影变换下空间圆中心在摄像机像平面上的畸变误差模型,进而消除透视投影畸变引起的偏心误差,计算标定参考点平面圆圆心精确的图像位置。
步骤15:基于相机成像理论模型以及各转动位置下光学参考点圆心的精确图像位置,通过非线性优化方法得到目标函数最优解,采用LM非线性优化方法求解精密二轴转台到相机的旋转矩阵、平移向量和相机内方位元素在最大似然准则下的最优解。
具体步骤如下:
由公式(5)可知,跟踪仪坐标系到转台坐标系的旋转矩阵是RiRTR,平移向量RitTR,设转台坐标系到相机坐标系的变换关系由旋转矩阵RRC和平移向量tRC表示,则此时标定参考点在相机坐标系下的坐标为:
对应的相机坐标系(OCXCYCZC)下的点,进行归一化(xn,vn)=(xc/zc,yc/zc),像点到主点距离无畸变条件下像点在以像素为单位的图像坐标系下的点为(ui,vi),根据相机成像模型,有:
其中,ax是u轴上的归一化焦距;同理ay称为v轴上的归一化焦距。(u0,v0)为像平面上的主点在像素坐标系下的坐标。上述模型是建立在理想的小孔成像线性模型基础上的,考虑径向畸变参数(k1,k2,k3)、切向畸变参数(p1,p2)、仿射和非正交变形参数(b1,b2)
因此,世界坐标系下的点在相机平面上的像素投影可以被修正为:
根据图像坐标系下检测值等于计算值,内方位元素、外方位元素为未知量,理论上差值等于零,由于误差的存在,实际上需要建立最小化目标函数:
其中xins=[ax,ay,u0,v0,k1,k2,k3,p1,p2,b1,b2]T;
采用LM非线性优化方法,选用合理的初值保证优化过程的计算速度和收敛性,进而求解旋转矩阵RTR、平移向量tTR和相机内参xins在最大似然准则下的最优解,完成标定。
以上所述,仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。
Claims (2)
1.一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法,其特征在于:实现步骤如下:
a、精密二轴转台与激光跟踪仪靶球一同做旋转运动,激光跟踪仪测量多个位置的靶球球心坐标,解算激光跟踪仪坐标系到精密二轴转台坐标系的转换矩阵;
b、跟踪仪靶球放置于精密二轴转台和激光跟踪仪前方,激光跟踪仪测量跟踪仪靶球坐标;光学反射球与跟踪仪靶球能实现精度无损互换,作为相机标定的光学参考点,则根据激光跟踪仪坐标系到精密二轴转台坐标系的转换矩阵,建立光学参考点在精密二轴转台坐标系下的位置;具体实现如下:
(1)将跟踪仪靶球置于精密二轴转台与激光跟踪仪前方,距离和位置可以根据应用要求自由调整,激光跟踪仪测量跟踪仪靶球球心的三维坐标;
(2)将跟踪仪靶球替换为直径相同,球心尺寸和位置精度相同,并带有高反光圆形特征平面的光学反射球,实现无精度损失的替换,作为相机标定的光学参考点,通过激光跟踪仪坐标系到精密二轴转台坐标系的转换矩阵,建立光学参考点中心位置在转台坐标系下的三维坐标;
c、相机与精密二轴转台一同做二维转动,拍摄光学参考点,并同步记录精密二轴转台各转动位置的角度值;根据光学参考点在精密二轴转台初始位置下的位置以及转台转动角度,建立起虚拟标定控制场;具体过程为:
(1)按照设定的角度间隔,相机与精密二轴转台一同做二维转动,在各转动位置下拍摄光学参考点,并同步记录各转动位置的角度值;
(2)根据标定光学参考点在精密二轴转台初始位置下的三维坐标,以及记录的转台转动角度,建立起虚拟的大尺度三维标定控制场;
d、在虚拟标定控制场各转动位置下采集的光学参考点图像中,拟合光学参考点的椭圆轮廓,获取椭圆中心,消除透视投影畸变,计算标定光学参考点圆心的图像位置;具体过程为:
(1)在虚拟标定控制场各转动位置下采集的光学参考点图像中,采用Canny算子和Steger方法相结合的方式来提取光学参考点的边缘点;
(2)剔除无法形成封闭圆环的噪声点之后,使用椭圆拟合算法,拟合光学参考点子像素级的椭圆轮廓,进而计算椭圆中心位置在图像坐标系下的坐标;
(3)根据空间圆中心在相机平面上的畸变误差模型,消除相机透视投影畸变引起的偏心误差,计算标定参考点平面圆圆心精确的图像位置;
e、步骤e中基于相机成像理论模型以及各转动位置下光学参考点圆心的精确图像位置,通过非线性优化方法得到目标函数最优解,采用Levenberg-Marquardt(LM)非线性优化方法求解精密二轴转台到相机的旋转矩阵、平移向量和相机内方位元素在最大似然准则下的最优解最小化目标函数如下:
其中xins=[ax,ay,u0,v0,k1,k2,k3,p1,p2,b1,b2]T代表摄像机内参;[RRC tRC]为转台坐标系到相机坐标系的转换矩阵,RRC为转台坐标系到相机坐标系的旋转矩阵,tRC为转台坐标系到相机坐标系的平移向量;
2.根据权利要求1所述的基于精密二轴转台和激光跟踪仪的大视场相机标定方法,其特征在于:步骤a中精密二轴转台与激光跟踪仪靶球一同做旋转运动,激光跟踪仪测量多个位置的靶球球心坐标,解算激光跟踪仪坐标系到精密二轴转台坐标系的转换矩阵;实现步骤如下:
(1)将跟踪仪靶球置于精密二轴转台任意位置,控制转台以设定角度间隔做旋转运动,激光跟踪仪采集并存储精密二轴转台的各转动位置下靶球球心坐标,得到一系列球心坐标,拟合空间圆,建立转台坐标系;
(2)根据转台坐标系原点在跟踪仪坐标系中的坐标位置,解算跟踪仪坐标系到转台坐标系的转换矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810420919.4A CN108921901B (zh) | 2018-05-04 | 2018-05-04 | 一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810420919.4A CN108921901B (zh) | 2018-05-04 | 2018-05-04 | 一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108921901A CN108921901A (zh) | 2018-11-30 |
CN108921901B true CN108921901B (zh) | 2022-03-22 |
Family
ID=64403099
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810420919.4A Active CN108921901B (zh) | 2018-05-04 | 2018-05-04 | 一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108921901B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110248148B (zh) | 2018-09-25 | 2022-04-15 | 浙江大华技术股份有限公司 | 一种确定定位参数的方法及装置 |
CN109712201B (zh) * | 2019-02-27 | 2023-09-15 | 中国科学院上海技术物理研究所 | 一种广域相机的定位能力标定装置及标定方法 |
CN110211175B (zh) * | 2019-05-21 | 2023-04-21 | 天津大学 | 准直激光器光束空间位姿标定方法 |
CN110246191B (zh) * | 2019-06-13 | 2021-03-16 | 易思维(杭州)科技有限公司 | 相机非参数模型标定方法及标定精度评估方法 |
CN110458894B (zh) * | 2019-07-31 | 2023-05-09 | 天津大学 | 一种测量机的相机与接触式测头的标定方法 |
CN110926366A (zh) * | 2019-12-13 | 2020-03-27 | 浙江省计量科学研究院 | 一种基于激光跟踪仪多站布局的曲面轮廓测量方法 |
CN111337908A (zh) * | 2020-03-30 | 2020-06-26 | 苏州华兴源创科技股份有限公司 | 激光雷达检测系统及其检测方法 |
CN112393882B (zh) * | 2020-04-21 | 2022-08-23 | 哈尔滨工业大学 | 一种基于显微成像微透镜参数检测的复眼成像调校方法 |
CN111551193B (zh) * | 2020-05-25 | 2023-09-05 | 天津科技大学 | 一种基于精密二轴转台的激光靶标定方法 |
CN111649668B (zh) * | 2020-06-01 | 2022-02-15 | 中国工程物理研究院激光聚变研究中心 | 诊断设备位置匹配方法、装置、激光跟踪仪和诊断设备标定系统 |
CN112435302B (zh) * | 2020-12-09 | 2024-05-31 | 北京理工大学 | 一种基于高精度转台和平行光管的远距离大视场鱼眼相机标定方法 |
CN112595231A (zh) * | 2020-12-10 | 2021-04-02 | 上海镭隆科技发展有限公司 | 一种基于智能图像识别二维随动系统动态监测装置及方法 |
CN112802144B (zh) * | 2021-01-27 | 2023-12-22 | 中国人民解放军火箭军工程大学 | 广角visar物面重构方法 |
CN113177987B (zh) * | 2021-04-30 | 2022-03-25 | 北京航空航天大学 | 一种视觉跟踪测量系统外场全局标定方法及系统 |
CN113188444B (zh) * | 2021-05-06 | 2023-01-13 | 上海航天测控通信研究所 | 激光测量系统不确定度测试实验与计算方法 |
CN113701992A (zh) * | 2021-09-07 | 2021-11-26 | 哈尔滨理工大学 | 一种三自由度运动的自动标定系统 |
CN113917477A (zh) * | 2021-10-08 | 2022-01-11 | 安徽创谱仪器科技有限公司 | 光路的搭建方法 |
CN114295051B (zh) * | 2021-12-29 | 2022-08-23 | 中国科学院近代物理研究所 | 一种谐波线圈磁场测量的自动化定位装置及方法 |
CN114509002B (zh) * | 2022-01-20 | 2023-06-06 | 渤海造船厂集团有限公司 | 一种大尺寸高精度环状t型结构精密检测装置及方法 |
CN114893669B (zh) * | 2022-05-26 | 2024-03-26 | 江苏省建筑工程集团有限公司 | 一种应用于塔机附着装置安装的监测装置及监测方法 |
CN115657734B (zh) * | 2022-11-07 | 2023-08-01 | 哈尔滨工业大学 | 超精密光学组件调整环姿态控制系统及调节方法 |
CN116152357B (zh) * | 2023-04-04 | 2023-07-28 | 国科天成科技股份有限公司 | 一种无限远对焦相机的参数标定系统和方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663763B (zh) * | 2012-04-25 | 2014-08-13 | 南京航空航天大学 | 一种面向大视场高精度视觉测量的摄像机标定方法 |
CN105913439A (zh) * | 2016-04-22 | 2016-08-31 | 清华大学 | 一种基于激光跟踪仪的大视场摄像机标定方法 |
CN107883870A (zh) * | 2017-10-24 | 2018-04-06 | 四川雷得兴业信息科技有限公司 | 基于双目视觉系统和激光跟踪仪测量系统的全局标定方法 |
-
2018
- 2018-05-04 CN CN201810420919.4A patent/CN108921901B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663763B (zh) * | 2012-04-25 | 2014-08-13 | 南京航空航天大学 | 一种面向大视场高精度视觉测量的摄像机标定方法 |
CN105913439A (zh) * | 2016-04-22 | 2016-08-31 | 清华大学 | 一种基于激光跟踪仪的大视场摄像机标定方法 |
CN107883870A (zh) * | 2017-10-24 | 2018-04-06 | 四川雷得兴业信息科技有限公司 | 基于双目视觉系统和激光跟踪仪测量系统的全局标定方法 |
Non-Patent Citations (2)
Title |
---|
《二维旋转平台下的相机参数标定》;顾国华等;《光学精密工程》;20170715;第25卷(第7期);第1891-1898页 * |
《基于虚拟标定场的数码相机内参数标定方法》;吴少平等;《系统仿真学报》;20061230;第18卷;第424-426页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108921901A (zh) | 2018-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108921901B (zh) | 一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法 | |
CN108648232B (zh) | 一种基于精密二轴转台的双目立体视觉传感器一体式标定方法 | |
CN111536902B (zh) | 一种基于双棋盘格的振镜扫描系统标定方法 | |
CN102376089B (zh) | 一种标靶校正方法及系统 | |
CN109859272B (zh) | 一种自动对焦双目摄像头标定方法及装置 | |
CN102364299B (zh) | 一种多个结构光投影三维型面测量头的标定技术 | |
CN109272574B (zh) | 基于投影变换的线阵旋转扫描相机成像模型构建方法和标定方法 | |
CN105716542B (zh) | 一种基于柔性特征点的三维数据拼接方法 | |
CN103267491B (zh) | 自动获取物体表面完整三维数据的方法及系统 | |
CN106871787B (zh) | 大空间线扫描成像三维测量方法 | |
CN110378969B (zh) | 一种基于3d几何约束的汇聚式双目相机标定方法 | |
CN108253939B (zh) | 可变视轴单目立体视觉测量方法 | |
Zhang et al. | A universal and flexible theodolite-camera system for making accurate measurements over large volumes | |
CN105698699A (zh) | 一种基于时间转轴约束的双目视觉测量方法 | |
CN110146038A (zh) | 筒形件装配转角的分布式单目相机激光测量装置及方法 | |
CN109579695B (zh) | 一种基于异构立体视觉的零件测量方法 | |
CN113205603A (zh) | 一种基于旋转台的三维点云拼接重建方法 | |
CN111768451A (zh) | 一种基于移动显示屏的大尺寸双目视觉离焦标定方法 | |
CN111879354A (zh) | 一种无人机精细化测量系统 | |
CN113119129A (zh) | 一种基于标准球的单目测距定位方法 | |
CN108537849A (zh) | 基于同心圆环的三维直角靶标的线阵相机的标定方法 | |
CN105374067A (zh) | 一种基于pal相机的三维重建方法及其重建系统 | |
JP2005003463A (ja) | キャリブレーションチャート画像表示装置、キャリブレーション装置、キャリブレーション方法 | |
CN113947638A (zh) | 鱼眼相机影像正射纠正方法 | |
CN113962853B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |