CN103438810A - 一种宽幅柔性板厚度实时检测装置及方法 - Google Patents

一种宽幅柔性板厚度实时检测装置及方法 Download PDF

Info

Publication number
CN103438810A
CN103438810A CN2013103487289A CN201310348728A CN103438810A CN 103438810 A CN103438810 A CN 103438810A CN 2013103487289 A CN2013103487289 A CN 2013103487289A CN 201310348728 A CN201310348728 A CN 201310348728A CN 103438810 A CN103438810 A CN 103438810A
Authority
CN
China
Prior art keywords
centerdot
coordinate
video camera
range cells
surface range
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
Application number
CN2013103487289A
Other languages
English (en)
Other versions
CN103438810B (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.)
JIANGSU XIECHENG SCIENCE AND TECHNOLOGY DEVELOPMENT Co Ltd
Tsinghua University
Original Assignee
JIANGSU XIECHENG SCIENCE AND TECHNOLOGY DEVELOPMENT Co Ltd
Tsinghua University
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 JIANGSU XIECHENG SCIENCE AND TECHNOLOGY DEVELOPMENT Co Ltd, Tsinghua University filed Critical JIANGSU XIECHENG SCIENCE AND TECHNOLOGY DEVELOPMENT Co Ltd
Priority to CN201310348728.9A priority Critical patent/CN103438810B/zh
Publication of CN103438810A publication Critical patent/CN103438810A/zh
Application granted granted Critical
Publication of CN103438810B publication Critical patent/CN103438810B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明涉及一种宽幅柔性板厚度实时检测装置及方法,该装置包括上、下表面测距单元,上、下导轨,上、下滑块;上、下导轨平行放置,上、下滑块分别沿上、下导轨滑动,上、下表面测距单元分别与上、下滑块固连,上滑块和所述下滑块同步运动;所述上、下表面测距单元均包括激光位移传感器、激光结构光源、摄像机和半透半反镜;半透半反镜法线方向与激光位移传感器的光线出射方向成135°,与激光结构光源出射方向成45°,且三者在同一平面内;摄像机的镜头朝向半透半反镜。本发明的方法包括以下步骤:1)对上表面测距单元进行标定;2)对下表面测距单元进行标定;3)对厚度检测装置进行标定;4)对厚度进行检测。本发明可以广泛应用在宽幅柔性板厚度实时检测中。

Description

一种宽幅柔性板厚度实时检测装置及方法
技术领域
本发明涉及一种厚度实时检测装置及方法,特别是关于一种宽幅柔性板厚度实时检测装置及方法。
背景技术
一般A2级防火材料板的生产需要首先将无机阻烧材料与有机胶混合,经振动加工辊挤后加热固化,形成厚度均匀的3mm柔性板材,成型后,再将芯板与上下两面铝板分别粘接。为提高生产效率并稳定产品质量,A2级防火材料板的生产需要进行全线自动化。为实现对防火板的厚度、平整度、缺陷、粘接强度等在线质量控制,其中经热固得到的芯板厚度及表面状态是反映生产过程状态的重要指标,对其进行实时检测是实现整个生产线自动化及质量控制的重要前提,为质量检验、工艺改进和生产线管理提供基础。
工业生产中使用的非接触式厚度检测系统多采用O形或C形结构,将两个传感器分别固连于机构的上下悬臂,条形或带形材料在两者之间连续地通过,由上下两个传感器分别测得的到对应表面的距离与两个传感器彼此之间的距离之差即为被测材料的厚度。非接触式传感器通常采用激光位移传感器、电涡流传感器、电磁传感器、超声波传感器、辐射传感器等。其中激光位移传感器是利用激光三角测距原理,由半导体激光器发射出的激光束在经被测物体表面反射后,再经透镜投射到线阵CCD上,并将线阵CCD上的光点位置经过三角函数计算得到位移传感器到被测物体表面的距离。该方法精度可达到微米级,广泛应用于位置、位移、厚度等工业测量领域。
美国Moduloc公司的MLC/MLM/MLS激光测量系统采用一对激光位移传感器,在待测板材上下分别进行到板材表面距离的测量,再通过d=H-1h-2h计算得到厚度的测量值d,其中d表示测量值,H表示两个激光位移传感器之间的距离,1h表示上方位移传感器得到的测量值,2h表示下方位移传感器得到的测量值。该装置固定使用时只能测量固定位置上的厚度值,而无法获得待测板材方向上厚度值的变化。德国MicroEpsilon公司的MTS8202.LLT CFK Cclamp宽幅板厚度测量系统采用了C形框架结构,将一对激光位移传感器分别固定于框架的上下悬臂,C形框架在导轨上作沿待测板材宽度方向的往复运动,以实现板材宽度方向全程的厚度测量。采用激光点或激光条三角测距原理分别获得两个位移传感器到待测板材上下表面的距离,再利用位移传感器之间的距离与标准厚度板的标定结果得到厚度测量结果。
现有方法存在如下的问题:对于在板材送进方向上存在倾斜的待测板材,采用激光点三角测距方法无法对其倾斜进行修正,而激光条三角测距法也只能修正板材送进方向的倾斜角带来的测量误差,而对于A2级防火材料芯板来说,板材宽度较大且材料具有一定柔性,在送进方向与宽度方向上都可能存在一定的弯曲或倾斜,激光测距传感器的激光指向方向与柔性板材的局部平面可能是非垂直的,若采用现有方法进行测量和计算,将有较大的厚度测量误差。
采用激光点三角测距原理进行厚度测量时,计算厚度的传统方法如下式所示:
d=D+1h0+2h0-1h-2h
但在待测柔性板材相对于测量激光入射方向非垂直时存在测量误差,测量时的系统误差如下式所示:
Figure BDA00003651612300021
其中,α接近π/2,l为小量,可知式中的第一项为误差的主要来源,θ角超过1.1°时厚度测量误差就超过2%,即当测量厚度为5mm时,系统误差超过0.1mm,此时无法满足生产要求;而采用激光条三角测距原理进行厚度测量的方法仅能对沿柔性板材送进方向的倾斜进行修正。
发明内容
针对上述问题,本发明的目的是提供一种宽幅柔性板厚度实时检测装置及方法,以实现宽幅柔性板厚度的实时准确检测与厚度信号特征值的提取,为产品质量控制提供基础。
为实现上述目的,本发明采取以下技术方案:一种宽幅柔性板厚度实时检测装置,其特征在于:它包括上表面测距单元、下表面测距单元、上导轨、下导轨、上滑块和下滑块;
所述上导轨与所述下导轨平行放置,所述上滑块沿所述上导轨滑动,所述下滑块沿所述下导轨滑动,且所述上表面测距单元与所述上滑块固连,所述下表面测距单元与所述下滑块固连,且所述上滑块和所述下滑块同步运动;
所述上表面测距单元和下表面测距单元均包括固定在盒体内的激光位移传感器、激光结构光源、摄像机和半透半反镜;所述半透半反镜法线方向与所述激光位移传感器的光线出射方向成135°,与所述激光结构光源出射方向成45°,且三者在同一平面内;所述摄像机的镜头朝向所述半透半反镜。
所述激光位移传感器接收端前方放置有一激光位移传感器滤光片,所述摄像机镜头前方放置有一摄像机滤光片。
所述激光位移传感器的激光波长为λ1,所述激光结构光源采用单圆、同心多圆、、网格和网点中的一种,激光波长为λ2,并满足下式:
12|≥20nm。
所述激光位移传感器滤光片为窄带通滤光片,中心波长为λ1,半带宽不大于10nm;所述摄像机滤光片为窄带通滤光片,中心波长为λ2,半带宽不大于10nm。
所述摄像机包括镜头与图像传感器,其中所述图像传感器采用CCD、CMOS和光电转换器件中的一种。
一种采用上述装置实现的宽幅柔性板厚度实时检测方法,其包括以下步骤:
1)对上表面测距单元进行标定;
2)对下表面测距单元进行标定;
3)对厚度检测装置进行标定;
4)对厚度进行检测。
所述步骤1)中对上表面测距单元进行标定包括以下步骤:
①记摄像机光心坐标系为{C1},摄像机对应的图像坐标系为{Im1},CCD摄像机标定过程中第n个标定板空间位置对应的标定板坐标系为{S1n};
②对摄像机进行标定,得到摄像机内部参数{1fx,1fy,1u0,1v0},其中1fx,1fy为镜头光心坐标系X,Y方向上的等效焦距,(1u0,1v0)为镜头光轴与图像平面交点在图像坐标系中的坐标;
③对激光结构光源与摄像机相对位置进行标定:将标定板放置于摄像机前方,将激光结构光投射于标定板上,变换标定板位置,拍摄一系列图像
Figure BDA00003651612300031
其中n=1,2,...,N,N≥4,右上标calib表示该图像用于标定,提取图像
Figure BDA00003651612300032
激光光条上不共线的五个特征点,分别记为
Figure BDA00003651612300033
其位置关系满足:
Figure BDA00003651612300034
共线、
Figure BDA00003651612300035
共线;特征点在图像坐标系中的坐标记为
Figure BDA00003651612300036
将摄像机光心坐标系记为1OC1XC1YC1ZC,第n个标定板坐标系记为1On1Xn1Yn1Zn,将两坐标系之间的转换矩阵记为[1Rn  1Tn],标定板的各个角点在1OC1XC1YC1ZC中的坐标记为
Figure BDA00003651612300037
由标定板图案角点位置关系知其在1On1Xn1Yn1Zn中的坐标为
Figure BDA00003651612300038
图像
Figure BDA00003651612300039
上角点的像素坐标为
Figure BDA000036516123000310
三者关系满足以下两式:
R n 1 T n 1 x n k 1 y n k 1 0 1 = x C k 1 y C k 1 z C k 1
f x 1 0 u 0 1 0 f y 1 v 0 1 0 0 1 x C k 1 z C k 1 y C k 1 z C k 1 1 = u n k 1 v n k 1 1
联立求解得到转换矩阵[1Rn  1Tn]
记点
Figure BDA00003651612300042
1OC1XC1YC1ZC中的坐标为
Figure BDA00003651612300043
由以下两式联立求解得到:
R n 1 T n 1 x n i 1 y n i 1 0 1 = x Cn i 1 y Cn i 1 z Cn i 1
f x 1 0 u 0 1 0 f y 1 v 0 1 0 0 1 x Cn i 1 z Cn i 1 y Cn i 1 z Cn i 1 1 = u n i 1 v n i 1 1
Figure BDA00003651612300046
在激光结构光的两个光平面上,光平面在1OC1XC1YC1ZC中的方程满足以下两式:
A 1 1 · x Cn i 1 + B 1 1 · y Cn i 1 + C 1 1 · z Cn i 1 = 1 , i ∈ { 0,1,3 } , n = { 1,2 , . . . , N }
A 2 1 · x Cn i 1 + B 2 1 · y Cn i 1 + C 2 1 · z Cn i 1 = 1 , i ∈ { 0,1,3 } , n = { 1,2 , . . . , N }
联立求解得到(1A1,1B1,1C1,1A2,1B2,1C2)。
所述步骤2)中对下表面测距单元进行标定包括以下步骤:
①记摄像机光心坐标系为{C2},摄像机对应的图像坐标系为{Im2},摄像机标定过程中第n个标定板空间位置对应的标定板坐标系为{S2n};
②对摄像机进行标定,得到摄像机内部参数{2fx,2fy,2u0,2v0},其中2fx,2fy为镜头光心坐标系X,Y方向上的等效焦距,(2u0,2v0)为镜头光轴与图像平面交点在图像坐标系中的坐标;
③对激光结构光源与摄像机的相对位置进行标定:将标定板放置于摄像机前方,将激光结构光投射于标定板上,变换标定板位置,拍摄一系列图像
Figure BDA00003651612300049
其中n=1,2,...,N,N≥4,右上标calib表示该图像用于标定,提取图像
Figure BDA000036516123000410
激光光条上不共线的五个特征点,分别记为
Figure BDA00003651612300051
其位置关系满足:
Figure BDA00003651612300052
共线,
Figure BDA00003651612300053
共线;特征点在图像坐标系中的坐标记为
Figure BDA00003651612300054
将摄像机光心坐标系记为2OC2XC2YC2ZC,第n个标定板坐标系记为2On2Xn2Yn2Zn,将两坐标系之间的转换矩阵记为[2Rn  2Tn],标定板的各个角点在2OC2XC2YC2ZC中的坐标记为
Figure BDA00003651612300055
由标定板图案角点位置关系知其在2On2Xn2Yn2Zn中的坐标为
Figure BDA00003651612300056
图像
Figure BDA00003651612300057
上角点的像素坐标为三者关系满足以下两式:
R n 2 T n 2 x n k 2 y n k 2 0 1 = x C k 2 y C k 2 z C k 2
f x 2 0 u 0 2 0 f y 2 v 0 2 0 0 1 x C k 2 z C k 2 y C k 2 z C k 2 1 = u n k 2 v n k 2 1
联立求解得到转换矩阵[2Rn  2Tn];记点
Figure BDA000036516123000511
2OC2XC2YC2ZC中的坐标为由以下两式联立求解得到:
R n 2 T n 2 x n i 2 y n i 2 0 1 = x Cn i 2 y Cn i 2 z Cn i 2
f x 2 0 u 0 2 0 f y 2 v 0 2 0 0 1 x Cn i 2 z Cn i 2 y Cn i 2 z Cn i 2 1 = u n i 2 v n i 2 1
Figure BDA000036516123000515
在激光结构光的两个光平面上,光平面在2OC2XC2YC2ZC中的方程满足以下两式:
A 1 2 · x Cn i 2 + B 1 2 · y Cn i 2 + C 1 2 · z Cn i 2 = 1 , i ∈ { 0,1,3 } , n = { 1,2 , . . . , N }
A 2 2 · x Cn i 2 + B 2 2 · y Cn i 2 + C 2 2 · z Cn i 2 = 1 , i ∈ { 0 , 2 , 4 } , n = { 1,2 , . . . , N }
联立求解得到(2A1,2B1,2C1,2A2,2B2,2C2)。
所述步骤3)中对厚度检测装置进行标定包括以下步骤:
将厚度为D的标准板置于上下正对安装的上表面测距单元与下表面测距单元之间,变化标准板的角度,从上表面测距单元与下表面测距单元中分别得到一系列测量值{1hm}与{2hm};上表面测距单元中的激光结构光源投射在标准板上表面的五个特征点
Figure BDA00003651612300061
在上表面测距单元中的摄像机拍摄的图像中坐标分别为下表面测距单元中的激光结构光源投射在标准板下表面的五个特征点
Figure BDA00003651612300063
在下表面测距单元中的摄像机拍摄的图像中坐标分别为
Figure BDA00003651612300064
i=0,1,2,3,4,m=1,2,...,M,M为标定次数;
Figure BDA00003651612300065
在上表面测距单元中摄像机光心坐标系1OC1XC1YC1ZC中的坐标
Figure BDA00003651612300066
满足以下三式:
f x 1 0 u 0 1 0 f y 1 v 0 1 0 0 1 x Cm i 1 z Cm i 1 y Cm i 1 z Cm i 1 1 = u m i 1 v m i 1 1
A 1 1 · x Cm i 1 + B 1 1 · y Cm i 1 + C 1 1 · z Cm i 1 = 1 , i ∈ { 0 , 1 , 3 }
A 2 1 · x Cm i 1 + B 2 1 · y Cm i 1 + C 2 1 · z Cm i 1 = 1 , i ∈ { 0 , 2 , 4 }
联立求解得到坐标
Figure BDA000036516123000610
满足下式:
A m top 1 · x Cm i 1 + B m top 1 · y Cm i 1 + C m top 1 · z Cm i 1 = 1 , i ∈ { 0,1,2,3,4 } , m ∈ { 1,2 , . . . , M }
解得
Figure BDA000036516123000612
在下表面测距单元中摄像机光心坐标系2OC2XC2YC2ZC中的坐标
Figure BDA000036516123000614
满足以下三式:
f x 2 0 u 0 2 0 f y 2 v 0 2 0 0 1 x Cm i 2 z Cm i 2 y Cm i 2 z Cm i 2 1 = u m i 2 v m i 2 1
A 1 2 · x Cm i 2 + B 1 2 · y Cm i 2 + C 1 2 · z Cm i 2 = 1 , i ∈ { 0 , 1 , 3 }
A 2 2 · x Cm i 2 + B 2 2 · y Cm i 2 + C 2 2 · z Cm i 2 = 1 , i ∈ { 0 , 2 , 4 }
联立求解得到坐标
Figure BDA000036516123000618
满足下式:
A m bottom 2 · x Cm i 2 + B m bottom 2 · y Cm i 2 + C m bottom 2 · z Cm i 2 = 1 , i ∈ { 0,1,2,3,4 } , m ∈ { 1,2 , . . . , M }
解得
Figure BDA00003651612300071
记标准板上表面法向的单位向量为
Figure BDA00003651612300072
其在1OC1XC1YC1ZC坐标系的表达式满足下式:
a m 1 b m 1 c m 1 = A m top 1 B m top 1 C m top 1 ( A m top 1 ) 2 + ( B m top 1 ) 2 + ( C m top 1 ) 2
记标准板下表面法向的单位向量为
Figure BDA00003651612300074
其在1OC1XC1YC1ZC坐标系的表达式满足下式:
a m 2 b m 2 c m 2 = A m top 2 B m top 2 C m top 2 ( A m top 2 ) 2 + ( B m top 2 ) 2 + ( C m top 2 ) 2
2OC2XC2YC2ZC坐标系到1OC1XC1YC1ZC坐标系的转换矩阵为 R 2 1 T 2 1 , 标准板的上下表面平行,上述两单位向量的表达式满足下式;
R 2 1 a m 2 b m 2 c m 2 = a m 1 b m 1 c m 1 , m = 1,2 , . . . , M
求解得到
Figure BDA00003651612300078
以上表面测距单元的激光位移传感器的激光出射方向为Z轴负向,以下表面测距单元的激光位移传感器的位移计算起始点到Z轴的垂向为X轴负向,建立坐标系OXYZ,上表面测距单元的激光位移传感器与下表面测距单元的激光位移传感器的激光出射方向所在直线L1与L2的方程满足以下两式:
x=0,y=0
x - l cos α cos β = y cos α sin β = z sin α
记上表面测距单元的激光位移传感器的位移计算起始点坐标为(0,0,H),记L1与S1的交点为P1(0,0,H-1hm),记L2与S2的交点为P2(l+2hmcosαcosβ,2hmcosαsinβ,2hmsinα);P1与P2分别在标准板上下表面所在平面S1与S2上,两平面在OXYZ坐标系中的方程满足以下两式:
Figure BDA000036516123000710
Figure BDA000036516123000711
两平面间的距离为D,满足下式:
D=k1m-k2m
联立上述三式,得到下式:
θm为OXYZ坐标系的Z轴与标准板上表面所在平面的夹角,近似得到下式:
sin θ m = D H - h m 1 - h m 2
记坐标系OXYZ的Z轴方向在1OC1XC1YC1ZC坐标系中表示为向量 z → c = a c b c c c , 满足下式:
&theta; m = &pi; 2 - < n &RightArrow; m 1 , z &RightArrow; c > = &pi; 2 - arccos ( n &RightArrow; m 1 &CenterDot; z &RightArrow; c ) = &pi; 2 - arccos [ a m 1 b m 1 c m 1 &CenterDot; a c b c c c T ]
联立上述两式得到下式:
a m 1 b m 1 c m 1 a c b c c c = D H - h m 1 - h m 2 , m = 1,2 , . . . , M
求解得到(acbccc)与H。
所述步骤4)中对厚度进行检测包括以下步骤:
将待测板放置于上表面测距单元与下表面测距单元之间,得到上表面测距单元中的激光位移传感器的测量值为1h,上表面测距单元中的激光位移传感的测量值为2h,下表面测距单元中的摄像机拍摄的图像为1I,下表面测距单元中的摄像机拍摄的图像为2I;
图像1I中投射在待测板上表面的激光的五个特征点分别为点1P0与四个端点1P1,1P2,1P3,1P4,通过图像处理,得到1P0,1P1,1P2,1P3,1P4的图像坐标分别为(1ui,1vi),i=0,1,2,3,4;
图像2I中投射在待测板下表面的激光的五个特征点分别为点2P0与四个端点2P1,2P2,2P3,2P4,通过图像处理,得到2P0,2P1,2P2,2P3,2P4的图像坐标分别为(2ui,2vi),i=0,1,2,3,4;
1Pi在上表面测距单元中的摄像机光心坐标系1OC1XC1YC1ZC中的坐标
Figure BDA00003651612300085
满足以下三式:
f x 1 0 u 0 1 0 f y 1 v 0 1 0 0 1 x C i 1 z C i 1 y C i 1 z C i 1 1 = u i 1 v i 1 1
A 1 1 &CenterDot; x C i 1 + B 1 1 &CenterDot; y C i 1 + C 1 1 &CenterDot; z C i 1 = 1 , i &Element; { 0 , 1 , 3 }
A 2 1 &CenterDot; x C i 1 + B 2 1 &CenterDot; y C i 1 + C 2 1 &CenterDot; z C i 1 = 1 , i &Element; { 0 , 2 , 4 }
联立解得坐标
Figure BDA00003651612300089
坐标满足下式;
A top 1 &CenterDot; x C i 1 + B top 1 &CenterDot; y C i 1 + C top 1 &CenterDot; z C i 1 = 1 , i &Element; { 0,1,2,3,4 }
联立解得
2Pi在下表面测距单元中的摄像机光心坐标系2OC2XC2YC2ZC中的坐标
Figure BDA00003651612300093
满足以下三式:
f x 2 0 u 0 2 0 f y 2 v 0 2 0 0 1 x C i 2 z C i 2 y C i 2 z C i 2 1 = u i 2 v i 2 1
A 1 2 &CenterDot; x C i 2 + B 1 2 &CenterDot; y C i 2 + C 1 2 &CenterDot; z C i 2 = 1 , i &Element; { 0 , 1 , 3 }
A 2 2 &CenterDot; x C i 2 + B 2 2 &CenterDot; y C i 2 + C 2 2 &CenterDot; z C i 2 = 1 , i &Element; { 0 , 2 , 4 }
联立解得坐标
Figure BDA00003651612300097
坐标满足下式:
A bottom 2 &CenterDot; x C i 2 + B bottom 2 &CenterDot; y C i 2 + C bottom 2 &CenterDot; z C i 2 = 1 , i &Element; { 0,1,2,3,4 }
联立解得(2Abottom,2Bbottom,2Cbottom);
坐标系1OC1XC1YC1ZC中待测板上下表面局部所在平面的法向量满足下式:
a 1 b 1 c 1 = a top 1 + a bottom 1 b top 1 + b bottom c top 1 + c bottom 1 ( a top 1 + a bottom 1 ) 2 + ( b top 1 + b bottom 1 ) 2 + ( c top 1 + c bottom 1 ) 2
待测板上下表面局部所在平面与坐标系OXYZ的Z轴之间夹角θ满足下式:
&theta; = &pi; 2 - arccos [ a 1 b 1 c 1 &CenterDot; a c b c c c T ]
待测板厚度由下式得到:
d=(H-1h-2h)sinθ。
本发明由于采取以上技术方案,其具有以下优点:1、本发明由于采用具有二维结构的激光结构光,比如:单圆、同心多圆、十字交叉、网格、网点等投射到柔性板材表面,再通过图像处理的方法,根据摄像机成像原理,计算得到柔性板材的测量局部平面相对于装置的倾斜角变化,据此可对厚度测量结果进行修正,提高测量精度。2、因为用于倾角检测的激光与用于位移传感的激光采用不同的波长,通过滤光可以使摄像机拍摄到的画面中仅出现用于倾角检测的二维激光模式、同时激光位移传感器的线阵仅接收到激光位移传感器所发射激光的反射光,因此两者互相不产生干扰。
附图说明
图1是本发明的应用状态示意图
图2是本发明的结构示意图
图3是本发明CCD摄像机拍摄的五个特征点示意图
图4是本发明两出射激光在待测板上的位置示意图
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明的装置包括上表面测距单元1、下表面测距单元2、上导轨3、下导轨4、上滑块5和下滑块6。
如图1、图2所示,本发明的上表面测距单元1与上滑块5固连,上表面测距单元1包括一固连激光位移传感器11、一激光结构光源12、一CCD摄像机13、一激光位移传感器滤光片14、一CCD摄像机滤光片15、一半透半反镜16和一上盒体17。激光位移传感器11顶部连接在上盒体17内顶部,光线出射方向垂直向下。激光结构光源12的底部固定在上盒体17的一侧,出射方向为平行上导轨3方向,并指向半透半反镜16。CCD摄像机13固定在与激光结构光源12同侧的上盒体17内,且镜头朝向半透半反镜16。激光位移传感器滤光片14放置于激光位移传感器11接收端前方,CCD摄像机滤光片15放置于CCD摄像机13镜头前方。半透半反镜16法线方向与激光位移传感器11的光线出射方向成135°、与激光结构光源12出射方向成45°,且三者在同一平面内。
在上述实施例中,本发明下表面测距单元2与下滑块6固连,下表面测距单元2与上表面测距单元1类似,包括一固连激光位移传感器21、一激光结构光源22、一CCD摄像机23、一位移传感器滤光片24、一CCD摄像机滤光片25、一半透半反镜26和一下盒体27。激光位移传感器21底部固定在下盒体27内底部,光线出射方向垂直向上。激光结构光源22的底部固定在上盒体27的一侧,出射方向为平行下导轨4方向,并指向半透半反镜26。CCD摄像机23固定在与激光结构光源22同侧的下盒体27内,且镜头朝向半透半反镜26。激光位移传感器滤光片24放置于激光位移传感器21接收端前方,CCD摄像机滤光片25放置于CCD摄像机23镜头前方。半透半反镜26法线方向与激光位移传感器21的光线出射方向成135°、与激光结构光源22出射方向成45°,且三者在同一平面内。
在一个优选的实施例中,激光位移传感器11、21的激光波长为λ1;激光结构光源12、22可以采用单圆、同心多圆、十字交叉、网格和网点中的一种,激光波长为λ2,并满足下式:
12|≥20nm。
在一个优选的实施例中,激光位移传感器滤光片14、24为窄带通滤光片,中心波长为λ1,半带宽不大于10nm;CCD摄像机滤光片15、25为窄带通滤光片,中心波长为λ2,半带宽不大于10nm。
在一个优选的实施例中,激光位移传感器11、21的激光波长λ1为650nm,激光结构光源12、22的激光波长λ2为405nm。
在一个优选的实施例中,CCD摄像机13、23包括镜头与图像传感器,其中图像传感器可以采用CCD、CMOS和光电转换器件中的一种。
本发明上导轨3与下导轨4平行放置,上滑块5沿上导轨3滑动,下滑块6沿下导轨4滑动,且上滑块5与下滑块6同步运动。待测板7从上导轨3与下导轨4之间通过,待测板7表面与导轨3、4平行,送进方向为上滑块5运动路径所在直线与下滑块6运动路径所在直线组成的平面的法线方向。
工作时,激光位移传感器11(21)的出射激光与激光结构光源12(22)的出射激光在经半透半反镜16(26)后近似垂直于待测板7表面的方向,指向待测板7。激光位移传感器11(21)的出射激光经待测板7反射后,经半透半反镜16(26)后返回激光位移传感器11(21)的接收器,由激光结构光源12(22)投射的入射光经待测板7反射后,反射光经半透半反镜16(26)进入CCD摄像机13(23)。
根据上述实施例中提供的实时检测装置,本发明还提出了一种宽幅柔性板厚度实时检测方法,其包括以下步骤:
1)对上表面测距单元1进行标定
①记CCD摄像机13光心坐标系为{C1},CCD摄像机13对应的图像坐标系为{Im1},CCD摄像机13标定过程中第n个标定板空间位置对应的标定板坐标系为{S1n};
②对CCD摄像机13进行标定,得到CCD摄像机内部参数{1fx,1fy,1u0,1v0},其中1fx,1fy为CCD摄像机13的镜头在{C1}的X,Y方向上的等效焦距,1u0,1v0为CCD摄像机13的镜头光轴与图像平面交点在{Im1}中的坐标;采用棋盘格形式的标定板对CCD摄像机13进行标定,CCD摄像机标定方法为已有技术,此处不再详述其操作;
③对激光结构光源12(以十字交叉激光结构为例)与CCD摄像机13的相对位置进行标定:将标定板放置于CCD摄像机13前方,将十字交叉激光结构光投射于标定板上,变换标定板的空间位置,每个位置拍摄一幅图像,共N幅,分别记为
Figure BDA00003651612300115
其中右下标n表示标定板空间位置的序号,n=1,2,...,N,N≥4,右上标calib表示该图像用于标定;如图2所示,将图像中激光光条的十字交叉点与四个端点分别记为
Figure BDA00003651612300112
并将其在{Im1}中的坐标记为
Figure BDA00003651612300113
将坐标系{C1}与{S1n}之间的转换矩阵记为[1Rn  1Tn],标定板棋盘格的各个角点在{C1}中的坐标记为
Figure BDA00003651612300114
其中k表示标定板棋盘格各角点的序号,由标定板棋盘格实际尺寸可知其在{S1n}中的坐标为
Figure BDA00003651612300121
上述角点在图像
Figure BDA00003651612300122
中像素坐标为
Figure BDA00003651612300123
其中k表示标定板棋盘格各角点的序号,联立式(1)、式(2),求解得到转换矩阵[1Rn  1Tn];记点
Figure BDA00003651612300124
在{C1}中的坐标依次为
Figure BDA00003651612300125
i=0,1,2,3,4,由式(3)与式(4)联立求解,得到
Figure BDA00003651612300126
十字交叉激光结构光的两个光平面在{C1}中的方程为式(5)与式(6),联立求解得到方程式(5)与式(6)的参数(1A1,1B1,1C1,1A2,1B2,1C2);
R n 1 T n 1 x n k 1 y n k 1 0 1 = x C k 1 y C k 1 z C k 1 - - - ( 1 )
f x 1 0 u 0 1 0 f y 1 v 0 1 0 0 1 x C k 1 z C k 1 y C k 1 z C k 1 1 = u n k 1 v n k 1 1 - - - ( 2 )
R n 1 T n 1 x n i 1 y n i 1 0 1 = x Cn i 1 y Cn i 1 z Cn i 1 - - - ( 3 )
f x 1 0 u 0 1 0 f y 1 v 0 1 0 0 1 x Cn i 1 z Cn i 1 y Cn i 1 z Cn i 1 1 = u n i 1 v n i 1 1 - - - ( 4 )
A 1 1 &CenterDot; x Cn i 1 + B 1 1 &CenterDot; y Cn i 1 + C 1 1 &CenterDot; z Cn i 1 = 1 , i &Element; { 0,1,3 } , n = { 1,2 , . . . , N } - - - ( 5 )
A 2 1 &CenterDot; x Cn i 1 + B 2 1 &CenterDot; y Cn i 1 + C 2 1 &CenterDot; z Cn i 1 = 1 , i &Element; { 0,1,3 } , n = { 1,2 , . . . , N } - - - ( 6 )
2)对下表面测距单元2进行标定
①记CCD摄像机23光心坐标系为{C2},CCD摄像机23对应的图像坐标系为{Im2},CCD摄像机23标定过程中第n个标定板空间位置对应的标定板坐标系为{S2n};
②对CCD摄像机23进行标定,得到CCD摄像机内部参数{2fx,2fy,2u0,2v0},其中2fx,2fy为CCD摄像机23的镜头在{C2}的X,Y方向上的等效焦距,(2u0,2v0)为CCD摄像机23的镜头光轴与图像平面交点在{Im2}中的坐标;采用棋盘格形式的标定板对CCD摄像机进行标定,CCD摄像机标定方法为已有技术,此处不再详述其操作;
③对激光结构光源22(以十字交叉激光结构为例)与CCD摄像机23的相对位置进行标定:将标定板放置于CCD摄像机23前方,将十字交叉激光结构光投射于标定板上,变换标定板的空间位置,每个位置拍摄一幅图像,共N幅,分别标记为
Figure BDA00003651612300131
其中右下标n表示标定板空间位置的序号,n=1,2,...,N,N≥4,右上标calib表示该图像用于标定;将图像
Figure BDA00003651612300132
上激光光条的十字交叉点与四个端点分别记为
Figure BDA00003651612300133
并将其在{Im2}中的坐标记为
Figure BDA00003651612300134
将坐标系{C2}与{Im2}之间的转换矩阵记为[2Rn  2Tn],标定板棋盘格的各个角点在{C2}中的坐标记为其中k表示标定板棋盘格各角点的序号,由标定板棋盘格实际尺寸可知其在{S2n}中的坐标为
Figure BDA00003651612300136
上述交点在图像
Figure BDA00003651612300137
中的坐标为
Figure BDA00003651612300138
其中k表示标定板棋盘格各角点的序号,联立式(7)、式(8),求解得到转换矩阵[2Rn  2Tn];记点
Figure BDA00003651612300139
在{C2}中的坐标依次为
Figure BDA000036516123001310
i=0,1,2,3,4,由式(9)与式(10)联立求解,得到
Figure BDA000036516123001311
十字交叉激光结构光的两个光平面在{C2}中的方程为式(11)与式(12),联立求解得到方程式(11)与式(12)的参数(2A1,2B1,2C1,2A2,2B2,2C2);
R n 2 T n 2 x n k 2 y n k 2 0 1 = x C k 2 y C k 2 z C k 2 - - - ( 7 )
f x 2 0 u 0 2 0 f y 2 v 0 2 0 0 1 x C k 2 z C k 2 y C k 2 z C k 2 1 = u n k 2 v n k 2 1 - - - ( 8 )
R n 2 T n 2 x n i 2 y n i 2 0 1 = x Cn i 2 y Cn i 2 z Cn i 2 - - - ( 9 )
f x 2 0 u 0 2 0 f y 2 v 0 2 0 0 1 x Cn i 2 z Cn i 2 y Cn i 2 z Cn i 2 1 = u n i 2 v n i 2 1 - - - ( 10 )
A 1 2 &CenterDot; x Cn i 2 + B 1 2 &CenterDot; y Cn i 2 + C 1 2 &CenterDot; z Cn i 2 = 1 , i &Element; { 0,1,3 } , n = { 1,2 , . . . , N } - - - ( 11 )
A 2 2 &CenterDot; x Cn i 2 + B 2 2 &CenterDot; y Cn i 2 + C 2 2 &CenterDot; z Cn i 2 = 1 , i &Element; { 0 , 2 , 4 } , n = { 1,2 , . . . , N } - - - ( 12 )
3)对厚度检测装置进行标定
将厚度为D的标准板置于上下正对安装的上表面测距单元1与下表面测距单元2之间,变化标准板角度,从上表面测距单元1得到一系列上表面测距单元中激光位移传感器的测量值1hm,从下表面测距单元2得到一系列下表面测距单元中激光位移传感器的测量值2hm,其中m=1,2,...,M;激光结构光源12投射在标准板上表面的十字交叉点与四个端点在{Im1}中坐标分别为
Figure BDA00003651612300143
其中i=0,1,2,3,4,m=1,2,...,M;激光结构光源22投射在标准板下表面的十字交叉点与四个端点在{Im2}中坐标分别为
Figure BDA00003651612300145
其中i=0,1,2,3,4,m=1,2,...,M,M为标定次数;
Figure BDA00003651612300146
在{C1}中的坐标为
Figure BDA00003651612300147
联立式(13)、式(14)与式(15),解得坐标
Figure BDA00003651612300148
标准板的上表面所在平面在{C1}中的方程如式(16)所示,求解得到方程式(16)的参数
Figure BDA00003651612300149
f x 1 0 u 0 1 0 f y 1 v 0 1 0 0 1 x Cm i 1 z Cm i 1 y Cm i 1 z Cm i 1 1 = u m i 1 v m i 1 1 - - - ( 13 )
A 1 1 &CenterDot; x Cm i 1 + B 1 1 &CenterDot; y Cm i 1 + C 1 1 &CenterDot; z Cm i 1 = 1 , i &Element; { 0 , 1 , 3 } - - - ( 14 )
A 2 1 &CenterDot; x Cm i 1 + B 2 1 &CenterDot; y Cm i 1 + C 2 1 &CenterDot; z Cm i 1 = 1 , i &Element; { 0 , 2 , 4 } - - - ( 15 )
A m top 1 &CenterDot; x Cm i 1 + B m top 1 &CenterDot; y Cm i 1 + C m top 1 &CenterDot; z Cm i 1 = 1 , i &Element; { 0,1,2,3,4 } , m &Element; { 1,2 , . . . , M } - - - ( 16 )
i在{C2}中的坐标为
Figure BDA000036516123001414
联立式(17)、式(18)与式(19),解得坐标
Figure BDA000036516123001415
标准板的下表面所在平面在{C2}中的方程如式(20),求解得到方程式(20)的参数
f x 2 0 u 0 2 0 f y 2 v 0 2 0 0 1 x Cm i 2 z Cm i 2 y Cm i 2 z Cm i 2 1 = u m i 2 v m i 2 1 - - - ( 17 )
A 1 2 &CenterDot; x Cm i 2 + B 1 2 &CenterDot; y Cm i 2 + C 1 2 &CenterDot; z Cm i 2 = 1 , i &Element; { 0 , 1 , 3 } - - - ( 18 )
A 2 2 &CenterDot; x Cm i 2 + B 2 2 &CenterDot; y Cm i 2 + C 2 2 &CenterDot; z Cm i 2 = 1 , i &Element; { 0 , 2 , 4 } - - - ( 19 )
A m bottom 2 &CenterDot; x Cm i 2 + B m bottom 2 &CenterDot; y Cm i 2 + C m bottom 2 &CenterDot; z Cm i 2 = 1 , i &Element; { 0,1,2,3,4 } , m &Element; { 1,2 , . . . , M } - - - ( 20 )
记标准板上表面的单位法向量为其在{C1}中的表达式为式(21),记标准板下表面法向的单位向量为
Figure BDA00003651612300153
其在{C2}中的表达式满足式(22),记{C2}到{C1}的转换矩阵为 R 2 1 T 2 1 , 标准板的上下表面平行,联立式(23),,
Figure BDA000036516123001513
M,求解得到
Figure BDA00003651612300155
a m 1 b m 1 c m 1 = A m top 1 B m top 1 C m top 1 ( A m top 1 ) 2 + ( B m top 1 ) 2 + ( C m top 1 ) 2 - - - ( 21 )
a m 2 b m 2 c m 2 = A m top 2 B m top 2 C m top 2 ( A m top 2 ) 2 + ( B m top 2 ) 2 + ( C m top 2 ) 2 - - - ( 22 )
R 2 1 a m 2 b m 2 c m 2 = a m 1 b m 1 c m 1 - - - ( 23 )
以激光位移传感器11的激光出射方向为Z轴负向,以激光位移传感器21的位移计算起始点到Z轴的垂向为X轴负向,建立坐标系{O},如图3所示,激光位移传感器11与激光位移传感器21的激光出射方向所在直线L1与L2的方程如式(24)与式(25)所示,记激光位移传感器11的位移计算起始点在坐标系中{O}坐标为(0,0,H);
记L1与S1的交点为P1,P1在{O}中的坐标为(0,0,H-1hm),记L2与S2的交点为P2,P2在{O}中的坐标为(l+2hmcosαcosβ,2hmcosαsinβ,2hmsinα);P1与P2分别在标准板上下表面所在平面S1与S2上,两平面在{O}中的方程如式(26)与式(27)所示;两平面间的距离为D,可得式(28);将P1与P2坐标代入式(26)与式(27),联立式(26)、式(27)、式(28),得到式(29);
x=0,y=0        (24)
x - l cos &alpha; cos &beta; = y cos &alpha; sin &beta; = z sin &alpha; - - - ( 25 )
Figure BDA000036516123001510
Figure BDA000036516123001511
D=k1m-k2m        (28)
Figure BDA000036516123001512
式(29)中的后三项为小量,可近似得到式(30),其中θm为{O}Z轴与标准板上表面所在平面的夹角;
sin &theta; m = D H - h m 1 - h m 2 - - - ( 30 )
记{O}的Z轴方向在{C1}中表示为向量 z &RightArrow; c = a c b c c c , 结合式(16)可得到与标准板上表面所在平面夹角θm,如式(31)所示;联立式(30)与式(31)得到式(32)形式的方程组,求解得到(acbccc)与H;
&theta; m = &pi; 2 - < n &RightArrow; m 1 , z &RightArrow; c > = &pi; 2 - arccos ( n &RightArrow; m 1 &CenterDot; z &RightArrow; c ) = &pi; 2 - arccos [ a m 1 b m 1 c m 1 &CenterDot; a c b c c c T ] - - - ( 31 )
a m 1 b m 1 c m 1 a c b c c c = D H - h m 1 - h m 2 , m = 1,2 , . . . , M - - - ( 32 )
4)对厚度进行检测
将待测板7放置于上表面测距单元1与下表面测距单元2之间,得到激光位移传感器11的测量值为1h,激光位移传感器21的测量值为2h,CCD摄像机13拍摄的图像为1I,CCD摄像机23拍摄的图像为2I;
图像1I中投射在被测板上表面的十字交叉激光的五个特征点分别为十字交叉点1P0与四个端点1P1,1P2,1P3,1P4,通过图像处理,得到1P0,1P1,1P2,1P3,1P4在{Im1}中的坐标分别为(1ui,1vi),其中i=0,1,2,3,4;
图像2I中投射在被测板下表面的十字交叉激光的五个特征点分别为十字交叉点2P0与四个端点2P1,2P2,2P3,2P4,通过图像处理,得到2P0,2P1,2P2,2P3,2P4在{Im1}中的坐标分别为(2ui,2vi),其中i=0,1,2,3,4;
1Pi在{C1}中的坐标为
Figure BDA00003651612300164
联立式(33)、式(34)与式(35),解得
Figure BDA00003651612300165
其中i=0,1,2,3,4;待测板7上表面所在平面在{C1}中的方程如式(36),联立解得(1Atop,1Btop,1Ctop);
f x 1 0 u 0 1 0 f y 1 v 0 1 0 0 1 x C i 1 z C i 1 y C i 1 z C i 1 1 = u i 1 v i 1 1 - - - ( 33 )
A 1 1 &CenterDot; x C i 1 + B 1 1 &CenterDot; y C i 1 + C 1 1 &CenterDot; z C i 1 = 1 , i &Element; { 0 , 1 , 3 } - - - ( 34 )
A 2 1 &CenterDot; x C i 1 + B 2 1 &CenterDot; y C i 1 + C 2 1 &CenterDot; z C i 1 = 1 , i &Element; { 0 , 2 , 4 } - - - ( 35 )
A top 1 &CenterDot; x C i 1 + B top 1 &CenterDot; y C i 1 + C top 1 &CenterDot; z C i 1 = 1 , i &Element; { 0,1,2,3,4 } - - - ( 36 )
2Pi在{C2}中的坐标为联立式(37)、式(38)与式(39),解得
Figure BDA00003651612300171
其中i=0,1,2,3,4;待测板7下表面所在平面在{C2}中的方程,如式(40),联立解得(2Abottom,2Bbottom,2Cbottom);
f x 2 0 u 0 2 0 f y 2 v 0 2 0 0 1 x C i 2 z C i 2 y C i 2 z C i 2 1 = u i 2 v i 2 1 - - - ( 37 )
A 1 2 &CenterDot; x C i 2 + B 1 2 &CenterDot; y C i 2 + C 1 2 &CenterDot; z C i 2 = 1 , i &Element; { 0 , 1 , 3 } - - - ( 38 )
A 2 2 &CenterDot; x C i 2 + B 2 2 &CenterDot; y C i 2 + C 2 2 &CenterDot; z C i 2 = 1 , i &Element; { 0 , 2 , 4 } - - - ( 39 )
A bottom 2 &CenterDot; x C i 2 + B bottom 2 &CenterDot; y C i 2 + C bottom 2 &CenterDot; z C i 2 = 1 , i &Element; { 0,1,2,3,4 } - - - ( 40 )
式(41)为{C1}中被测板上下表面法向量的平均值;由式(42)得到被测板上下表面所在平面与{O}的Z轴之间夹角θm
a 1 b 1 c 1 = a top 1 + a bottom 1 b top 1 + b bottom c top 1 + c bottom 1 ( a top 1 + a bottom 1 ) 2 + ( b top 1 + b bottom 1 ) 2 + ( c top 1 + c bottom 1 ) 2 - - - ( 41 )
&theta; = &pi; 2 - arccos [ a 1 b 1 c 1 &CenterDot; a c b c c c T ] - - - ( 42 )
被测板厚度由式(43)得到;
d=(H-1h-2h)sinθ        (43)
上述实施例仅用于说明本发明,其中激光结构光源波长、光条二维结构、上表面测距单元与下表面测距单元中的分光器件、CCD摄像机的图像传感器类型、以及上表面测距单元与下表面测距单元中分光器件-激光位移传感器-CCD摄像机-激光结构光源的相对位置等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。

Claims (10)

1.一种宽幅柔性板厚度实时检测装置,其特征在于:它包括上表面测距单元、下表面测距单元、上导轨、下导轨、上滑块和下滑块;
所述上导轨与所述下导轨平行放置,所述上滑块沿所述上导轨滑动,所述下滑块沿所述下导轨滑动,且所述上表面测距单元与所述上滑块固连,所述下表面测距单元与所述下滑块固连,且所述上滑块和所述下滑块同步运动;
所述上表面测距单元和下表面测距单元均包括固定在盒体内的激光位移传感器、激光结构光源、摄像机和半透半反镜;所述半透半反镜法线方向与所述激光位移传感器的光线出射方向成135°,与所述激光结构光源出射方向成45°,且三者在同一平面内;所述摄像机的镜头朝向所述半透半反镜。
2.如权利要求1所述的一种宽幅柔性板厚度实时检测装置,其特征在于:所述激光位移传感器接收端前方放置有一激光位移传感器滤光片,所述摄像机镜头前方放置有一摄像机滤光片。
3.如权利要求1或2所述的一种宽幅柔性板厚度实时检测装置,其特征在于:所述激光位移传感器的激光波长为λ1,所述激光结构光源采用单圆、同心多圆、十字交叉、网格和网点中的一种,激光波长为λ2,并满足下式:
12|≥20nm。
4.如权利要求3所述的一种宽幅柔性板厚度检测装置,其特征在于:所述激光位移传感器滤光片为窄带通滤光片,中心波长为λ1,半带宽不大于10nm;所述摄像机滤光片为窄带通滤光片,中心波长为λ2,半带宽不大于10nm。
5.如权利要求1到4任一项所述的一种宽幅柔性板厚度检测装置,其特征在于:所述摄像机包括镜头与图像传感器,其中所述图像传感器采用CCD、CMOS和光电转换器件中的一种。
6.一种采用如权利要求1到5任一项所述装置实现的宽幅柔性板厚度实时检测方法,其包括以下步骤:
1)对上表面测距单元进行标定;
2)对下表面测距单元进行标定;
3)对厚度检测装置进行标定;
4)对厚度进行检测。
7.如权利要求6所述的一种宽幅柔性板厚度实时检测方法,其特征在于:所述步骤1)中对上表面测距单元进行标定包括以下步骤:
①记摄像机光心坐标系为{C1},摄像机对应的图像坐标系为{Im1},CCD摄像机标定过程中第n个标定板空间位置对应的标定板坐标系为{S1n};
②对摄像机进行标定,得到摄像机内部参数{1fx,1fy,1u0,1v0},其中1fx,1fy为镜头光心坐标系X,Y方向上的等效焦距,(1u0,1v0)为镜头光轴与图像平面交点在图像坐标系中的坐标;
③对激光结构光源与摄像机相对位置进行标定:将标定板放置于摄像机前方,将激光结构光投射于标定板上,变换标定板位置,拍摄一系列图像
Figure FDA00003651612200021
其中n=1,2,...,N,N≥4,右上标calib表示该图像用于标定,提取图像
Figure FDA00003651612200022
激光光条上不共线的五个特征点,分别记为其位置关系满足:
Figure FDA00003651612200024
共线、
Figure FDA00003651612200025
共线;特征点在图像坐标系中的坐标记为将摄像机光心坐标系记为1OC1XC1YC1ZC,第n个标定板坐标系记为1On1Xn1Yn1Zn,将两坐标系之间的转换矩阵记为[1Rn  1Tn],标定板的各个角点在1OC1XC1YC1ZC中的坐标记为
Figure FDA00003651612200027
由标定板图案角点位置关系知其在1On1Xn1Yn1Zn中的坐标为
Figure FDA00003651612200028
图像
Figure FDA00003651612200029
上角点的像素坐标为三者关系满足以下两式:
R n 1 T n 1 x n k 1 y n k 1 0 1 = x C k 1 y C k 1 z C k 1
f x 1 0 u 0 1 0 f y 1 v 0 1 0 0 1 x C k 1 z C k 1 y C k 1 z C k 1 1 = u n k 1 v n k 1 1
联立求解得到转换矩阵[1Rn  1Tn]
记点
Figure FDA000036516122000213
1OC1XC1YC1ZC中的坐标为
Figure FDA000036516122000214
由以下两式联立求解得到:
R n 1 T n 1 x n i 1 y n i 1 0 1 = x Cn i 1 y Cn i 1 z Cn i 1
f x 1 0 u 0 1 0 f y 1 v 0 1 0 0 1 x Cn i 1 z Cn i 1 y Cn i 1 z Cn i 1 1 = u n i 1 v n i 1 1
在激光结构光的两个光平面上,光平面在1OC1XC1YC1ZC中的方程满足以下两式:
A 1 1 &CenterDot; x Cn i 1 + B 1 1 &CenterDot; y Cn i 1 + C 1 1 &CenterDot; z Cn i 1 = 1 , i &Element; { 0,1,3 } , n = { 1,2 , . . . , N }
A 2 1 &CenterDot; x Cn i 1 + B 2 1 &CenterDot; y Cn i 1 + C 2 1 &CenterDot; z Cn i 1 = 1 , i &Element; { 0,1,3 } , n = { 1,2 , . . . , N }
联立求解得到(1A1,1B1,1C1,1A2,1B2,1C2)。
8.如权利要求6所述的一种用于宽幅柔性板厚度实时检测方法,其特征在于:所述步骤2)中对下表面测距单元进行标定包括以下步骤:
①记摄像机光心坐标系为{C2},摄像机对应的图像坐标系为{Im2},摄像机标定过程中第n个标定板空间位置对应的标定板坐标系为{S2n};
②对摄像机进行标定,得到摄像机内部参数{2fx,2fy,2u0,2v0},其中2fx,2fy为镜头光心坐标系X,Y方向上的等效焦距,(2u0,2v0)为镜头光轴与图像平面交点在图像坐标系中的坐标;
③对激光结构光源与摄像机的相对位置进行标定:将标定板放置于摄像机前方,将激光结构光投射于标定板上,变换标定板位置,拍摄一系列图像其中n=1,2,...,N,N≥4,右上标calib表示该图像用于标定,提取图像
Figure FDA00003651612200036
激光光条上不共线的五个特征点,分别记为其位置关系满足:
Figure FDA00003651612200038
共线,
Figure FDA00003651612200039
共线;特征点在图像坐标系中的坐标记为
Figure FDA000036516122000310
将摄像机光心坐标系记为2OC2XC2YC2ZC,第n个标定板坐标系记为2On2Xn2Yn2Zn,将两坐标系之间的转换矩阵记为[2Rn  2Tn],标定板的各个角点在2OC2XC2YC2ZC中的坐标记为
Figure FDA000036516122000311
由标定板图案角点位置关系知其在2On2Xn2Yn2Zn中的坐标为
Figure FDA000036516122000312
图像
Figure FDA000036516122000313
上角点的像素坐标为
Figure FDA000036516122000314
三者关系满足以下两式:
R n 2 T n 2 x n k 2 y n k 2 0 1 = x C k 2 y C k 2 z C k 2
f x 2 0 u 0 2 0 f y 2 v 0 2 0 0 1 x C k 2 z C k 2 y C k 2 z C k 2 1 = u n k 2 v n k 2 1
联立求解得到转换矩阵[2Rn  2Tn];记点2OC2XC2YC2ZC中的坐标为由以下两式联立求解得到:
R n 2 T n 2 x n i 2 y n i 2 0 1 = x Cn i 2 y Cn i 2 z Cn i 2
f x 2 0 u 0 2 0 f y 2 v 0 2 0 0 1 x Cn i 2 z Cn i 2 y Cn i 2 z Cn i 2 1 = u n i 2 v n i 2 1
Figure FDA00003651612200046
在激光结构光的两个光平面上,光平面在2OC2XC2YC2ZC中的方程满足以下两式:
A 1 2 &CenterDot; x Cn i 2 + B 1 2 &CenterDot; y Cn i 2 + C 1 2 &CenterDot; z Cn i 2 = 1 , i &Element; { 0,1,3 } , n = { 1,2 , . . . , N }
A 2 2 &CenterDot; x Cn i 2 + B 2 2 &CenterDot; y Cn i 2 + C 2 2 &CenterDot; z Cn i 2 = 1 , i &Element; { 0 , 2 , 4 } , n = { 1,2 , . . . , N }
联立求解得到(2A1,2B1,2C1,2A2,2B2,2C2)。
9.如权利要求6所述的一种用于宽幅柔性板厚度实时检测方法,其特征在于:所述步骤3)中对厚度检测装置进行标定包括以下步骤:
将厚度为D的标准板置于上下正对安装的上表面测距单元与下表面测距单元之间,变化标准板的角度,从上表面测距单元与下表面测距单元中分别得到一系列测量值{1hm}与{2hm};上表面测距单元中的激光结构光源投射在标准板上表面的五个特征点在上表面测距单元中的摄像机拍摄的图像中坐标分别为
Figure FDA000036516122000410
下表面测距单元中的激光结构光源投射在标准板下表面的五个特征点
Figure FDA000036516122000411
在下表面测距单元中的摄像机拍摄的图像中坐标分别为
Figure FDA000036516122000412
i=0,1,2,3,4,m=1,2,...,M,M为标定次数;
Figure FDA000036516122000413
在上表面测距单元中摄像机光心坐标系1OC1XC1YC1ZC中的坐标
Figure FDA000036516122000414
满足以下三式:
f x 1 0 u 0 1 0 f y 1 v 0 1 0 0 1 x Cm i 1 z Cm i 1 y Cm i 1 z Cm i 1 1 = u m i 1 v m i 1 1
A 1 1 &CenterDot; x Cm i 1 + B 1 1 &CenterDot; y Cm i 1 + C 1 1 &CenterDot; z Cm i 1 = 1 , i &Element; { 0 , 1 , 3 }
A 2 1 &CenterDot; x Cm i 1 + B 2 1 &CenterDot; y Cm i 1 + C 2 1 &CenterDot; z Cm i 1 = 1 , i &Element; { 0 , 2 , 4 }
联立求解得到坐标满足下式:
A m top 1 &CenterDot; x Cm i 1 + B m top 1 &CenterDot; y Cm i 1 + C m top 1 &CenterDot; z Cm i 1 = 1 , i &Element; { 0,1,2,3,4 } , m &Element; { 1,2 , . . . , M }
解得
Figure FDA00003651612200056
Figure FDA00003651612200057
在下表面测距单元中摄像机光心坐标系2OC2XC2YC2ZC中的坐标
Figure FDA00003651612200058
满足以下三式:
f x 2 0 u 0 2 0 f y 2 v 0 2 0 0 1 x Cm i 2 z Cm i 2 y Cm i 2 z Cm i 2 1 = u m i 2 v m i 2 1
A 1 2 &CenterDot; x Cm i 2 + B 1 2 &CenterDot; y Cm i 2 + C 1 2 &CenterDot; z Cm i 2 = 1 , i &Element; { 0 , 1 , 3 }
A 2 2 &CenterDot; x Cm i 2 + B 2 2 &CenterDot; y Cm i 2 + C 2 2 &CenterDot; z Cm i 2 = 1 , i &Element; { 0 , 2 , 4 }
联立求解得到坐标
Figure FDA000036516122000512
满足下式:
A m bottom 2 &CenterDot; x Cm i 2 + B m bottom 2 &CenterDot; y Cm i 2 + C m bottom 2 &CenterDot; z Cm i 2 = 1 , i &Element; { 0,1,2,3,4 } , m &Element; { 1,2 , . . . , M }
解得
Figure FDA000036516122000514
记标准板上表面法向的单位向量为
Figure FDA000036516122000515
其在1OC1XC1YC1ZC坐标系的表达式满足下式:
a m 1 b m 1 c m 1 = A m top 1 B m top 1 C m top 1 ( A m top 1 ) 2 + ( B m top 1 ) 2 + ( C m top 1 ) 2
记标准板下表面法向的单位向量为
Figure FDA000036516122000517
其在1OC1XC1YC1ZC坐标系的表达式满足下式:
a m 2 b m 2 c m 2 = A m top 2 B m top 2 C m top 2 ( A m top 2 ) 2 + ( B m top 2 ) 2 + ( C m top 2 ) 2
2OC2XC2YC2ZC坐标系到1OC1XC1YC1ZC坐标系的转换矩阵为 R 2 1 T 2 1 , 标准板的上下表面平行,上述两单位向量的表达式满足下式;
R 2 1 a m 2 b m 2 c m 2 = a m 1 b m 1 c m 1 , m = 1,2 , . . . , M
求解得到
Figure FDA00003651612200062
以上表面测距单元的激光位移传感器的激光出射方向为Z轴负向,以下表面测距单元的激光位移传感器的位移计算起始点到Z轴的垂向为X轴负向,建立坐标系OXYZ,上表面测距单元的激光位移传感器与下表面测距单元的激光位移传感器的激光出射方向所在直线L1与L2的方程满足以下两式:
x=0,y=0
x - l cos &alpha; cos &beta; = y cos &alpha; sin &beta; = z sin &alpha;
记上表面测距单元的激光位移传感器的位移计算起始点坐标为(0,0,H),记L1与S1的交点为P1(0,0,H-1hm),记L2与S2的交点为P2(l+2hmcosαcosβ,2hmcosαsinβ,2hmsinα);P1与P2分别在标准板上下表面所在平面S1与S2上,两平面在OXYZ坐标系中的方程满足以下两式:
Figure FDA00003651612200064
Figure FDA00003651612200065
两平面间的距离为D,满足下式:
D=k1m-k2m
联立上述三式,得到下式:
Figure FDA00003651612200066
θm为OXYZ坐标系的Z轴与标准板上表面所在平面的夹角,近似得到下式:
sin &theta; m = D H - h m 1 - h m 2
记坐标系OXYZ的Z轴方向在1OC1XC1YC1ZC坐标系中表示为向量 z &RightArrow; c = a c b c c c , 满足下式:
&theta; m = &pi; 2 - < n &RightArrow; m 1 , z &RightArrow; c > = &pi; 2 - arccos ( n &RightArrow; m 1 &CenterDot; z &RightArrow; c ) = &pi; 2 - arccos [ a m 1 b m 1 c m 1 &CenterDot; a c b c c c T ]
联立上述两式得到下式:
a m 1 b m 1 c m 1 a c b c c c = D H - h m 1 - h m 2 , m = 1,2 , . . . , M
求解得到(acbccc)与H。
10.如权利要求6所述的一种用于宽幅柔性板厚度实时检测方法,其特征在于:所述步骤4)中对厚度进行检测包括以下步骤:
将待测板放置于上表面测距单元与下表面测距单元之间,得到上表面测距单元中的激光位移传感器的测量值为1h,上表面测距单元中的激光位移传感的测量值为2h,下表面测距单元中的摄像机拍摄的图像为1I,下表面测距单元中的摄像机拍摄的图像为2I;
图像1I中投射在待测板上表面的激光的五个特征点分别为点1P0与四个端点1P1,1P2,1P3,1P4,通过图像处理,得到1P0,1P1,1P2,1P3,1P4的图像坐标分别为(1ui,1vi),i=0,1,2,3,4;
图像2I中投射在待测板下表面的激光的五个特征点分别为点2P0与四个端点2P1,2P2,2P3,2P4,通过图像处理,得到2P0,2P1,2P2,2P3,2P4的图像坐标分别为(2ui,2vi),i=0,1,2,3,4;
1Pi在上表面测距单元中的摄像机光心坐标系1OC1XC1YC1ZC中的坐标
Figure FDA00003651612200071
满足以下三式:
f x 1 0 u 0 1 0 f y 1 v 0 1 0 0 1 x C i 1 z C i 1 y C i 1 z C i 1 1 = u i 1 v i 1 1
A 1 1 &CenterDot; x C i 1 + B 1 1 &CenterDot; y C i 1 + C 1 1 &CenterDot; z C i 1 = 1 , i &Element; { 0 , 1 , 3 }
A 2 1 &CenterDot; x C i 1 + B 2 1 &CenterDot; y C i 1 + C 2 1 &CenterDot; z C i 1 = 1 , i &Element; { 0 , 2 , 4 }
联立解得坐标坐标满足下式;
A top 1 &CenterDot; x C i 1 + B top 1 &CenterDot; y C i 1 + C top 1 &CenterDot; z C i 1 = 1 , i &Element; { 0,1,2,3,4 }
联立解得
Figure FDA00003651612200077
2Pi在下表面测距单元中的摄像机光心坐标系2OC2XC2YC2ZC中的坐标
Figure FDA00003651612200078
满足以下三式:
f x 2 0 u 0 2 0 f y 2 v 0 2 0 0 1 x C i 2 z C i 2 y C i 2 z C i 2 1 = u i 2 v i 2 1
A 1 2 &CenterDot; x C i 2 + B 1 2 &CenterDot; y C i 2 + C 1 2 &CenterDot; z C i 2 = 1 , i &Element; { 0 , 1 , 3 }
A 2 2 &CenterDot; x C i 2 + B 2 2 &CenterDot; y C i 2 + C 2 2 &CenterDot; z C i 2 = 1 , i &Element; { 0 , 2 , 4 }
联立解得坐标
Figure FDA00003651612200083
坐标满足下式:
A bottom 2 &CenterDot; x C i 2 + B bottom 2 &CenterDot; y C i 2 + C bottom 2 &CenterDot; z C i 2 = 1 , i &Element; { 0,1,2,3,4 }
联立解得(2Abottom,2Bbottom,2Cbottom);
坐标系1OC1XC1YC1ZC中待测板上下表面局部所在平面的法向量满足下式:
a 1 b 1 c 1 = a top 1 + a bottom 1 b top 1 + b bottom c top 1 + c bottom 1 ( a top 1 + a bottom 1 ) 2 + ( b top 1 + b bottom 1 ) 2 + ( c top 1 + c bottom 1 ) 2
待测板上下表面局部所在平面与坐标系OXYZ的Z轴之间夹角θ满足下式:
&theta; = &pi; 2 - arccos [ a 1 b 1 c 1 &CenterDot; a c b c c c T ]
待测板厚度由下式得到:
d=(H-1h-2h)sinθ。
CN201310348728.9A 2013-08-12 2013-08-12 一种宽幅柔性板厚度实时检测装置及方法 Expired - Fee Related CN103438810B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310348728.9A CN103438810B (zh) 2013-08-12 2013-08-12 一种宽幅柔性板厚度实时检测装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310348728.9A CN103438810B (zh) 2013-08-12 2013-08-12 一种宽幅柔性板厚度实时检测装置及方法

Publications (2)

Publication Number Publication Date
CN103438810A true CN103438810A (zh) 2013-12-11
CN103438810B CN103438810B (zh) 2016-01-06

Family

ID=49692506

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310348728.9A Expired - Fee Related CN103438810B (zh) 2013-08-12 2013-08-12 一种宽幅柔性板厚度实时检测装置及方法

Country Status (1)

Country Link
CN (1) CN103438810B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107816942A (zh) * 2017-09-28 2018-03-20 中国东方电气集团有限公司 一种基于十字结构光视觉系统的平面二维尺寸测量方法
CN108801914A (zh) * 2018-05-29 2018-11-13 华中科技大学 一种对多沟槽型面板材成形缺陷的检测方法及检测系统
CN110332896A (zh) * 2019-04-25 2019-10-15 太仓北新建材有限公司 一种石膏板厚度监测装置
CN110411358A (zh) * 2019-07-30 2019-11-05 歌尔股份有限公司 深度相机测量方法、装置、设备及计算机可读存储介质
CN110455238A (zh) * 2019-08-12 2019-11-15 佛山市景旭机电设备有限公司 橡塑卷材在线自动检测厚度的装置及方法
CN111735487A (zh) * 2020-05-18 2020-10-02 清华大学深圳国际研究生院 传感器、传感器标定方法与设备、存储介质
CN111886475A (zh) * 2018-03-16 2020-11-03 明产株式会社 片材厚度测量装置
CN112013784A (zh) * 2019-05-28 2020-12-01 深圳中科飞测科技有限公司 一种检测方法和检测系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2911739A1 (de) * 1979-03-26 1980-10-02 Wilhelm Becker Punktplanimeterartiges geraet
JPH08304032A (ja) * 1994-06-27 1996-11-22 Honshu Paper Co Ltd コンデンサ用金属化フィルムのマージン部幅および金属蒸着膜厚さ測定方法
JP2000321546A (ja) * 1999-05-07 2000-11-24 Sharp Corp 液晶表示装置、そのセル厚測定装置及び測定法並びに位相差板
CN1800776A (zh) * 2004-12-30 2006-07-12 中南大学 微米级精度激光测厚装置
CN202267459U (zh) * 2011-09-28 2012-06-06 中电装备东芝(常州)变压器有限公司 硅钢片叠厚测量装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2911739A1 (de) * 1979-03-26 1980-10-02 Wilhelm Becker Punktplanimeterartiges geraet
JPH08304032A (ja) * 1994-06-27 1996-11-22 Honshu Paper Co Ltd コンデンサ用金属化フィルムのマージン部幅および金属蒸着膜厚さ測定方法
JP2000321546A (ja) * 1999-05-07 2000-11-24 Sharp Corp 液晶表示装置、そのセル厚測定装置及び測定法並びに位相差板
CN1800776A (zh) * 2004-12-30 2006-07-12 中南大学 微米级精度激光测厚装置
CN202267459U (zh) * 2011-09-28 2012-06-06 中电装备东芝(常州)变压器有限公司 硅钢片叠厚测量装置

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107816942A (zh) * 2017-09-28 2018-03-20 中国东方电气集团有限公司 一种基于十字结构光视觉系统的平面二维尺寸测量方法
CN111886475A (zh) * 2018-03-16 2020-11-03 明产株式会社 片材厚度测量装置
CN111886475B (zh) * 2018-03-16 2022-06-03 明产株式会社 片材厚度测量装置
CN108801914A (zh) * 2018-05-29 2018-11-13 华中科技大学 一种对多沟槽型面板材成形缺陷的检测方法及检测系统
CN110332896A (zh) * 2019-04-25 2019-10-15 太仓北新建材有限公司 一种石膏板厚度监测装置
CN112013784A (zh) * 2019-05-28 2020-12-01 深圳中科飞测科技有限公司 一种检测方法和检测系统
CN110411358A (zh) * 2019-07-30 2019-11-05 歌尔股份有限公司 深度相机测量方法、装置、设备及计算机可读存储介质
CN110411358B (zh) * 2019-07-30 2021-08-10 歌尔光学科技有限公司 深度相机测量方法、装置、设备及计算机可读存储介质
CN110455238A (zh) * 2019-08-12 2019-11-15 佛山市景旭机电设备有限公司 橡塑卷材在线自动检测厚度的装置及方法
CN111735487A (zh) * 2020-05-18 2020-10-02 清华大学深圳国际研究生院 传感器、传感器标定方法与设备、存储介质

Also Published As

Publication number Publication date
CN103438810B (zh) 2016-01-06

Similar Documents

Publication Publication Date Title
CN103438810B (zh) 一种宽幅柔性板厚度实时检测装置及方法
CN104567679B (zh) 一种涡轮叶片视觉检测的系统
US9019351B2 (en) Three-dimensional image measuring apparatus
CN102941410B (zh) 一种点扫描三维测量系统振镜标定方法
CN107192340A (zh) 一种测量孔、槽的位置与几何尺寸的机器视觉系统
CN103115580B (zh) 基于光学相干层析扫描的三维孔形检测方法及系统
CN101373134B (zh) 一种线状激光窄带滤光psd热辐射板厚测量方法
CN101910782B (zh) 检测移动衬底的形状变化的方法和装置
CN102500627A (zh) 一种板材的宽度、边部形状、边部缺陷测量仪及测量方法
CN101900744A (zh) 一种用于粒子图像测速仪的三维激光测准定位仪
CN113394141B (zh) 一种芯片结构缺陷的质量评估系统及方法
CN104613883A (zh) 金属薄板印刷涂层的湿膜厚度差分测量及均匀性评估方法
CN109632828A (zh) 一种平板玻璃缺陷复检系统及复检方法
CN203069151U (zh) 基于光学相干层析扫描的三维孔形检测系统
CN108775863A (zh) 光谱共焦位移传感器定位装置及定位方法
CN105519264B (zh) 非接触式综合测量大型平板裂缝天线表面三维形貌和缝隙尺寸的方法及其系统
CN105783743B (zh) 基于红外反射法的金属薄板印涂湿膜厚度在线检测系统
JP5535095B2 (ja) ワーク寸法測定装置
CN102095381A (zh) 一种面向刚性透明材质的激光测厚方法
CN103822580B (zh) 超长框架变形与姿态的多点实时测量系统与方法
CN105180827A (zh) 一种检测起重机主梁上拱曲线的装置
TWI306149B (en) Optical device for sensing distance
CN109186503A (zh) 基于四面锥形直角等腰反射棱镜的直线度误差测量系统
CN110455181B (zh) 一种位姿快速测量系统及方法
CN103256893B (zh) 测量太赫兹时域光谱系统光斑位置和大小的测量仪及方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160106