CN107024186A - 波长移相中基于fast‑ica算法的超薄平板干涉信号分离处理方法 - Google Patents
波长移相中基于fast‑ica算法的超薄平板干涉信号分离处理方法 Download PDFInfo
- Publication number
- CN107024186A CN107024186A CN201710180047.4A CN201710180047A CN107024186A CN 107024186 A CN107024186 A CN 107024186A CN 201710180047 A CN201710180047 A CN 201710180047A CN 107024186 A CN107024186 A CN 107024186A
- Authority
- CN
- China
- Prior art keywords
- ultra
- interference signal
- rear surfaces
- interference
- signal
- 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
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/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/2441—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures using interferometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02083—Interferometers characterised by particular signal processing and presentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Signal Processing (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
本发明涉及一种波长移相中基于FAST‑ICA算法的超薄平板干涉信号分离处理方法。本发明在波长移相干涉测量获得的超薄平行平板干涉信号后利用盲信号分离技术,引入FAST‑ICA算法,通过快速迭代算法分离超薄平板前后表面干涉条纹,然后通过时域移相算法处理分离后的干涉条纹,从而得到前后表面的面形信息。该方法主要针对波长移相中超薄平行平板的条纹图处理,为这类超薄平板前后表面信息分离提供了一种方法,可用于实现超薄平板前后表面面形的分离测量。
Description
技术领域
本发明涉及光学检测领域,具体涉及一种波长移相中基于FAST-ICA算法的超薄平板干涉信号分离处理方法。
背景技术
移相干涉技术从实现方式上可以分为以压电陶瓷为代表的硬件移相技术和波长移相技术,波长移相技术又称为准软件移相技术。当测量平行平板时,硬件移相技术不能克服平行平板前后表面干涉形成的寄生条纹的影响,需要进行涂抹凡士林等措施消除寄生条纹,但由此会对测量引入误差。而波长移相技术可以克服上述缺点,利用算法实现前后表面干涉信号分离,从而以较高的精度实现平行平板前后表面面形的测量。
在波长移相干涉仪中,激光器既作为光源,其波长又可以连续改变,起到移相器的作用。当测量平行平板时,根据波长移相原理,前后表面产生的干涉信号虽然空域的干涉图上是混叠在一起的,但在频域上,前后表面产生的干涉信息频率是不同的,因此,原理上前后表面的干涉信息是可分离的,从而可以通过设计合适的算法实现前后表面信息的分离并获得各自的轮廓信息。
波长移相测量技术的优点总结如下:
(1)波长移相技术不需要推动大口径的参考镜,可以简化参考镜的机械结构。
(2)波长移相干涉仪在测量中,系统的机械部分保持不动,消除了由于硬件移动而引起的误差,进一步提高了测量精度。
(3)可以通过软件算法消除寄生条纹的影响。
在利用波长移相干涉仪测量平行平板时,如果平板太薄,受限于可调谐半导体激光器的波长调谐分辨率,前表面干涉光强信号与后表面干涉光强频率信息在频域内是会产生叠加的,利用传统的波长移相算法无法进行有效分离,或分离误差较大,从而无法实现这些超薄平板的测量。
发明内容
本发明的目的在于针对现有波长移相算法无法分离超薄平行平板前后表面干涉条纹信息的现状,提供一种波长移相中基于FAST-ICA算法的超薄平板干涉信号分离处理方法,在波长移相干涉测量获得的超薄平行平板干涉信号后引入FAST-ICA算法,通过快速迭代算法分离超薄平板前后表面干涉信号,然后通过时域移相算法处理分离后的干涉信号,从而得到前后表面的面形信息。实现波长移相对超薄平行平板的测量。
为达到上述发明目的,本发明采用下述技术方案:
一种波长移相中基于FAST-ICA算法的超薄平板干涉信号分离处理方法,由以下步骤实现:
步骤一:波长移相干涉仪通过波长调谐测量超薄平行平板,得到包含前后表面面形信息的一系列波长移相干涉图;
步骤二:根据像素列出时域序列干涉信号,对序列干涉信号去均值化:首先求出干涉信号的均值,然后将干涉信号减去均值,得到去均值化后的干涉信号;
步骤三:利用白化矩阵将均值化后的干涉信号进行白化;
步骤四:根据独立分量提取算法,利用分离矩阵处理白化后的去均值信号;
步骤五:判断是否收敛:如果收敛,则分离矩阵已最优,即得到待分离的前后表面干涉信号,执行第六步;如果不收敛,则分离矩阵不是最优,需要重新优化分离矩阵,再返回第四步;
步骤六:将分离后的前后表面干涉信号利用时域移相算法,计算求解得到前后表面面形。
与现有技术相比,本发明具有如下的优点:
本发明针对传统波长移相干涉测量方法的局限性,在波长移相干涉测量获得的超薄平行平板干涉信号后采用盲信号分离技术。盲分离技术是指在单传感器采集条件下,同时采集到目标信号和一个或者多个干扰信号,在缺乏目标信号和干扰信号的先验信息的情况下,利用接收到的混合信号,采用盲信号分离方法,分离出目标信号和干扰信号,或者提取出感兴趣的信号。在波长移相干涉测量获得的超薄平行平板干涉信号后引入FAST-ICA算法,通过快速迭代算法分离超薄平板前后表面干涉条纹,然后通过时域移相算法处理分离后的干涉条纹,从而得到前后表面的面形信息。该方法主要针对波长移相中超薄平行平板的条纹图处理,为这类超薄平板前后表面信息分离提供了一种方法,可用于实现超薄平板前后表面面形的分离测量。解决了传统波长移相测量超薄平板前后表面干涉光强信号频谱叠加的问题,为实现波长移相对超薄平板表面面形的测量提供了一种方法。
附图说明
图1是本发明分离处理方法流程图。
图2是干涉图处理示例。
具体实施方式
下面结合附图和优选实施例对本发明进行说明。
参见图1,本波长移相中基于FAST-ICA算法的超薄平板干涉信号分离处理方法,首先利用波长移相干涉仪通过波长调谐测量超薄平行平板,得到包含前后表面面形信息的一系列波长移相干涉图。然后对采集到的这一系列波长移相干涉信号去均值化,接着利用白化矩阵将均值化后的干涉信号进行白化。下一步根据独立分量提取算法,利用分离矩阵处理白化后的去均值信号。判断分离矩阵处理白化后的去均值信号是否收敛。如果收敛,则分离矩阵已最优;如果不收敛,则分离矩阵不是最优,需要重新优化分离矩阵,然后再利用优化后的分离矩阵处理白化后的去均值信号,直到优化后的分离矩阵处理白化后的去均值信号收敛,即可得到待分离的前后表面干涉信号。最后将分离后的前后表面干涉信号利用时域移相算法,计算求解得到前后表面面形。
利用FAST-ICA算法分离超薄平板前后表面干涉光强信息,具体方法如下:
假设前表面干涉调制信号I1为式(1)表示,如图2中(a)所示:
I1=sin(2*π*f1*t) (1)
式中f1为频率,t代表时间,取f1=60。
假设后表面干涉调制信号I2为式(2)表示,如图2中(b)所示:
I2=cos(2*π*f2*t); (2)
式中f2为频率,取f2=66。
假设寄生条纹信号I3为随机值,用式(3)表示,如图2中(c)所示:
I3=randn(size(t)) (3)
则前后表面干涉混合产生的调制信号I(包含寄生条纹)用式(4)表示,如图2中(d)所示:
I=A*sin(2*π*f1*t)+B* cos(2*π*f2*t)+C* randn(size(t)) (4)
式中A,B,C为任意值。
利用matlab模拟FAST-ICA算法处理调制信号I,即式(4),分离后的信号如图2中(e),(f),(g)所示。根据仿真结果显示,FAST-ICA可以较精确的实现超薄平板前后表面干涉信号分离。
Claims (1)
1.一种波长移相中基于FAST-ICA算法的超薄平板干涉信号分离处理方法,其特征在于,由以下步骤实现:
步骤一:波长移相干涉仪通过波长调谐测量超薄平行平板,得到包含前后表面面形信息的一系列波长移相干涉图;
步骤二:根据像素列出时域序列干涉信号,对序列干涉信号去均值化:首先求出干涉信号的均值,然后将干涉信号减去均值,得到去均值化后的干涉信号;
步骤三:利用白化矩阵将均值化后的干涉信号进行白化;
步骤四:根据独立分量提取算法,利用分离矩阵处理白化后的去均值信号;
步骤五:判断是否收敛:如果收敛,则分离矩阵已最优,即得到待分离的前后表面干涉信号,执行第六步;如果不收敛,则分离矩阵不是最优,需要重新优化分离矩阵,再返回第四步;
步骤六:将分离后的前后表面干涉信号利用时域移相算法,计算求解得到前后表面面形。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710180047.4A CN107024186A (zh) | 2017-03-24 | 2017-03-24 | 波长移相中基于fast‑ica算法的超薄平板干涉信号分离处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710180047.4A CN107024186A (zh) | 2017-03-24 | 2017-03-24 | 波长移相中基于fast‑ica算法的超薄平板干涉信号分离处理方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107024186A true CN107024186A (zh) | 2017-08-08 |
Family
ID=59525864
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710180047.4A Pending CN107024186A (zh) | 2017-03-24 | 2017-03-24 | 波长移相中基于fast‑ica算法的超薄平板干涉信号分离处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107024186A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108489422A (zh) * | 2018-03-12 | 2018-09-04 | 四川大学 | 一种变频相移最小二乘迭代叠加面形分离的方法 |
CN111798380A (zh) * | 2019-04-09 | 2020-10-20 | Tcl集团股份有限公司 | 一种图像处理方法、图像处理装置及终端设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103424134A (zh) * | 2013-07-30 | 2013-12-04 | 北京交通大学 | 一种多参量同时测量的光纤传感方法 |
CN103536282A (zh) * | 2013-11-06 | 2014-01-29 | 中国人民解放军第三军医大学 | 基于Fast-ICA法的磁感应心肺活动信号分离方法 |
CN103763049A (zh) * | 2014-01-23 | 2014-04-30 | 中国科学院半导体研究所 | 一种基于固定点独立分量分析算法的合作频谱感知方法 |
CN105910987A (zh) * | 2016-04-27 | 2016-08-31 | 广东工业大学 | 基于盲分离技术的多层复杂结构材料高精度透视检测方法 |
-
2017
- 2017-03-24 CN CN201710180047.4A patent/CN107024186A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103424134A (zh) * | 2013-07-30 | 2013-12-04 | 北京交通大学 | 一种多参量同时测量的光纤传感方法 |
CN103536282A (zh) * | 2013-11-06 | 2014-01-29 | 中国人民解放军第三军医大学 | 基于Fast-ICA法的磁感应心肺活动信号分离方法 |
CN103763049A (zh) * | 2014-01-23 | 2014-04-30 | 中国科学院半导体研究所 | 一种基于固定点独立分量分析算法的合作频谱感知方法 |
CN105910987A (zh) * | 2016-04-27 | 2016-08-31 | 广东工业大学 | 基于盲分离技术的多层复杂结构材料高精度透视检测方法 |
Non-Patent Citations (4)
Title |
---|
于瀛洁 等: "加权多步波长移相算法研究", 《仪器仪表学报》 * |
刘永坚 等: "《现代电子战支援侦察系统分析与设计》", 31 December 2016, 国防工业出版社 * |
庞中华 等: "《系统辨识与自适应控制MATLAB仿真 修订版》", 31 August 2013, 北京航空航天大学出版社 * |
詹科利: "《大学物理 下》", 31 May 2015, 东华大学出版社 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108489422A (zh) * | 2018-03-12 | 2018-09-04 | 四川大学 | 一种变频相移最小二乘迭代叠加面形分离的方法 |
CN111798380A (zh) * | 2019-04-09 | 2020-10-20 | Tcl集团股份有限公司 | 一种图像处理方法、图像处理装置及终端设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4885154B2 (ja) | 複数波長による表面形状の測定方法およびこれを用いた装置 | |
CN113091634B (zh) | 一种适用于白光扫描干涉的微观形貌快速测量方法 | |
CN102322822B (zh) | 一种三频彩色条纹投影三维测量方法 | |
CN102155924A (zh) | 基于绝对相位恢复的四步相移方法 | |
CN101986098A (zh) | 基于三色光栅投影的傅里叶变换三维测量法 | |
CN204116710U (zh) | 一种两步衍射相位成像系统 | |
US20190137265A1 (en) | Method and system for the optical inspection and measurement of a face of an object | |
CN107024186A (zh) | 波长移相中基于fast‑ica算法的超薄平板干涉信号分离处理方法 | |
CN108053392A (zh) | 结合深度测量的车轮踏面擦伤视觉检测方法 | |
EP3062065A1 (en) | Imaging device and phase difference detection method | |
García-Isáis et al. | One shot profilometry using a composite fringe pattern | |
US20160301855A1 (en) | Imaging device and phase difference detection method | |
CN104697454A (zh) | 一种基于双光栅的细丝直径测量方法和装置 | |
CN106247980A (zh) | 基于白光干涉彩色条纹处理的多波长相移干涉测量方法 | |
CN110779454B (zh) | 基于双通道结构散斑互相关算法的改进数字图像相关方法 | |
KR101637252B1 (ko) | 백색광 위상천이 간섭계의 간섭신호 차수 오차 보정방법 | |
CN102519375B (zh) | 基于光循环与谱域载频的超大量程间距测量系统 | |
CN111121661B (zh) | 用于光滑表面形貌测量的窄带非单色光n+1幅相移测试算法 | |
CN101718533B (zh) | 镜像莫尔测量装置及方法 | |
CN106556340A (zh) | 基于调制度的宽光谱干涉的零级条纹的寻找方法 | |
CN112212806A (zh) | 一种基于相位信息导向的立体相位展开方法 | |
CN106949846B (zh) | 处理波长移相干涉中超薄平行平板交叉式干涉条纹方法 | |
CN102645174A (zh) | 用于高信噪比散斑干涉的实时相移方法 | |
US7304745B2 (en) | Phase measuring method and apparatus for multi-frequency interferometry | |
CN108775861A (zh) | 一种基于等效波长π/(2k)相移的双波长时域相位解调方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170808 |