CN102645280B - 高效的光谱还原方法 - Google Patents
高效的光谱还原方法 Download PDFInfo
- Publication number
- CN102645280B CN102645280B CN201210128157.3A CN201210128157A CN102645280B CN 102645280 B CN102645280 B CN 102645280B CN 201210128157 A CN201210128157 A CN 201210128157A CN 102645280 B CN102645280 B CN 102645280B
- Authority
- CN
- China
- Prior art keywords
- sequence
- complex numbers
- interferogram
- light distribution
- zero
- 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.)
- Expired - Fee Related
Links
Landscapes
- Spectrometry And Color Measurement (AREA)
Abstract
本发明涉及一种高效的光谱还原方法,该还原方法步骤为首先对干涉图进行过零单边采样并单边补零,干涉图短双边采样并双边补零;其次,把补零后的过零单边数据和短双边数据组成复数序列,其中过零单边数据为实部,短双边数据为虚部;然后采用下陷式三角切趾函数对复数序列切趾,再进行快速傅里叶变换;最后根据共轭对称复数序列的性质,进行光谱相位校正,得到准确的光谱分布。
Description
技术领域
本发明涉及一种新的干涉图切趾函数及其高效的光谱还原方法。
背景技术
在傅里叶变换光谱仪干涉图的实际测量中,一般采样到有限的极大光程差L处,实际上意味着强制干涉图函数从这一点开始及其后突然下降到零,这导致干涉图在该处出现尖锐的不连续性,从而引起变换光谱在颇大波数范围内的“微扰”,即旁瓣。为了抑制旁瓣对光谱失真造成的影响,可以采用三角切趾函数对干涉图加权,以缓和在极大光程差L处出现的不连续性。
另外,理想情况下,干涉图关于零光程差点是对称的,但是由于检测干涉信号很难从零光程差开始、分束板存在吸收损耗和不均匀性、电子线路误差等,导致干涉图出现一定程度的不对称性和相位误差,降低探测灵敏度,并导致光谱失真。因此,相位误差的校正是傅里叶变换光谱分析仪计算的重要环节。
在乘积法相位误差校正中,当采用过零单边干涉图进行光谱复原时,短双边干涉图数据被重复使用两次,导致光谱曲线变得平滑、分辨率下降。为了减小过零单边干涉图中心条纹领域内的数据被计算两次带来的误差,过零单边干涉图需要乘以一个权重函数,即采用Mertz提出的切趾函数,对短双边干涉图各点乘以不同的系数,使具有相同光程差的两点的干涉光强之和保持单点的幅值。
目前,能够抑制旁瓣并且减小短双边干涉图数据被使用两次带来的误差,获得光谱的最佳方案如下:
步骤(一):通过测量得到过零单边干涉图光强分布Is(x)和短双边干涉图光强分布Id(x),最大光程差分别为L1、L2(L2<L1),采样间隔为T;
步骤(二):在光程差-L1~-L2之间对过零单边采样数据补零,补零间隔为T,然后采用Mertz切趾函数对补零后的过零单边干涉图进行切趾,切趾之后的光强分布I1(x)为:
其中,x是光程差,A1(x)是Mertz切趾函数,表示如下
步骤(三):在光程差-L1~-L2和L2~L1之间对短双边采样数据进行补零,补零间隔为T,补零后得到新的光强分布I2(x)为:
步骤(四):根据实序列傅里叶变换的性质,把光强分布序列I1(x)和I2(x)分别作为一组复数序列的实部和虚部,组成一组复数序列f0(x),表示如下
f0(x)=I1(x)+i·I2(x) (4)
然后采用三角切趾函数对此复数序列进行切趾,切趾后的复数序列f(x)为:
f(x)=f0(x)·A2(x) (5)
其中,A2(x)是三角切趾函数,表示如下:
步骤(五):对复数序列f(x)进行快速傅里叶变换,得到变换后的复数序列 F(k),根据共轭对称复数序列的性质,把复数序列F(k)分成复数序列C(k)和D(k)两部分,其中
其中,a1(k)、b1(k)分别是C(k)中点k的实部和虚部,a2(k)、b2(k)分别是D(k)中点k的实部和虚部,N是光强分布I1(x)的总点数,对复数序列F(k)进行相位校正,就可以得到相位校正后点k的光强了,校正方法以及校正后点k的光强分布如下:
步骤(六):对光强进行归一化。假设光强分布B0(k)的最大值是a,最小值是b,对B0(k)采用如下方法进行归一化处理,就可以得到归一化的光强分布B(k):
根据式(11)中点数k与波数v的关系,就可以把光强分布B(k)转换为光强关于波数v的分布;
假设单边采样点数为N1,过零采样点数为N2,则此方案中步骤(一)至步骤(四)的计算量为
加法:4N1+2N2-1 (12)
乘法:4N1+6N2+1 (13)
采用这种方案获得还原光谱,为了减小过零单边干涉图中心条纹领域内的 数据被计算两次带来的误差,过零单边干涉图需要乘以Mertz切趾函数,使短双边干涉图具有相同光程差的两点的干涉光强之和保持单点的幅值。为了抑制旁瓣对光谱失真的影响,还需要采用三角切趾函数对干涉图切趾处理。由此可见,在已有的乘积法相位误差校正方法中,针对适用于三角切趾函数加权的干涉图,为了减小短双边干涉图数据被使用两次带来的误差和抑制旁瓣,需要分别采用Mertz切趾函数和三角切趾函数进行两次切趾处理,计算量大,效率低;
另外,使用Mertz切趾函数对过零单边干涉图加权,减小短双边干涉图数据被使用两次带来的误差,该方法具有明显的缺点:当相位误差不满足随波数缓变的条件时,干涉图的非对称性比较严重,此时这种加权方法不但不能减弱非对称性,反而会加重不对称性,达不到预期的加权效果,导致光谱失真,计算精度低。而且随着过零点数据量的减少这种现象变得更加明显。
因此,采用何种切趾函数,在抑制旁瓣的同时,能够减小短双边干涉图数据被使用两次带来的误差,减少校正相位误差的计算量;以及当相位误差不满足随波数缓变的条件时,采用何种切趾函数,对非对称性比较严重的干涉图加权,减小干涉图的非对称性带来的计算误差,已经成为亟待解决的问题。
发明内容
针对上述缺点,本发明的目的在于克服上述技术问题的缺陷,提出一种傅里叶变换光谱仪下陷式三角切趾函数及其高效的光谱还原方法,研究结果表明,本发明的下陷式三角切趾函数在抑制旁瓣的同时,能够使短双边干涉图中具有相同光程差的两个对应点的光强与各自对应的旋转因子的乘积之和为加权之前这两个对应点的光强与各自对应的旋转因子的乘积之和的平均值,减小短双边干涉图数据被使用两次带来的误差。当相位误差不满足随波数缓变的条件时, 此切趾函数也能够减弱干涉图的非对称性,达到预期的加权效果。
本发明的目的之一是通过以下技术方案来实现的:
步骤(一):首先通过测量得到干涉图光强分布I(x);
步骤(二):提出下陷式三角切趾函数A(x):
其中,x是光程差,L1、L2分别是过零单边干涉图和短双边干涉图的最大光程差;
步骤(三):在光程差-L2~L1范围内,对干涉图进行过零单边采样,采样间隔为T,在光程差-L1~-L2之间对过零单边采样数据补零,补零间隔为T,补零后的光强分布I1(x)为
其中,Is(x)是过零单边采样的光强分布;
步骤(四):在光程差-L2~L2范围内,对干涉图进行短双边采样,采样间隔为T,然后对短双边采样数据在光程差-L1~-L2和L2~L1之间进行补零,补零间隔为T,补零后得到新的光强分布I2(x)为:
其中,Id(x)是短双边采样的光强分布;
步骤(五):根据实序列傅里叶变换的性质,把光强分布序列I1(x)和I2(x)分别作为一组复数序列的实部和虚部,组成一组复数序列f0(x),表示如下:
f0(x)=I1(x)+i·I2(x) (4)
然后采用下陷式三角切趾函数对此复数序列进行切趾,切趾后的复数序列f(x)为:
f(x)=f0(x)·A(x) (5)
根据本发明的优选实施例,其中切趾之后,使短双边干涉图中具有相同光程差的两个对应点的光强与各自对应的旋转因子的乘积之和为加权之前这两个对应点的光强与各自对应的旋转因子的乘积之和的平均值。
步骤(六):对复数序列f(x)进行快速傅里叶变换,得到变换后的复数序列F(k),根据共轭对称复数序列的性质,把复数序列F(k)分成复数序列C(k)和D(k)两部分,其中
其中,a1(k)、b1(k)分别是C(k)中点k的实部和虚部,a2(k)、b2(k)分别是D(k)中点k的实部和虚部,N是光强分布I1(x)的总点数,对复数序列F(k)进行相位校正,就可以得到相位校正后点k的光强了,校正方法以及校正后点k的光强分布如下:
步骤(七):对光强进行归一化。假设光强分布B0(k)的最大值是a,最小值是b,对B0(k)采用如下方法进行归一化处理,就可以得到归一化的光强分布B(k):
根据式(10)中点数k与波数v的关系,就可以把光强分布B(k)转换为光强关于波数v的分布;
利用本发明的下陷式三角切趾函数及光谱还原方法只需对干涉图进行一次切趾处理,就可以获得高精度的光谱,并且具有高效率、高精度优点。
虽然在下文中将结合一些示例性实施及使用方法来描述本发明,但本领域技术人员应当理解,并不旨在将本发明限制于这些实施例。反之,旨在覆盖包含在所附的权利要求书所定义的本发明的精神与范围内的所有替代品、修正及等效物。
本发明的其它优点、目的和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目的和其它优点可以通过下面的说明书,权利要求书,以及附图中所特别指出的结构来实现和获得。
附图说明
下面结合附图对本发明的具体实施例作进一步详细的说明。
图1本发明实现反演光谱流程图
图2本发明计算效率提高曲线图
图3相位误差不满足随波数缓变的条件时,采集的短双边干涉图
图4相位误差不满足随波数缓变时的相位谱图
图5采用传统方法和本发明方法反演的光谱图
具体实施方式
下面将结合附图,对本发明的优选实施例进行详细的描述;应当理解,优 选实施例仅为了说明本发明,而不是为了限制本发明的保护范围。
图1示出了本发明的技术方案:首先对干涉图进行过零单边采样并单边补零,干涉图短双边采样并双边补零;其次,把补零后的过零单边采样数据和短双边采样数据组成复数序列,其中过零单边数据为实部,短双边数据为虚部;然后采用下陷式三角切趾函数对复数序列切趾,再进行快速傅里叶变换;最后根据共轭对称复数序列的性质,进行光谱相位校正,得到准确的光谱分布。
下陷式三角切趾函数对复数序列切趾获得高精度还原光谱的过程为:
步骤(一):首先通过测量得到干涉图光强分布I(x);
步骤(二):提出下陷式三角切趾函数A(x):
其中,x是光程差,L1、L2分别是过零单边干涉图和短双边干涉图的最大光程差;
步骤(三):在光程差-L2~L1范围内,对干涉图进行过零单边采样,采样间隔为T,在光程差-L1~-L2之间对过零单边采样数据补零,补零间隔为T,补零后的光强分布I1(x)为
其中,Is(x)是过零单边采样的光强分布;
步骤(四):在光程差-L2~L2范围内,对干涉图进行短双边采样,采样间隔为T,然后对短双边采样数据在光程差-L1~-L2和L2~L1之间进行补零,补零间隔为T,补零后的光强分布I2(x)为:
其中,Id(x)是短双边采样的光强分布;
步骤(五):根据实序列傅里叶变换的性质,把光强分布序列I1(x)和I2(x)分别作为一组复数序列的实部和虚部,组成一组复数序列f0(x),表示如下:
f0(x)=I1(x)+i·I2(x) (4)
然后采用下陷式三角切趾函数对此复数序列进行切趾,在抑制旁瓣的同时,使短双边中具有相同光程差的两个对应点的光强与各自对应的旋转因子的乘积之和为加权之前这两个对应点的光强与各自对应的旋转因子的乘积之和的平均值。举例来说,比如具有相同光程差的两个对应点分别是A、B,那么加权后A点的光强与A点对应的旋转因子的乘积+加权后B点的光强与B点对应的旋转因子的乘积=(加权前A点的光强与A点对应的旋转因子的乘积+加权前B点的光强与B点对应的旋转因子的乘积)/2。其中,不同点对应的旋转因子是不同的。切趾后的复数序列f(x)为:
f(x)=f0(x)·A(x) (5)
步骤(六):对复数序列f(x)进行快速傅里叶变换,得到变换后的复数序列F(k),根据共轭对称复数序列的性质,把复数序列F(k)分成复数序列C(k)和D(k)两部分,其中
其中,a1(k)、b1(k)分别是C(k)中点k的实部和虚部,a2(k)、b2(k)分别是D(k)中 点k的实部和虚部,N是光强分布I1(x)的总点数,对复数序列F(k)进行相位校正,就可以得到相位校正后点k的光强了,校正方法以及校正后点k的光强分布如下:
步骤(七):对光强进行归一化。假设光强分布B0(k)的最大值是a,最小值是b,对B0(k)采用如下方法进行归一化处理,就可以得到归一化的光强分布B(k):
根据式(10)中点数k与波数v的关系,就可以把光强分布B(k)转换为光强关于波数v的分布;
假设单边采样点数为N1,过零采样点数为N2,则此方案中步骤(一)至步骤(五)的计算量为:
加法:4N1-2 (11);
乘法:4N1+2N2-1 (12)。
由于本发明中步骤(六)、步骤(七)与现有的最佳方案中步骤(五)、步骤(六)的计算量相同,所以本发明可以节约的计算量为:
加法:2N2+1 (13)
乘法:4N2+2 (14)
本发明提出的下陷式三角切趾函数及光谱还原方法,在抑制旁瓣的同时,使短双边干涉图中具有相同光程差的两个对应点的光强与各自对应的旋转因子的乘积之和为加权之前这两个对应点的光强与各自对应的旋转因子的乘积之和 的平均值,减小了短双边干涉图数据被使用两次带来的误差。只需要一次切趾就可以获得光谱,计算效率更高。从背景技术现有最佳方案中公式(12)、(13)以及具体实施方式中公式(13)、(14)可以计算出,随着过零采样点数N2的增加,本发明的计算效率也将逐步得到提高,当单边采样点数N1=8200时,本发明与最好的现有技术相比,计算效率提高曲线如图2所示。
当相位误差不满足随波数缓变的条件时,干涉图的非对称性比较严重,通过以上措施,能够获得比Mertz乘积法更高精度的光谱。在仿真试验中,过零单边采样的最大光程差L1=0.4502cm,短双边采样的最大光程差L2=0.0136cm,采样间隔T=632.8×10-7cm,图3就是此时采集到的短双边干涉图,相位误差如图4所示。从图3和图4中可以看出,此时的相位误差不再满足随波数缓变的条件,采集的干涉图非对称性也比较严重,用Mertz切趾函数对此干涉图加权,就达不到使具有相同光程差的两点的干涉光强之和保持单点的幅值的效果。而采用本发明的下陷式三角切趾函数对此干涉图加权,在抑制旁瓣的同时,使短双边干涉图中具有相同光程差的两个对应点的光强与各自对应的旋转因子的乘积之和为加权之前这两个对应点的光强与各自对应的旋转因子的乘积之和的平均值,减小了短双边干涉图数据被使用两次带来的误差。与采用两次切趾获得的光谱相比,通过以上措施,能够获得更高精度的还原光谱,图5是分别采用两次切趾和本发明方法获得的光谱图,从图中可以看出,与现有的最好技术相比,本发明获得的光谱与实际光谱更加吻合,进一步减小了光谱失真,从而证明了本发明方法的可行性及有效性。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范 围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (2)
1.一种高效光谱还原方法,步骤为首先对干涉图过零单边采样并单边补零,干涉图短双边采样并双边补零;其次把补零后过零单边数据和短双边数据组成复数序列,其中过零单边数据为实部,短边数据为虚部;对复数序列进行快速傅里叶变换;最后根据共轭对称复数序列的性质,进行光谱相位校正,得到准确的光谱分布,其特征在于:在对复数序列进行快速傅里叶变换前,采用下陷式三角切趾函数对复数序列切趾。
2.如权利要求1所述的高效光谱还原方法,其特征在于:下陷式三角切趾函数对复数序列切趾获得光谱的过程为:
步骤(一):首先通过测量得到干涉图光强分布I(x);
步骤(二):提出下陷式三角切趾函数A(x):
步骤(三):在光程差-L2~L1范围内,对干涉图进行过零单边采样,采样间隔为T,在光程差-L1~-L2之间对过零单边采样数据补零,补零间隔为T,补零后的光强分布I1(x)为
步骤(四):在光程差-L2~L2范围内,对干涉图进行短双边采样,采样间隔为T,然后对短双边采样数据在光程差-L1~-L2和L2~L1之间进行补零,补零间隔为T,补零后得到新的光强分布I2(x)为:
步骤(五):根据实序列傅里叶变换的性质,把光强分布序列I1(x)和I2(x)分别作为一组复数序列的实部和虚部,组成一组复数序列f0(x),表示如下:
f0(x)=I1(x)+i·I2(x) (4)
然后采用下陷式三角切趾函数对此复数序列进行切趾,切趾后的复数序列f(x)为:
f(x)=f0(x)·A(x) (5)
切趾之后,使短双边中具有相同光程差的两个对应点的光强与各自对应的旋转因子的乘积之和为加权之前这两个对应点的光强与各自对应的旋转因子的乘积之和的平均值;
步骤(六):对复数序列f(x)进行快速傅里叶变换,得到变换后的复数序列F(k),根据共轭对称复数序列的性质,把复数序列F(k)分成复数序列C(k)和D(k)两部分,其中
其中,a1(k)、b1(k)分别是C(k)中点k的实部和虚部,a2(k)、b2(k)分别是D(k)中点k的实部和虚部,N是光强分布I1(x)的总点数,对复数序列F(k)进行相位校正,就可以得到相位校正后点k的光强了,校正方法以及校正后点k的光强分布如下:
步骤(七):对光强进行归一化;假设光强分布B0(k)的最大值是a,最小值是b,对B0(k)采用如下方法进行归一化处理,就可以得到归一化的光强分布B(k):
根据式(10)中点数k与波数v的关系,就可以把光强分布B(k)转换为光强关于波数v的分布:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210128157.3A CN102645280B (zh) | 2012-04-27 | 2012-04-27 | 高效的光谱还原方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210128157.3A CN102645280B (zh) | 2012-04-27 | 2012-04-27 | 高效的光谱还原方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102645280A CN102645280A (zh) | 2012-08-22 |
CN102645280B true CN102645280B (zh) | 2014-12-03 |
Family
ID=46658248
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210128157.3A Expired - Fee Related CN102645280B (zh) | 2012-04-27 | 2012-04-27 | 高效的光谱还原方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102645280B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109767544B (zh) * | 2017-01-10 | 2022-02-15 | 中钞印制技术研究院有限公司 | 有价证券的图像分析方法及图像分析系统 |
CN109840548A (zh) * | 2019-01-02 | 2019-06-04 | 广东工业大学 | 一种基于bp神经网络变压器故障诊断方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2372747B1 (en) * | 2010-03-31 | 2018-08-01 | Thermo Fisher Scientific (Bremen) GmbH | Methods and apparatus for producing a mass spectrum |
-
2012
- 2012-04-27 CN CN201210128157.3A patent/CN102645280B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102645280A (zh) | 2012-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102519597B (zh) | 一种傅里叶变换光谱仪相位校正切趾方法 | |
CN107490397B (zh) | 高精度自适应滤波fbg光谱快速寻峰方法 | |
CN103308804B (zh) | 基于快速k-s变换电能质量扰动信号时频参数提取方法 | |
CN101105446B (zh) | 差分吸收光谱空气质量检测系统 | |
CN102680020B (zh) | 一种基于波长调制光谱技术的气体参数在线测量方法 | |
CN103033481A (zh) | 基于fft的激光分析仪二次谐波滤波方法 | |
CN105954735B (zh) | 一种用于fmcw绝对距离测量技术中改进的高速色散失配校正方法 | |
CN106483374A (zh) | 一种基于Nuttall双窗全相位FFT的谐波间谐波检测方法 | |
CN102043091B (zh) | 数字化高精度相位检测器 | |
CN106771709A (zh) | 一种多端口网络的s参数去嵌入方法 | |
CN103645154A (zh) | 一种利用太赫兹光谱信号提取材料光学常数的方法 | |
CN105784672A (zh) | 一种基于双树复小波算法的毒品检测仪标准化方法 | |
CN104198416A (zh) | 一种光谱仪波长漂移引起的测量误差实时补偿方法 | |
CN108918446A (zh) | 一种超低浓度二氧化硫紫外差分特征提取算法 | |
CN103941090A (zh) | 基于谱线能量插值的谐波测量方法 | |
CN103644969B (zh) | 一种弹光调制干涉信号预处理的方法 | |
CN102645280B (zh) | 高效的光谱还原方法 | |
CN107036720A (zh) | 一种超短啁啾脉冲时域相位与频域相位测量方法及系统 | |
CN102818630B (zh) | 一种干涉型成像光谱仪的光谱定标方法 | |
CN101718816B (zh) | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 | |
CN103197143A (zh) | 基于汉宁窗fft算法与遍历滤波的谐波、间谐波检测方法 | |
CN103163078A (zh) | 一种提高椭偏仪测量精度的方法 | |
CN104215833A (zh) | 电力系统频率测量方法及装置 | |
CN104914439A (zh) | 一种超声波测距的双相位测量方法 | |
CN103630336B (zh) | 基于随机快轴方位角延迟阵列的动态干涉测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141203 Termination date: 20200427 |