CN1401990A - 序列散斑场强度扫描位移测量方法 - Google Patents
序列散斑场强度扫描位移测量方法 Download PDFInfo
- Publication number
- CN1401990A CN1401990A CN 02130912 CN02130912A CN1401990A CN 1401990 A CN1401990 A CN 1401990A CN 02130912 CN02130912 CN 02130912 CN 02130912 A CN02130912 A CN 02130912A CN 1401990 A CN1401990 A CN 1401990A
- Authority
- CN
- China
- Prior art keywords
- interval
- intensity
- sequence
- phase
- speckle
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 26
- 230000011218 segmentation Effects 0.000 claims description 4
- 230000007423 decrease Effects 0.000 claims description 3
- 230000000630 rising effect Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 description 8
- 238000005259 measurement Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005305 interferometry Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 239000004744 fabric Substances 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 229920003229 poly(methyl methacrylate) Polymers 0.000 description 2
- 239000004926 polymethyl methacrylate Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241000931526 Acer campestre Species 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005538 encapsulation Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 239000003973 paint Substances 0.000 description 1
- 230000005622 photoelectricity Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000007921 spray Substances 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
一种序列散斑场强度扫描位移测量方法,属光电无损检测领域。本方法由下列步骤完成:对被测物通过双光束干涉获得单调连续变化的序列散斑图;对所得的序列散斑图在时间域内用相关法划分出不同的强度统计和扫描区间;对划分好区间的序列散斑图进行强度统计平均和扫描获得该区间相位函数的主值;对获得主值相位函数区间内的序列散斑图所对应的余弦强度序列再进行六区域分割并对其强度值逐点比较,获得余弦序列所对应的相位变化周期;将主值相位和相位变化周期迭加,获得由相关划分所确定的区间中的连续位相分布;将获得的所有区间的连续相位值叠加即可完成位移场检测。采用本方法可以在不加外载波和相移***的条件下,对时变位移场进行定量测量。
Description
技术领域
本发明涉及序列散斑场强度扫描位移测量方法,属光电无损检测领域。
背景技术
散斑干涉或电子散斑干涉(ESPI)是目前被广泛应用于漫射体表面位移或变形检测的现代光学计量方法。它们基于双光束散斑干涉而形成双曝光散斑图,然后通过逐点或全场分析(对散斑干涉)或不同运动与变形状态下的数字化散斑干涉图,通过数字图像处理进行分析(对电子散斑干涉)。在以往的计量中它们主要是以强度相关条纹的形式来表征所需要测量的位移或变形量的。但是由于条纹本身的复杂性与随机性,很难形成统一的条纹处理模式和定量解释方式。因此,对条纹判读主要是依靠人工进行,如确定条纹中心、为条纹级数赋值等。这种传统的图象处理方法速度慢而且精度低,条纹的稀疏性和条纹的规则性都会影响判断的准确性。随着计算机图像处理技术的发展,用计算机可以完成图像的平滑,二值化,细化,条纹跟踪,条纹级数赋值等,使条纹图的分析由人工向自动化迈进了一步,但由于受条纹场的复杂性及噪声水平的影响,要获得高精度的测量结果需要付出相当的工作量。因此,以相关条纹场为计量表征模式的散斑干涉或ESPI主要被应用于定性分析。
为了解决散斑计量中的定量检测问题,在ESPI技术中引入了相移方法,即相移电子散斑干涉术(PSESPI)。该方法能够很容易地对变形场进行定量检测。但其也存在着明显的不足,即它需要同一变形状态下的至少三幅干涉图进行相位求解,同时相移装置的引入增加了实验装置的复杂性和对环境噪声的敏感性。因此,相移方法很难适用于时变或较恶劣环境下的动态位移或变形检测。另一种方法是在散斑干涉和ESPI检测中引入载波条纹(通过在检测***中引入载波装置或在检测中旋转参考光的反射镜引入空间条纹而形成载波),也称作载波电子散斑干涉术(CESPI)。它利用傅立叶变换技术,在只需要一幅干涉图的情况下通过应用傅立叶空间相位解调方法,就可以解调出全场相位。所以CESPI可以方便地应用于无法用相移技术的场合,例如时变场或瞬态变化场等。
然而,不论是从传统的散斑法、ESPI到PSESPI和CESPI,位移和变形场的测量都是基于两个状态,即变形前后的两个状态,没有考虑到时间参数。但时变场是普遍存在的,如热变形场,由材料或结构的疲劳破坏引起的蠕变等。这些场通常是基于时间变化并且变形量是比较大的。因此非常需要能高精度定量测量时变场和大变形的方法。九十年代中后期,由于高速固体摄像器件(CCD)的出现,使短时间隔的数字序列图像的采集成为可能。因而以时间参量为载波的序列场分析被提到日程上来,并进而对高速变形场的数字化检测产生了推动作用。德国斯图加特大学的Charles Joenathan等人提出了一种基于傅立叶空间滤波原理的时变散斑干涉方法。在这一方法中,他们利用计算机控制的旋转台使物体产生连续的运动,从而引起散斑场强度的时间调制变化。这样在物体变形的过程中采集大量的散斑图,并通过傅立叶变换技术提取出全场的相位信息。因为要使用傅立叶变换解调技术来提取相位信息,所以在Joenathan方法中需要大量的摄取散斑干涉图,并且要求被测物体的变形或***能提供近似线性的时间载波,以便用傅立叶方法分析时能统一进行旁瓣滤波。因此该方法不适应被测物体的任意变形,需要附加线性载波,所需设备复杂操作较困难。
发明内容
本发明的目的是提供一种不需附加线性载波和相移的时变场相位检测方法。
本发明序列散斑场强度扫描位移测量方法由下列步骤完成:1、对被测物体首先通过双光束散斑干涉获得单调连续变化的序列散斑图;2、对所得的序列散斑图在时间域内应用相关法划分出不同的强度统计和扫描区间;3、对划分好区间的序列散斑图进行时间域强度统计平均和强度扫描获得该区间内相位函数的主值;4、对获得主值相位函数的区间内的序列散斑图所对应的余弦强度序列再进行六个区域分割,即第一采样点、上升区、波峰、下降区、波谷和最后采样点;5、对所分割的六区域强度值逐点比较,获得余弦序列所对应的相位变化周期;6、将主值相位和相位变化周期合成,获得由相关划分所确定的该区间中的连续位相分布;7、重复步骤3到步骤6的过程,获得所有区间的连续相位分布,并将所有区间对应的连续相位叠加,即可完成位移场检测。
具体处理方法如下:
1、双光束干涉获得序列散斑场:
I(x,y,ti)=I0(x,y)+Im(x,y)cos(x,y,ti) (1)
其中,I0(x,y)和Im(x,y)分别为干涉场的平均强度和调制函数,(x,y,t)是与物体表面位移或变形有关的时变相位函数。
2、时间域强度相关划分确定统计和扫描区间nk:
第一个相关区间: 当Clj≤C0,n1=|j-l| (2)
··· ···
第k个相关区间 当Ckj≤C0,nk=|j-k| (3)其中,j,k为正整数,Ckj为归一化相关系数C0为相关阈值,可取C0≥0.6·Ckk。
3、区段nk中时间域强度统计平均和强度扫描求主值相位:首先,应用时间域统计平均运算获得平均强度I0(x,y),即
然后,在选定的时间序列区间,进行序列场强度扫描,获得强度极值Imax(x,y),于是可得调制函数: 最后由(1)式可以获得主值相位函数(x,y,t)为:
4、区段nk中余弦强度序列六个区域分割:
对区段nk中余弦强度序列按图1分别分割成第一采样点、上升区、波峰、下降区、波谷和最后采样点等六个区域。
5、6、对所分割的六区域强度值逐点比较,并将主值相位和相位变化周期合成,获得由相关划分所确定的区间中的连续位相分布,具体过程为:
首先,设I(x,y,ti)表示时刻i时所对应的时变散斑场,p为余弦函数周期数,c(x,y,ti)表示相位去包裹后所获得的i时刻的连续相位函数,则(x,y,ti)的时间域去包裹过程为:
若I1(x,y,t1)≥I2(x,y,t2),则ck(x,y,t1)=(x,y,t1)
(7)
若I1(x,y,t1)<I2(x,y,t2),则ck(x,y,t1)=2π-(x,y,t1)
当在第i时刻到第nk-1时刻时:
i=2,p=0
若Ii-1(x,y,ti-1)>Ii(x,y,ti)>Ii+1(x,y,ti+1),则ck(x,y,ti)=2pπ+(x,y,ti) (8)
若Ii-1(x,y,ti-1)<Ii(x,y,ti)<Ii+1(x,y,ti+1),则ck(x,y,ti)=2(p+1)π-(x,y,ti)
当Ii-1(x,y,ti-1)>Ii(x,y,ti)且Ii(x,y,ti)<Ii+1(x,y,ti+1)时:
若Ii-1(x,y,ti-1)>Ii+1(x,y,ti+1),则ck(x,y,ti)=2pπ+(x,y,ti) (9)
若Ii-1(x,y,ti-1)<Ii+1(x,y,ti+1),则ck(x,y,ti)=2(p+1)π-(x,y,ti)
当Ii-1(x,y,ti-1)<Ii(x,y,ti)且Ii(x,y,ti)>Ii+1(x,y,ti+1)时:
则p=p+1
(10)
若Ii-1(x,y,ti-1)<Ii+1(x,y,ti+1),ck(x,y,ti)=2pπ-(x,y,ti)
若Ii-1(x,y,ti-1)>Ii+1(x,y,ti+1),ck(x,y,ti)=2pπ+(x,y,ti)
最后,在nk时刻:若
则
(11)若
则
7、重复步骤3到步骤6获得所有区间的连续相位分布和,并结合检测***参数完成位移检测,即:
通过方程(7)到(11)的运算,即可将主值相位函数(x,y,ti)转换为连续相位函数ck(x,y,ti)。如设ck(x,y,t1)表示初始时刻的相位值,则由位移和变形引起的相位变化为:Δk(x,y,ti)=ck(x,y,ti)-ck(x,y,t1) (12)
于是,所有区间所对应的连续相位为: 其中K为相关划分所产生的总区间个数。进而其位移场为:离面位移:W(x,y,ti)=λ(x,y,ti)/[2π(1+cosθ)] (14)面内位移:U(x,y,ti)=λ(x,y,ti)/(4πsinθ)V(x,y,ti)=λ(x,y,ti)/(4πsinθ) (15)
其中,θ为光轴和照明物光轴间的夹角,U(x,y,t),V(x,y,t)和W(x,y,t)分别为ti时刻物体沿x,y及z轴的位移分量,λ为所用相干光波长。
对于上文的主值相位和强度变化周期合成过程(对应步骤4、5、6),也可以用图2所示的流程图表示。这样,通过上文七个步骤,即完成时变场位移测量。本发明序列散斑场强度扫描位移测量方法可以不加载波和相移而对时变位移或变形场进行测量。
附图说明
图1是相关区段nk段中余弦强度序列的六个区域分割图
图2是时间域强度扫描相位去包裹流程图
图3是实施例中要检测的轿车前灯配光镜
图4是实施例中车灯配光镜热变形产生的序列散斑干涉图
图5是对散斑图实现区间划分后截出的一个序列区间图
图6是对划分好区间中的序列散斑图进行主值相位和周期数求解结果图
图7,图8分别是用灰度和三维图表示的离面位移场分布
具体实施方式
实施例:车灯配光镜热变形检测
如图3为实际要检测的轿车前灯配光镜,其中试件长宽限度分别为120mm×50mm,材料为聚甲基丙烯酸甲酯(PMMA)。检测量为车灯配光镜表面受热变形后,产生的离面位移。根据上文提供的测量方法,检测过程为:1)、检测前先在待检测的试件表面均匀地喷一层白漆,以增强物体表面对激光的反射强度;2)、将配光镜放在双光束干涉光路中,并牢固地固定在试件夹持架上;3)、通过干涉***记录车灯配光镜热变形产生的序列散斑干涉图,见图4(对应上文中步骤1);4)、对记录的序列散斑图作相关运算实现区间划分(对应上文中步骤2),图5为划分中截出的一个序列区间;5)、对划分好区间中的序列散斑图进行主值相位和周期数求解(见上文步骤3、4、5、6),结果见图6;7)将所有区间连续相位叠加,并结合检测***的激光波长、***几何参数,应用公式(14)计算车灯配光镜热变形产生的离面位移。图7和8分别为灰度和三维图表示的离面位移。
Claims (1)
1.一种序列散斑场强度扫描位移测量方法,其特征是本方法由下列步骤完成:
(1)对被测物体首先通过双光束散斑干涉获得单调连续变化的序列散斑图;
(2)对所得的序列散斑图在时间域应用相关法划分出不同的强度统计和扫描区间;
(3)对划分好区间的序列散斑图进行时间域强度统计平均和强度扫描获得该区间内相位函数的主值;
(4)对获得主值相位函数的区间内的序列散斑图所对应的余弦强度序列再进行六个区域分割,即第一采样点、上升区、波峰、下降区、波谷和最后采样点;
(5)对所分割的六区域强度值逐点比较,获得余弦序列所对应的相位变化周期;
(6)将主值相位和相位变化周期合成,获得由相关划分所确定的该区间中的连续位相分布;
(7)重复步骤(3)到步骤(6)的过程,获得所有区间的连续相位分布,并将所有区间对应的连续相位叠加,即可完成位移场检测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 02130912 CN1401990A (zh) | 2002-09-13 | 2002-09-13 | 序列散斑场强度扫描位移测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 02130912 CN1401990A (zh) | 2002-09-13 | 2002-09-13 | 序列散斑场强度扫描位移测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN1401990A true CN1401990A (zh) | 2003-03-12 |
Family
ID=4746503
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 02130912 Pending CN1401990A (zh) | 2002-09-13 | 2002-09-13 | 序列散斑场强度扫描位移测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1401990A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100552431C (zh) * | 2004-05-12 | 2009-10-21 | 中国科学院长春光学精密机械与物理研究所 | 激光散斑干涉测量方法及装置 |
CN102967261A (zh) * | 2012-11-10 | 2013-03-13 | 华中科技大学 | 一种基于数字散斑相关法的激光位移测量方法 |
-
2002
- 2002-09-13 CN CN 02130912 patent/CN1401990A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100552431C (zh) * | 2004-05-12 | 2009-10-21 | 中国科学院长春光学精密机械与物理研究所 | 激光散斑干涉测量方法及装置 |
CN102967261A (zh) * | 2012-11-10 | 2013-03-13 | 华中科技大学 | 一种基于数字散斑相关法的激光位移测量方法 |
CN102967261B (zh) * | 2012-11-10 | 2015-04-15 | 华中科技大学 | 一种基于数字散斑相关法的激光位移测量方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5069548A (en) | Field shift moire system | |
Morimoto et al. | Fringe pattern analysis by a phase-shifting method using Fourier transform | |
Judge et al. | Holographic deformation measurements by Fourier transform technique with automatic phase unwrapping | |
CN113091634B (zh) | 一种适用于白光扫描干涉的微观形貌快速测量方法 | |
CN1093935C (zh) | 一种快速投影结构光的三维轮廓相位测量方法及装置 | |
EP0647310B1 (en) | Method of an apparatus for interferometrically inspecting a surface of an object | |
CN1932433A (zh) | 载频电子散斑位移场的分离方法 | |
US5239364A (en) | Light phase difference measuring method using an interferometer | |
Aziz | Interferometric measurement of surface roughness in engine cylinder walls | |
CN101718533B (zh) | 镜像莫尔测量装置及方法 | |
Wang | Development and application of computer-aided fringe analysis | |
CN1401990A (zh) | 序列散斑场强度扫描位移测量方法 | |
Zhou et al. | Eliminating the influence of source spectrum of white light scanning interferometry through time-delay estimation algorithm | |
US6738511B1 (en) | Reduced noise sensitivity method and apparatus for converting an interferogram phase map to a surface profile map | |
Davies et al. | Television holography and its applications | |
CN1401991A (zh) | 时间序列散斑场脉冲计数位移测量方法及其装置 | |
Morimoto et al. | Theory and application of sampling moiré method | |
CN1112567C (zh) | 用二元编码模板产生正弦结构光场的方法 | |
KR19980084166A (ko) | 타이어의 접지압 분포 측정장치 및 그 측정방법 | |
Yatagai | Interferometric testing technology developments and applications in Japan | |
Castellini et al. | System for measuring the coordinates of tire surfaces in transient conditions when rolling over obstacles: description of the system and performance analysis | |
Carlsson et al. | Automatic measurement of surface height with the interferoscope | |
Pedraza-Ortega et al. | Three-Dimensional Reconstruction System based on a Segmentation Algorithms and a Modified Fourier Transform Profilometry | |
Nukala et al. | Novel composite coded pattern for small angle measurement using imaging method | |
Yoshizawa et al. | Interpolation of contour intervals in Moire topography by phase modulation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |