CN114397262A - 一种傅里叶变换红外光谱仪波数漂移的校正方法及系统 - Google Patents

一种傅里叶变换红外光谱仪波数漂移的校正方法及系统 Download PDF

Info

Publication number
CN114397262A
CN114397262A CN202111595873.8A CN202111595873A CN114397262A CN 114397262 A CN114397262 A CN 114397262A CN 202111595873 A CN202111595873 A CN 202111595873A CN 114397262 A CN114397262 A CN 114397262A
Authority
CN
China
Prior art keywords
correction
wave number
interference signal
fourier transform
data
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
Application number
CN202111595873.8A
Other languages
English (en)
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.)
Hangzhou Chunlai Technology Co Ltd
Original Assignee
Hangzhou Chunlai Technology Co Ltd
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 Hangzhou Chunlai Technology Co Ltd filed Critical Hangzhou Chunlai Technology Co Ltd
Priority to CN202111595873.8A priority Critical patent/CN114397262A/zh
Publication of CN114397262A publication Critical patent/CN114397262A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/3504Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/27Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
    • G01N21/274Calibration, base line adjustment, drift correction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N2021/3595Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using FTIR

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Pathology (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Immunology (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Spectrometry And Color Measurement (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明属于光谱校正技术领域,具体涉及一种傅里叶变换红外光谱仪波数漂移的校正方法及系统,其方法包括以下步骤:S1、开启校正光源,校正光源发出的红外光束经过校正气体室后进入光谱仪,光谱仪得到检测干涉信号;其中,校正气体室内充满氟化氢气体;S2、对检测干涉信号以及光谱仪预设的校正干涉信号的数据列分别进行补零处理,之后再进行傅里叶变换,得到检测吸收光谱和校正吸收光谱;S3、根据氟化氢气体的吸收峰对应的波数以及校正干涉信号的采样间距,分别确定检测吸收光谱和校正吸收光谱中吸收峰对应的第一波数和第二波数;S4、根据第一波数和第二波数对待校正的检测干涉信号进行波数漂移校正。本发明的校正精度高。

Description

一种傅里叶变换红外光谱仪波数漂移的校正方法及系统
技术领域
本发明属于光谱校正技术领域,具体涉及一种傅里叶变换红外光谱仪波数漂移的校正方法及系统。
背景技术
傅里叶变换红外光谱技术作为一种综合性探测技术,近年来得以快速发展。大气中大多数的微量、痕量气体都是红外活性气体,在2~30μm波段范围内具有吸收和发射红外特征光谱的能力,上述波段称为中红外区指纹波段,对于光谱测量非常有利。
傅里叶变换红外光谱技术利用分光束干涉原理,将迈克尔逊干涉仪、调制技术与计算机技术结合,通过傅里叶变换的方法,实现由干涉图到光谱的还原;再根据得到的光谱学反演出待测气体的浓度。该技术在空气污染物监测中应用前景非常广泛。
基于傅里叶变换红外光谱技术的开路式气体分析仪广泛应用于工业园区周边,能实现工业园区环境质量的在线监测。工业园区周边的污染物成分复杂,各成分间的吸收光谱存在交叉干扰。虽然可以通过智能识别算法去除污染物之间的交叉干扰,达到精确定量的效果;但是,智能识别算法对光谱仪精度的依赖性很高,如果光谱仪波数稍有漂移就会导致测量结果偏差巨大。
因此,对光谱仪波数漂移定期自动校正对保证数据质量有着重要的意义。
发明内容
基于现有技术中存在的上述缺点和不足,本发明的目的是提供满足前述需求的一种傅里叶变换红外光谱仪波数漂移的校正方法及系统。
为了达到上述发明目的,本发明采用以下技术方案:
一种傅里叶变换红外光谱仪波数漂移的校正方法,包括以下步骤:
S1、开启校正光源,校正光源发出的红外光束经过校正气体室后进入光谱仪,光谱仪得到检测干涉信号;其中,校正气体室内充满氟化氢气体;
S2、对检测干涉信号以及光谱仪预设的校正干涉信号的数据列分别进行补零处理,之后再进行傅里叶变换,得到检测吸收光谱和校正吸收光谱;
S3、根据氟化氢气体的吸收峰对应的波数以及校正干涉信号的采样间距,分别确定检测吸收光谱和校正吸收光谱中吸收峰对应的第一波数和第二波数;
S4、根据第一波数和第二波数对待校正的检测干涉信号进行波数漂移校正。
作为优选方案,所述步骤S2中,对检测干涉信号以及校正干涉信号的数据列分别进行补零处理,具体包括:对检测干涉信号X(n)和校正干涉信号Y(n)的数据列的两端分别填入m位0数据,得到数据列X(n+2m)和Y(n+2m);其中,m为正整数,n为采样点数;
之后再对数据列X(n+2m)和Y(n+2m)进行傅里叶变换,分别得到检测吸收光谱S(n+2m)和校正吸收光谱T(n+2m)。
作为优选方案,所述步骤S3中,氟化氢气体的吸收峰对应的波数设为4000cm-1,校正干涉信号的采样间距A nm,则4000cm-1对应光谱的第
Figure BDA0003431101780000021
个数据;其中,<·>为取整符号。
作为优选方案,所述步骤S3中,分别确定检测吸收光谱和校正吸收光谱中吸收峰对应的第一波数和第二波数,包括:
在4000cm-1对应光谱的第
Figure BDA0003431101780000022
个波数附近取透光率最小值对应的波数:
Figure BDA0003431101780000031
Figure BDA0003431101780000032
其中,p取值为2-20之间的整数,i为第一波数,j为第二波数。
作为优选方案,所述步骤S4,具体包括以下步骤:
S41、对待校正的检测干涉信号X’(n)进行傅里叶变换,得到数据列χ(n)=FFT[X’(n)];
S42、将数据列χ(n)分为两个数列,即χ(n)=[χ1(n/2),χ2(n/2)],并在中间补上k个0,得到数据列α(n+k)=[χ1(n/2),0,……,0,χ2(n/2)];k取值为正整数;
S43、基于数据列α(n+k)筛选得到波数漂移校正后的检测干涉信号数列B(<n*j/i>)=[α<(n+k)*j/(n*i)>,α<2(n+k)*j/(n*i)>,α<3(n+k)*j/(n*i)>,……,α<(n-2)*(n+k)*j/(n*i)>,α<(n-1)*(n+k)*j/(n*i)>,α<n*(n+k)*j/(n*i)>];其中,<·>为取整符号;
S44、对波数漂移校正后的检测干涉信号数列进行傅里叶变换得到波数漂移校正后的红外光谱。
本发明还提供一种傅里叶变换红外光谱仪波数漂移的校正系统,包括:
光源控制模块,用于控制校正光源的开启及关闭;
采集模块,用于采集校正光源发出的红外光束经过校正气体室后进入光谱仪并得到的检测干涉信号以及光谱仪预设的校正干涉信号;其中,校正气体室内充满氟化氢气体;
数据补零处理模块,用于对检测干涉信号以及光谱仪预设的校正干涉信号的数据列分别进行补零处理;
傅里叶变换模块,用于对补零处理之后的检测干涉信号和校正干涉信号的数据列进行傅里叶变换,得到检测吸收光谱和校正吸收光谱;
波数计算模块,用于根据氟化氢气体的吸收峰对应的波数以及校正干涉信号的采样间距,分别确定检测吸收光谱和校正吸收光谱中吸收峰对应的第一波数和第二波数;
校正模块,用于根据第一波数和第二波数对待校正的检测干涉信号进行波数漂移校正。
作为优选方案,所述数据补零处理模块对检测干涉信号以及光谱仪预设的校正干涉信号的数据列分别进行补零处理,包括:
对检测干涉信号X(n)和校正干涉信号Y(n)的数据列的两端分别填入m位0数据,得到数据列X(n+2m)和Y(n+2m);其中,m为正整数,n为采样点数;
相应地,所述傅里叶变换模块用于对数据列X(n+2m)和Y(n+2m)进行傅里叶变换,分别得到检测吸收光谱S(n+2m)和校正吸收光谱T(n+2m)。
作为优选方案,所述波数计算模块确定第一波数和第二波数的过程,包括:
氟化氢气体的吸收峰对应的波数设为4000cm-1,校正干涉信号的采样间距A nm,则4000cm-1对应光谱的第
Figure BDA0003431101780000041
个数据;其中,<·>为取整符号;
在4000cm-1对应光谱的第
Figure BDA0003431101780000042
个波数附近取透光率最小值对应的波数:
Figure BDA0003431101780000043
Figure BDA0003431101780000044
其中,p取值为2-20之间的整数,i为第一波数,j为第二波数。
作为优选方案,所述校正模块包括拆分模块和筛选模块;
根据第一波数和第二波数对待校正的检测干涉信号进行波数漂移校正的过程,包括:
首先,利用傅里叶变换模块对待校正的检测干涉信号X’(n)进行傅里叶变换,得到数据列χ(n)=FFT[X’(n)];
然后,利用拆分模块将数据列χ(n)分为两个数列,即χ(n)=[χ1(n/2),χ2(n/2)],并利用数据补零处理模块在中间补上k个0,得到数据列α(n+k)=[χ1(n/2),0,……,0,χ2(n/2)];k取值为正整数;
接着,利用筛选模块基于数据列α(n+k)筛选得到波数漂移校正后的检测干涉信号数列B(<n*j/i>)=[α<(n+k)*j/(n*i)>,α<2(n+k)*j/(n*i)>,α<3(n+k)*j/(n*i)>,……,α<(n-2)*(n+k)*j/(n*i)>,α<(n-1)*(n+k)*j/(n*i)>,α<n*(n+k)*j/(n*i)>];其中,<·>为取整符号;
最后,利用傅里叶变换模块对波数漂移校正后的检测干涉信号数列进行傅里叶变换得到波数漂移校正后的红外光谱。
本发明与现有技术相比,有益效果是:
(1)本发明利用氟化氢气体作为校正气体,基于氟化氢气体的吸收峰线宽较窄,且不与其他气体产生干扰,有利于提升校正的精度;
(2)对信号数据列进行补零处理,提高傅里叶变换后的分辨率,进一步提升后续校正的精度;
(3)根据检测吸收光谱和校正吸收光谱中吸收峰对应的第一波数和第二波数,并对应数据列α(n+k)中的数据匹配筛选得到波数漂移校正后的检测干涉信号数列,实现波数漂移校正,校正的精度高。
附图说明
图1是本发明实施例的傅里叶变换红外光谱仪波数漂移的校正方法的流程图;
图2是本发明实施例的傅里叶变换红外光谱仪波数漂移的校正系统的模块构架图。
具体实施方式
为了更清楚地说明本发明实施例,下面将对照附图说明本发明的具体实施方式。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图,并获得其他的实施方式。
如图1所示,本发明实施例的傅里叶变换红外光谱仪波数漂移的校正方法,包括以下步骤:
S1、开启校正光源,校正光源发出的红外光束经过校正气体室后进入光谱仪,光谱仪得到检测干涉信号;其中,校正气体室内充满氟化氢气体;
具体地,在开放光路设备自动监测过程中,定期进入光谱仪校正模式,此时切断外界光源,开启光谱仪内部的校正光源进行校正;
校正光源发出的红外光束,经过校正气体室后进入光谱仪,光谱仪检测到氟化氢气体的干涉信号。由于气体吸收峰位置是由气体分子振动频率决定的,不会发生改变,因此选用气体吸收峰作为光谱仪波数的基准。又因为氟化氢气体的吸收峰线宽较窄,且不与其他气体产生干扰,因此选择氟化氢作为校正气体。
S2、对检测干涉信号以及光谱仪预设的校正干涉信号的数据列分别进行补零处理,之后再进行傅里叶变换,得到检测吸收光谱和校正吸收光谱;
其中,对检测干涉信号以及校正干涉信号的数据列分别进行补零处理,提高傅里叶变换后的分辨率,具体过程包括:对检测干涉信号X(n)和校正干涉信号Y(n)的数据列的两端分别填入m位0数据,得到数据列X(n+2m)和Y(n+2m);其中,m为正整数,n为采样点数;
之后再对数据列X(n+2m)和Y(n+2m)进行傅里叶变换,分别得到检测吸收光谱S(n+2m)和校正吸收光谱T(n+2m):
S(n+2m)=FFT[X(n+2m)]
T(n+2m)=FFT[Y(n+2m)]。
S3、根据氟化氢气体的吸收峰对应的波数以及校正干涉信号的采样间距,分别确定检测吸收光谱和校正吸收光谱中吸收峰对应的第一波数和第二波数;
其中,氟化氢气体的吸收峰对应的波数设为4000cm-1,校正干涉信号的采样间距Anm,则4000cm-1对应光谱的第
Figure BDA0003431101780000071
个数据;其中,<·>为取整符号。
另外,分别确定检测吸收光谱和校正吸收光谱中吸收峰对应的第一波数和第二波数,具体包括:
在4000cm-1对应光谱的第
Figure BDA0003431101780000072
个波数附近取透光率最小值(即吸收峰)对应的波数:
Figure BDA0003431101780000073
Figure BDA0003431101780000074
其中,p取值为2-20之间的整数,例如,p取值为5;i为第一波数,j为第二波数。
i与j的差值为波数漂移数。
将i和j储存于光谱仪的内存中,退出校正模式,关闭光谱仪内部的校正光源,重新检测外界光源信号。
S4、获取待校正的检测干涉信号,并根据第一波数和第二波数对待校正的检测干涉信号进行波数漂移校正。
其中,在自动监测过程中获取待校正的检测干涉信号之后,需要对其进行重构,使其采样间距A’小于校正干涉信号的采样间距A,具体地,上述步骤S4包括以下步骤:
S41、对待校正的检测干涉信号X’(n)进行傅里叶变换,得到数据列χ(n)=FFT[X’(n)];
S42、将数据列χ(n)分为两个数列,即χ(n)=[χ1(n/2),χ2(n/2)],并在中间补上k个0,得到数据列α(n+k)=[χ1(n/2),0,……,0,χ2(n/2)];k取值为正整数;
S43、基于数据列α(n+k)筛选得到波数漂移校正后的检测干涉信号数列B(<n*j/i>)=[α<(n+k)*j/(n*i)>,α<2(n+k)*j/(n*i)>,α<3(n+k)*j/(n*i)>,……,α<(n-2)*(n+k)*j/(n*i)>,α<(n-1)*(n+k)*j/(n*i)>,α<n*(n+k)*j/(n*i)>];其中,<·>为取整符号;
S44、对波数漂移校正后的检测干涉信号数列进行傅里叶变换以及信号处理即可得到波数漂移校正后的红外光谱,实现波数漂移校正。
基于本发明实施例的上述校正方法,本发明实施例还提供一种傅里叶变换红外光谱仪波数漂移的校正系统,如图2所示,包括光源控制模块100、采集模块200、数据补零处理模块300、傅里叶变换模块400、波数计算模块500和校正模块600,其中,校正模块600还包括拆分模块6001和筛选模块6002。
具体地,光源控制模块100用于控制校正光源的开启及关闭,以便光谱仪进入校正模式以及退出校正模式。
采集模块200用于采集校正光源发出的红外光束经过校正气体室后进入光谱仪并得到的检测干涉信号以及光谱仪预设的校正干涉信号;其中,校正气体室内充满氟化氢气体;还用于采集自动监测过程中获取的待校正的检测干涉信号;
数据补零处理模块300用于对检测干涉信号以及光谱仪预设的校正干涉信号的数据列分别进行补零处理;
傅里叶变换模块400用于对补零处理之后的检测干涉信号和校正干涉信号的数据列进行傅里叶变换,得到检测吸收光谱和校正吸收光谱;
波数计算模块500用于根据氟化氢气体的吸收峰对应的波数以及校正干涉信号的采样间距,分别确定检测吸收光谱和校正吸收光谱中吸收峰对应的第一波数和第二波数;
校正模块600用于根据第一波数和第二波数对待校正的检测干涉信号进行波数漂移校正。
其中,数据补零处理模块对检测干涉信号以及光谱仪预设的校正干涉信号的数据列分别进行补零处理,包括:
对检测干涉信号X(n)和校正干涉信号Y(n)的数据列的两端分别填入m位0数据,得到数据列X(n+2m)和Y(n+2m);其中,m为正整数,n为采样点数;
相应地,傅里叶变换模块用于对数据列X(n+2m)和Y(n+2m)进行傅里叶变换,分别得到检测吸收光谱S(n+2m)和校正吸收光谱T(n+2m):
S(n+2m)=FFT[X(n+2m)]
T(n+2m)=FFT[Y(n+2m)]
上述波数计算模块确定第一波数和第二波数的过程,包括:
氟化氢气体的吸收峰对应的波数设为4000cm-1,校正干涉信号的采样间距A nm,则4000cm-1对应光谱的第
Figure BDA0003431101780000091
个数据;其中,<·>为取整符号;
在4000cm-1对应光谱的第
Figure BDA0003431101780000092
个波数附近取透光率最小值对应的波数:
Figure BDA0003431101780000093
Figure BDA0003431101780000094
其中,p取值为2-20之间的整数,例如,p取值为5;i为第一波数,j为第二波数。
校正模块根据第一波数和第二波数对待校正的检测干涉信号进行波数漂移校正的过程,包括:
首先,利用傅里叶变换模块400对待校正的检测干涉信号X’(n)进行傅里叶变换,得到数据列χ(n)=FFT[X’(n)];
然后,利用拆分模块6001将数据列χ(n)分为两个数列,即χ(n)=[χ1(n/2),χ2(n/2)],并利用数据补零处理模块300在中间补上k个0,得到数据列α(n+k)=[χ1(n/2),0,……,0,χ2(n/2)];k取值为正整数;
接着,利用筛选模块6002基于数据列α(n+k)筛选得到波数漂移校正后的检测干涉信号数列B(<n*j/i>)=[α<(n+k)*j/(n*i)>,α<2(n+k)*j/(n*i)>,α<3(n+k)*j/(n*i)>,……,α<(n-2)*(n+k)*j/(n*i)>,α<(n-1)*(n+k)*j/(n*i)>,α<n*(n+k)*j/(n*i)>];其中,<·>为取整符号;
最后,利用傅里叶变换模块400对波数漂移校正后的检测干涉信号数列进行傅里叶变换得到波数漂移校正后的红外光谱。
以上所述仅是对本发明的优选实施例及原理进行了详细说明,对本领域的普通技术人员而言,依据本发明提供的思想,在具体实施方式上会有改变之处,而这些改变也应视为本发明的保护范围。

Claims (9)

1.一种傅里叶变换红外光谱仪波数漂移的校正方法,其特征在于,包括以下步骤:
S1、开启校正光源,校正光源发出的红外光束经过校正气体室后进入光谱仪,光谱仪得到检测干涉信号;其中,校正气体室内充满氟化氢气体;
S2、对检测干涉信号以及光谱仪预设的校正干涉信号的数据列分别进行补零处理,之后再进行傅里叶变换,得到检测吸收光谱和校正吸收光谱;
S3、根据氟化氢气体的吸收峰对应的波数以及校正干涉信号的采样间距,分别确定检测吸收光谱和校正吸收光谱中吸收峰对应的第一波数和第二波数;
S4、根据第一波数和第二波数对待校正的检测干涉信号进行波数漂移校正。
2.根据权利要求1所述的一种傅里叶变换红外光谱仪波数漂移的校正方法,其特征在于,所述步骤S2中,对检测干涉信号以及校正干涉信号的数据列分别进行补零处理,具体包括:对检测干涉信号X(n)和校正干涉信号Y(n)的数据列的两端分别填入m位0数据,得到数据列X(n+2m)和Y(n+2m);其中,m为正整数,n为采样点数;
之后再对数据列X(n+2m)和Y(n+2m)进行傅里叶变换,分别得到检测吸收光谱S(n+2m)和校正吸收光谱T(n+2m)。
3.根据权利要求2所述的一种傅里叶变换红外光谱仪波数漂移的校正方法,其特征在于,所述步骤S3中,氟化氢气体的吸收峰对应的波数设为4000cm-1,校正干涉信号的采样间距A nm,则4000cm-1对应光谱的第
Figure FDA0003431101770000011
个数据;其中,<·>为取整符号。
4.根据权利要求3所述的一种傅里叶变换红外光谱仪波数漂移的校正方法,其特征在于,所述步骤S3中,分别确定检测吸收光谱和校正吸收光谱中吸收峰对应的第一波数和第二波数,包括:
在4000cm-1对应光谱的第
Figure FDA0003431101770000021
个波数附近取透光率最小值对应的波数:
Figure FDA0003431101770000022
Figure FDA0003431101770000023
其中,p取值为2-20之间的整数,i为第一波数,j为第二波数。
5.根据权利要求4所述的一种傅里叶变换红外光谱仪波数漂移的校正方法,其特征在于,所述步骤S4,具体包括以下步骤:
S41、对待校正的检测干涉信号X’(n)进行傅里叶变换,得到数据列χ(n)=FFT[X’(n)];
S42、将数据列χ(n)分为两个数列,即χ(n)=[χ1(n/2),χ2(n/2)],并在中间补上k个0,得到数据列α(n+k)=[χ1(n/2),0,……,0,χ2(n/2)];k取值为正整数;
S43、基于数据列α(n+k)筛选得到波数漂移校正后的检测干涉信号数列B(<n*j/i>)=[α<(n+k)*j/(n*i)>,α<2(n+k)*j/(n*i)>,α<3(n+k)*j/(n*i)>,……,α<(n-2)*(n+k)*j/(n*i)>,α<(n-1)*(n+k)*j/n*i>,α<n*(n+k)*j/(n*i)>];其中,<·>为取整符号;
S44、对波数漂移校正后的检测干涉信号数列进行傅里叶变换得到波数漂移校正后的红外光谱。
6.一种傅里叶变换红外光谱仪波数漂移的校正系统,其特征在于,包括:
光源控制模块,用于控制校正光源的开启及关闭;
采集模块,用于采集校正光源发出的红外光束经过校正气体室后进入光谱仪并得到的检测干涉信号以及光谱仪预设的校正干涉信号;其中,校正气体室内充满氟化氢气体;
数据补零处理模块,用于对检测干涉信号以及光谱仪预设的校正干涉信号的数据列分别进行补零处理;
傅里叶变换模块,用于对补零处理之后的检测干涉信号和校正干涉信号的数据列进行傅里叶变换,得到检测吸收光谱和校正吸收光谱;
波数计算模块,用于根据氟化氢气体的吸收峰对应的波数以及校正干涉信号的采样间距,分别确定检测吸收光谱和校正吸收光谱中吸收峰对应的第一波数和第二波数;
校正模块,用于根据第一波数和第二波数对待校正的检测干涉信号进行波数漂移校正。
7.根据权利要求6所述的一种傅里叶变换红外光谱仪波数漂移的校正系统,其特征在于,所述数据补零处理模块对检测干涉信号以及光谱仪预设的校正干涉信号的数据列分别进行补零处理,包括:
对检测干涉信号X(n)和校正干涉信号Y(n)的数据列的两端分别填入m位0数据,得到数据列X(n+2m)和Y(n+2m);其中,m为正整数,n为采样点数;
相应地,所述傅里叶变换模块用于对数据列X(n+2m)和Y(n+2m)进行傅里叶变换,分别得到检测吸收光谱S(n+2m)和校正吸收光谱T(n+2m)。
8.根据权利要求7所述的一种傅里叶变换红外光谱仪波数漂移的校正系统,其特征在于,所述波数计算模块确定第一波数和第二波数的过程,包括:
氟化氢气体的吸收峰对应的波数设为4000cm-1,校正干涉信号的采样间距A nm,则4000cm-1对应光谱的第
Figure FDA0003431101770000031
个数据;其中,<·>为取整符号;
在4000cm-1对应光谱的第
Figure FDA0003431101770000032
个波数附近取透光率最小值对应的波数:
Figure FDA0003431101770000041
Figure FDA0003431101770000042
其中,p取值为2-20之间的整数,i为第一波数,j为第二波数。
9.根据权利要求8所述的一种傅里叶变换红外光谱仪波数漂移的校正系统,其特征在于,所述校正模块包括拆分模块和筛选模块;
根据第一波数和第二波数对待校正的检测干涉信号进行波数漂移校正的过程,包括:
首先,利用傅里叶变换模块对待校正的检测干涉信号X’(n)进行傅里叶变换,得到数据列χ(n)=FFT[X’(n)];
然后,利用拆分模块将数据列χ(n)分为两个数列,即χ(n)=[χ1(n/2),χ2(n/2)],并利用数据补零处理模块在中间补上k个0,得到数据列α(n+k)=[χ1(n/2),0,……,0,χ2(n/2)];k取值为正整数;
接着,利用筛选模块基于数据列α(n+k)筛选得到波数漂移校正后的检测干涉信号数列B(<n*j/i>)=[α<(n+k)*j/(n*i)>,α<2(n+k)*j/(n*i)>,α<3(n+k)*j/(n*i)>,……,α<(n-2)*(n+k)*j/(n*i)>,α<(n-1)*(n+k)*j/(n*i)>,α<n*(n+k)*j/(n*i)>];其中,<·>为取整符号;
最后,利用傅里叶变换模块对波数漂移校正后的检测干涉信号数列进行傅里叶变换得到波数漂移校正后的红外光谱。
CN202111595873.8A 2021-12-24 2021-12-24 一种傅里叶变换红外光谱仪波数漂移的校正方法及系统 Pending CN114397262A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111595873.8A CN114397262A (zh) 2021-12-24 2021-12-24 一种傅里叶变换红外光谱仪波数漂移的校正方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111595873.8A CN114397262A (zh) 2021-12-24 2021-12-24 一种傅里叶变换红外光谱仪波数漂移的校正方法及系统

Publications (1)

Publication Number Publication Date
CN114397262A true CN114397262A (zh) 2022-04-26

Family

ID=81227541

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111595873.8A Pending CN114397262A (zh) 2021-12-24 2021-12-24 一种傅里叶变换红外光谱仪波数漂移的校正方法及系统

Country Status (1)

Country Link
CN (1) CN114397262A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115015120A (zh) * 2022-06-20 2022-09-06 厦门宇昊软件有限公司 一种傅里叶红外光谱仪及其温漂在线校正方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070145258A1 (en) * 2005-12-16 2007-06-28 Nelson Matthew P Method and apparatus for automated spectral calibration
JP2008197120A (ja) * 2008-05-19 2008-08-28 Seiko Epson Corp 赤外吸収分光器を用いた温室効果ガス測定方法
JP2014137328A (ja) * 2013-01-18 2014-07-28 Konica Minolta Inc フーリエ変換型分光計およびフーリエ変換型分光計の波長校正方法
US20190025125A1 (en) * 2017-07-20 2019-01-24 The Charles Stark Draper Laboratory, Inc. Grouped Molecular Absorption Line Wavelength Calibration Apparatus And Method
CN110017898A (zh) * 2019-03-29 2019-07-16 南京理工大学 一种基于干涉成像光谱仪的高精度波长定标方法
US20200393368A1 (en) * 2018-03-12 2020-12-17 Kanto Denka Kogyo Co., Ltd. Method and device for analyzing gas

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070145258A1 (en) * 2005-12-16 2007-06-28 Nelson Matthew P Method and apparatus for automated spectral calibration
JP2008197120A (ja) * 2008-05-19 2008-08-28 Seiko Epson Corp 赤外吸収分光器を用いた温室効果ガス測定方法
JP2014137328A (ja) * 2013-01-18 2014-07-28 Konica Minolta Inc フーリエ変換型分光計およびフーリエ変換型分光計の波長校正方法
US20190025125A1 (en) * 2017-07-20 2019-01-24 The Charles Stark Draper Laboratory, Inc. Grouped Molecular Absorption Line Wavelength Calibration Apparatus And Method
US20200393368A1 (en) * 2018-03-12 2020-12-17 Kanto Denka Kogyo Co., Ltd. Method and device for analyzing gas
CN110017898A (zh) * 2019-03-29 2019-07-16 南京理工大学 一种基于干涉成像光谱仪的高精度波长定标方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
邹曜璞 等: "傅里叶光谱仪高精度光谱定标研究", 《光谱学与光谱分析》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115015120A (zh) * 2022-06-20 2022-09-06 厦门宇昊软件有限公司 一种傅里叶红外光谱仪及其温漂在线校正方法
CN115015120B (zh) * 2022-06-20 2022-12-20 厦门宇昊软件有限公司 一种傅里叶红外光谱仪及其温漂在线校正方法

Similar Documents

Publication Publication Date Title
CN104990886B (zh) 一种傅里叶变换红外光谱仪波数漂移的校正方法
CN112326565A (zh) 紫外-可见光谱法检测水质cod时浊度影响的校正方法
Xie et al. Using FT-NIR spectroscopy technique to determine arginine content in fermented Cordyceps sinensis mycelium
CN114923873A (zh) 一种气体浓度检测方法、装置、电子设备及存储介质
CN114397262A (zh) 一种傅里叶变换红外光谱仪波数漂移的校正方法及系统
CN113607686B (zh) 一种基于修正tdlas背景信号漂移的气体浓度检测方法、系统及计算机可读存储介质
CN109520941B (zh) 在线光谱测量仪器的响应函数校正方法
CN113390825B (zh) 基于tdlas的时频域联合的气体浓度反演方法及装置
CN102495014A (zh) 一种差分吸收光谱测量中的光谱漂移修正方法
CN103558182A (zh) 一种用于激光气体在线分析仪确定气体浓度的方法
JP5424636B2 (ja) Ftir法を用いたガス分析装置及びこれに用いるプログラム
CN211347925U (zh) 一种气体浓度测量装置
CN112666104A (zh) 一种基于doas的气体浓度反演方法
CN115950827A (zh) 一种降低开放光路光谱检测信号噪声影响的方法
CN116183537A (zh) 基于差分消元的抗干扰ndir混合气体检测方法和系统
US20060238768A1 (en) Fourier transform infrared spectrometer
Zhu et al. Improved Savitzky-Golay filtering algorithm for measuring a pharmaceutical vial’s oxygen content based on wavelength modulation spectroscopy
Gao et al. A system for gas sensing employing correlation spectroscopy and wavelength modulation techniques with a multimode diode laser
JP3004747B2 (ja) フーリエ変換赤外分光計を用いた定量分析方法
CN110346346B (zh) 一种基于压缩感知关联算法的拉曼气体检测方法
CN114166798A (zh) 一种基于曲率的tdlas露点测量信号增强方法
CN108982401B (zh) 一种从混合气体的红外吸收光谱中解析单组分流量的方法
JP2649667B2 (ja) 分光分析における多成分分析方法
CN112729543B (zh) 一种光学检测中窄带吸收的背景校正方法、设备及介质
CN118130400A (zh) 一种基于解广义逆矩阵的doas系统气体浓度测量方法

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: 20220426