一种基于射影几何的线结构光标定方法与流程

本发明涉及一种视觉结构光测量领域，具体涉及一种基于射影几何的线结构光标定方法。
背景技术：
：随着工业相机与计算机的发展，视觉检测在工业检测中被大量使用，其中结构光视觉系统由于无需粘贴或定制标志点，包含信息多且精度高而应用更为广泛，最常见的即为线结构光与单目相机结合的测量系统，在使用系统进行检测前，必须标定光平面在相机坐标系下的表达式，进而由相机获取承载光条信息的被测物图像，透视成像模型结合光平面方程，实现三维测量。常见的线结构光标定方法有拉丝法、锯齿靶法、基于交比不变的标定方法、基于神经网络的标定方法。拉丝法是根据光平面投射到空间分布不共面的细线上会形成亮点，采用其他仪器测出亮点在空间中的三维坐标，并结合亮点在相机中成像坐标得到光平面方程，该法需借助其他仪器且细丝上光点产生的漫反射会带来较大误差；锯齿靶法是采用一个锯齿状的靶标与一个一维工作平台实现较高精度的标定，但该法需要调节工作台使得光平面与靶标棱线垂直，操作较复杂，且对于靶标要求较高；基于交比不变的方法是利用透视投影共线点满足交比不变的性质，计算激光线上点在靶标坐标系下的三维坐标，再通过相机与每幅图像的相对外参转换到相机坐标系下，最后进行平面拟合，此法操作简单但外参矩阵的误差会传递导致标志点在相机坐标系下误差，并且获取的光平面点较少，基于神经网络的标定方法是借助光点瞄准装置以及高精度移动平台获取特征点的三维坐标，利用样本数据训练bp神经网络，优化寻找光平面方程系数，该法精度高，但需要其他高精度仪器，且操作较复杂不适合现场标定。技术实现要素：本发明的目的在于针对现有标定线结构光技术需借助其他仪器或精度方面的不足，提出一种基于射影几何的线结构光标定方法，采用普通的棋盘格靶标，借助射影几何中空间直线方向向量与成像消隐点的对应关系以及交比不变的性质，结合透视成像模型完成标定，实现简单易行、精度较高的线结构光平面标定。本发明采用如下技术方案：一种基于射影几何的线结构光标定方法，包括相机、线激光器、计算机及棋盘格，具体步骤如下：第一步：固定相机与线激光器位置，将棋盘格作为靶标置于相机视场内，激光器打在靶标形成条纹，相机采集包含条纹的靶标图像，并输入计算机；第二步：对采集的靶标图像进行内参与畸变系数的标定，得到内参矩阵k、畸变系数向量kc；第三步：利用内参矩阵与畸变系数校正图像，经高斯滤波去除噪点和阈值化后，利用canny算子提取激光条纹边缘，从而利用最小二乘法拟合得到激光条纹在图像坐标系下的方程l0，用[a0b0c0]表示；第四步：利用亚像素角点检测的方法在截取的roi区域中确定第一排棋盘格角点a、b、c以及第四排棋盘格角点aa、bb、cc和与a、b点同列的某点m、n的坐标，进而得到直线ab方程l1：[a1b1c1]、直线aabb方程l2：[a2b2c2]、直线am方程l3：[a3b3c3]、直线bn方程l4：[a4b4c4]，并求l0与l1的交点p坐标(u1,v1)、l0与l2的交点pp点坐标(u2,v2)，p、pp点相机坐标系下坐标满足：其中，i＝1,2；第五步：由直线l1与l2方程得到对应棋盘格平面的消隐点vanish_point1，同样由l3与l4方程得到另一消隐点vanish_point2，进一步得到消隐线vanish_l方程，结合激光条纹l0方程确定对应激光条纹的消隐点vanish_point0齐次坐标，进而得到光条在相机坐标系下的三维方向向量，并进一步得到p、pp点相机坐标系下坐标满足：第六步：利用射影几何交比不变的性质通过已知点的图像距离与靶标距离得到p、pp点欧式距离，进一步确认p、pp点相机坐标系下坐标满足：第七步：根据方程(4)、(6)、(8)确定点p、pp在相机坐标系下的三维坐标；第八步：重复第一步到第七步，改变靶标位置，采集n张靶标图像对采集的n张图像重复上述操作，得到2n个光平面上的特征点，由此进行拟合，通过优化目标函数，确定光平面方程系数，确定光平面方程，完成标定。所述靶标为8×9的黑白相间棋盘格，每一格边长尺寸为dx＝25mm，采集图像过程中保证靶标不与相机成像平面平行。所述第二步中，采用张氏标定法对采集的靶标图像进行内参与畸变系数的标定，得到内参矩阵k及畸变系数向量kckc＝[r1r2p1p20](2)其中，fx为x轴归一化焦距，fy为y轴归一化焦距，(u0,v0)为相机光轴中心的图像坐标，r1、r2为一阶、二阶径向畸变系数，p1、p2为一阶、二阶切向畸变系数；另外，根据成像模型：其中，(u,v)为某点图像坐标，(xc,yc,zc)为对应点在相机坐标系下坐标。在射影几何中，空间中一条直线的消隐点由平行于该直线并过相机中心的射线与图像平面交点得到，因此消隐点仅依赖于直线的方向，从而直线l1与l2映射为同一消隐点，l3与l4映射为同一消隐点，由各自消隐点的齐次坐标叉乘即可得到消隐线方程，在提高精度的目的下，提取多条平行直线图像方程，采用最小二乘法以点到各直线距离和为目标函数优化消隐点坐标。所述p、pp点欧式距离：其中，dis为p、pp点欧式距离，dx为靶标方格边长的标称尺寸，λ为直线l1与直线l2间方格列数，距离为欧式不变量。本发明的有益效果：(1)本发明方法仅需借助一个棋盘格靶标，计算内参与标定光平面可采用同一组图像，且对靶标与光平面的相对位姿无特殊要求，操作简单易行，适用于现场标定。(2)本发明通过消隐点与交比构造等式约束求解特征点三维坐标，无需借助外参进行多次坐标系间变换，减小误差的累积，并可获得较多的特征点坐标，提高光平面拟合精度。附图说明图1为本发明的总体流程图；图2为本发明射影几何消隐点原理图；图3为本发明射影几何交比不变原理图；图4为基于射影几何的线结构光标定方法靶标示意图；图5为本发明拟合平面图。具体实施方式下面结合实施例及附图，对本发明作进一步地详细说明，但本发明的实施方式不限于此。实施例如图1所示，一种基于射影几何的线结构光标定方法，包括第一步，固定相机与线激光器位置后，将棋盘格靶标置于相机视场中，棋盘格作为标定板，线激光器打在标定板上形成条纹，相机采集图像并保存入计算机，调整靶标位置姿态，重复操作，采集n(20～30)张图像；第二步采用张氏标定法对采集的靶标图像进行内参与畸变系数的标定，得到内参矩阵k、畸变系数向量kc：kc＝[r1r2p1p20](2)其中，fx为x轴归一化焦距，fy为y轴归一化焦距，(u0,v0)为相机光轴中心的图像坐标，r1、r2为一阶、二阶径向畸变系数，p1、p2为一阶、二阶切向畸变系数；另外，根据成像模型：其中，(u,v)为某点图像坐标，(xc,yc,zc)为对应点在相机坐标系下坐标；第三步利用内参矩阵与畸变系数校正图像，经高斯滤波去除噪点并阈值化后，利用canny算子提取激光条纹边缘，从而利用最小二乘法拟合得到激光条纹在图像坐标系下的方程l0：[a0b0c0](表示a0u+b0v+c0＝0,下同)，利用亚像素角点检测的方法在截取的roi区域中确定第一排棋盘格角点a、b、c以及第四排棋盘格角点aa、bb、cc和与a、b点同列的某点m、n的坐标，进而得到直线ab方程l1：[a1b1c1]、直线aabb方程l2：[a2b2c2]、直线am方程l3：[a3b3c3]、直线bn方程l4：[a4b4c4]，并求交点得到点p坐标(u1,v1)、pp点坐标(u2,v2)，由式(3)可知p、pp点相机坐标系下坐标满足(i＝1,2)：roi(regionofinterest感兴趣区域)可通过先粗略读取角点位置(画图软件鼠标移到相应位置即可读取)后给定区域长宽得到。通过在roi内的图像采用亚像素角点检测即可检测到角点的准确位置，通过计算灰度梯度最大值及方向即可。第四步由直线l1与l2方程得到对应棋盘格平面的消隐点vanish_point1，同样由l3与l4方程得到另一消隐点vanish_point2，进一步得到消隐线vanish_l方程，结合激光条纹l0方程确定对应激光条纹的消隐点vanish_point0齐次坐标，进而由下式得到光条在相机坐标系下的三维方向向量：d＝k-1·vanish_point0＝[d1d2d3]t(5)从而，p、pp点相机坐标系下坐标满足：第五步利用射影几何交比不变的性质通过已知点的图像距离与靶标距离得到p、pp点欧式距离，采用下式计算：其中，dis为p、pp点欧式距离，dx为靶标方格边长的标称尺寸，λ为直线l1与直线l2间方格列数，距离为欧式不变量，从而p、pp点相机坐标系下坐标满足：第六步联立方程(4)、(6)、(8)可确定点p、pp在相机坐标系下的三维坐标，对采集的n张图像重复上述操作，得到2n个光平面上的特征点，由此进行拟合，通过优化目标函数，确定光平面方程系数：确定光平面方程为ax+by+cz+d＝0，完成标定。如图2所示，在射影几何中，空间中一条直线的消隐点由平行于该直线并过相机中心的射线与图像平面交点得到，因此消隐点仅依赖于直线的方向，从而直线l1与l2映射为同一消隐点vanish_point1，l3与l4映射为同一消隐点vanish_point2，由各自消隐点的齐次坐标叉乘即可得到消隐线方程vanish_l，在提高精度的目的下，可提取多条平行直线图像方程，采用最小二乘法以点到各直线距离和为目标函数优化消隐点vanish_point0坐标。如图3及图4所示，在射影几何中，对于共线上的4点满足交比不变的特性，即如图3所示的共线上的点a、b、c、p及其成像在图像平面的点a、b、c、p满足式(7)关系，通过图像平面点的坐标计算交比，进而可以得到bp、bbpp的欧式距离，进一步由已知的棋盘格尺寸可以计算得到p与pp间欧式距离。在本实施例中，采用靶标为8×9的黑白相间棋盘格，每一格边长尺寸为dx＝25mm，λ＝4，采集图像过程中保证靶标不与相机成像平面平行。相机采用德国balser工业ccd相机，分辨率为1600×1200，采用张氏标定法得到的相机内参如下表1所示：表1相机内参与畸变系数fxfyu0(pixel)v0(pixel)k1k2p1p22687.512687.29749.27568.50-0.06690.12170.00110.0003通过所述方法求得光平面特征点(p、pp)坐标如下表2所示，表2特征点坐标pointxpointypointzpointxpointypointz-23.3089-108.987859.7911-21.2956-92.3189864.1578-6.43272-0.03872894.1116-12.94426.75672880.302-29.3592-79.2536850.244-3.96281-74.4572904.2096-11.550713.24586885.7533-5.9071-11.8413898.810418.16498-81.4566950.2754-74.2995-72.0592751.575424.8643617.70306961.902-49.010316.4036803.99180.833039-93.7765910.5384-109.396-76.1727675.333122.84769-3.64605956.6095-81.755410.52019733.047882.39167-43.88351087.478-59.3109-86.6249780.246254.7272130.388511024.975-40.4145.037658819.734645.25209-41.96961006.506-60.7632-74.1962780.92115.760826.89107940.2586-44.864820.10983812.739637.04282-102.419991.2281-73.4591-91.1901754.592418.21101-3.04706947.6638-51.748911.06049798.2864-21.5025-71.9344865.8015-57.9157-67.8432783.9389-7.4904923.08991893.6505-41.652626.78855817.1374-64.3855-33.8113768.7468-98.1294-98.4791698.731-43.063453.23104813.125-79.1065-3.97422737.9853-73.4186-53.2433750.9705-97.5141-80.1515702.1413-53.102343.0153792.8189-86.846217.32126721.888最小二乘法拟合得到的平面方程为：2.3686x-0.0275y-1.1005z+1000＝0(10)特征点分布及拟合平面示意如图5所示。本发明通过图像处理得到相关直线与交点的信息，利用平行直线的齐次坐标表达得到对应靶标平面的消隐线方程，进而得到对应激光条纹的消隐点坐标，进一步转换得到激光条纹在相机坐标系下的方向向量；利用交比不变的射影性质得到激光条纹上两特征点的欧式距离，结合方向向量与成像模型信息，联立方程组得到特征点在相机坐标系下三维坐标，对所有图像特征点坐标利用最小二乘进行平面拟合，得到光平面方程，完成线结构光标定。上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受所述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。当前第1页12