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
- signal
- filter
- white
- 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
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 (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104792270A (zh) * | 2015-04-16 | 2015-07-22 | 石家庄铁路职业技术学院 | 一种基于白光数字图像频域分析法的深基坑检测方法 |
CN108645431A (zh) * | 2018-05-30 | 2018-10-12 | 西安工业大学 | 一种用于光纤法珀传感器腔长相关解调的拟合寻峰方法 |
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 (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 | 西安工业大学 | 一种用于光纤法珀传感器腔长相关解调的拟合寻峰方法 |
CN108645431B (zh) * | 2018-05-30 | 2020-08-07 | 西安工业大学 | 一种用于光纤法珀传感器腔长相关解调的拟合寻峰方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104006962A (zh) | 一种齿轮故障特征提取方法及系统 | |
CN103455813A (zh) | 一种ccd图像测量系统光斑中心定位的方法 | |
CN108225764A (zh) | 一种基于包络提取的高精度无键相信号阶次跟踪方法及系统 | |
CN103292792A (zh) | 一种适用海底探测与假地形处理的实测svp重构方法 | |
CN101929859B (zh) | 一种基于图像全帧扫描的空间碎片检测方法 | |
CN112069962B (zh) | 一种基于图像识别强噪声背景下振动频谱的方法 | |
CN112150606A (zh) | 一种基于点云数据的螺纹表面三维重构方法 | |
CN103808268A (zh) | 一种白光干涉传感器解调信号处理方法 | |
CN104809702A (zh) | 基于频率域处理的脉诊曲线图像网格消除方法 | |
CN111693283B (zh) | 一种轴承非稳态故障特征无转速提取方法 | |
Yang et al. | Multi-scale edge detection of wood defect images based on the dyadic wavelet transform | |
CN103679790A (zh) | 一种基于图像处理技术的井下煤仓煤位检测方法 | |
CN103487804B (zh) | 基于冰雷达技术的南极冰下湖半自动识别方法 | |
CN202794336U (zh) | 一种基于欠采样技术的射电天文窄带信号数字分析终端 | |
CN117974733A (zh) | 适用于环境变化的高精度点云匹配方法和系统及存储介质 | |
CN104215339B (zh) | 一种基于连续远场的波前复原系统及方法 | |
CN115061136B (zh) | 一种基于sar图像的河湖岸线变化点检测方法及系统 | |
CN103956013A (zh) | 一种风雨扰动信号实时判定方法 | |
CN114594503B (zh) | 一种浅海地形反演方法、计算机设备及存储介质 | |
CN106761692B (zh) | 一种压力传感器测温方法 | |
CN113985276B (zh) | 风力发电机组的故障诊断方法和装置 | |
CN114088077B (zh) | 一种改进的半球谐振陀螺信号去噪方法 | |
CN112149626B (zh) | 一种发动机爆震特征频率提取方法 | |
CN103903272B (zh) | 一种基于Hadoop的StaMPS算法并行化处理方法 | |
Huang et al. | Image Stitching Algorithm Based on Region Division for Underwater Dam Crack Image |
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 |