CN108592791B - 一种凹坑检查方法 - Google Patents

一种凹坑检查方法 Download PDF

Info

Publication number
CN108592791B
CN108592791B CN201810390814.9A CN201810390814A CN108592791B CN 108592791 B CN108592791 B CN 108592791B CN 201810390814 A CN201810390814 A CN 201810390814A CN 108592791 B CN108592791 B CN 108592791B
Authority
CN
China
Prior art keywords
pit
array
ccd
coordinate system
laser array
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810390814.9A
Other languages
English (en)
Other versions
CN108592791A (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.)
Hon Hai Precision Industry Co Ltd
Original Assignee
Hon Hai Precision Industry Co Ltd
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 Hon Hai Precision Industry Co Ltd filed Critical Hon Hai Precision Industry Co Ltd
Priority to CN201810390814.9A priority Critical patent/CN108592791B/zh
Publication of CN108592791A publication Critical patent/CN108592791A/zh
Application granted granted Critical
Publication of CN108592791B publication Critical patent/CN108592791B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/002Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/02Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
    • G01B11/03Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness by measuring coordinates of points
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/22Measuring arrangements characterised by the use of optical techniques for measuring depth

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种凹坑检查系统,它由机架、CCD、激光器阵列、线路板和计算机组成。激光器阵列安装在机架下面的中央,两侧分别安装CCD。激光器阵列由个排成阵列的激光器组成。开启激光器阵列的所有激光器,在待检体的表面出现个阵列式排列的激光光斑。两个CCD同时拍摄光斑阵列的图像。利用CCD成像的特点可以确定照射在待检体上的光斑中心在测量坐标系中的位置及其最大深度,由此可确定待检体上有无凹坑,并在有凹坑的情况下确定凹坑中的点。通过寻求凹坑中点的边界,可以进一步确定凹坑的最小宽度以及波浪度。本发明的有益效果在于:1.无接触式测量,不损伤待检体表面;2.测量系统结构简单,造价低廉;3.引入三维测量坐标系,因此能够确定凹坑中点的位置,为获取其边界并求出凹坑最小宽度提供了可能。

Description

一种凹坑检查方法
技术领域
本发明涉及一种凹坑检查方法。
背景技术
飞机在起飞和着陆过程中,蒙皮、壁板和整流包皮因受外来物撞击,例如空中的飞鸟,或者其它意外事故,会形成局部凹坑。凹坑将增加飞机的飞行阻力,影响飞机的飞行性能。在加载条件下,凹坑处产生的次级弯曲应力很大,将导致疲劳强度大大下降。波浪度是凹坑的最大深度处的最小宽度与凹坑的最大深度的比值,是凹坑损伤非常重要的一个参数,也是凹坑超过允许损伤与否的重要判据。波浪度的大小直接反映了凹坑的尖锐程度,波浪度越小,凹坑越尖锐,在加载条件下,凹坑处产生的次级弯曲应力就越大。此外,凹坑的波浪度越小,对飞机产生的飞行阻力越大,对飞机飞行性能的影响也越大。因此,如何探测凹坑并求出凹坑的波浪度,是对飞机蒙皮进行及时修理以保证飞行安全的极其重要的环节。
发明内容
为了探测凹坑,并在凹坑存在的情况下确定其波浪度,本发明公开了一种凹坑检查系统,它由机架、CCD 1、CCD 2、激光器阵列、线路板和计算机组成。激光器阵列安装在机架下面的中央,两侧分别安装CCD 1和CCD 2。激光器阵列由n×n个排成阵列的激光器组成。CCD 1和CCD 2的像平面位于同一平面内,且对应的水平像素线在同一直线上。
进行凹坑检查时,凹坑检查系统放在待检体上,开启激光器阵列的所有激光器,在待检体的表面出现n×n个阵列式排列的激光光斑。CCD 1和CCD 2同时拍摄光斑阵列的图像。
建立两个像平面坐标系和一个三维测量坐标系之后,利用CCD成像的特点可以确定照射在待检体上的光斑中心在测量坐标系中的位置及其最大深度,由此可确定待检体上有无凹坑,并在有凹坑的情况下确定凹坑中的点。通过寻求凹坑中的点边界,可以进一步确定凹坑的最小宽度以及波浪度。
本发明的有益效果在于:
1.无接触式测量,不损伤待检体表面;
2.测量系统结构简单,造价低廉;
3.引入三维测量坐标系,因此能够确定凹坑中点的位置,为获取其边界并求出凹坑最小宽度提供了可能。
附图说明
图1是凹坑检查系统工作原理示意图,图2是凹坑边界点及最小宽度示意图。
标号说明
1机架,2CCD 1,3激光器阵列,4CCD 2,5待检体,6激光器,7激光束,8凹坑。
具体实施方式
如图1所示,凹坑检查系统由机架(1)、CCD 1(2)、激光器阵列(3)、CCD 2(4)、线路板和计算机组成。激光器阵列(3)安装在机架(1)下面的中央,两侧分别安装CCD 1(2)和CCD2(4)。激光器阵列(3)由n×n个排成阵列的激光器组成。CCD 1(2)和CCD 2(4)的像平面位于同一平面内,且对应的水平像素线在同一直线上。
进行凹坑检查时,凹坑检查系统放在待检体(5)上,开启激光器阵列(3)的所有激光器,在待检体(5)的表面出现n×n个阵列式排列的激光光斑。
CCD 1(2)和CCD 2(4)同时拍摄光斑阵列的图像,计算机通过图像处理确定有无凹坑,并在有凹坑的情况下确定凹坑的最大深度Y、最小宽度A以及波浪度λ=A/Y。
为此,需要建立两个像平面坐标系和一个三维测量坐标系,确定光斑在三维测量坐标系中的位置、处在凹坑中的光斑、凹坑最大深度Y、凹坑中光斑的边界,并由凹坑中光斑的边界确定凹坑最小宽度A。
1.坐标系的建立及光斑在三维测量坐标系中的位置
以o1和o2分别表示CCD 1(2)和CCD 2(4)像面中点,2d表示o1,o2之间的距离,以o1,o2的连线作为x轴,o1,o2连线的中点作为坐标原点o,建立三维测量坐标系o-xyz.在两个像面内分别建立像平面坐标系o1-x1y1以及o2-x2y2,x1和x2轴与测量坐标系的x轴重合且方向相同。于是两个CCD像平面的中心在测量坐标系中的坐标分别为(-d,0,0)和(d,0,0)。
以f,l分别表示两个CCD的焦距和像素长度。设点目标P在测量坐标系中的坐标为A(x,y,z),在两个像平面坐标系中的坐标分别为P(m,n)和P(m',n'),其中|m|、 |n|、|m'|、|n'|分别为两个CCD横、纵方向上的像素个数。则由小孔成像原理知
Figure GDA0002417896620000021
Figure GDA0002417896620000031
Figure GDA0002417896620000032
Figure GDA0002417896620000033
因此在测量坐标系中
Figure GDA0002417896620000034
Figure GDA0002417896620000035
Figure GDA0002417896620000036
对两幅图像进行图像处理,获取阵列光斑的中心在两个像平面坐标系中的像素位置pij(mij,nij)和pij'(mij',nij'),并由此获取阵列光斑在三维测量坐标系中的位置
Figure GDA0002417896620000037
Figure GDA0002417896620000038
Figure GDA0002417896620000039
(i,j=1,2,…,n)
2.处在凹坑中的光斑及凹坑最大深度Y
取一个很小的正数δ作为判别阈值,令
Figure GDA00024178966200000310
Figure GDA00024178966200000311
则当M-m>δ令Y=M-m,否则Y=0;而当zij>m+δ时判定点(xij,yij,zij) (i,j=1,2,…,n)位于凹坑之内。以N表示位于凹坑内点的个数,则所有凹坑之内点的前两个坐标可表示为S={(Xi,Yi)|i=1,2,…,N}.
3.凹坑中光斑的边界
当N>2时,如图2所示,在S中取
Figure GDA00024178966200000312
使
Figure 12
Figure GDA00024178966200000314
是边界上的第一个点。以L表示依进入次序排列的边界上点的集合,则此时
Figure DEST_PATH_10
Figure GDA0002417896620000042
中取
Figure GDA0002417896620000043
使
Figure 100002_5
其中“·”表示向量内积运算。令
Figure 100002_8
这里“
Figure GDA00024178966200000415
”为赋值号。继续执行如下程序:
1)k=2;
2)在S中取
Figure GDA0002417896620000046
使
Figure 100002_9
Figure GDA0002417896620000048
转3);否则令
Figure 100002_10
转2);
3)停止。
此时光斑的边界由L中的k个点确定。
4.凹坑最小宽度以及波浪度
当N≤2时,令A=x11-x12;当N>2时,以lj表示
Figure 13
中的点到直线
Figure GDA00024178966200000411
之间的距离的最大值,j=1,2,…,k-1,lk表示
Figure 14
中的点到直线
Figure GDA00024178966200000413
之间的距离的最大值,令
Figure 590450DEST_PATH_5
由此可以得到凹坑探测结果:当Y=0时无凹坑;当Y≠0时有凹坑,其波浪度为λ=A/Y。

Claims (1)

1.一种凹坑检查方法,其特征在于:凹坑检查系统包括机架(1)、CCD 1(2)、激光器阵列(3)、CCD 2(4)、线路板和计算机;激光器阵列(3)安装在机架(1)下面的中央,两侧分别安装CCD 1(2)和CCD 2(4);激光器阵列(3)由n×n个排成阵列的激光器组成,CCD 1(2)和CCD 2(4)的像平面位于同一平面内,且对应的水平像素线在同一直线上;凹坑检查系统放在待检体(5)上,开启激光器阵列(3)的所有激光器,CCD 1(2)和CCD 2(4)同时拍摄待检体(5)上光斑阵列的图像,并送计算机处理,处理过程分如下步骤:
1)建立坐标系并确定光斑在三维测量坐标系中的位置
以o1和o2分别表示CCD 1(2)和CCD 2(4)像面中点,2d表示o1,o2之间的距离,以o1,o2的连线作为x轴,o1,o2连线的中点作为坐标原点o,建立三维测量坐标系o-xyz,在两个像面内分别建立像平面坐标系o1-x1y1以及o2-x2y2,x1和x2轴与测量坐标系的x轴重合且方向相同;以f,l分别表示两个CCD的焦距和像素长度,对i,j=1,2,…,n分别以pij(mij,nij)和pij'(mij',nij')表示阵列光斑的中心在两个像平面坐标系中的像素位置,则在三维测量坐标系中的坐标为
Figure RE-FDA0002480686750000011
Figure RE-FDA0002480686750000012
Figure RE-FDA0002480686750000013
2)确定凹坑中的光斑及凹坑最大深度Y
取一个很小的正数δ作为判别阈值,令
Figure RE-FDA0002480686750000014
Figure RE-FDA0002480686750000015
当M-m>δ令Y=M-m,否则Y=0;而对于i,j=1,2,…,n,当zij>m+δ时判定点(xij,yij,zij)位于凹坑之内,以N表示位于凹坑内点的个数,则所有凹坑之内点的前两个坐标可表示为S={(Xi,Yi)|i=1,2,…,N};
3)确定凹坑中光斑的边界
当N>2时,在S中取
Figure RE-FDA0002480686750000021
使
Figure RE-FDA0002480686750000022
Figure RE-FDA0002480686750000023
是边界上的第一个点,以L表示依次选入的边界上点的集合,则此时
Figure 10
Figure RE-FDA0002480686750000025
中取
Figure RE-FDA0002480686750000026
使
Figure RE-FDA0002480686750000027
其中“·”表示向量内积运算;令
Figure 9
这里
Figure RE-FDA0002480686750000029
为赋值号;继续执行如下步骤:
i)k=2;
ii)在S中取
Figure RE-FDA00024806867500000210
使
Figure RE-FDA00024806867500000211
Figure RE-FDA00024806867500000212
转iii);否则令
Figure 8
转ii);
iii)停止;
此时光斑的边界由L中的k个点确定;
4)确定凹坑最小宽度以及波浪度
当N≤2时,令A=x11-x12;当N>2时,以lj表示
Figure 7
中的点到直线
Figure RE-FDA0002480686750000031
之间的距离的最大值,j=1,2,…,k-1,lk表示
Figure 3
中的点到直线
Figure RE-FDA0002480686750000033
之间的距离的最大值,令
Figure 5
由此可以得到凹坑探测结果:当Y=0时无凹坑;当Y≠0时有凹坑,其波浪度为λ=A/Y。
CN201810390814.9A 2018-04-27 2018-04-27 一种凹坑检查方法 Active CN108592791B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810390814.9A CN108592791B (zh) 2018-04-27 2018-04-27 一种凹坑检查方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810390814.9A CN108592791B (zh) 2018-04-27 2018-04-27 一种凹坑检查方法

Publications (2)

Publication Number Publication Date
CN108592791A CN108592791A (zh) 2018-09-28
CN108592791B true CN108592791B (zh) 2020-06-16

Family

ID=63610529

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810390814.9A Active CN108592791B (zh) 2018-04-27 2018-04-27 一种凹坑检查方法

Country Status (1)

Country Link
CN (1) CN108592791B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114441549B (zh) * 2021-12-31 2024-02-23 武汉钢铁有限公司 钢卷质量检测系统及方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5589942A (en) * 1990-04-05 1996-12-31 Intelligent Automation Systems Real time three dimensional sensing system
CN1587900A (zh) * 2004-07-09 2005-03-02 中国科学院计算技术研究所 三维表面测量方法及装置
CN102494609B (zh) * 2011-11-18 2013-09-18 李志扬 一种基于激光探针阵列的三维摄影方法及装置
CN104296681B (zh) * 2014-10-16 2016-12-07 浙江大学 基于激光点阵标识的三维地形传感方法
CN104634276B (zh) * 2015-02-12 2018-08-07 上海图漾信息科技有限公司 三维测量系统、拍摄设备和方法、深度计算方法和设备
CN105300316B (zh) * 2015-09-22 2017-10-13 大连理工大学 基于灰度重心法的光条中心快速提取方法

Also Published As

Publication number Publication date
CN108592791A (zh) 2018-09-28

Similar Documents

Publication Publication Date Title
CN112669318B (zh) 表面缺陷检测方法、装置、设备及存储介质
CN108876744B (zh) 一种基于区域分割的大尺度点云噪声去噪方法
WO2019150872A1 (ja) 画像処理装置及び画像処理方法
CN109829883B (zh) 产品质量检测方法及装置
EP2430588A1 (en) Object recognition method, object recognition apparatus, and autonomous mobile robot
JP5086970B2 (ja) 木材の外観検査装置、木材の外観検査方法
KR102027986B1 (ko) 비전 카메라를 이용한 비드 인식 장치 및 그 방법
JPH06137828A (ja) 障害物位置検出方法
CN111352123A (zh) 用于车辆检查的机器人和方法、确定车辆底盘中线的方法
CN110910382A (zh) 一种集装箱检测系统
CN108592791B (zh) 一种凹坑检查方法
JP6555211B2 (ja) 二次元画像のエッジ抽出方法
JP6566903B2 (ja) 表面欠陥検出方法および表面欠陥検出装置
CN112565615A (zh) 飞拍拍摄触发点的确定方法和装置
JP2004219154A (ja) 物体の表面形状計測方法及び自動溶接装置
JP2015129751A (ja) 検査方法及びその装置
JP5403367B2 (ja) 物体形状評価装置
JP2010009236A (ja) 平面領域推定装置及びプログラム
CN114663882B (zh) 基于深度学习的电动汽车底盘划痕三维检测方法
WO2019138930A1 (ja) 錠剤検査方法および錠剤検査装置
CN116309882A (zh) 一种面向无人叉车应用的托盘检测和定位方法及系统
Sansoni et al. Design and development of a 3D system for the measurement of tube eccentricity
CN115248187A (zh) 一种偏振分束同时成像的线结构光视觉传感器标定方法
CN115239801A (zh) 一种对象定位方法及装置
CN111950351B (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
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Pit inspection method

Effective date of registration: 20211122

Granted publication date: 20200616

Pledgee: Shandong Nanshan Science and Technology Research Institute Co.,Ltd.

Pledgor: Yantai Nanshan University

Registration number: Y2021980012898

PE01 Entry into force of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20230612

Granted publication date: 20200616

Pledgee: Shandong Nanshan Science and Technology Research Institute Co.,Ltd.

Pledgor: Yantai Nanshan University

Registration number: Y2021980012898

PC01 Cancellation of the registration of the contract for pledge of patent right