CN110276808A - 一种单相机结合二维码测量玻璃板不平度的方法 - Google Patents
一种单相机结合二维码测量玻璃板不平度的方法 Download PDFInfo
- Publication number
- CN110276808A CN110276808A CN201910501700.1A CN201910501700A CN110276808A CN 110276808 A CN110276808 A CN 110276808A CN 201910501700 A CN201910501700 A CN 201910501700A CN 110276808 A CN110276808 A CN 110276808A
- Authority
- CN
- China
- Prior art keywords
- camera
- point
- image
- matrix
- dimensional code
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/30—Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06K—GRAPHICAL DATA READING; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
- G06K7/00—Methods or arrangements for sensing record carriers, e.g. for reading patterns
- G06K7/10—Methods or arrangements for sensing record carriers, e.g. for reading patterns by electromagnetic radiation, e.g. optical sensing; by corpuscular radiation
- G06K7/14—Methods or arrangements for sensing record carriers, e.g. for reading patterns by electromagnetic radiation, e.g. optical sensing; by corpuscular radiation using light without selection of wavelength, e.g. sensing reflected white light
- G06K7/1404—Methods for optical code recognition
- G06K7/1408—Methods for optical code recognition the method being specifically adapted for the type of code
- G06K7/1417—2D bar codes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06K—GRAPHICAL DATA READING; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
- G06K7/00—Methods or arrangements for sensing record carriers, e.g. for reading patterns
- G06K7/10—Methods or arrangements for sensing record carriers, e.g. for reading patterns by electromagnetic radiation, e.g. optical sensing; by corpuscular radiation
- G06K7/14—Methods or arrangements for sensing record carriers, e.g. for reading patterns by electromagnetic radiation, e.g. optical sensing; by corpuscular radiation using light without selection of wavelength, e.g. sensing reflected white light
- G06K7/1404—Methods for optical code recognition
- G06K7/1439—Methods for optical code recognition including a method step for retrieval of the optical code
- G06K7/1443—Methods for optical code recognition including a method step for retrieval of the optical code locating of the code in an image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06K—GRAPHICAL DATA READING; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
- G06K7/00—Methods or arrangements for sensing record carriers, e.g. for reading patterns
- G06K7/10—Methods or arrangements for sensing record carriers, e.g. for reading patterns by electromagnetic radiation, e.g. optical sensing; by corpuscular radiation
- G06K7/14—Methods or arrangements for sensing record carriers, e.g. for reading patterns by electromagnetic radiation, e.g. optical sensing; by corpuscular radiation using light without selection of wavelength, e.g. sensing reflected white light
- G06K7/1404—Methods for optical code recognition
- G06K7/1439—Methods for optical code recognition including a method step for retrieval of the optical code
- G06K7/1452—Methods for optical code recognition including a method step for retrieval of the optical code detecting bar code edges
-
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Electromagnetism (AREA)
- General Health & Medical Sciences (AREA)
- Toxicology (AREA)
- Artificial Intelligence (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种单相机结合二维码测量玻璃板不平度的方法,利用单相机获得玻璃板的图像,通过精确识别玻璃板上角点坐标计算的二维码(如QR码)阵列,从而建立图像点和空间三维点的关系,来确定物体三维信息。整个测量过程包括相机标定,求取相机内外参数;基于玻璃板反射的相机外参数求解;对二维码先计算***轮廓中心距,作为定位点,然后根据轮廓顶点关系确定二维码的四个顶点坐标;将采集图像上的变形角点坐标与已知角点坐标拟合求得反射对应关系,根据其对应关系来求玻璃板的每点深度表达式,恢复表面形貌,分析其不平度。本发明采用单相机结合二维码来实现具有反射特性的物体重建,可减小测量***体积,降低测量成本,操作方便。
Description
技术领域
本发明涉及视觉测量技术领域,尤其涉及一种单相机结合二维码测量玻璃板不平度的方法。
背景技术
物体表面的三维形貌是物体表面最基本的几何性质,对物体表面形貌信息的获取在工业生产制造中具有重要的意义,所以如何高效、准确地获取物体表面点的三维数据是现代工业发展中一个需要解决的难题。近年来,随着传感技术,图像处理以及计算机技术的发展,具有测量快速、准确,且可测量范围很大等优点的非接触式三维重建方法,在物体表面外形检测领域中得到了非常广泛的应用。
目前常用的三维测量方式包括使用双目相机或者单相机加一定的光学条件。在使用双目相机测量时,对于两个相机的摆放位置有一定要求,这会使得整个测量***的体积较大。并且大多数测量物体都是刚性物体,且多为漫反射表面。
发明内容
本发明目的就是为了弥补已有技术的缺陷,提供一种单相机结合二维码测量玻璃板不平度的方法。
本发明是通过以下技术方案实现的:
本发明利用单相机获得玻璃板的图像,通过精确识别玻璃板上角点坐标可计算的二维码(如QR码)阵列,从而建立图像点和空间三维点的关系,来确定物体三维信息。整个测量过程包括相机标定,求取相机内外参数;基于玻璃板反射的相机外参数求解;对二维码先计算***轮廓中心距,作为定位点,然后根据轮廓顶点关系确定二维码的四个顶点坐标;将采集图像上的变形角点坐标与已知角点坐标拟合求得反射对应关系,根据其对应关系来求玻璃板的每点深度表达式,恢复表面形貌,从而可分析其不平度。
一种单相机结合二维码测量玻璃板不平度的方法,包括以下步骤:
S1:相机标定,得到相机内参数和外参数初始值,建立图像像素的位置坐标与场景点所处空间位置之间的关系。所述相机的内参数包含相机的有效焦距、光心及描述相机像平面u轴与v轴的相互垂直程度的参数。所述外参数是指旋转矩阵和平移向量。
S2:基于玻璃板反射的相机外参数求解,得到相机与投影设备的位置关系。
S3:玻璃板反射角点坐标已知的二维码阵列,用相机采集图像。对二维码先计算***轮廓中心距,作为定位点,然后根据轮廓顶点关系确定二维码的四个顶点坐标,提取图像中所有二维码中的角点。
S4:将检测到的特征点与标记点进行拟合,根据所得表达式求取待测面每点深度,从而求得玻璃板的表面形貌。
所述的S1的具体过程如下:
S1.1:建立空间三维点与相机平面上二维点映射关系,空间任意一点P在图像上的成像位置可以用针孔模型近似表示。
空间点P在世界坐标系和相机坐标系下的齐次坐标分别为(xw,yw,zw,1)和(x,y,z,1),可以写为如下表达式:
其中M2表示世界坐标系和相机坐标系之间的转换关系,为4×4的矩阵;
将相机坐标系内物点P以理想小孔模型映射到图像坐标系,表达式如下:
其中(X,Y)为P点对应的图像坐标,f代表相机的焦距;
相机CMOS靶面上,像元水平方向尺寸用dx表示,垂直方向尺寸用dy表示,则理想小孔模型下的像点p,在图像坐标系下物理坐标到像素坐标之间的转化可用如下表达式表示:
由公式(1)至(3)可得到,一物点P从世界坐标系经过忽略畸变的小孔模型,投影到像面坐标系下的数学表达式如下:
其中内参数矩阵fx、fy是单目相机u轴、v轴上的两个有效焦距;γ是用来描述相机像平面u轴与v轴的相互垂直程度的参数,通常令γ=0,(u0,v0)为光心坐标,为空间点P在世界坐标系下的齐次坐标。
相机外参数中R是旋转矩阵,T是平移向量。
S1.2:对相机进行标定,获得相机内外参数。通过计算单应性矩阵Hi求解出相机内参数矩阵A和相机相对于不同位置靶标的旋转矩阵Ri,以及平移向量Ti(i=1,2…..N),利用径向畸变系数与切向畸变系数k1,k2,p1,p2校正理想模型,最后使用L-M优化算法对所有参数的初始值进行优化,提高精度。
平面靶标上角点P在世界坐标系下的坐标记为P=[xw,yw,zw]T,经过相机数学模型成像后的像点坐标记为p=[u,v]T,将二者的坐标写作齐次坐标为和根据公式(4)可得出标志角点P与其像点p在相机成像数学模型中的转化关系:
上式中s为一任意的非零尺度因子,A为相机内参数矩阵,R是靶标坐标系与相机坐标系之间的旋转矩阵,T是靶标坐标系到相机坐标系的平移向量。
为计算方便,令标志点坐标的z值为零,即zw=0,令用单应性矩阵H描述靶标上点P到其像点p的变换,其中H是一个3×3的矩阵:
如果标定所用的平面靶标含有n个标定点,每个标定点可以列出一个关系式,所有关系式联立可以写成其中L是2n×9的系数矩阵。当所制备的平面靶标上标定特征点的个数n>4时,是一个超定方程组,利用SVD分解算法能够求解其最优解由能够解出H。
标定过程中相机总共采集了N个靶标空间不同位置的图像,总共有N个单应性矩阵Hi(i=1,2…..N)与每个空间位置的靶标相对应,通过计算N个单应性矩阵可以求出摄像机的内参数矩阵A,以及相机相对于每个靶标的旋转矩阵Ri以及平移向量Ti(i=1,2…..N)。根据旋转矩阵R的正交性可以得到内参数矩阵A的两个基本约束条件:
由式(8)可以计算A-TA-1,并令记单应性矩阵H第i行为hi=[hi1hi2hi3]T,根据计算可以得到:
通过SVD分解法求出最优的值,由于未知数的数量为6,张正友标定法至少需要3幅标定图像,即N≥3。根据的值可以很方便的求解出相机的内参矩阵A:
求得单应性阵H=[h1h2h3]与内参矩阵A之后,由(11)可求得相机相对靶标位置的外参数:
由于实验中存在误差,由公式(11)计算得到的一般并不具有旋转矩阵的正交性,因此,Q只能作为相机参数的初始估计值求解最佳旋转矩阵R,优化的目标函数为
为提高精度,对相机参数进行优化,利用最优化方法中的L-M优化方法,建立目标函数如下:
其中mij表示相机图像坐标系下实验获得的坐标,表示相机图像坐标系下经过理论计算的坐标。
所述的S2的具体过程如下:
S2.1:建立镜面反射数学模型。
镜面到屏幕参考点P的距离与镜面到屏幕参考点镜像的距离是相等的,由此可得公式:
在相机坐标系下,nj是镜面π在j位置下的表面法向量,dj是镜面在位置j下到相机坐标系原点的距离,是屏幕上第i个点到镜面位置j的距离。
到相机坐标系原点的距离可以表达为故可得以下公式:
由此,屏幕参考点与参考点镜像之间的关系可以写成:
其中S是描述镜面反射的矩阵:
S2.2:基于镜面反射的相机外参数求解。
用R表示镜像屏幕坐标系与相机坐标系之间的旋转矩阵,T表示二者之间的平移矩阵,则将投影模型与镜面反射模型相结合可以得到屏幕坐标系到相机坐标系的变换关系:
t=2djRnj+T (19)
其中与t分别为相机相对屏幕坐标系的旋转矩阵与平移矩阵,这种经过反射的旋转矩阵被称之为非正规旋转矩阵。
所述的S3的具体过程如下:
S3.1在玻璃板上投射二维码阵列;
S3.2对拍摄的玻璃板图像进行处理,转为灰度图像;
S3.3对图像进行校正、去噪、图像增强处理,使特征信息更明显;
S3.4二维码中心点提取,包含行列信息;
S3.5计算得到所有二维码四个角点的坐标信息,提取图像中所有二维码中的角点。
所述的S4的具体过程如下:
S4.1根据之前标定求得的相机内外参数,以及基于玻璃板的相机外参,得到目标位置与姿态的六个参数;
S4.2根据之前提到的小孔成像模型,由入射光线和反射光线的关系求得待测物上点的法向量;
S4.3由投影点在相机坐标系下的坐标和图像点在相机坐标系下的坐标拟合求得反射对应关系,而后求得待测物上点的深度函数g(x,y);
S4.4对于图像中的每点求取深度g,再与图像点相乘得到待测物空间三维点信息;
S4.5得到被测物三维点坐标信息后,可用积分方式得到面型,而后用重建面中的高度差来衡量其不平度。
本发明的优点是:本发明采用单目相机来实现三维测量,降低了成本和测量***的体积,并且便于操作。玻璃板反射角点坐标可求得的二维码阵列,所以可由单张图像进行处理得到物体三维信息,降低了测量***的复杂度。该方法可以通过对待测面的拟合,由表面最高点和最低点建立两个平面,从而用它们的高度差来表示玻璃板不平度。
附图说明
图1本发明的方法流程图。
图2本发明中相机成像模型。
图3本发明中相机标定示意图。
图4本发明中基于镜面的标定示意图。
图5本发明中二维码阵列示意图。
具体实施方式
如图1所示,一种单相机结合二维码测量玻璃板不平度的方法,包含了以下几个细分步骤:
步骤1:相机标定
本发明的相机标定采用张正友二维靶标标定法。设定整个工作范围为1~3m,为了满足测量要求,在考虑视场、焦距等各参数之后,选用USB2.0接口CMOS工业数字相机和焦距为16mm的固定焦距镜头。
步骤1.1参见图2所示,建立二维坐标和三维坐标之间的映射关系,空间任意一点P在图像上的成像位置可以用针孔模型近似表示。
空间点P在世界坐标系和相机坐标系下的齐次坐标分别为(xw,yw,zw,1)和(x,y,z,1),可以写为如下表达式:
其中M2表示世界坐标系和相机坐标系之间的转换关系,为4×4的矩阵;
将相机坐标系内物点P以理想小孔模型映射到图像坐标系,表达式如下:
其中(X,Y)为P点对应的图像坐标,f代表相机的焦距;
相机CMOS靶面上,像元水平方向尺寸用dx表示,垂直方向尺寸用dy表示,则理想小孔模型下的像点p,在图像坐标系下物理坐标到像素坐标之间的转化可用如下表达式表示:
由公式(1)至(3)可得到,一物点P从世界坐标系经过忽略畸变的小孔模型,投影到像面坐标系下的数学表达式如下:
其中内参数矩阵fx、fy是单目相机u轴、v轴上的两个有效焦距;γ是用来描述相机像平面u轴与v轴的相互垂直程度的参数,通常令γ=0,(u0,v0)为光心坐标,为空间点P在世界坐标系下的齐次坐标。
相机外参数中R是旋转矩阵,T是平移向量。
步骤1.2对相机进行标定,获得相机内外参数。如图3所示,具体按以下步骤进行:
(1)已靶标左上第一个特征角点为坐标原点,靶标所在平面为XY平面,Z轴垂直于XY平面,使其符合右手坐标系。
(2)调节相机焦距及光圈,使其拍摄图像尽可能清晰,避免出现曝光。
(3)拍摄靶标,手持靶标在相机拍摄视野内变换位置,尽量将物***于图像中间,注意要有旋转和平移,不能只在同一平面内移动。
(4)采集15幅靶标图像,利用MATLAB中程序对其标定,初次标定发现有些参与标定的靶标图片重投影误差过大,这可能是角点提取有误,或者手持靶标进行拍摄过程中出现抖动,将重投影误差过大的靶标图片逐一剔除,进行再次标定,最终留下了10幅图像作为双目相机标定的图像,求得相机内外参数。
通过计算单应性矩阵Hi求解出相机内参数矩阵A和相机相对于不同位置靶标的旋转矩阵Ri,以及平移向量Ti(i=1,2…..N),利用径向畸变系数与切向畸变系数k1,k2,p1,p2校正理想模型,最后使用L-M优化算法对所有参数的初始值进行优化,提高精度。
步骤2:玻璃板标定
步骤2.1参见图4所示,建立镜面反射数学模型。具体按以下步骤进行:
(1)在相机标定好后,接着进行玻璃板在相机下的标定。首先投射编程产生的与相机分辨率相同棋盘格靶标图像,确定棋盘格靶标的宽度和高度。
(2)完成在显示靶标图像后,利用相机拍摄玻璃板上屏幕靶标在三个位置下的镜像,然后根据之前标定所得的相机的内参数,利用OpenCV的相关函数计算相机相对于这三个镜像靶标的外参数,继而使用基于镜面的相机外参标定算法编写程序求解相机相对于真实屏幕靶标位置的外参数。
(3)利用基于镜面的相机外参数标定算法求得了相机相对于真实屏幕靶标位置的外参数。
镜面到屏幕参考点P的距离与镜面到屏幕参考点镜像的距离是相等的,由此可得公式:
在相机坐标系下,nj是镜面π在j位置下的表面法向量,dj是镜面在位置j下到相机坐标系原点的距离,是屏幕上第i个点到镜面位置j的距离。
到相机坐标系原点的距离可以表达为故可得以下公式:
由此,屏幕参考点与参考点镜像之间的关系可以写成:
其中S是描述镜面反射的矩阵:
步骤2.2基于镜面反射的相机外参数求解。
用R表示镜像屏幕坐标系与相机坐标系之间的旋转矩阵,T表示二者之间的平移矩阵,则将投影模型与镜面反射模型相结合可以得到屏幕坐标系到相机坐标系的变换关系:
t=2djRnj+T (11)
其中与t分别为相机相对屏幕坐标系的旋转矩阵与平移矩阵,这种经过反射的旋转矩阵被称之为非正规旋转矩阵。
步骤3:投射编码产生的二维码阵列并采集玻璃板上镜像。所用二维码阵列如图5所示,具体按以下步骤进行:
(1)利用行列坐标作为二维码编码所存储的信息,产生阵列码,投射到墙面上,相机拍摄得到图像。
(2)对拍摄的玻璃板图像进行处理,转为灰度图像。
(3)对图像进行校正、去噪、图像增强处理,使特征信息更明显。
(4)二维码中心点提取,包含行列信息。
(5)计算得到所有二维码四个角点的坐标信息,提取图像中所有二维码中的角点。
步骤4:检测特征点,根据图像与与参考面之间的关系,求得待测物表面信息,具体按以下步骤进行:
(1)根据之前标定求得的相机内外参数,以及基于玻璃板的相机外参,得到目标位置与姿态的六个参数;
(2)根据之前提到的小孔成像模型,由入射光线和反射光线的关系求得待测物上点的法向量
(3)由投影点在相机坐标系下的坐标和图像点在相机坐标系下的坐标拟合求得反射对应关系m(x,y),而后求得待测物上点的深度函数g(x,y);
(4)对于图像中的每点求取深度g,再根据p(x,y)=g(x,y)v得到待测物空间三维点信息;
(5)得到被测物三维点坐标信息后,可用积分方式得到面型,而后用重建面中的高度差来衡量其不平度。
Claims (5)
1.一种单相机结合二维码测量玻璃板不平度的方法,其特征在于:包括以下步骤:
S1:相机标定,得到相机内参数和外参数初始值,建立图像像素的位置坐标与场景点所处空间位置之间的关系;所述的相机的内参数包含相机的有效焦距、光心及描述相机像平面u轴与v轴的相互垂直程度的参数,所述外参数是指旋转矩阵和平移向量;
S2:基于玻璃板反射的相机外参数求解,得到相机与投影设备的位置关系;
S3:玻璃板反射角点坐标已知的二维码阵列,用相机采集图像,对二维码先计算***轮廓中心距,作为定位点,然后根据轮廓顶点关系确定二维码的四个顶点坐标,提取图像中所有二维码中的角点;
S4:将检测到的特征点与标记点进行拟合,根据所得表达式求取待测面每点深度,从而求得玻璃板的表面形貌。
2.根据权利要求1所述的一种单相机结合二维码测量玻璃板不平度的方法,其特征在于:步骤S1所述的相机标定,得到相机内参数和外参数初始值,建立图像像素的位置坐标与场景点所处空间位置之间的关系,具体过程如下:
S1.1:建立空间三维点与相机平面上二维点映射关系,空间任意一点P在图像上的成像位置用针孔模型近似表示;
空间点P在世界坐标系和相机坐标系下的齐次坐标分别为(xw,yw,zw,1)和(x,y,z,1),将P从世界坐标系转至相机坐标系的刚体变换,写为如下表达式:
其中M2表示世界坐标系和相机坐标系之间的转换关系,为4×4的矩阵;
将相机坐标系内物点P以理想小孔模型映射到图像坐标系,表达式如下:
其中(X,Y)为P点对应的图像坐标,f代表相机的焦距;
相机CMOS靶面上,像元水平方向尺寸用dx表示,垂直方向尺寸用dy表示,则理想小孔模型下的像点p,在图像坐标系下物理坐标到像素坐标之间的转化用如下表达式表示:
由公式(1)至(3)得到,物点P从世界坐标系经过忽略畸变的小孔模型,投影到像面坐标系下的数学表达式如下:
其中内参数矩阵fx、fy是单目相机u轴、v轴上的两个有效焦距;γ是用来描述相机像平面u轴与v轴的相互垂直程度的参数,令γ=0,(u0,v0)为光心坐标,为空间点P在世界坐标系下的齐次坐标;
相机外参数中R是旋转矩阵,T是平移向量;
S1.2:对相机进行标定,获得相机内外参数;通过计算单应性矩阵Hi求解出相机内参数矩阵A和相机相对于不同位置靶标的旋转矩阵Ri,以及平移向量Ti(i=1,2…..N),利用径向畸变系数与切向畸变系数k1,k2,p1,p2校正理想模型,最后使用L-M优化算法对所有参数的初始值进行优化;
平面靶标上角点P在世界坐标系下的坐标记为P=[xw,yw,zw]T,经过相机数学模型成像后的像点坐标记为p=[u,v]T,将二者的坐标写作齐次坐标为和根据公式(4)得出标志角点P与其像点p在相机成像数学模型中的转化关系:
上式中s为一任意的非零尺度因子,A为相机内参数矩阵,R是靶标坐标系与相机坐标系之间的旋转矩阵,T是靶标坐标系到相机坐标系的平移向量;
令标志点坐标的z值为零,即zw=0,令用单应性矩阵H描述靶标上点P到其像点p的变换,其中H是一个3×3的矩阵:
如果标定所用的平面靶标含有n个标定点,每个标定点列出一个关系式,所有关系式联立写成其中L是2n×9的系数矩阵;当所制备的平面靶标上标定特征点的个数n>4时,是一个超定方程组,利用SVD分解算法求解其最优解由解出H;
标定过程中相机总共采集了N个靶标空间不同位置的图像,总共有N个单应性矩阵Hi(i=1,2…..N)与每个空间位置的靶标相对应,通过计算N个单应性矩阵求出摄像机的内参数矩阵A,以及相机相对于每个靶标的旋转矩阵Ri以及平移向量Ti(i=1,2…..N);根据旋转矩阵R的正交性得到内参数矩阵A的两个约束条件:
由式(8)计算A-TA-1,并令记单应性矩阵H第i行为hi=[hi1 hi2 hi3]T,根据计算得到:
通过SVD分解法求出最优的值,由于未知数的数量为6,张正友标定法至少需要3幅标定图像,即N≥3;根据的值求解出相机的内参矩阵A:
求得单应性阵H=[h1 h2 h3]与内参矩阵A之后,由(11)求得相机相对靶标位置的外参数:
由于实验中存在误差,由公式(11)计算得到的并不具有旋转矩阵的正交性,因此,Q只能作为相机参数的初始估计值求解最佳旋转矩阵R,优化的目标函数为
对相机参数进行优化,利用最优化方法中的L-M优化方法,建立目标函数如下:
其中mij表示相机图像坐标系下实验获得的坐标,表示相机图像坐标系下经过理论计算的坐标。
3.根据权利要求2所述的一种单相机结合二维码测量玻璃板不平度的方法,其特征在于:步骤S2所述的基于玻璃板反射的相机外参数求解,得到相机与投影设备的位置关系,具体过程如下:
S2.1:建立镜面反射数学模型
镜面到屏幕参考点P的距离与镜面到屏幕参考点镜像的距离是相等的,由此可得公式:
在相机坐标系下,nj是镜面π在j位置下的表面法向量,dj是镜面在位置j下到相机坐标系原点的距离,是屏幕上第i个点到镜面位置j的距离;
到相机坐标系原点的距离表达为故得以下公式:
由此,屏幕参考点与参考点镜像之间的关系写成:
其中S是描述镜面反射的矩阵:
公式(16)中I表示3×3的单位矩阵;
S2.2:基于镜面反射的相机外参数求解;
用R表示镜像屏幕坐标系与相机坐标系之间的旋转矩阵,T表示二者之间的平移矩阵,则将投影模型与镜面反射模型相结合得到屏幕坐标系到相机坐标系的变换关系:
t=2djRnj+T (19)
其中与t分别为相机相对屏幕坐标系的旋转矩阵与平移矩阵,这种经过反射的旋转矩阵被称之为非正规旋转矩阵。
4.根据权利要求3所述的一种单相机结合二维码测量玻璃板不平度的方法,其特征在于:步骤S3所述的玻璃板反射角点坐标已知的二维码阵列,用相机采集图像,对二维码先计算***轮廓中心距,作为定位点,然后根据轮廓顶点关系确定二维码的四个顶点坐标,提取图像中所有二维码中的角点,具体过程如下:
S3.1 在玻璃板上投射二维码阵列;
S3.2 对拍摄的玻璃板图像进行处理,转为灰度图像;
S3.3 对图像进行校正、去噪、图像增强处理;
S3.4 二维码中心点提取,包含行列信息;
S3.5 计算得到所有二维码四个角点的坐标信息,提取图像中所有二维码中的角点。
5.根据权利要求4所述的一种单相机结合二维码测量玻璃板不平度的方法,其特征在于:步骤S4所述的将检测到的特征点与标记点进行拟合,根据所得表达式求取待测面每点深度,从而求得玻璃板的表面形貌,具体过程如下:
S4.1 根据之前标定求得的相机内外参数,以及基于玻璃板的相机外参,得到目标位置与姿态的六个参数;
S4.2 根据之前提到的小孔成像模型,由入射光线和反射光线的关系求得待测物上点的法向量;
S4.3 由投影点在相机坐标系下的坐标和图像点在相机坐标系下的坐标拟合求得反射对应关系,而后求得待测物上点的深度函数g(x,y);
S4.4 对于图像中的每点求取深度g,再与图像点相乘得到待测物空间三维点信息;
S4.5 得到被测物三维点坐标信息后,用积分方式得到面型,而后用重建面中的高度差来衡量其不平度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910501700.1A CN110276808B (zh) | 2019-06-11 | 2019-06-11 | 一种单相机结合二维码测量玻璃板不平度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910501700.1A CN110276808B (zh) | 2019-06-11 | 2019-06-11 | 一种单相机结合二维码测量玻璃板不平度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110276808A true CN110276808A (zh) | 2019-09-24 |
CN110276808B CN110276808B (zh) | 2022-10-18 |
Family
ID=67962129
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910501700.1A Active CN110276808B (zh) | 2019-06-11 | 2019-06-11 | 一种单相机结合二维码测量玻璃板不平度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110276808B (zh) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110823106A (zh) * | 2019-10-16 | 2020-02-21 | 江苏大学 | 一种基于激光连续波调制原理的平板玻璃质量检测方法 |
CN110853103A (zh) * | 2019-11-12 | 2020-02-28 | 佛山智能装备技术研究院 | 一种用于深度学习姿态估计的数据集制作方法 |
CN110930458A (zh) * | 2019-10-22 | 2020-03-27 | 同济大学 | 一种简易的Nao机器人摄像头外参数标定方法 |
CN111121818A (zh) * | 2019-12-11 | 2020-05-08 | 中山大学 | 一种无人车中相机与二维码的标定方法 |
CN111141212A (zh) * | 2019-12-28 | 2020-05-12 | 西安交通大学 | 一种大规模R-LATs测量场光信号可达性分析方法 |
CN111311671A (zh) * | 2020-05-12 | 2020-06-19 | 创新奇智(南京)科技有限公司 | 一种工件测量方法、装置、电子设备及存储介质 |
CN111487043A (zh) * | 2020-05-07 | 2020-08-04 | 北京的卢深视科技有限公司 | 单目散斑结构光***的散斑投射器标定参数确定方法 |
CN111504182A (zh) * | 2020-04-10 | 2020-08-07 | 四川超影科技有限公司 | 一种利用二维码辅助校准的物体体积视觉测量方法 |
CN111912387A (zh) * | 2020-06-17 | 2020-11-10 | 成都飞机工业(集团)有限责任公司 | 基于光学反射原理的飞机全机水平测量点集中测量点法 |
CN111951336A (zh) * | 2020-08-18 | 2020-11-17 | 成都经纬快图科技有限公司 | 基于gpu的三维大坐标处理方法 |
CN112308899A (zh) * | 2020-11-09 | 2021-02-02 | 北京经纬恒润科技股份有限公司 | 一种挂车角度识别方法及装置 |
CN112344876A (zh) * | 2020-11-03 | 2021-02-09 | 上海港旺新能源科技有限公司 | 检测槽式太阳光热装置聚光反射板安装误差的***及方法 |
CN112507755A (zh) * | 2020-12-22 | 2021-03-16 | 芜湖英视迈智能科技有限公司 | 最小化二维码角点重投影误差的目标物体六自由度定位方法及*** |
CN112611343A (zh) * | 2020-10-30 | 2021-04-06 | 江苏天艾美自动化科技有限公司 | 一种用于5g环形器的芯片平整度检测装置及其检测方法 |
CN112634375A (zh) * | 2020-12-21 | 2021-04-09 | 杭州东信北邮信息技术有限公司 | 一种ai智能检测中的平面标定和三维重建方法 |
CN112950720A (zh) * | 2021-02-06 | 2021-06-11 | 南京理工大学智能计算成像研究院有限公司 | 一种远心相机参数超定方程求解方法 |
CN113822942A (zh) * | 2021-09-09 | 2021-12-21 | 南京中科逆熵科技有限公司 | 一种基于二维码的单目摄像头测量物体尺寸的方法 |
CN114509048A (zh) * | 2022-01-20 | 2022-05-17 | 中科视捷(南京)科技有限公司 | 基于单目相机的架空输电线路空间三维信息获取方法及*** |
CN114663517A (zh) * | 2022-02-07 | 2022-06-24 | 哈尔滨工业大学 | 仿真目标位姿获取方法、基于MBDyn仿真的目标捕获方法及装置、航天在轨目标捕获方法 |
CN114913243A (zh) * | 2022-05-30 | 2022-08-16 | 歌尔股份有限公司 | 一种光学材料的畸变检测方法、装置及介质 |
CN116067283A (zh) * | 2023-04-03 | 2023-05-05 | 成都飞机工业(集团)有限责任公司 | 一种深腔测量方法、装置、设备及介质 |
US12010466B2 (en) | 2021-06-22 | 2024-06-11 | Industrial Technology Research Institute | Visual recognition based method and system for projecting patterned light, method and system applied to oral inspection, and machining system |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011141776A (ja) * | 2010-01-08 | 2011-07-21 | Pegasus Net Kk | 光学的情報読取装置及びホログラム型2次元コード |
CN109506589A (zh) * | 2018-12-25 | 2019-03-22 | 东南大学苏州医疗器械研究院 | 一种基于结构光场成像的三维轮廓测量方法 |
-
2019
- 2019-06-11 CN CN201910501700.1A patent/CN110276808B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011141776A (ja) * | 2010-01-08 | 2011-07-21 | Pegasus Net Kk | 光学的情報読取装置及びホログラム型2次元コード |
CN109506589A (zh) * | 2018-12-25 | 2019-03-22 | 东南大学苏州医疗器械研究院 | 一种基于结构光场成像的三维轮廓测量方法 |
Non-Patent Citations (1)
Title |
---|
张勇等: "用单数码相机实现物体表面的三维重建", 《计算机工程与设计》 * |
Cited By (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110823106A (zh) * | 2019-10-16 | 2020-02-21 | 江苏大学 | 一种基于激光连续波调制原理的平板玻璃质量检测方法 |
CN110823106B (zh) * | 2019-10-16 | 2021-09-10 | 江苏大学 | 一种基于激光连续波调制原理的平板玻璃质量检测方法 |
CN110930458A (zh) * | 2019-10-22 | 2020-03-27 | 同济大学 | 一种简易的Nao机器人摄像头外参数标定方法 |
CN110930458B (zh) * | 2019-10-22 | 2023-05-02 | 同济大学 | 一种简易的Nao机器人摄像头外参数标定方法 |
CN110853103A (zh) * | 2019-11-12 | 2020-02-28 | 佛山智能装备技术研究院 | 一种用于深度学习姿态估计的数据集制作方法 |
CN110853103B (zh) * | 2019-11-12 | 2022-06-10 | 佛山智能装备技术研究院 | 一种用于深度学习姿态估计的数据集制作方法 |
CN111121818A (zh) * | 2019-12-11 | 2020-05-08 | 中山大学 | 一种无人车中相机与二维码的标定方法 |
CN111141212A (zh) * | 2019-12-28 | 2020-05-12 | 西安交通大学 | 一种大规模R-LATs测量场光信号可达性分析方法 |
CN111504182A (zh) * | 2020-04-10 | 2020-08-07 | 四川超影科技有限公司 | 一种利用二维码辅助校准的物体体积视觉测量方法 |
CN111487043A (zh) * | 2020-05-07 | 2020-08-04 | 北京的卢深视科技有限公司 | 单目散斑结构光***的散斑投射器标定参数确定方法 |
CN111311671B (zh) * | 2020-05-12 | 2020-08-07 | 创新奇智(南京)科技有限公司 | 一种工件测量方法、装置、电子设备及存储介质 |
CN111311671A (zh) * | 2020-05-12 | 2020-06-19 | 创新奇智(南京)科技有限公司 | 一种工件测量方法、装置、电子设备及存储介质 |
CN111912387A (zh) * | 2020-06-17 | 2020-11-10 | 成都飞机工业(集团)有限责任公司 | 基于光学反射原理的飞机全机水平测量点集中测量点法 |
CN111951336A (zh) * | 2020-08-18 | 2020-11-17 | 成都经纬快图科技有限公司 | 基于gpu的三维大坐标处理方法 |
CN111951336B (zh) * | 2020-08-18 | 2024-01-19 | 成都经纬快图科技有限公司 | 基于gpu的三维大坐标处理方法 |
CN112611343A (zh) * | 2020-10-30 | 2021-04-06 | 江苏天艾美自动化科技有限公司 | 一种用于5g环形器的芯片平整度检测装置及其检测方法 |
CN112344876A (zh) * | 2020-11-03 | 2021-02-09 | 上海港旺新能源科技有限公司 | 检测槽式太阳光热装置聚光反射板安装误差的***及方法 |
CN112308899A (zh) * | 2020-11-09 | 2021-02-02 | 北京经纬恒润科技股份有限公司 | 一种挂车角度识别方法及装置 |
CN112308899B (zh) * | 2020-11-09 | 2024-05-07 | 北京经纬恒润科技股份有限公司 | 一种挂车角度识别方法及装置 |
CN112634375B (zh) * | 2020-12-21 | 2022-08-05 | 杭州东信北邮信息技术有限公司 | 一种ai智能检测中的平面标定和三维重建方法 |
CN112634375A (zh) * | 2020-12-21 | 2021-04-09 | 杭州东信北邮信息技术有限公司 | 一种ai智能检测中的平面标定和三维重建方法 |
CN112507755A (zh) * | 2020-12-22 | 2021-03-16 | 芜湖英视迈智能科技有限公司 | 最小化二维码角点重投影误差的目标物体六自由度定位方法及*** |
CN112507755B (zh) * | 2020-12-22 | 2024-05-24 | 芜湖英视迈智能科技有限公司 | 最小化二维码角点重投影误差的目标物体六自由度定位方法及*** |
CN112950720A (zh) * | 2021-02-06 | 2021-06-11 | 南京理工大学智能计算成像研究院有限公司 | 一种远心相机参数超定方程求解方法 |
US12010466B2 (en) | 2021-06-22 | 2024-06-11 | Industrial Technology Research Institute | Visual recognition based method and system for projecting patterned light, method and system applied to oral inspection, and machining system |
CN113822942A (zh) * | 2021-09-09 | 2021-12-21 | 南京中科逆熵科技有限公司 | 一种基于二维码的单目摄像头测量物体尺寸的方法 |
CN113822942B (zh) * | 2021-09-09 | 2023-11-17 | 南京中科逆熵科技有限公司 | 一种基于二维码的单目摄像头测量物体尺寸的方法 |
CN114509048B (zh) * | 2022-01-20 | 2023-11-07 | 中科视捷(南京)科技有限公司 | 基于单目相机的架空输电线路空间三维信息获取方法及*** |
CN114509048A (zh) * | 2022-01-20 | 2022-05-17 | 中科视捷(南京)科技有限公司 | 基于单目相机的架空输电线路空间三维信息获取方法及*** |
CN114663517A (zh) * | 2022-02-07 | 2022-06-24 | 哈尔滨工业大学 | 仿真目标位姿获取方法、基于MBDyn仿真的目标捕获方法及装置、航天在轨目标捕获方法 |
CN114913243A (zh) * | 2022-05-30 | 2022-08-16 | 歌尔股份有限公司 | 一种光学材料的畸变检测方法、装置及介质 |
CN116067283B (zh) * | 2023-04-03 | 2023-08-04 | 成都飞机工业(集团)有限责任公司 | 一种深腔测量方法、装置、设备及介质 |
CN116067283A (zh) * | 2023-04-03 | 2023-05-05 | 成都飞机工业(集团)有限责任公司 | 一种深腔测量方法、装置、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN110276808B (zh) | 2022-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110276808A (zh) | 一种单相机结合二维码测量玻璃板不平度的方法 | |
CN109949899B (zh) | 图像三维测量方法、电子设备、存储介质及程序产品 | |
CN103559735B (zh) | 一种三维重建方法及*** | |
CN104335005B (zh) | 3d扫描以及定位*** | |
CN102410832B (zh) | 位置姿势测量设备和位置姿势测量方法 | |
JP4245963B2 (ja) | 較正物体を用いて複数のカメラを較正するための方法およびシステム | |
EP2104365A1 (en) | Method and apparatus for rapid three-dimensional restoration | |
Willi et al. | Robust geometric self-calibration of generic multi-projector camera systems | |
Aliaga et al. | A self-calibrating method for photogeometric acquisition of 3D objects | |
CN106500596A (zh) | 结构光全景测量***的测量方法 | |
CN109272555B (zh) | 一种rgb-d相机的外部参数获得及标定方法 | |
CN110827360B (zh) | 一种光度立体式测量***及其标定光源方向的方法 | |
CN110940295B (zh) | 基于激光散斑极限约束投影的高反射物体测量方法及*** | |
CN114283203B (zh) | 一种多相机***的标定方法及*** | |
CN107941153B (zh) | 一种激光测距优化标定的视觉*** | |
CN112489109B (zh) | 一种三维成像***方法、装置及三维成像*** | |
CN110360930A (zh) | 一种激光位移法线传感器及其测量方法 | |
WO2018120168A1 (zh) | 一种视觉检测方法及*** | |
CN113091646A (zh) | 基于条纹标定的三维形貌测量方法 | |
Xie et al. | Mesh-based computation for solving photometric stereo with near point lighting | |
CN116295113A (zh) | 一种融合条纹投影的偏振三维成像方法 | |
CN114459384A (zh) | 一种基于多角度正弦条纹光场融合的相移轮廓术 | |
Yamauchi et al. | Calibration of a structured light system by observing planar object from unknown viewpoints | |
Draréni et al. | Geometric video projector auto-calibration | |
CN113658267B (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 |