CN103808268A - 一种白光干涉传感器解调信号处理方法 - Google Patents
一种白光干涉传感器解调信号处理方法 Download PDFInfo
- Publication number
- CN103808268A CN103808268A CN201210452137.1A CN201210452137A CN103808268A CN 103808268 A CN103808268 A CN 103808268A CN 201210452137 A CN201210452137 A CN 201210452137A CN 103808268 A CN103808268 A CN 103808268A
- Authority
- CN
- China
- Prior art keywords
- white light
- light interference
- filter
- signal
- sensor
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 7
- 230000003287 optical effect Effects 0.000 claims abstract description 21
- 238000006073 displacement reaction Methods 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 7
- 238000005070 sampling Methods 0.000 claims description 14
- 238000001914 filtration Methods 0.000 claims description 12
- 238000013507 mapping Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 4
- 239000013307 optical fiber Substances 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- DMYOHQBLOZMDLP-UHFFFAOYSA-N 1-[2-(2-hydroxy-3-piperidin-1-ylpropoxy)phenyl]-3-phenylpropan-1-one Chemical compound C1CCCCN1CC(O)COC1=CC=CC=C1C(=O)CCC1=CC=CC=C1 DMYOHQBLOZMDLP-UHFFFAOYSA-N 0.000 description 1
- 229920001079 Thiokol (polymer) Polymers 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明属于信号处理领域,特别是涉及一种白光干涉传感器解调信号处理方法。本发明采取的方案为:对于白光干涉原理的光学传感器,把CCD读取的含有白光干涉条纹的原始信号用差值处理方式去除干涉信号中的基底噪声,再利用自适应滤波器去除白光干涉信号中的环境噪声。最后,从去噪后的信号中提取出干涉条纹,根据干涉条纹的分布确定传感器的位移量。本发明有效解决了从CCD中读取的白光干涉信号噪声大,干涉条纹不易提取,且既有高频噪声又有低频噪声,不易进行处理的问题。从而保证了利用白光干涉条纹计算光传感器的位移量。
Description
技术领域
本发明属于信号处理领域,特别是涉及一种白光干涉传感器解调信号处理方法。
背景技术
白光干涉型位移传感器是光位移传感器研制领域的一个重要分支。国外早在上世纪90年代开展了白光干涉型位移传感器研究。ATK Aerospace公司的Thiokol Propulsion喷气推进部分在火箭发动机上已经验证了使用楔形介质膜的白光干涉型线位移传感器,并得到了满意的结果,美国Davidson公司也正在NAVY的先进战舰SC-21上试验了这种新型的传感器。
白光干涉信号解调是白光干涉型位移传感器的重要组成部分。采用光楔为解调器件的白光干涉信号光楔解调系统其最大特点是没有运动部件,从理论上讲具有很好的长期可靠性,另外,由于该解调系统不需要昂贵的宽带光源和光谱接收器件,成本非常低,将是今后白光干涉解调的发展方向。白光干涉信号经过光楔解调系统后被CCD接收,CCD接收到的光信号包括干涉信号和背景噪声,通常情况下传感器有用信号会被淹没在噪音中,如何将有用的干涉信号从噪声中滤出,是本发明的主要内容。
发明内容
本发明的目的:提出一种能自适应的对光学传感器的白光干涉传感器解调信号处理办法。
本发明采取的技术方案为:一种白光干涉传感器解调信号处理方法,包括以下步骤:
步骤一、数据采集:读取光学传感器CCD原始采样数据;
步骤二、差值计算:将包含白光干涉条纹的采样数据I(t)减去未包含白光干涉条纹的采样数据I(t0),从而得到去除了基底噪声后的采样数据D(t);
步骤三、对滤波器的截止频率下限ωp及截止频率上限ωs进行参数选取:ωp从0~15π/16之间逐渐变化,由于ωs=ωp+π/16,则ωs从π/16~π之间随之逐渐变化,同时用该截止频率构造的滤波器滤波得到的结果进行统计,得出白光干涉条纹积分能量最大值时对应的ωp和ωs;
步骤四、滤波处理:利用步骤三得到的滤波器参数ωp和ωs,对去除了基底噪声后的采样数据D(t)进行滤波处理,采样数据D(t)与滤波器Filter(ωp,ωs)进行卷积运算得到F(t),F(t)=D(t)*Filter(ωp,ωs);
步骤五、信号提取:遍历滤波处理后的信号F(t),找出其中绝对值最大点对应的CCD位置;
步骤六:位移映射:根据映射关系,得到传感器的位移量。
优选地,步骤四采用巴特沃茨滤波器进行卷积运算。
本发明具有的优点和有益效果:本发明有效解决了从CCD中读取的白光干涉信号噪声大,干涉条纹不易提取,且既有高频噪声又有低频噪声,不易进行处理的问题。从而保证了利用白光干涉条纹计算光传感器的位移量。
附图说明
图1是光学传感器干涉信号处理实验平台示意图;
图2是白光干涉传感器解调信号处理流程图;
图3是包含噪声的白光干涉信号示意图;
图4是差值处理后的白光干涉信号示意图;
图5是自适应滤波后的白光干涉条纹示意图。
其中,1:传感器干涉信号,2:光纤,3:扩束透镜,4:柱透镜,5:光楔,6:CCD感光器件;
具体实施方式
对于白光干涉原理的光学传感器,把CCD读取的含有白光干涉条纹的原始信号用差值处理方式,去除干涉信号中的基底噪声。随后,利用本发明提出的自适应滤波器,根据白光干涉信号的特征,计算滤波器的截止频率参数,用构造好的滤波器去除白光干涉信号中的环境噪声。最后,从去噪后的信号中提取出干涉条纹,根据干涉条纹的分布确定传感器的位移量。
本发明对白光干涉传感器解调信号进行处理。首先将含有干涉条纹的采样信号减去上一时刻未含有干涉条纹的采样信号,从而消除CCD读取的原始信号中的基底噪声。然后根据所得到的信号的特点进行计算,得到自适应调整滤波器的参数,进行滤波处理,进一步消除掉CCD信号中的环境噪声。将白光干涉条纹很好的从CCD源信号中提取出来。从而解决了从CCD中读取的白光干涉信号噪声大,且既有高频噪声又有低频噪声,不易进行处理的问题。从而保证了利用白光干涉条纹计算光传感器的位移量。
下面结合附图对本发明做详细说明,本发明的光学传感器干涉信号处理实验平台如图1所示,光学传感器的干涉信号从光纤中发出,到达扩束透镜的表面。经过透镜折射后,光信号穿过柱透镜。通过柱透镜后,光信号照射在光楔上。经过光楔的反射,白光解调信号照射在CCD表面。CCD采集的信号由计算机收集并存储下来。
一种白光干涉传感器解调信号处理方法,包括以下步骤:
步骤一、数据采集:利用光学传感器干涉信号处理实验平台读取光学传感器CCD原始采样数据,数据以数值的形式存储在计算机中;
步骤二、差值计算:将没有干涉条纹时的CCD采样信号存储在计算机中,用I(t0)表示此时的信号。在信号源加载干涉条纹之后,再次记录CCD感光后的采样信号,用I(t)表示此时的信号。随后,用包含白光干涉条纹的采样数据I(t)减去未包含白光干涉条纹的采样数据I(t0),从而得到去除了基底噪声后的采样数据D(t);
D(t)=I(t)-I(t0);
步骤三、对滤波器的截止频率下限ωp及截止频率上限ωs进行参数选取:ωp从0~15π/16之间逐渐变化,由于ωs=ωp+π/16,则ωs从π/16~π之间随之逐渐变化,同时用该截止频率构造的滤波器滤波得到的结果进行统计。对于每次滤波处理后的信号进行积分运算,记录下白光干涉条纹积分能量值最大时对应的ωp和ωs,作为构造滤波器的截止频率下限和截止频率上限;
步骤四、滤波处理:利用步骤三得到的滤波器参数ωp和ωs,对去除了基底噪声后的采样数据D(t)进行滤波处理,采样数据D(t)与滤波器Filter(ωp,ωs)都是一维数据矩阵,进行卷积运算得到F(t),F(t)=D(t)*Filter(ωp,ωs);
步骤五、信号提取:遍历滤波处理后的信号F(t),找出其中绝对值最大点对应的CCD位置;
步骤六:位移映射:根据映射关系,得到传感器的位移量。
实例
以某一具体事例为例,对本发明做进一步详细描述。
步骤一、数据采集:光学传感器干涉信号处理实验平台如图1所示,光学传感器的干涉信号从光纤中发出,到达扩束透镜的表面。经过透镜折射后,光信号穿过柱透镜。通过柱透镜后,光信号照射在光楔上。经过光楔的反射,白光解调信号照射在CCD表面。CCD采集的信号由计算机收集并存储下来,在计算机上显示的结果如图3所示;
步骤二、差值计算:将没有干涉条纹时的CCD采样信号存储在计算机中,用I(t0)表示此时的信号。在信号源加载干涉条纹之后,再次记录CCD感光后的采样信号,用I(t)表示此时的信号。随后,用包含白光干涉条纹的采样数据I(t)减去未包含白光干涉条纹的采样数据I(t0),从而得到去除了基底噪声后的采样数据D(t),得到的结果如图4所示,计算公式如下;
D(t)=I(t)-I(t0);
步骤三、对滤波器的截止频率下限ωp及截止频率上限ωs进行参数选取:ωp从0~15π/16之间逐渐变化,由于ωs=ωp+π/16,则ωs从π/16~π之间随之逐渐变化,同时用该截止频率构造的滤波器滤波得到的结果进行统计,得出白光干涉条纹积分能量最大值时对应的ωp和ωs,此时的ωp=5π/16,而ωs=6π/16;
步骤四、滤波处理:利用步骤三得到的滤波器参数ωp=5π/16和ωs=6π/16,对去除了基底噪声后的采样数据D(t)进行滤波处理,滤波器为巴特沃茨滤波器。采样数据D(t)与滤波器Filter(ωp,ωs)进行卷积运算得到F(t),F(t)=D(t)*Filter(ωp,ωs),如图5所示;
步骤五、信号提取:遍历滤波处理后的信号F(t),找出其中绝对值最大点对应的CCD位置为第1476个像素点;
步骤六:位移映射:根据映射关系,得到传感器的位移量。
Claims (2)
1.一种白光干涉传感器解调信号处理方法,其特征在于,包括以下步骤:
步骤一、数据采集:读取光学传感器CCD原始采样数据;
步骤二、差值计算:将包含白光干涉条纹的采样数据I(t)减去未包含白光干涉条纹的采样数据I(t0),从而得到去除了基底噪声后的采样数据D(t);
步骤三、对滤波器的截止频率下限ωp及截止频率上限ωs进行参数选取:ωp从0~15π/16之间逐渐变化,由于ωs=ωp+π/16,则ωs从π/16~π之间随之逐渐变化,同时用该截止频率构造的滤波器滤波得到的结果进行统计,得出白光干涉条纹积分能量最大值时对应的ωp和ωs;
步骤四、滤波处理:利用步骤三得到的滤波器参数ωp和ωs,对去除了基底噪声后的采样数据D(t)进行滤波处理,采样数据D(t)与滤波器Filter(ωp,ωs)进行卷积运算得到F(t),F(t)=D(t)*Filter(ωp,ωs);
步骤五、信号提取:遍历滤波处理后的信号F(t),找出其中绝对值最大点对应的CCD位置;
步骤六:位移映射:根据映射关系,得到传感器的位移量。
2.根据权利要求1所述的白光干涉传感器解调信号处理方法,其特征在于,步骤四采用巴特沃茨滤波器进行卷积运算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210452137.1A CN103808268A (zh) | 2012-11-12 | 2012-11-12 | 一种白光干涉传感器解调信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210452137.1A CN103808268A (zh) | 2012-11-12 | 2012-11-12 | 一种白光干涉传感器解调信号处理方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103808268A true CN103808268A (zh) | 2014-05-21 |
Family
ID=50705387
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210452137.1A Pending CN103808268A (zh) | 2012-11-12 | 2012-11-12 | 一种白光干涉传感器解调信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103808268A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104792270A (zh) * | 2015-04-16 | 2015-07-22 | 石家庄铁路职业技术学院 | 一种基于白光数字图像频域分析法的深基坑检测方法 |
CN108645431A (zh) * | 2018-05-30 | 2018-10-12 | 西安工业大学 | 一种用于光纤法珀传感器腔长相关解调的拟合寻峰方法 |
CN119223175A (zh) * | 2024-12-04 | 2024-12-31 | 上海拜安传感技术有限公司 | 一种白光干涉式mems光纤位移传感器 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1556914A (zh) * | 2001-09-21 | 2004-12-22 | Kmac株式会社 | 利用二维反射计测量多层薄膜的厚度轮廓和折射率分布的装置及其测量方法 |
US20050046862A1 (en) * | 2003-08-25 | 2005-03-03 | Ivan Melnyk | Fiber optic sensing device for measuring a physical parameter |
CN101315301A (zh) * | 2008-06-18 | 2008-12-03 | 中国科学院上海光学精密机械研究所 | 双光路去噪声载波-包络相位测量装置 |
CN101561536A (zh) * | 2008-04-16 | 2009-10-21 | 中国科学院半导体研究所 | 一种抗振动与噪声干扰的光纤传感波长解调系统 |
CN101608932A (zh) * | 2008-06-19 | 2009-12-23 | 上海前所光电科技有限公司 | 一种分组同步式光纤传感分析仪 |
-
2012
- 2012-11-12 CN CN201210452137.1A patent/CN103808268A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1556914A (zh) * | 2001-09-21 | 2004-12-22 | Kmac株式会社 | 利用二维反射计测量多层薄膜的厚度轮廓和折射率分布的装置及其测量方法 |
US20050046862A1 (en) * | 2003-08-25 | 2005-03-03 | Ivan Melnyk | Fiber optic sensing device for measuring a physical parameter |
CN101561536A (zh) * | 2008-04-16 | 2009-10-21 | 中国科学院半导体研究所 | 一种抗振动与噪声干扰的光纤传感波长解调系统 |
CN101315301A (zh) * | 2008-06-18 | 2008-12-03 | 中国科学院上海光学精密机械研究所 | 双光路去噪声载波-包络相位测量装置 |
CN101608932A (zh) * | 2008-06-19 | 2009-12-23 | 上海前所光电科技有限公司 | 一种分组同步式光纤传感分析仪 |
Non-Patent Citations (2)
Title |
---|
杨玉孝等: "光纤白光干涉法与膜厚纳米测量新技术研究", 《光子学报》 * |
王选择等: "应用光栅刻线与CCD 测量运动位移的研究", 《光电工程》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104792270A (zh) * | 2015-04-16 | 2015-07-22 | 石家庄铁路职业技术学院 | 一种基于白光数字图像频域分析法的深基坑检测方法 |
CN108645431A (zh) * | 2018-05-30 | 2018-10-12 | 西安工业大学 | 一种用于光纤法珀传感器腔长相关解调的拟合寻峰方法 |
CN108645431B (zh) * | 2018-05-30 | 2020-08-07 | 西安工业大学 | 一种用于光纤法珀传感器腔长相关解调的拟合寻峰方法 |
CN119223175A (zh) * | 2024-12-04 | 2024-12-31 | 上海拜安传感技术有限公司 | 一种白光干涉式mems光纤位移传感器 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102692216B (zh) | 基于机器视觉技术的实时光纤绕制缺陷检测方法 | |
CN109726642A (zh) | 一种基于变分模态分解的分布式光纤振动信号的降噪方法 | |
CN103940371A (zh) | 一种用于跃变物体的高精度三维面形测量的方法 | |
CN114646381B (zh) | 一种旋转机械振动测量方法、系统、设备及存储介质 | |
CN112669265B (zh) | 基于傅里叶变换和图像梯度特征实现表面缺陷检测的方法 | |
CN105761231B (zh) | 一种用于去除高分辨率天文图像中条纹噪声的方法 | |
CN107506710A (zh) | 一种滚动轴承复合故障提取方法 | |
CN105133471A (zh) | 一种基于线结构光路面检测系统的路面深度图像生产方法 | |
CN103808268A (zh) | 一种白光干涉传感器解调信号处理方法 | |
CN107167087A (zh) | 一种基于经验小波变换的物体变形测量方法 | |
CN104809702A (zh) | 基于频率域处理的脉诊曲线图像网格消除方法 | |
CN103679790B (zh) | 一种基于图像处理技术的井下煤仓煤位检测方法 | |
CN109425897B (zh) | 消除地震资料野值干扰的方法及系统 | |
CN112903296B (zh) | 一种滚动轴承故障检测方法及系统 | |
Yang et al. | Multi-scale edge detection of wood defect images based on the dyadic wavelet transform | |
Ma et al. | Edge detection on pantograph slide image | |
CN105844593B (zh) | 一种单幅干涉圆条纹预处理的自动化处理方法 | |
CN108053379A (zh) | 一种基于改进的变分模态分解的dspi相位提取方法 | |
CN118067294A (zh) | 一种基于计算机视觉的斜拉桥拉索索力检测方法 | |
CN103345727A (zh) | 一种二元光学图像光谱的重建方法 | |
CN104077772A (zh) | 基于模糊图像相关和分形小波插值的mems面内微运动测量方法 | |
CN115294286A (zh) | 一种消除杂散频谱的数字全息物像频谱自适应提取方法 | |
CN116152716A (zh) | 一种双目视觉动力学模态参数识别中丢失模态的辨识方法 | |
Zhang et al. | An EMD-based denoising method for lidar signal | |
CN103440661B (zh) | 基于hsv空间与区域选择的微米木纤维直径检测算法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140521 |