CN108917641B - 基于激光器波数合成的样件内部轮廓检测方法及系统 - Google Patents

基于激光器波数合成的样件内部轮廓检测方法及系统 Download PDF

Info

Publication number
CN108917641B
CN108917641B CN201810461489.0A CN201810461489A CN108917641B CN 108917641 B CN108917641 B CN 108917641B CN 201810461489 A CN201810461489 A CN 201810461489A CN 108917641 B CN108917641 B CN 108917641B
Authority
CN
China
Prior art keywords
wave number
frequency
sequence
determining
tunable single
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
CN201810461489.0A
Other languages
English (en)
Other versions
CN108917641A (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.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
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 Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201810461489.0A priority Critical patent/CN108917641B/zh
Publication of CN108917641A publication Critical patent/CN108917641A/zh
Application granted granted Critical
Publication of CN108917641B publication Critical patent/CN108917641B/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/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/2441Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures using interferometry

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种基于激光器波数合成的样件内部轮廓检测方法及系统。所述检测方法包括:获取多个可调谐单频激光器的波数;对多个所述波数在时间轴上依次进行线性调频,确定波数序列;根据所述波数序列确定时间序列干涉条纹图像信号;对所述时间序列干涉条纹图像信号沿时间轴进行希尔伯特变换,构造解析信号;利用傅里叶变换对所述解析信号进行相位补偿,确定补偿后的解析信号;根据所述补偿后的解析信号确定所述被测样件内部的各个界面层的轮廓。采用本发明所提供的检测方法及系统能够提高轴向分辨率以及计算效率。

Description

基于激光器波数合成的样件内部轮廓检测方法及系统
技术领域
本发明涉及激光波数调频干涉领域,特别是涉及一种基于激光器波数合成的样件内部轮廓检测方法及系统。
背景技术
样件的三维轮廓特征是样件重要特征之一。通过高精度检测样件的三维轮廓分布,可以实现对样件缺陷、模具三维特征检测等相关操作。随着工业的发展,现今对轮廓测量提出了以下两点要求:1)高精度,即测量精度应达到微纳米级别;2)透视性,即不仅可对样件表面进行测量,更能对样件内部结构进行测量。近几年来,人们提出了激光波数扫频干涉检测技术,它可以将被测样件内部各个界面层的轮廓结构编码为干涉频谱,通过求解干涉频谱的相位,从而重构出轮廓。由于该技术是以光的波长作为计量单位,同时建立了被测样件内部各个界面层轮廓与干涉频谱相位之间的一一映射关系,它可同时满足上述两点测量要求,具有广阔的应用前景。
目前,制约着激光波数扫频干涉检测技术应用的瓶颈问题是:轴向分辨率不高。该问题的产生主要是激光器波数扫频范围有限。当被测材料内部两界面层在深度方向上的位置较为接近时,干涉频谱将会发生严重的频谱混叠现象,导致无法分离该两界面层的干涉频谱。如何能在激光扫频范围有限条件下大幅度提升轴向分辨率,是目前激光波数扫频干涉检测技术迫切需要解决的技术问题之一。
通常情况下,采用盲分离信息处理技术解决轴向分辨率过低的问题。但是盲分离技术中的数学最优化是基于非线性最小二乘法计算得到的,因此,需要较为准确的数值迭代初始值,算法才能收敛;而在盲分离技术方案应用中,需要事先获取准确的被测材料结构先验知识,这在实际中往往较难实现,因此,现有的盲分离信息处理技术不能很好的解决轴向分辨率过低的问题;此外,这种数值迭代的方式计算效率极低。
发明内容
本发明的目的是提供一种基于激光器波数合成的样件内部轮廓检测方法及系统,以解决现有技术中样件内部轮廓轴向分辨率过低以及计算效率低的问题。
为实现上述目的,本发明提供了如下方案:
一种基于激光器波数合成的样件内部轮廓检测方法,包括:
获取多个可调谐单频激光器的波数;多个所述可调谐单频激光器扫描被测样件内部;
对多个所述波数在时间轴上依次进行线性调频,确定波数序列;所述波数序列分别照射到反射镜和所述被测样件内部的各个界面层;
根据所述波数序列确定时间序列干涉条纹图像信号;
对所述时间序列干涉条纹图像信号沿时间轴进行希尔伯特变换,构造解析信号;
利用傅里叶变换对所述解析信号进行相位补偿,确定补偿后的解析信号;
根据所述补偿后的解析信号确定所述被测样件内部的各个界面层的轮廓。
可选的,所述对多个所述波数在时间轴上依次进行线性调频,确定波数序列,具体包括:
根据公式
Figure GDA0002534941710000021
对多个所述波数在时间轴上依次进行线性调频,确定波数序列;
其中,t为时间轴;kl为可调谐单频激光器l的波数,l=1,2,…,L;k0为可调谐单频激光器1的起始波数;Δkl为可调谐单频激光器l的波数调频范围;Dkl为可调谐单频激光器l与可调谐单频激光器l-1之间的波数间隔;区间[t2l-1,t2l]为可调谐单频激光器l调频的时间范围。
可选的,所述根据所述波数序列确定时间序列干涉条纹图像信号,具体包括:
将所述波数序列分别照射到反射镜和所述被测样件内部的各个界面层上;
获取所述反射镜反射的第一光强以及所述各个界面层反射的第二光强;
根据所述第一光强以及所述第二光强确定时间序列干涉条纹图像信号。
可选的,所述对所述时间序列干涉条纹图像信号沿时间轴进行希尔伯特变换,构造解析信号之前,还包括:
对所述时间序列干涉条纹图像信号进行滤波处理,滤除所述时间序列干涉条纹图像信号的直流分量以及低频分量。
可选的,所述对所述时间序列干涉条纹图像信号沿时间轴进行希尔伯特变换,构造解析信号,具体包括:
根据公式
Figure GDA0002534941710000031
对所述时间序列干涉条纹图像信号沿时间轴进行希尔伯特变换,构造解析信号;
其中,(x,y)为空间坐标点;I为反射光强;p为所述被测样件内部的第p个界面层;q为所述被测样件内部的第q个界面层;R为反射镜;ΛRp(x,y)为光程差;
Figure GDA0002534941710000032
,fRp(x,y)=ΛRp(x,y)/π,fRp(x,y)为干涉频率;
Figure GDA0002534941710000033
Figure GDA0002534941710000034
为干涉相位。
可选的,所述利用傅里叶变换对所述解析信号进行相位补偿,确定补偿后的解析信号,具体包括:
根据公式
Figure GDA0002534941710000041
利用傅里叶变换对所述解析信号进行相位补偿,确定补偿后的解析信号;
其中,j为虚数单位;ξl对应第l个可调谐单频激光器与第l-1个可调谐单频激光器之间的相位补偿值,
Figure GDA0002534941710000042
可选的,所述根据所述补偿后的解析信号确定所述被测样件内部的各个界面层的轮廓,具体包括:
根据公式
Figure GDA0002534941710000043
确定所述被测样件内部的各个界面层的轮廓;
其中,dRp(x,y)为被测件内部各个界面层的轮廓信息;p=1,2,…,M,符号“Im”表示对复数求虚部,“Re”表示对复数求实部;fRp为各个干涉峰值对应的傅里叶变换频率;ρRp为折射率。
一种基于激光器波数合成的样件内部轮廓检测系统,包括:
波数获取模块,用于获取多个可调谐单频激光器的波数;多个所述可调谐单频激光器扫描被测样件内部;
波数序列确定模块,用于对多个所述波数在时间轴上依次进行线性调频,确定波数序列;所述波数序列分别照射到反射镜和所述被测样件内部的各个界面层;
时间序列干涉条纹图像信号确定模块,用于根据所述波数序列确定时间序列干涉条纹图像信号;
解析信号构造模块,用于对所述时间序列干涉条纹图像信号沿时间轴进行希尔伯特变换,构造解析信号;
相位补偿模块,用于利用傅里叶变换对所述解析信号进行相位补偿,确定补偿后的解析信号;
轮廓确定模块,用于根据所述补偿后的解析信号确定所述被测样件内部的各个界面层的轮廓。
可选的,所述波数序列确定模块具体包括:
波数序列确定单元,用于根据公式
Figure GDA0002534941710000051
对多个所述波数在时间轴上依次进行线性调频,确定波数序列;
其中,t为时间轴;kl为可调谐单频激光器l的波数,l=1,2,…,L;k0为可调谐单频激光器1的起始波数;Δkl为可调谐单频激光器l的波数调频范围;Dkl为可调谐单频激光器l与可调谐单频激光器l-1之间的波数间隔;区间[t2l-1,t2l]为可调谐单频激光器l调频的时间范围。
可选的,所述时间序列干涉条纹图像信号确定模块具体包括:
照射单元,用于将所述波数序列分别照射到反射镜和所述被测样件内部的各个界面层上;
光强获取模块,用于获取所述反射镜反射的第一光强以及所述各个界面层反射的第二光强;
时间序列干涉条纹图像信号确定单元,用于根据所述第一光强以及所述第二光强确定时间序列干涉条纹图像信号。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明提供了一种基于激光器波数合成的样件内部轮廓检测方法及系统,将多个可调谐单频激光器的波数序列合成一段较宽的波数输出,大幅度扩展激光的波数调频范围,从而提高了样件内部各个界面层的轮廓轴向分辨率的精度。
在此情况下,使用傅里叶变换计算的干涉频谱分辨率更高。本发明基于傅里叶变换对解析信号进行相位补偿,因此,不存在盲分离信息处理技术中的初始值问题,从而能够提高样件内部各个界面层的轮廓轴向分辨率的精度。
另外,傅里叶变换在计算机上可通过快速傅里叶变换实现,计算效率相比盲分离信息处理技术中的数值迭代优化得到显著提升,具备动态在线测量的应用前景,提高了计算效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明所提供的基于激光器波数合成的样件内部轮廓检测方法流程图;
图2为本发明所提供的可协调单频激光器1和可协调单频激光器2的波数序列图;
图3为本发明所提供的干涉条纹示意图;
图4为本发明所提供的干涉峰值图;
图5为本发明所提供的被测件内部一个界面层的轮廓场示意图;
图6为本发明所提供的被测件内部另一个界面层的轮廓场示意图;
图7为本发明所提供的基于激光器波数合成的样件内部轮廓检测系统结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种基于激光器波数合成的样件内部轮廓检测方法及系统,能够提高轮廓轴向分辨率的精度及计算精度。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明所提供的基于激光器波数合成的样件内部轮廓检测方法流程图,如图1所示,一种基于激光器波数合成的样件内部轮廓检测方法,包括:
步骤101:获取多个可调谐单频激光器的波数;多个所述可调谐单频激光器扫描被测样件内部。
步骤102:对多个所述波数在时间轴上依次进行线性调频,确定波数序列;所述波数序列分别照射到反射镜和所述被测样件内部的各个界面层。
设可调谐单频激光器1,可调谐单频激光器2,可调谐单频激光器3,…,可调谐单频激光器l的波数分别为k1,k2,…,kL;它们在时间轴t上被依次线性调频如下:
Figure GDA0002534941710000071
式中,l=1,2,…,L;k0为可调谐单频激光器1的起始波数;Δkl为可调谐单频激光器l的波数调频范围;Dkl为可调谐单频激光器l与可调谐单频激光器l-1之间的波数间隔;区间[t2l-1,t2l]为可调谐单频激光器l调频的时间范围。
步骤103:根据所述波数序列确定时间序列干涉条纹图像信号。
波数序列kl(t)分别照射到反射镜和被测材料内部的各个界面层上,反射镜反射的光强和被测材料内部各个界面层反射的光强相互叠加形成如图3所示的时间序列干涉条纹图像信号:
Figure GDA0002534941710000081
式中,第一项为时间序列干涉条纹图像的直流分量,第二项为时间序列干涉图像的低频分量,第三项为时间序列干涉图像的高频分量;由于第三项包含着被测材料内部各个界面层的轮廓场信息,其为待求解的信号;(x,y)为空间坐标点;I为反射光强,下标p,q,R分别为被测材料内部第p,q个界面层以及反射镜;ΛRp(x,y)为光程差,即为待求解的轮廓场信息。
步骤104:对所述时间序列干涉条纹图像信号沿时间轴进行希尔伯特变换,构造解析信号。
滤除直流分量和低频分量后,对时间序列干涉条纹图像信号沿时间轴作希尔伯特变换,并构造解析信号如下:
Figure GDA0002534941710000082
Figure GDA0002534941710000083
,式中,fRp(x,y)=ΛRp(x,y)/π为干涉频率;
Figure GDA0002534941710000084
为干涉相位。
步骤105:利用傅里叶变换对所述解析信号进行相位补偿,确定补偿后的解析信号。
当t∈[t2l-1,t2l],l=2,3,…,L时,解析信号s(x,y,t)的初始相位会因可调谐单频激光器之间的波数间隔产生相位偏置值:
Figure GDA0002534941710000091
导致干涉频谱能量不集中,无法完成数据解调;为消除该相位偏置值从而实现波数序列k1,k2,…,kL的合成,在傅里叶变换域中对解析信号作相位补偿如下:
Figure GDA0002534941710000092
Figure GDA0002534941710000093
上式中,j为虚数单位;ξl对应第l个可调谐单频激光器与第l-1个可调谐单频激光器之间的相位补偿值。
步骤106:根据所述补偿后的解析信号确定所述被测样件内部的各个界面层的轮廓。
被测材料内部各个界面层的轮廓信息被编码为频域上的各个干涉峰值,对干涉峰值处的频谱信息解码得到被测件内部各个界面层的轮廓信息:
Figure GDA0002534941710000094
上式中,dRp(x,y)为被测件内部各个界面层的轮廓信息;符号“Im”表示对复数求虚部,“Re”表示对复数求实部;fRp为各个干涉峰值对应的傅里叶变换频率;ρRp为折射率。
下面结合两个可调谐单频激光器构成的具体实施例说明本发明的技术方案。
1)如图2所示,可调谐单频激光器1、2的波数k1、k2在时间轴上被依次线性调频如下:
Figure GDA0002534941710000101
本实施例中,可调谐单频激光器1的起始波数为k0=7.30×106m-1,波数调频范围为Δk1=2.57×103m-1;可调谐单频激光器1、2之间的波数间隔为Dk1=9.88×103m-1;可调谐单频激光器2的波数调频范围为Δk2=4.54×103m-1;区间[t1,t2]离散化为1,2,…,456;区间[t3,t4]离散化为457,458,…,1257;
2)波数序列kl(t)分别照射到反射镜和被测材料内部的各个界面层上(本实施例中,被测件内部的界面层为2个,M=2;两个界面层相距148.19μm),反射镜反射的光强和被测材料内部各个界面层反射的光强相互叠加形成如图2所示的时间序列干涉条纹图像信号:
Figure GDA0002534941710000102
式中,第一项为时间序列干涉条纹图像的直流分量,第二项为时间序列干涉图像的低频分量,第三项为时间序列干涉图像的高频分量;由于第三项包含着被测材料内部各个界面层的轮廓场信息,其为待求解的信号;(x,y)为空间坐标点;I为反射光强,下标p,q,R分别为被测材料内部第p,q个界面层以及反射镜;ΛRp(x,y)为光程差,即为待求解的轮廓场信息;
3)滤除直流分量和低频分量后,对时间序列干涉条纹图像信号沿时间轴作希尔伯特变换,并构造解析信号如下:
Figure GDA0002534941710000111
Figure GDA0002534941710000112
,式中,fRp(x,y)=ΛRp(x,y)/π为干涉频率;
Figure GDA0002534941710000113
为干涉相位;
4)当t∈[t3,t4],l=2,3,…,L时,解析信号s(x,y,t)的初始相位会因可调谐单频激光器之间的波数间隔产生相位偏置值:
Figure GDA0002534941710000114
导致干涉频谱能量不集中,无法完成数据解调;为消除该相位偏置值从而实现波数序列k1,k2的合成,在傅里叶变换域中对解析信号作相位补偿如下:
Figure GDA0002534941710000115
ξ1=-2π·f·(Δk1+Dk1);
上式中,j为虚数单位;ξ1对应第2个可调谐单频激光器与第1个可调谐单频激光器之间的相位补偿值;
5)如图4所示,经步骤3)-4)后,被测材料内部2个界面层的轮廓信息被编码为频域上的2个干涉峰值;如图5-图6所示,对干涉峰值处的频谱信息解码得到被测件内部2个界面层的轮廓信息:
Figure GDA0002534941710000121
图5为本发明所提供的被测件内部一个界面层的轮廓场示意图,图6为本发明所提供的被测件内部另一个界面层的轮廓场示意图,如图5-图6所示,本发明针对激光波数扫频干涉检测系统,采用相位补偿的方式,将多个可调谐单频激光器的波数序列合成一段较宽的波数输出,大幅度提升光源的波数扫频范围,使得轴向分辨率提高至微米级别。
图7为本发明所提供的基于激光器波数合成的样件内部轮廓检测系统结构图,如图7所示,一种基于激光器波数合成的样件内部轮廓检测系统,包括:
波数获取模块701,用于获取多个可调谐单频激光器的波数;多个所述可调谐单频激光器扫描被测样件内部。
波数序列确定模块702,用于对多个所述波数在时间轴上依次进行线性调频,确定波数序列;所述波数序列分别照射到反射镜和所述被测样件内部的各个界面层。
所述波数序列确定模块702具体包括:波数序列确定单元,用于根据公式
Figure GDA0002534941710000131
对多个所述波数在时间轴上依次进行线性调频,确定波数序列;其中,t为时间轴;kl为可调谐单频激光器l的波数,l=1,2,…,L;k0为可调谐单频激光器1的起始波数;Δkl为可调谐单频激光器l的波数调频范围;Dkl为可调谐单频激光器l与可调谐单频激光器l-1之间的波数间隔;区间[t2l-1,t2l]为可调谐单频激光器l调频的时间范围。
时间序列干涉条纹图像信号确定模块703,用于根据所述波数序列确定时间序列干涉条纹图像信号。
所述时间序列干涉条纹图像信号确定模块703具体包括:照射单元,用于将所述波数序列分别照射到反射镜和所述被测样件内部的各个界面层上;光强获取模块,用于获取所述反射镜反射的第一光强以及所述各个界面层反射的第二光强;时间序列干涉条纹图像信号确定单元,用于根据所述第一光强以及所述第二光强确定时间序列干涉条纹图像信号。
解析信号构造模块704,用于对所述时间序列干涉条纹图像信号沿时间轴进行希尔伯特变换,构造解析信号。
相位补偿模块705,用于利用傅里叶变换对所述解析信号进行相位补偿,确定补偿后的解析信号。
轮廓确定模块706,用于根据所述补偿后的解析信号确定所述被测样件内部的各个界面层的轮廓。
本发明所提供的样件内部轮廓检测方法及系统能够将多个可调谐单频激光器的波数序列合成一段较宽的波数输出,大幅度扩展激光的波数调频范围,从而提高了轴向分辨率精度;且采用傅里叶变换对所述解析信号进行相位补偿,不需要考虑数值迭代初始值的精度,从而提高了计算效率。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (9)

1.一种基于激光器波数合成的样件内部轮廓检测方法,其特征在于,包括:
获取多个可调谐单频激光器的波数;多个所述可调谐单频激光器扫描被测样件内部;
对多个所述波数在时间轴上依次进行线性调频,确定波数序列;所述波数序列分别照射到反射镜和所述被测样件内部的各个界面层;将多个可调谐单频激光器的波数序列合成一段较宽的波数输出,扩展激光的波数调频范围;
根据所述波数序列确定时间序列干涉条纹图像信号;
对所述时间序列干涉条纹图像信号沿时间轴进行希尔伯特变换,构造解析信号;
利用傅里叶变换对所述解析信号进行相位补偿,确定补偿后的解析信号;
根据所述补偿后的解析信号确定所述被测样件内部的各个界面层的轮廓;被测材料内部各个界面层的轮廓信息被编码为频域上的各个干涉峰值,对干涉峰值处的频谱信息解码得到被测件内部各个界面层的轮廓信息:
Figure FDA0002534941700000011
其中,dRp(x,y)为被测件内部各个界面层的轮廓信息;p=1,2,…,M,符号“Im”表示对复数求虚部,“Re”表示对复数求实部;fRp为各个干涉峰值对应的傅里叶变换频率;ρRp为折射率。
2.根据权利要求1所述的样件内部轮廓检测方法,其特征在于,所述对多个所述波数在时间轴上依次进行线性调频,确定波数序列,具体包括:
根据公式
Figure FDA0002534941700000012
对多个所述波数在时间轴上依次进行线性调频,确定波数序列;
其中,t为时间轴;kl为可调谐单频激光器l的波数,l=1,2,…,L;k0为可调谐单频激光器1的起始波数;Δkl为可调谐单频激光器l的波数调频范围;Dkl为可调谐单频激光器l与可调谐单频激光器l-1之间的波数间隔;区间[t2l-1,t2l]为可调谐单频激光器l调频的时间范围。
3.根据权利要求2所述的样件内部轮廓检测方法,其特征在于,所述根据所述波数序列确定时间序列干涉条纹图像信号,具体包括:
将所述波数序列分别照射到反射镜和所述被测样件内部的各个界面层上;
获取所述反射镜反射的第一光强以及所述各个界面层反射的第二光强;
根据所述第一光强以及所述第二光强确定时间序列干涉条纹图像信号。
4.根据权利要求3所述的样件内部轮廓检测方法,其特征在于,所述对所述时间序列干涉条纹图像信号沿时间轴进行希尔伯特变换,构造解析信号之前,还包括:
对所述时间序列干涉条纹图像信号进行滤波处理,滤除所述时间序列干涉条纹图像信号的直流分量以及低频分量。
5.根据权利要求4所述的样件内部轮廓检测方法,其特征在于,所述对所述时间序列干涉条纹图像信号沿时间轴进行希尔伯特变换,构造解析信号,具体包括:
根据公式
Figure FDA0002534941700000021
对所述时间序列干涉条纹图像信号沿时间轴进行希尔伯特变换,构造解析信号;
其中,(x,y)为空间坐标点;I为反射光强;p为所述被测样件内部的第p个界面层;q为所述被测样件内部的第q个界面层;R为反射镜;ΛRp(x,y)为光程差;
Figure FDA0002534941700000031
Figure FDA0002534941700000032
fRp(x,y)=ΛRp(x,y)/π,fRp(x,y)为干涉频率;
Figure FDA0002534941700000033
Figure FDA0002534941700000034
为干涉相位。
6.根据权利要求5所述的样件内部轮廓检测方法,其特征在于,所述利用傅里叶变换对所述解析信号进行相位补偿,确定补偿后的解析信号,具体包括:
根据公式
Figure FDA0002534941700000035
利用傅里叶变换对所述解析信号进行相位补偿,确定补偿后的解析信号;
其中,j为虚数单位;ξl对应第l个可调谐单频激光器与第l-1个可调谐单频激光器之间的相位补偿值,
Figure FDA0002534941700000036
7.一种基于激光器波数合成的样件内部轮廓检测系统,其特征在于,包括:
波数获取模块,用于获取多个可调谐单频激光器的波数;多个所述可调谐单频激光器扫描被测样件内部;
波数序列确定模块,用于对多个所述波数在时间轴上依次进行线性调频,确定波数序列;所述波数序列分别照射到反射镜和所述被测样件内部的各个界面层;将多个可调谐单频激光器的波数序列合成一段较宽的波数输出,扩展激光的波数调频范围;
时间序列干涉条纹图像信号确定模块,用于根据所述波数序列确定时间序列干涉条纹图像信号;
解析信号构造模块,用于对所述时间序列干涉条纹图像信号沿时间轴进行希尔伯特变换,构造解析信号;
相位补偿模块,用于利用傅里叶变换对所述解析信号进行相位补偿,确定补偿后的解析信号;
轮廓确定模块,用于根据补偿后的所述解析信号确定所述被测样件内部的各个界面层的轮廓;被测材料内部各个界面层的轮廓信息被编码为频域上的各个干涉峰值,对干涉峰值处的频谱信息解码得到被测件内部各个界面层的轮廓信息:
Figure FDA0002534941700000041
dRp(x,y)为被测件内部各个界面层的轮廓信息;其中,p=1,2,…,M,符号“Im”表示对复数求虚部,“Re”表示对复数求实部;fRp为各个干涉峰值对应的傅里叶变换频率;ρRp为折射率。
8.根据权利要求7所述的样件内部轮廓检测系统,其特征在于,所述波数序列确定模块具体包括:
波数序列确定单元,用于根据公式
Figure FDA0002534941700000042
对多个所述波数在时间轴上依次进行线性调频,确定波数序列;
其中,t为时间轴;kl为可调谐单频激光器l的波数,l=1,2,…,L;k0为可调谐单频激光器1的起始波数;Δkl为可调谐单频激光器l的波数调频范围;Dkl为可调谐单频激光器l与可调谐单频激光器l-1之间的波数间隔;区间[t2l-1,t2l]为可调谐单频激光器l调频的时间范围。
9.根据权利要求8所述的样件内部轮廓检测系统,其特征在于,所述时间序列干涉条纹图像信号确定模块具体包括:
照射单元,用于将所述波数序列分别照射到反射镜和所述被测样件内部的各个界面层上;
光强获取模块,用于获取所述反射镜反射的第一光强以及所述各个界面层反射的第二光强;
时间序列干涉条纹图像信号确定单元,用于根据所述第一光强以及所述第二光强确定时间序列干涉条纹图像信号。
CN201810461489.0A 2018-05-15 2018-05-15 基于激光器波数合成的样件内部轮廓检测方法及系统 Active CN108917641B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810461489.0A CN108917641B (zh) 2018-05-15 2018-05-15 基于激光器波数合成的样件内部轮廓检测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810461489.0A CN108917641B (zh) 2018-05-15 2018-05-15 基于激光器波数合成的样件内部轮廓检测方法及系统

Publications (2)

Publication Number Publication Date
CN108917641A CN108917641A (zh) 2018-11-30
CN108917641B true CN108917641B (zh) 2020-08-04

Family

ID=64403910

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810461489.0A Active CN108917641B (zh) 2018-05-15 2018-05-15 基于激光器波数合成的样件内部轮廓检测方法及系统

Country Status (1)

Country Link
CN (1) CN108917641B (zh)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10196214B4 (de) * 2000-05-19 2010-01-28 Zygo Corp., Middlefield Höhenabtast-Interferometer zum Bestimmen der Absolutposition und des Oberflächenprofils eines Objekts im Hinblick auf ein Datum
JPWO2003036229A1 (ja) * 2001-10-25 2005-02-17 東レエンジニアリング株式会社 表面形状測定方法およびその装置
CN104034277B (zh) * 2014-04-08 2016-08-17 华南师范大学 一种双波长微纳结构相位测量方法
CN105890538A (zh) * 2014-12-30 2016-08-24 广东工业大学 三表面干涉式高精度曲面轮廓测量系统及方法
CN105928877B (zh) * 2016-04-27 2017-03-15 广东工业大学 强脉冲噪声环境下基于方差比盲分离的复杂结构材料高精度透视检测方法
CN106500589B (zh) * 2016-12-05 2020-09-25 苏州大学 一种多波长可调谐显微干涉的测量方法及其装置

Also Published As

Publication number Publication date
CN108917641A (zh) 2018-11-30

Similar Documents

Publication Publication Date Title
CN109163672B (zh) 一种基于白光干涉零光程差位置拾取算法的微观形貌测量方法
CN1784588B (zh) 使用扫描干涉测量形成复杂表面结构的轮廓以及对其表征
US7388676B2 (en) Image processing apparatus and refractive index distribution measuring apparatus
CN1784587B (zh) 使用扫描干涉测量形成复杂表面结构的轮廓
KR20120113653A (ko) 간섭을 이용한 막 두께 계측 장치 및 간섭을 이용한 막 두께 계측 방법
EP2677271B1 (en) Broadband interferometer for determining a property of a thin film
CN113418469B (zh) 光谱共焦扫描共光路数字全息测量系统及测量方法
Kim et al. Fast and reliable measurement of thin film thickness profile based on wavelet transform in spectrally resolved white-light interferometry
CN107490464B (zh) 基于波形叠加原理的非线性波浪的反射波分离方法
JP2003344024A (ja) 膜厚測定方法及び膜厚測定装置
CN108917641B (zh) 基于激光器波数合成的样件内部轮廓检测方法及系统
CN117128877B (zh) 一种薄膜厚度检测方法、计算机及系统
CN117109450B (zh) 一种自适应参数谱估计的高分辨光谱共焦干涉薄膜测厚法
CN103196387B (zh) 柱面面型检测系统及方法
CN105928877A (zh) 基于方差比盲分离的复杂结构材料高精度透视检测方法
CN111157115B (zh) 一种水下布里渊散射光谱获取方法及装置
Kotlikov et al. Comparative analysis of the stability criteria of interference coatings
CN116086338A (zh) 一种数字全息变形检测方法
CN105115940A (zh) 光学材料折射率曲线测量方法及装置
CN113418470B (zh) 光谱扫描共焦单次曝光数字全息测量系统及测量方法
CN114998469A (zh) 一种基于神经卷积网络的区域红外数字全息方法
Zhu et al. Photoacoustic speckle pattern interferometry for detecting cracks of different sizes
Gao et al. Characterization of the topography fidelity of 3D optical microscopy
CN118190878B (zh) 光纤型光谱共焦相干层析光学系统及应用
Wei et al. Determining a bandpass filter for interfering fringe spectrum extraction using clustering analysis

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