[go: up one dir, main page]

CN117611684A - Structural parameter optimization calibration method for biprism virtual binocular vision system - Google Patents

Structural parameter optimization calibration method for biprism virtual binocular vision system Download PDF

Info

Publication number
CN117611684A
CN117611684A CN202311641278.2A CN202311641278A CN117611684A CN 117611684 A CN117611684 A CN 117611684A CN 202311641278 A CN202311641278 A CN 202311641278A CN 117611684 A CN117611684 A CN 117611684A
Authority
CN
China
Prior art keywords
virtual
biprism
calibration
vision system
camera
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
CN202311641278.2A
Other languages
Chinese (zh)
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.)
Tianjin University of Technology
Original Assignee
Tianjin University of Technology
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 Tianjin University of Technology filed Critical Tianjin University of Technology
Priority to CN202311641278.2A priority Critical patent/CN117611684A/en
Publication of CN117611684A publication Critical patent/CN117611684A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • G06T2207/10012Stereo images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence

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

The invention relates to a structural parameter optimization calibration method for a biprism virtual binocular vision system, which comprises the following steps: s1: adjusting the relative position of the camera and the biprism; s2: acquiring three-dimensional images of a plurality of groups of calibration plates; s3: acquiring initial values of internal parameters and structural parameters of a virtual binocular camera; s4: acquiring a stereoscopic image sequence, and establishing a space measurement field consisting of characteristic points with regular distribution; s5: performing distortion correction on sub-images corresponding to the left virtual camera and the right virtual camera by using the internal parameters of the virtual binocular camera obtained in the step S3; s6: acquiring feature point matching point pairs; s7: measuring three-dimensional coordinates of feature points in a field; s8: and constructing an objective function based on the space distance error and epipolar constraint, performing nonlinear iterative optimization on the structural parameters, and updating the three-dimensional coordinates of the feature points in the step S7 by using the structural parameters obtained by each iterative calculation. The method can meet the comprehensive requirements of the biprism virtual binocular system in terms of three-dimensional measurement and two-dimensional epipolar geometric constraint.

Description

一种用于双棱镜虚拟双目视觉系统的结构参数优化标定方法A structural parameter optimization and calibration method for biprism virtual binocular vision system

技术领域Technical field

本发明属于图像数字化处理技术领域,涉及视觉标定方法,具体是一种用于双棱镜虚拟双目视觉系统的结构参数优化标定方法。The invention belongs to the technical field of image digital processing and relates to a visual calibration method, specifically a structural parameter optimization calibration method for a biprism virtual binocular vision system.

背景技术Background technique

依据双棱镜折射原理成像的虚拟双目视觉系统因其结构紧凑且无双目同步取像误差的特点被大量用于电弧熔池形貌测量、立体显微镜等领域。由双棱镜折射后并由单个相机获取的立体图像可等效为由两个虚拟相机拍到的双目立体视觉图像对,并称其为虚拟双目视觉系统。该系统的成像原理如图1所示。The virtual binocular vision system based on the principle of biprism refraction is widely used in arc molten pool morphology measurement, stereomicroscope and other fields because of its compact structure and no binocular synchronization imaging error. The stereoscopic image refracted by the biprism and acquired by a single camera can be equivalent to a pair of binocular stereoscopic vision images captured by two virtual cameras, and is called a virtual binocular vision system. The imaging principle of this system is shown in Figure 1.

众所周知,高精度的标定是三维视觉测量的关键。根据以往的研究可知,将传统的双目标定方法应用到双棱镜虚拟双目系统中是切实可行且有效的。传统的双目标定方法一般是两个相机获取同一靶标图像并提取特征点,通过线性方法得到内参与结构参数初值,最后结合误差约束对标定初值进行非线性优化。在此领域的研究一般着力于引入不同的误差约束,并采用不同的优化目标函数对参数进行非线性优化标定。此外,将基于双棱镜的单镜头立体视觉系统当作一种虚拟双目系统的设想较为理想化,比如:该虚拟双目系统不存在物理意义上的相机主点。因此,直接将传统的标定方法应用到虚拟双目系统中会相应地限制其测量精度。As we all know, high-precision calibration is the key to three-dimensional vision measurement. According to previous research, it is feasible and effective to apply the traditional binocular calibration method to the biprism virtual binocular system. The traditional dual-target calibration method generally uses two cameras to acquire the same target image and extract feature points, obtain the initial values of the internal participation structure parameters through a linear method, and finally perform nonlinear optimization of the initial calibration values using error constraints. Research in this field generally focuses on introducing different error constraints and using different optimization objective functions to perform nonlinear optimization calibration of parameters. In addition, it is ideal to regard the single-lens stereo vision system based on biprism as a virtual binocular system. For example, the virtual binocular system does not have a physical camera principal point. Therefore, directly applying traditional calibration methods to a virtual binocular system will limit its measurement accuracy accordingly.

相比于传统的双目立体视觉系统,双棱镜虚拟双目视觉系统更适合在一定视野范围内进行测量工作,被测物与视觉系统的相对位置比较固定。例如通过检索发现,公告号为CN108830906B的中国专利提供了一种基于虚拟双目视觉原理的摄像机参数自动标定方法,采用装置包括云台旋转支架、摄像机、外部触发电路、PC、绘有黑白格棋盘图案的标定板和两面成夹角设置的平面镜;两面平面镜设于云台旋转支架后方;标定板置于云台旋转支架处;摄像机设于云台旋转支架前方且摄像方面朝向两平面镜;当摄像机对标定板摄像时,其摄取的标定板图像包括标定板在两平面镜内形成的左右两个虚像;当标定板每转动一个角度时,外部触发电路控制摄像机对标定板摄像;PC对摄像机拍摄的标定板图像进行处理,以获取并标定摄像机的内部参数、畸变参数以及在世界坐标系下的外部参数。Compared with the traditional binocular stereo vision system, the biprism virtual binocular vision system is more suitable for measurement work within a certain field of view, and the relative position of the measured object and the visual system is relatively fixed. For example, through search, it was found that the Chinese patent with announcement number CN108830906B provides an automatic calibration method of camera parameters based on the principle of virtual binocular vision. The device includes a pan-tilt rotating bracket, a camera, an external trigger circuit, a PC, and a chessboard with black and white grids. The calibration plate of the pattern and two plane mirrors arranged at an angle; the two plane mirrors are arranged behind the pan/tilt rotating bracket; the calibration plate is placed at the pan/tilt rotating bracket; the camera is arranged in front of the pan/tilt rotating bracket and the camera side faces the two plane mirrors; when the camera When the calibration board is photographed, the image of the calibration board it captures includes the left and right virtual images formed by the calibration board in the two plane mirrors; when the calibration board rotates an angle, the external trigger circuit controls the camera to take pictures of the calibration board; the PC captures the image of the calibration board. The calibration plate image is processed to obtain and calibrate the internal parameters, distortion parameters and external parameters of the camera in the world coordinate system.

公开号为CN114111637A的中国专利提供了一种基于虚拟双目的条纹结构光三维重建方法,包括构建基于双棱镜的虚拟双目条纹结构光视觉系统;对虚拟双目视觉系统进行标定,获得虚拟双目两个相机的内外参数;向物体投射三频四步正弦条纹图像;使用多频外差解包裹的方法得到分别对应左右虚拟相机的展开相位;根据虚拟相机的内外参数,对展开相位进行极线校正;对左右展开相位进行立体匹配,获得左右展开相位对应匹配点的视差;结合虚拟双目的内外参数,利用视差原理将视差转换为深度信息;根据物体的三维空间坐标生成物体稠密的点云图,完成物体三维重建。上述两项发明专利旨在对用虚拟双目视觉系统代替传统的传统双目系统,对参数进行标定,但标定方法过于繁琐,又如CN107121109A一种基于前镀膜平面镜的结构光参数标定装置及方法,包括:将前镀膜平面镜与平面玻璃靶标置于摄像机前位置,摄像机同时拍摄平面玻璃靶标图像与其镜像图像,建立虚拟双目测量模型,通过非线性优化方法求取前镀膜平面镜坐标系到摄像机坐标系旋转矩阵和平移矢量的优化解,采用最小二乘法对备选特征点进行图像消影点求取;将白色打印纸置于前镀膜平面镜前,摄像机同时拍摄实际光条图像与镜像光条图像,提取光条图像中心点,计算匹配点的三维坐标求解光平面方程。本发明改善光条质量提高光条中心提取精度,并提供微米级别位置精度的特征点,获得更多数量的标定点,具有更高的标定精度,标定结果更稳定,其标定客体并非双棱镜或相似结构的虚拟双目视觉系统,在标定涉及的算法逻辑有所不同。因此,结合双棱镜虚拟双目的系统特点发明一种双棱镜虚拟双目标定方法对于提高该系统的三维测量精度起着关键作用。The Chinese patent with publication number CN114111637A provides a three-dimensional reconstruction method based on virtual binocular striped structured light, including constructing a virtual binocular striped structured light vision system based on a biprism; calibrating the virtual binocular vision system to obtain a virtual binocular The internal and external parameters of the two cameras are projected; a three-frequency four-step sinusoidal stripe image is projected to the object; the unfolded phase corresponding to the left and right virtual cameras is obtained using the multi-frequency heterodyne unwrapping method; the unfolded phase is polarized according to the internal and external parameters of the virtual camera. Line correction; perform stereo matching on the left and right unfolded phases to obtain the parallax of the matching points corresponding to the left and right unfolded phases; combine the internal and external parameters of the virtual binoculars and use the parallax principle to convert the parallax into depth information; generate dense points of the object based on the three-dimensional spatial coordinates of the object Cloud image to complete the three-dimensional reconstruction of objects. The above two invention patents aim to replace the traditional binocular system with a virtual binocular vision system to calibrate parameters, but the calibration method is too cumbersome. Another example is CN107121109A, a structured light parameter calibration device and method based on a front-coated plane mirror. , including: placing the front-coated plane mirror and the plane glass target in front of the camera, the camera simultaneously shooting the plane glass target image and its mirror image, establishing a virtual binocular measurement model, and obtaining the front-coated plane mirror coordinate system to the camera coordinates through nonlinear optimization methods is the optimized solution of the rotation matrix and translation vector. The least squares method is used to calculate the image elimination point of the candidate feature points. The white printing paper is placed in front of the front coated plane mirror, and the camera simultaneously captures the actual light strip image and the mirrored light strip image. , extract the center point of the light strip image, calculate the three-dimensional coordinates of the matching point, and solve the light plane equation. The present invention improves the quality of the light strip, improves the accuracy of extracting the center of the light strip, and provides feature points with micron-level position accuracy, obtains a greater number of calibration points, has higher calibration accuracy, and the calibration results are more stable. The calibration object is not a biprism or a prism. Virtual binocular vision systems with similar structures have different algorithm logic involved in calibration. Therefore, inventing a biprism virtual binocular calibration method based on the characteristics of the biprism virtual binocular system plays a key role in improving the three-dimensional measurement accuracy of the system.

发明内容Contents of the invention

基于现有技术的不足,本发明的目的在于,提供一种用于双棱镜虚拟双目视觉系统的结构参数优化标定方法,用以提高系统的三维测量精度。Based on the shortcomings of the existing technology, the purpose of the present invention is to provide a structural parameter optimization and calibration method for a biprism virtual binocular vision system to improve the three-dimensional measurement accuracy of the system.

本发明解决其技术问题是采取以下技术方案实现的:The present invention solves its technical problems by adopting the following technical solutions:

一种用于双棱镜虚拟双目视觉系统的结构参数优化标定方法,A structural parameter optimization and calibration method for biprism virtual binocular vision system,

S1:调整相机与双棱镜的相对位置,使双棱镜虚拟双目视觉系统获取分割均匀的左、右子图像;S1: Adjust the relative position of the camera and the biprism so that the biprism virtual binocular vision system obtains evenly divided left and right sub-images;

S2:将标定板放置在双棱镜虚拟双目视觉系统景深范围下的不同位置,获取多组标定板的立体图像;S2: Place the calibration plate at different positions under the depth of field range of the biprism virtual binocular vision system, and obtain three-dimensional images of multiple sets of calibration plates;

S3:对标定板的原始立体图像进行均匀分割,得到两个尺寸一致的子图像并分别提取左、右子图上的特征点,对双棱镜虚拟双目视觉系统进行粗标定,获取虚拟双目相机的内参与结构参数初值;S3: Evenly segment the original stereo image of the calibration plate to obtain two sub-images with the same size and extract the feature points on the left and right sub-images respectively. Perform rough calibration of the biprism virtual binocular vision system to obtain the virtual binocular Initial values of the camera’s intrinsic structural parameters;

S4:将测量板放置于视觉系统的前方,固定距离多次移动测量板并获取立体图像序列,建立由分布规则的特征点组成的空间测量场;S4: Place the measurement board in front of the vision system, move the measurement board multiple times at a fixed distance and obtain a stereoscopic image sequence, and establish a spatial measurement field composed of regularly distributed feature points;

S5:对空间测量场内的若干组立体图像进行均匀分割,利用步骤S3得到的虚拟双目相机的内参对左、右虚拟相机对应的子图像进行畸变矫正;S5: Evenly segment several sets of stereo images in the spatial measurement field, and use the internal parameters of the virtual binocular camera obtained in step S3 to perform distortion correction on the sub-images corresponding to the left and right virtual cameras;

S6:提取畸变矫正后的测量板图像上的圆心特征点,根据圆心特征点序列的排列次序,对左、右子图的圆心特征点进行匹配,获取特征点匹配点对;S6: Extract the center feature points on the distortion-corrected measurement plate image, match the center feature points of the left and right sub-images according to the order of the circle center feature point sequences, and obtain feature point matching point pairs;

S7:结合虚拟双目相机的内参与结构参数,通过线性三角法恢复出空间测量场内特征点的三维坐标;S7: Combined with the intrinsic structural parameters of the virtual binocular camera, the three-dimensional coordinates of the feature points in the spatial measurement field are restored through linear triangulation;

S8:构建基于空间距离误差与极线约束的目标函数,对结构参数进行非线性迭代优化,且每次迭代计算得到的结构参数将用于更新步骤S7中的特征点的三维坐标。S8: Construct an objective function based on spatial distance error and epipolar constraints, perform nonlinear iterative optimization of structural parameters, and the structural parameters calculated in each iteration will be used to update the three-dimensional coordinates of the feature points in step S7.

而且,所述双棱镜虚拟双目视觉系统是由双棱镜折射后并由单个相机获取的立体图像可等效为由两个虚拟相机拍到的双目立体视觉图像对。Moreover, in the biprism virtual binocular vision system, the stereoscopic image refracted by the biprism and acquired by a single camera can be equivalent to a binocular stereoscopic image pair captured by two virtual cameras.

而且,所述S1的具体步骤为:利用精密位移台调整相机与双棱镜的位置,使图像垂直方向上的中线与双棱镜的中心的棱边趋近地重合。Moreover, the specific steps of S1 are: using a precision shift stage to adjust the positions of the camera and the biprism so that the center line in the vertical direction of the image closely coincides with the edge of the center of the biprism.

而且,利用张正友标定法对所述双棱镜虚拟双目视觉系统进行粗标定。Moreover, the biprism virtual binocular vision system was roughly calibrated using Zhang Zhengyou's calibration method.

而且,所述S4的具体步骤为:将一个表面具有规则排列的圆形靶标的测量板固定在高精密的电动位移台的移动平台上,测量板的空间位置相对平行于相机的成像面,操作控制器使测量板多次且定向地移动固定距离,同时获取每次位移位置的测量板立体图像。Moreover, the specific steps of S4 are: fix a measurement plate with regularly arranged circular targets on the surface on the mobile platform of a high-precision electric displacement stage. The spatial position of the measurement plate is relatively parallel to the imaging plane of the camera. Operation The controller makes the measurement plate move a fixed distance multiple times and in a direction, and at the same time obtains a three-dimensional image of the measurement plate at each displacement position.

而且,所述S8的具体步骤为:Moreover, the specific steps of S8 are:

S8.1:对空间测量场体积内的任意一个控制点进行六个方向上的定点距长度测量,基于空间距离约束的误差函数设计为:S8.1: Measure the fixed-point distance length in six directions for any control point within the volume of the spatial measurement field. The error function based on the spatial distance constraint is designed as:

Lq(q=1,2,……,6)代表指定控制点与其周围第q个特征点之间的实测长度,其长度的真值为Dq(q=1,2,……,6),i与j分别代表测量板的位移次序与测量板上的特征点圆心次序;Lq(q=1,2,…,6) represents the measured length between the specified control point and the q-th feature point around it, and the true value of its length is Dq(q=1,2,…,6), i and j respectively represent the displacement sequence of the measurement plate and the order of the center points of the characteristic points on the measurement plate;

S8.2:定义为一组匹配点的齐次坐标,两点到其所对应的极线ll与lr的距离分别为disl、disr,构建基于极线约束的误差函数为:S8.2: Definition are the homogeneous coordinates of a set of matching points. The distances between the two points and their corresponding epipolar lines l l and l r are dis l and dis r respectively. The error function based on the epipolar constraint is constructed as:

其中,为/>对应的左极线,/>为/>对应的右极线,F代表虚拟双目系统的基础矩阵,即:F=R·S,R=[r11,r12,r13;r21,r22,r23;r31,r32,r33]代表虚拟双目的旋转矩阵,S代表虚拟双目平移矩阵T=[t1,t2,t3]的反对称矩阵:in, for/> The corresponding left polar line,/> for/> The corresponding right pole line, F represents the basic matrix of the virtual binocular system, that is: F=R·S, R=[r 11 , r 12 , r 13 ; r 21 , r 22 , r 23 ; r 31 , r 32 , r 33 ] represents the virtual binocular rotation matrix, and S represents the antisymmetric matrix of the virtual binocular translation matrix T = [t 1 , t 2 , t 3 ]:

S8.3:结合两种误差函数组成目标函数:S8.3: Combine two error functions to form an objective function:

Hobj(x)=arg min(EL+EP)H obj (x)=arg min(E L +E P )

目标函数的输入变量为结构参数,即:x=(r11,r12,r13,r21,r22,r23,r31,r32,r33,t1,t2,t3),结合SQP(Sequential Quadratic Programming)非线性最优化算子对结构参数进行迭代求解,迭代截止以目标函数输出局部最小值为准。 The input variables of the objective function are structural parameters , that is : , combined with the SQP (Sequential Quadratic Programming) nonlinear optimization operator to iteratively solve the structural parameters, and the iteration cutoff is based on the local minimum value of the objective function output.

而且,所述S3的具体步骤为:Moreover, the specific steps of S3 are:

S3.1:将标定板的原始立体图像沿着图像垂直方向的中线均匀分割,得到两个像素尺寸相同的标定板子图像IL、IR,左、右两个子图像被分别视为左、右虚拟相机VCL、VCR拍摄所得;S3.1: Divide the original stereo image of the calibration plate evenly along the center line of the vertical direction of the image to obtain two calibration plate sub-images I L and I R with the same pixel size. The left and right sub-images are regarded as left and right respectively. Captured by virtual cameras VCL and VCR ;

S3.2:规定标定板上特征点的世界坐标并提取该特征点所对应的像素坐标/> 建立各个标定板图像上特征点的2D像素坐标与3D世界坐标的映射关系,以左虚拟相机为例,得到特征点世界坐标与像素坐标的转换关系:S3.2: Specify the world coordinates of feature points on the calibration plate And extract the pixel coordinates corresponding to the feature point/> Establish the mapping relationship between the 2D pixel coordinates and the 3D world coordinates of the feature points on each calibration plate image. Taking the left virtual camera as an example, the conversion relationship between the world coordinates and the pixel coordinates of the feature points is obtained:

其中,i与j分别代表拍摄的标定板序列与标定板上的特征点序列,KL为左虚拟相机的内参矩阵,为/>代表VCL的光心相对于第i个标定板的旋转矩阵与平移矩阵;Among them, i and j respectively represent the captured calibration plate sequence and the feature point sequence on the calibration plate, K L is the internal parameter matrix of the left virtual camera, for/> Represents the rotation matrix and translation matrix of the optical center of VC L relative to the i-th calibration plate;

S3.3:采用张正友标定法获取VCL,VCR的内参矩阵、畸变系数以及每幅标定板图像所对应的VCL与VCR光心之间的坐标转换关系,即:虚拟双目的结构参数,将多组虚拟双目的结构参数取平均值的结果作为系统标定的虚拟双目的结构参数。S3.3: Use Zhang Zhengyou’s calibration method to obtain the internal parameter matrices and distortion coefficients of VCL and VCR , as well as the coordinate conversion relationship between the optical centers of VCL and VCR corresponding to each calibration plate image, that is: the virtual binocular structure Parameters, the average result of multiple sets of virtual binocular structural parameters is used as the system calibrated virtual binocular structural parameters.

而且,每次迭代得到的所述结构参数将结合粗标定时获取的虚拟相机内参更新空间测量场内的控制点的三维坐标,新的结构参数与特征点阵列的三维坐标将会用于下一次的迭代计算。Moreover, the structural parameters obtained in each iteration will be combined with the virtual camera internal parameters obtained during rough calibration to update the three-dimensional coordinates of the control points in the spatial measurement field. The new structural parameters and the three-dimensional coordinates of the feature point array will be used for the next time iterative calculation.

而且,为校正镜头畸变,引入传统的径向畸变与切向畸变于虚拟双目视觉系统中,以左虚拟相机为例,其畸变模型被定义为:Moreover, in order to correct lens distortion, traditional radial distortion and tangential distortion are introduced into the virtual binocular vision system. Taking the left virtual camera as an example, the distortion model is defined as:

其中与mnl(xlm,ylm)分别为畸变点与去畸变点在归一化平面上的坐标,kl(1-2)与pl(1-2)各代表径向与切向畸变系数,/> in and m nl (x lm ,y lm ) are the coordinates of the distortion point and the de-distortion point on the normalized plane respectively, k l(1-2) and p l(1-2) respectively represent the radial and tangential distortion Coefficient,/>

一种电子设备,其特征在于:包括处理器以及存储器,存储器所述处理器之间互相通信连接,所述存储器中存储有计算机指令,所述处理器通过执行所述计算机指令,从而执行上述任一项所述的方法。An electronic device, characterized in that: it includes a processor and a memory, the memory and the processor are connected to each other, the memory stores computer instructions, and the processor executes any of the above by executing the computer instructions. method described in one item.

本发明的优点和积极效果是:The advantages and positive effects of the present invention are:

本发明适用于基于双棱镜的虚拟双目视觉系统的标定工作。由于双棱镜虚拟双目系统是一种特殊的双目视觉系统,传统的双目标定方法在应用于双棱镜虚拟双目系统中时的标定效果相对有限。此外,基于目标函数设计的双目视觉标定方法虽可提高标定精度,但这些方法也存在明显的缺陷,比如:使用了复杂的约束无法确定对标定精度提升贡献较大的约束类型,优化参数过多容易造成耦合解。因此,本发明具体有以下几个优点:The invention is suitable for calibration work of a virtual binocular vision system based on biprisms. Since the biprism virtual binocular system is a special binocular vision system, the calibration effect of traditional binocular calibration methods is relatively limited when applied to the biprism virtual binocular system. In addition, although binocular vision calibration methods based on objective function design can improve calibration accuracy, these methods also have obvious shortcomings. For example, the use of complex constraints cannot determine the type of constraints that contribute significantly to the improvement of calibration accuracy. How easy it is to cause coupled solutions. Therefore, the present invention specifically has the following advantages:

(1)本发明方法提出的标定方法充分结合了双棱镜虚拟双目视觉系统的特点,在人为定义的有限空间范围内对视觉系统进行标定与优化,提高了标定质量。(1) The calibration method proposed by the method of the present invention fully combines the characteristics of the biprism virtual binocular vision system, performs calibration and optimization of the visual system within an artificially defined limited space, and improves the calibration quality.

(2)本发明方法提出的基于结构参数优化的目标函数设计更加简单。将空间距离约束与极线约束结合在一起,满足了双棱镜虚拟双目系统在三维测量与二维对极几何约束方面的综合需求。(2) The objective function design based on structural parameter optimization proposed by the method of the present invention is simpler. Combining spatial distance constraints and epipolar constraints meets the comprehensive needs of the biprism virtual binocular system in terms of three-dimensional measurement and two-dimensional epipolar geometric constraints.

附图说明Description of drawings

图1是双棱镜虚拟双目系统的成像原理示意图;Figure 1 is a schematic diagram of the imaging principle of the biprism virtual binocular system;

图2是本发明实施例中的双棱镜虚拟双目标定方法实现的总体流程图;Figure 2 is an overall flow chart of the implementation of the biprism virtual dual target determination method in the embodiment of the present invention;

图3是本发明实施例中的构建空间测量场的示意图;Figure 3 is a schematic diagram of constructing a spatial measurement field in an embodiment of the present invention;

图4是本发明实施例中的结构参数非线性优化方法的流程图;Figure 4 is a flow chart of the nonlinear optimization method of structural parameters in the embodiment of the present invention;

图5是本发明实施例中的空间距离误差测量示意图;Figure 5 is a schematic diagram of spatial distance error measurement in the embodiment of the present invention;

图6是本发明实施例中的极线约束误差示意图,(a)图为对极几何示意图,(b)图为极线距离定义图。Figure 6 is a schematic diagram of the epipolar constraint error in the embodiment of the present invention. (a) is a schematic diagram of epipolar geometry, and (b) is a diagram of epipolar distance definition.

具体实施方式Detailed ways

为了使本申请的目的、技术方案及优点更加清晰且易于理解,下面结合附图和具体实施方式,对本发明作进一步说明。应当理解,此处描述的具体实施例仅用以解释本申请,并不用于限定本申请。In order to make the purpose, technical solutions and advantages of the present application clearer and easier to understand, the present invention will be further described below in conjunction with the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described here are only used to explain the application and are not used to limit the application.

本发明提供了一种用于双棱镜虚拟双目视觉系统的结构参数优化标定方法,整体标定流程如图2所示。本发明的整体流程如下:The present invention provides a structural parameter optimization calibration method for a biprism virtual binocular vision system. The overall calibration process is shown in Figure 2. The overall process of the present invention is as follows:

步骤1:搭建基于双棱镜折射的虚拟双目立体视觉系统。通过精密位移台调整相机与双棱镜的相对位置,使双棱镜相对平行于相机成像面以及使视觉系统获取分割均匀的左、右子图像。双棱镜虚拟双目视觉系统的成像原理如图1所示。获取的立体图像沿着图像中线均匀分割,得到两个像素尺寸相同的子图像。其左、右两个子图可被视为两个不同位置的左(VCL)、右(VCR)虚拟相机拍摄所得。Step 1: Build a virtual binocular stereoscopic vision system based on biprism refraction. The relative position of the camera and the biprism is adjusted through a precision shift stage so that the biprism is relatively parallel to the imaging surface of the camera and the visual system obtains evenly divided left and right sub-images. The imaging principle of the biprism virtual binocular vision system is shown in Figure 1. The acquired stereoscopic image is evenly divided along the center line of the image to obtain two sub-images with the same pixel size. The left and right sub-images can be regarded as taken by two virtual cameras on the left (VC L ) and right (VC R ) at different positions.

步骤2:将表面带有实心圆点的标定板随机放置在虚拟双目的共同视野的不同位置,获取立体标定板图像。Step 2: Randomly place the calibration plate with solid dots on the surface at different positions in the common field of view of the virtual binoculars to obtain the stereoscopic calibration plate image.

步骤3:沿着图像中线均匀分割步骤2所得的标定板图像,提取各个标定板上的圆心特征点的三维坐标与二维像素坐标/>以左虚拟相机为例,建立标定板特征点的世界坐标与像素坐标的映射关系:Step 3: Evenly divide the calibration plate image obtained in Step 2 along the center line of the image, and extract the three-dimensional coordinates of the center feature points on each calibration plate. and two-dimensional pixel coordinates/> Taking the left virtual camera as an example, establish the mapping relationship between the world coordinates and pixel coordinates of the calibration plate feature points:

其中,i与j分别代表拍摄的标定板序列与标定板上的特征点序列。KL为左虚拟相机的内参矩阵,为/>代表VCL的光心相对于第i个标定板的旋转矩阵与平移矩阵。此外,为校正镜头畸变,本实例引入传统的径向畸变与切向畸变于虚拟双目视觉系统中,畸变模型为:Among them, i and j respectively represent the photographed calibration plate sequence and the feature point sequence on the calibration plate. K L is the internal parameter matrix of the left virtual camera, for/> Represents the rotation matrix and translation matrix of the optical center of VC L relative to the i-th calibration plate. In addition, in order to correct lens distortion, this example introduces traditional radial distortion and tangential distortion into the virtual binocular vision system. The distortion model is:

其中为畸变点在归一化平面上的坐标,mnl(xlm,ylm)为去畸变点在归一化平面上的坐标。kl(1-2)与pl(1-2)各代表径向与切向畸变系数,/> in is the coordinate of the distortion point on the normalized plane, m nl (x lm , y lm ) is the coordinate of the dedistorted point on the normalized plane. k l(1-2) and p l(1-2) respectively represent the radial and tangential distortion coefficients,/>

步骤4:根据步骤3所述的虚拟相机成像模型,采用张正友标定法对左、右两个虚拟相机进行粗标定。粗标定可以得到包括两个虚拟相机各自的内参矩阵、畸变系数以及每幅标定板图像所对应的VCL与VCR光心之间的坐标转换关系,即:虚拟双目的结构参数。本发明将多组结构参数取平均值的结果作为系统标定的结构参数初值。Step 4: Based on the virtual camera imaging model described in Step 3, use Zhang Zhengyou’s calibration method to roughly calibrate the left and right virtual cameras. Coarse calibration can obtain the internal parameter matrices, distortion coefficients of each of the two virtual cameras, and the coordinate transformation relationship between the VC L and VC R optical centers corresponding to each calibration plate image, that is, the virtual binocular structural parameters. In the present invention, the average value of multiple sets of structural parameters is used as the initial value of the structural parameters for system calibration.

步骤5:如图3所示,将一个表面数量有10×21且规则排布的测量板垂直固定在一个高精度位移台上,测量板的位移方向相对垂直于相机的成像面。操控位移台使测量板每次移动1mm,共移动51次。测量板上的每一个圆心特征点被视为一个控制点,根据测量板的线性位移在虚拟双目的共同视场中构建了规格排布且相邻控制点距离可知的空间测量场。Step 5: As shown in Figure 3, a measuring plate with a surface number of 10×21 and a regular arrangement is vertically fixed on a high-precision displacement stage. The displacement direction of the measuring plate is relatively perpendicular to the imaging plane of the camera. The displacement stage is controlled to move the measurement plate 1mm each time for a total of 51 times. Each circle center feature point on the measurement plate is regarded as a control point. According to the linear displacement of the measurement plate, a spatial measurement field with a standardized arrangement and a known distance between adjacent control points is constructed in the common field of view of the virtual binoculars.

步骤6:均匀分割用于构建空间测量场的各组立体图像,根据步骤4所得的内参结果分别对左、右视图的测量板图像进行畸变矫正并提取控制点的像素坐标。根据控制点在测量板上的排列次序与测量板在三维空间中的位移次序获得空间控制场中控制点的匹配坐标。而后构建目标函数对结构参数进行非线性优化,具体流程如图4所示。Step 6: Uniformly segment each group of stereo images used to construct the spatial measurement field. According to the internal reference results obtained in step 4, perform distortion correction on the left and right view measurement plate images and extract the pixel coordinates of the control points. The matching coordinates of the control points in the spatial control field are obtained according to the arrangement order of the control points on the measurement plate and the displacement order of the measurement plate in the three-dimensional space. Then the objective function is constructed to perform nonlinear optimization of the structural parameters. The specific process is shown in Figure 4.

步骤7:采用线性三角法恢复出各个控制点以左相机光心为坐标系的三维坐标,以此建立基于空间距离误差的约束函数。空间距离误差的定义如图5所示,选取测量空间中的某一个控制点,测量该控制点与其所在测量板平面上的左斜上方向M1点、右斜上方向M2点、左斜下方向M3点、右斜下方向M4点以及相对于其所在测量板前向M5点、后向M6点这六条线段的空间距离。基于空间距离约束的误差函数设计为:Step 7: Use the linear trigonometric method to recover the three-dimensional coordinates of each control point with the optical center of the left camera as the coordinate system, thereby establishing a constraint function based on the spatial distance error. The definition of spatial distance error is shown in Figure 5. Select a certain control point in the measurement space, and measure the control point and the point M1 in the left diagonal upward direction, the point M2 in the right diagonal upward direction, and the left diagonal downward direction on the plane of the measurement plate. The spatial distances between point M3, point M4 in the right and downward direction, and six line segments relative to the forward point M5 and the backward point M6 of the measurement board where they are located. The error function based on spatial distance constraints is designed as:

Lq(q=1,2,……,6)代表指定控制点与其周围第q个控制点之间的实测长度,其长度的真值为Dq(q=1,2,……,6)。i与j分别代表测量板的位移次序与测量板上的控制点次序。Lq (q=1,2,...,6) represents the measured length between the specified control point and the qth control point around it, and the true value of its length is Dq (q=1,2,...,6). i and j respectively represent the displacement sequence of the measurement plate and the sequence of control points on the measurement plate.

步骤7:设计基于极线约束的约束函数。图6(a)展示了双目视觉系统的对极几何模型。VCL、VCR的光心是Ocl与Ocr,连接Ocl与Ocr的线段定义为双目系统的基线,其与左虚拟相机的成像面与右虚拟相机成像面的交点分别为左、右极点:el、er为匹配点对在像素坐标系上的齐次坐标。匹配点对应的左极线与右极线可分别被表示为:/>F代表虚拟双目系统的基础矩阵,即:F=R·S。Step 7: Design the constraint function based on epipolar constraints. Figure 6(a) shows the epipolar geometric model of the binocular vision system. The optical centers of VCL and VCR are Ocl and Ocr . The line segment connecting Ocl and Ocr is defined as the baseline of the binocular system. Its intersection points with the imaging surface of the left virtual camera and the imaging surface of the right virtual camera are respectively , right pole: e l , e r . is the homogeneous coordinates of the matching point pair on the pixel coordinate system. The left polar line and right polar line corresponding to the matching point can be expressed as:/> F represents the basic matrix of the virtual binocular system, that is: F=R·S.

R=[r11,r12,r13;r21,r22,r23;r31,r32,r33]代表虚拟双目的旋转矩阵,S代表平移矩阵T=[t1,t2,t3]的反对称矩阵:R=[r 11 , r 12 , r 13 ; r 21 , r 22 , r 23 ; r 31 , r 32 , r 33 ] represents the virtual binocular rotation matrix, S represents the translation matrix T=[t 1 , t 2 ,t 3 ] antisymmetric matrix:

理想情况下,指定点的匹配点应该在由指定点像素坐标所计算得到的极线上,即匹配点到其对应的极线的距离为0。但实际的情况并不如此,匹配点距离其对应极线总是存在一定的垂直距离,该垂直距离越小说明结构参数越符合当前立体视觉系统的对极几何约束。图6(b)展示了极线距离误差示意图,定义到极线ll与lr的距离分别为disl、disr。基于极线约束的误差函数设计为:Ideally, the matching point of the specified point should be on the epipolar line calculated from the pixel coordinates of the specified point, that is, the distance from the matching point to its corresponding epipolar line is 0. But this is not the actual situation. There is always a certain vertical distance between the matching point and its corresponding epipolar line. The smaller the vertical distance, the more the structural parameters are in line with the epipolar geometric constraints of the current stereo vision system. Figure 6(b) shows the schematic diagram of the epipolar distance error, defined The distances to the polar lines l l and l r are dis l and dis r respectively. The error function based on epipolar constraints is designed as:

其中,i、j分别代表测量板的位移次序与测量板上控制点阵列的排列次序。Among them, i and j respectively represent the displacement sequence of the measurement plate and the arrangement order of the control point array on the measurement plate.

步骤8:上述的步骤6、7构建了基于空间距离误差与极线约束的误差约束函数。将两种约束函数结合并构建目标函数为:Step 8: The above steps 6 and 7 construct an error constraint function based on spatial distance error and epipolar constraints. Combine the two constraint functions and construct the objective function as:

Hobj(x)=arg min(EL+EP)。H obj (x) = arg min ( EL + E P ).

在此,需要优化的参数有12个变量,即:x=(r11,r12,r13,r21,r22,r23,r31,r32,r33,t1,t2,t3)。采用非线性优化算子顺序二次规划法(Sequential Quadratic Programming,SQP)对目标函数进行求解。结构参数的非线性优化是一个迭代的过程,每次迭代求解得到的误差输出将进行一次迭代截止条件的判断。若不符合迭代截止条件则会将本次迭代输出的结构参数作为下次迭代的初值,且会根据当前的结构参数更新空间测量场中控制点的三维坐标。当目标函数输出达到截止条件时,则当前结构参数则认为达到了非线性优化的局部最优值。 Here , there are 12 variables that need to be optimized , namely : t3 ). The objective function is solved using the nonlinear optimization operator Sequential Quadratic Programming (SQP) method. The nonlinear optimization of structural parameters is an iterative process, and the error output obtained from each iteration will be used to judge the iteration cutoff conditions. If the iteration cutoff conditions are not met, the structural parameters output from this iteration will be used as the initial values for the next iteration, and the three-dimensional coordinates of the control points in the spatial measurement field will be updated based on the current structural parameters. When the objective function output reaches the cut-off condition, the current structural parameters are considered to have reached the local optimal value of nonlinear optimization.

本发明还提供了一种电子设备,包括:The invention also provides an electronic device, including:

存储器和处理器,所述存储器和所述处理器之间互相通信连接,所述存储器中存储有计算机指令,所述处理器通过执行所述计算机指令,从而执行如上述任一种方法。A memory and a processor. The memory and the processor are communicatively connected to each other. Computer instructions are stored in the memory, and the processor executes any of the above methods by executing the computer instructions.

本发明还提供了一种计算机存储介质,所述计算机可读存储介质上存储有计算机指令,所述计算机指令用于使计算机执行如上述任一种方法。The present invention also provides a computer storage medium. Computer instructions are stored on the computer-readable storage medium, and the computer instructions are used to cause the computer to execute any of the above methods.

本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。本申请实施例中的方案可以采用各种计算机语言实现,例如,面向对象的程序设计语言Java和直译式脚本语言JavaScript等。Those skilled in the art will understand that embodiments of the present application may be provided as methods, systems, or computer program products. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment that combines software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein. The solutions in the embodiments of this application can be implemented using various computer languages, such as the object-oriented programming language Java and the literal scripting language JavaScript.

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each process and/or block in the flowchart illustrations and/or block diagrams, and combinations of processes and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing device to produce a machine, such that the instructions executed by the processor of the computer or other programmable data processing device produce a use A device for realizing the functions specified in one process or multiple processes of the flowchart and/or one block or multiple blocks of the block diagram.

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer-readable memory that causes a computer or other programmable data processing apparatus to operate in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including the instruction means, the instructions The device implements the functions specified in a process or processes of the flowchart and/or a block or blocks of the block diagram.

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions may also be loaded onto a computer or other programmable data processing device, causing a series of operating steps to be performed on the computer or other programmable device to produce computer-implemented processing, thereby executing on the computer or other programmable device. Instructions provide steps for implementing the functions specified in a process or processes of a flowchart diagram and/or a block or blocks of a block diagram.

以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。The above-mentioned embodiments only express several implementation modes of the present application. The descriptions are relatively specific and detailed, but should not be construed as limiting the patent scope of the present invention. It should be noted that, for those of ordinary skill in the art, several modifications and improvements can be made without departing from the concept of the present application, and these all fall within the protection scope of the present application. Therefore, the protection scope of this patent application should be determined by the appended claims.

Claims (10)

1. A structural parameter optimization calibration method for a biprism virtual binocular vision system is characterized by comprising the following steps of:
s1: adjusting the relative position of the camera and the biprism to enable the biprism virtual binocular vision system to acquire left and right sub-images which are evenly divided;
s2: the calibration plates are placed at different positions under the depth of field range of the biprism virtual binocular vision system, and three-dimensional images of a plurality of groups of calibration plates are obtained;
s3: uniformly dividing an original stereoscopic image of a calibration plate to obtain two sub-images with the same size, respectively extracting characteristic points on a left sub-image and a right sub-image, roughly calibrating a biprism virtual binocular vision system, and obtaining an internal reference and a structural parameter initial value of a virtual binocular camera;
s4: the measuring plate is placed in front of the vision system, the measuring plate is moved for a plurality of times at a fixed distance, a stereoscopic image sequence is obtained, and a space measuring field composed of characteristic points with regular distribution is established;
s5: uniformly dividing a plurality of groups of stereoscopic images in a space measurement field, and carrying out distortion correction on sub-images corresponding to the left and right virtual cameras by utilizing the internal parameters of the virtual binocular cameras obtained in the step S3;
s6: extracting circle center characteristic points on the distorted measurement plate image, and matching the circle center characteristic points of the left and right subgraphs according to the arrangement sequence of the circle center characteristic point sequences to obtain characteristic point matching point pairs;
s7: combining internal parameters and structural parameters of the virtual binocular camera, and recovering three-dimensional coordinates of characteristic points in the space measurement field by a linear trigonometry;
s8: and constructing an objective function based on the space distance error and epipolar constraint, performing nonlinear iterative optimization on the structural parameters, and updating the three-dimensional coordinates of the feature points in the step S7 by using the structural parameters obtained by each iterative calculation.
2. The structural parameter optimization calibration method for a biprism virtual binocular vision system of claim, wherein the method comprises the following steps: the biprism virtual binocular vision system is a binocular stereoscopic vision image pair which is obtained by a single camera after being refracted by a biprism and can be equivalent to binocular stereoscopic vision images shot by two virtual cameras.
3. The structural parameter optimization calibration method for a biprism virtual binocular vision system according to claim 1, wherein: the specific steps of the S1 are as follows: the positions of the camera and the biprism are adjusted by the precise displacement table, so that the center line in the vertical direction of the image is approximately overlapped with the edge of the center of the biprism.
4. The structural parameter optimization calibration method for a biprism virtual binocular vision system according to claim 1, wherein: and (5) performing coarse calibration on the biprism virtual binocular vision system by using a Zhang Zhengyou calibration method.
5. The structural parameter optimization calibration method for a biprism virtual binocular vision system according to claim 1, wherein: the specific steps of the S4 are as follows: a measuring plate with a circular target arranged regularly on the surface is fixed on a moving platform of a high-precision electric displacement table, the spatial position of the measuring plate is relatively parallel to an imaging surface of a camera, and an operation controller enables the measuring plate to move a fixed distance for a plurality of times and directionally, and meanwhile, a three-dimensional image of the measuring plate at each displacement position is acquired.
6. The structural parameter optimization calibration method for a biprism virtual binocular vision system according to claim 1, wherein: the specific steps of the S8 are as follows:
s8.1: the fixed point distance length measurement in six directions is carried out on any one control point in the volume of the space measurement field, and an error function based on space distance constraint is designed as follows:
lq (q=1, 2, … …, 6) represents the measured length between the designated control point and the q-th feature point around it, the true value of the length is Dq (q=1, 2, … …, 6), i and j represent the displacement order of the measuring plate and the center of circle order of the feature point on the measuring plate, respectively;
s8.2: definition of the definitionIs the homogeneous coordinates of a group of matching points, and the two points are connected to the corresponding polar line l l And/l r Is dis respectively l 、dis r Constructing an error function based on epipolar constraint as follows:
wherein,is->Corresponding left pole line,/->Is->The corresponding right epipolar line, F, represents the basis matrix of the virtual binocular system, namely: f=r·s, r= [ R ] 11 ,r 12 ,r 13 ;r 21 ,r 22 ,r 23 ;r 31 ,r 32 ,r 33 ]Representing a virtual binocular rotation matrix, S representing a virtual binocular translation matrix t= [ T ] 1 ,t 2 ,t 3 ]Is an antisymmetric matrix of (a):
s8.3: combining two error functions to form an objective function:
H obj (x)=arg min(E L +E P )
the input variables of the objective function are structural parameters, namely: x= (r 11 ,r 12 ,r 13 ,r 21 ,r 22 ,r 23 ,r 31 ,r 32 ,r 33 ,t 1 ,t 2 ,t 3 ) And carrying out iterative solution on the structural parameters by combining SQP (Sequential Quadratic Programming) nonlinear optimization operators, wherein the iterative cutoff is based on the output local minimum value of the objective function.
7. The structural parameter optimization calibration method for the biprism virtual binocular vision system according to claim 4, wherein the structural parameter optimization calibration method is characterized by comprising the following steps of: the specific steps of the S3 are as follows:
s3.1: uniformly dividing an original stereoscopic image of the calibration plate along a central line in the vertical direction of the image to obtain two calibration plate sub-images I with the same pixel size L 、I R The left and right sub-images are respectively regarded as left and right virtual cameras VC L 、VC R Shooting the obtained product;
s3.2: world coordinates defining feature points on calibration plateAnd extracting the pixel coordinates corresponding to the feature points Establishing a mapping relation between 2D pixel coordinates and 3D world coordinates of the feature points on each calibration plate image, and taking a left virtual camera as an example, obtaining a conversion relation between the feature point world coordinates and the pixel coordinates:
wherein i and j respectively represent the shot calibration plate sequence and the characteristic point sequence on the calibration plate, K L Is an internal reference matrix of the left virtual camera,is->Representing VC L A rotation matrix and a translation matrix of the optical center of the (c) relative to the ith calibration plate;
s3.3: VC is obtained by adopting Zhang Zhengyou calibration method L ,VC R Is characterized by comprising an internal reference matrix, distortion coefficients and VC corresponding to each calibration plate image L With VC R Coordinate conversion relation between optical centers, namely: and the virtual double-purpose structure parameters are used as virtual double-purpose structure parameters for system calibration by taking the average value of the multiple groups of virtual double-purpose structure parameters.
8. The structural parameter optimization calibration method for the biprism virtual binocular vision system according to claim 6, wherein: the structural parameters obtained by each iteration are combined with the virtual camera internal parameters obtained during coarse calibration to update the three-dimensional coordinates of the control points in the space measurement field, and the new structural parameters and the three-dimensional coordinates of the feature point array are used for the next iteration calculation.
9. The structural parameter optimization calibration method for a biprism virtual binocular vision system according to claim 7, wherein: to correct lens distortion, conventional radial distortion and tangential distortion are introduced into a virtual binocular vision system, taking a left virtual camera as an example, the distortion model is defined as:
wherein the method comprises the steps ofAnd m is equal to nl (x lm ,y lm ) The coordinates of the distortion point and the de-distortion point on the normalized plane are respectively k l(1-2) And p is as follows l(1-2) Each representing a radial and tangential distortion coefficient, +.>
10. An electronic device, characterized in that: comprising a processor and a memory, said processors being communicatively coupled to each other, said memory having stored therein computer instructions, said processor executing the method of any of claims 1 to 9 by executing said computer instructions.
CN202311641278.2A 2023-12-04 2023-12-04 Structural parameter optimization calibration method for biprism virtual binocular vision system Pending CN117611684A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311641278.2A CN117611684A (en) 2023-12-04 2023-12-04 Structural parameter optimization calibration method for biprism virtual binocular vision system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311641278.2A CN117611684A (en) 2023-12-04 2023-12-04 Structural parameter optimization calibration method for biprism virtual binocular vision system

Publications (1)

Publication Number Publication Date
CN117611684A true CN117611684A (en) 2024-02-27

Family

ID=89951332

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311641278.2A Pending CN117611684A (en) 2023-12-04 2023-12-04 Structural parameter optimization calibration method for biprism virtual binocular vision system

Country Status (1)

Country Link
CN (1) CN117611684A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118298033A (en) * 2024-06-06 2024-07-05 江苏魔视智能科技有限公司 Parameter calibration method, device, equipment and storage medium of binocular camera
CN118691687A (en) * 2024-08-22 2024-09-24 中国科学院深海科学与工程研究所 An epipolar correction algorithm for a binocular underwater camera system

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118298033A (en) * 2024-06-06 2024-07-05 江苏魔视智能科技有限公司 Parameter calibration method, device, equipment and storage medium of binocular camera
CN118691687A (en) * 2024-08-22 2024-09-24 中国科学院深海科学与工程研究所 An epipolar correction algorithm for a binocular underwater camera system

Similar Documents

Publication Publication Date Title
CN109272570B (en) Space point three-dimensional coordinate solving method based on stereoscopic vision mathematical model
CN107063129B (en) A kind of array parallel laser projection three-dimensional scan method
CN110288642B (en) Three-dimensional object rapid reconstruction method based on camera array
CN117611684A (en) Structural parameter optimization calibration method for biprism virtual binocular vision system
Zhao et al. Calibration for stereo vision system based on phase matching and bundle adjustment algorithm
CN109767476A (en) A kind of calibration of auto-focusing binocular camera and depth computing method
CN108489395A (en) Vision measurement system structural parameters calibration and affine coordinate system construction method and system
CN112288826B (en) Calibration method and device of binocular camera and terminal
CN105654476B (en) Bi-objective determination method based on chaotic particle swarm optimization algorithm
CN108288291A (en) Polyphaser calibration based on single-point calibration object
CN109712232B (en) Object surface contour three-dimensional imaging method based on light field
CN112465912A (en) Three-dimensional camera calibration method and device
CN112215880B (en) Image depth estimation method and device, electronic equipment and storage medium
CN113592721B (en) Photogrammetry method, device, equipment and storage medium
CN114359406A (en) Calibration of auto-focusing binocular camera, 3D vision and depth point cloud calculation method
CN110345921A (en) Stereoscopic fields of view vision measurement and vertical axial aberration and axial aberration bearing calibration and system
CN110044301A (en) Three-dimensional point cloud computing method based on monocular and binocular mixed measurement
CN103473758A (en) Secondary calibration method of binocular stereo vision system
CN113160393B (en) High-precision three-dimensional reconstruction method and device based on large depth of field and related components thereof
WO2023040095A1 (en) Camera calibration method and apparatus, electronic device, and storage medium
CN116957987A (en) Multi-eye polar line correction method, device, computer equipment and storage medium
CN115359127B (en) A polarization camera array calibration method suitable for multi-layer medium environment
CN114926538B (en) External parameter calibration method and device for monocular laser speckle projection system
CN114998449B (en) A high-precision calibration method for zoom binocular vision measurement system
CN112396663A (en) Visual calibration method, device, equipment and medium for multi-depth camera

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