CN106841086A - 一种提高大气探测傅里叶光谱仪信噪比的方法 - Google Patents
一种提高大气探测傅里叶光谱仪信噪比的方法 Download PDFInfo
- Publication number
- CN106841086A CN106841086A CN201611045931.9A CN201611045931A CN106841086A CN 106841086 A CN106841086 A CN 106841086A CN 201611045931 A CN201611045931 A CN 201611045931A CN 106841086 A CN106841086 A CN 106841086A
- Authority
- CN
- China
- Prior art keywords
- signal
- data
- zero crossing
- sampled
- resampling
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000005070 sampling Methods 0.000 claims abstract description 36
- 238000012952 Resampling Methods 0.000 claims abstract description 30
- 238000001914 filtration Methods 0.000 claims abstract description 20
- 230000003287 optical effect Effects 0.000 claims abstract description 16
- 210000001367 artery Anatomy 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 210000003462 vein Anatomy 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 description 26
- 230000003595 spectral effect Effects 0.000 description 11
- 238000013461 design Methods 0.000 description 8
- 238000007796 conventional method Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000005855 radiation Effects 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 101100113692 Caenorhabditis elegans clk-2 gene Proteins 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000013139 quantization Methods 0.000 description 3
- 230000000630 rising effect Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 101100328957 Caenorhabditis elegans clk-1 gene Proteins 0.000 description 2
- 238000000862 absorption spectrum Methods 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3504—Investigating 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/45—Interferometric spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N2021/3595—Investigating 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)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
一种提高大气探测傅里叶红外光谱仪信噪比的方法，根据傅里叶红外光谱仪的工作原理，通过高频时钟对傅里叶光谱仪激光干涉信号的过零脉冲进行采样，并把每一个过零点脉冲周期细分成若干等时间间隔的子脉冲，利用细分后的子脉冲作为触发信号，实现对红外干涉信号的过采样，将采集到的红外干涉信号采样点的时间坐标信息变换成光程差坐标信息，并对过零点之间的采样数据进行等光程差插值重采样计算，对重采样后数据进行数字降采样滤波。本发明方法提高了数据采样频率并实现了数字滤波，不仅可以消除采样中的混叠噪声和量化噪声，也避免了模拟滤波元器件本身引入噪声，从而提高了光谱仪的信噪比。
Description
技术领域
本发明涉及一种提高光谱仪信噪比的方法。特别是对高光谱分辨率、干涉信号微弱的大气探测傅里叶光谱仪提高光谱信噪比的实现方法。
背景技术
大气探测傅里叶红外光谱仪是一种以迈克耳逊干涉仪为核心的时间调制型红外干涉仪，依靠动镜的机械摆扫运动来对光进行调制。傅里叶光谱仪只能直接得到目标物的干涉曲线，某一谱段上的能量是通过干涉图的傅里叶变换获得，所以每一个光谱波数上的信噪比和整个干涉图的信噪比相关。在光谱仪数据实际反演分析中，吸收系数的计算结果准确度依赖于探测信号对于各种影响因素的敏感程度，当影响吸收特性的信号被噪声淹没时，吸收系数的计算精度也就受到了限制，对理想情况而言，1％的反演精度需要100∶1的光谱信噪比(SNR)。因此干涉图信噪比是傅里叶光谱仪的一个重要指标，是决定系统辐射分辨率的重要因素，对数据恢复精度有着密切关系。
实现较高的探测光谱SNR一般从两个方面入手，一是增强探测目标强度，二是减小系统噪声水平。大气成份探测利用大气对太阳光的吸收光谱，直接采用太阳光作为系统的输入辐射源，因此很难再依靠辐射源来提高探测光谱SNR。从降低系统噪声水平的角度，傅里叶光谱仪在接收测量光谱信息的同时也接收测量了噪声信号。这些可能的噪声包括探测器噪声、辐射源光子噪声，还包括光源强度微小变化引起的噪声、外界振动干扰引入的噪声、电子线路引起的噪声、数字量化噪声等，这些变化的噪声信号会加到被测光谱信号中，使输出的光谱既包含被测光谱信号也包含噪声信号。
使用宽频率响应的探测器引入不必要的噪声会干扰正常的红外干涉信号.使傅里叶变换得到的光谱产生偏差，同时强的环境噪声还会“掩盖”光谱信号中的弱吸收峰，从而限制光谱分析的检出限。因此，必须对光谱转换前的红外干涉信号作滤波处理，尽可能消除噪声对真实信号的干扰。
传统的傅里叶光谱仪在降噪处理上采用探测器预放大之后加装高阶模拟滤波器的方法，例如采用多片MAX274芯片级连构建模拟滤波器，以阻止噪声信号进入数据采集模块，达到提高信噪比和降噪目的。
采用硬件实现高阶模拟滤波器时存在以下不足之处：
(1)受模拟滤波芯片数量和功能的限制及专用设计软件的限制，模拟滤波器设计复杂，在参数调整上非常困难；
(2)模拟滤波器元器件本身也会引入噪声；
(3)大气探测光谱仪的光谱带宽窄，系统的采样频率较低，使带外噪声容易混叠到探测光谱带内，增加了探测光谱带内的噪声，较低了信噪比；
(4)模拟滤波器无法滤除模数转换器采样引入的量化噪声；
发明内容
本发明的技术解决问题是：克服现有技术的不足，提出了一种可有效滤除噪声，提高傅里叶光谱仪信噪比的方法。
本发明的技术解决方案是：一种提高大气探测傅里叶光谱仪信噪比的方法，步骤如下：
(1)将周期为T0的高频时钟作为参考时钟信号，对傅里叶光谱仪激光干涉信号的过零点脉冲进行采样，记录过零点脉冲当前周期内参考时钟的整周期个数计数值为n2，过零点脉冲当前周期的前一周期内参考时钟的整周期个数计数值为n1。所述过零点脉冲周期即为傅里叶光谱仪激光干涉信号的过零点脉冲周期；
(2)根据步骤(1)的n1，将当前周期的傅里叶光谱仪激光干涉信号过零点脉冲细分成若干等时间间隔的子脉冲，本发明中细分的子脉冲个数优选为4；
(3)对进入傅里叶光谱仪的红外干涉信号进行采样，红外干涉信号的采样触发信号为步骤(2)中细分后的子脉冲，采样后得到红外干涉数据；
(4)将步骤(3)中采集到的红外干涉数据的采样点时间坐标信息变换成光程差坐标信息；
(5)对步骤(4)中的红外干涉数据以每相邻的过零点之间为一数据段，进行多项式插值计算，并按等光程差重采样，得到重采样数据；
(6)对步骤(5)中重采样数据进行降采样数字滤波，得到去噪声的红外干涉数据。
所述步骤(2)中将当前周期的傅里叶光谱仪激光干涉信号过零点脉冲细分成若干等时间间隔的子脉冲的方法为：将前一周期内过零点脉冲的参考时钟计数值n1当作当前过零点脉冲周期的时间间隔，则对当前过零点脉冲周期τ0等分(本专利中τ0＝4)，新的子脉冲周期为int为取整。
所述步骤(4)中红外干涉数据的采样点时间坐标信息变换成光程差坐标信息的方法为：转换公式为OPD＝8Rωt，其中，OPD为光程差坐标信息，R为发出激光干涉信号的傅里叶光谱仪的摆臂长，ω为发出激光干涉信号的傅里叶光谱仪的摆臂转动角速度，t为时间坐标信息。
所述步骤(5)中多项式插值重采样的方法为：把每一个周期过零点脉冲作为一个插值区间，对当前周期过零点脉冲τ0等分，则插值区间内共有τ0个红外干涉采样数据，加上插值区间下边界的一个红外干涉采样数据，共τ0+1个红外干涉采样数据作为本区间的插值点进行多项式插值，计算插值函数L(x)，并按照等光程差对插值函数L(x)进行重采样，得到重采样数据，重采样点数为τ1；
所述的步骤(6)中降采样数字滤波的方法为：将每一个周期过零点脉冲内按照等光程差重采样的τ1个数据进行求均值，此均值作为当前周期过零点脉冲的输出值，则把采样数据由τ1个降采样为1个。
本发明与现有技术相比的优点在于：
(1)本发明方法完全是在数字处理器内实现，处理的均为数字量的数据，避免了复杂的模拟电路设计，简化了系统的设计难度；
(2)本发明方法将傅里叶光谱仪激光干涉信号过零点脉冲细分，提高了光谱仪系统的采样数据率，增加了系统的采样带宽。根据采样定理模数转换的信号频带必须带宽限到采样频率的一半，否则便会产生混叠，然而实际应用中，红外干涉信号不是窄带信号，且宽带的加性噪声占据了高频区域，这些噪声分量也会使采样结果发生混叠，本发明方法通过提高采样频率，使带外噪声混叠到探测光谱频带以外，降低了探测光谱频带内的噪声，提高了系统的信噪比；
(3)本发明方法对红外干涉信号进行采样时，利用激光干涉信号过零点位置作为触发采样的采样点，即在每一个激光干涉信号过零点位置都会采集到红外干涉数据，避免了各过零点数据的相位漂移，提高探测光谱数据恢复精度；
(4)本发明方法的降采样滤波器有二个作用，一是滤除过采样后的随机噪声，减小量化噪声的功率；二是相对于最终的采样率它必须起到抗混叠滤波器的作用。本发明所采用算数平均滤波器，结构简单易于实现，滤波器可选择更大阶数，且滤波的过渡带可以很窄，不存在模拟滤波器过渡带过宽的弊端，并且不会引入额外的噪声；
(5)量化噪声是模拟信号到数字信号转换过程中不可消除的噪声，是系统中模数转换器固有的噪声，本发明方法是在红外干涉信号模数转换之后实现，可以滤除红外干涉信号量化所固有的随机噪声，提高系统的信噪比；
附图说明
图1为本发明一种提高大气探测傅里叶光谱仪信噪比的方法实现框图；
图2为本发明激光干涉过零点脉冲周期计数时序图；
图3为本发明激光干涉过零点脉冲周期细分流程框图；
图4的(a)为光谱仪采用传统方法设计采集到的系统噪声测试曲线；(b)为光谱仪采用本发明方法采集到的系统的噪声测试曲线
图5(a)为光谱仪采用传统方法设计，傅里叶变换后的光谱曲线，(b)为光谱仪采用本发明方法采集红外干涉数据经傅里叶变换得到的光谱曲线；
图6(a)为光谱仪采用传统方法设计得到的红外干涉信号光谱信噪比曲线，(b)为光谱仪采用本发明方法设计得到的红外干涉信号光谱信噪比曲线。
具体实施方式
本发明的基本思路为：提出一种提高大气探测傅里叶红外光谱仪信噪比的方法，根据傅里叶红外光谱仪的工作原理，通过高频时钟对傅里叶光谱仪激光干涉信号的过零脉冲进行采样，并把每一个过零点脉冲周期细分成若干等时间间隔的子脉冲，利用细分后的子脉冲作为触发信号，实现对红外干涉信号的过采样，将采集到的红外干涉信号采样点的时间坐标信息变换成光程差坐标信息，并对过零点之间的采样数据进行等光程差插值重采样计算，对重采样后数据进行数字降采样滤波。本发明方法提高了数据采样频率并实现了数字滤波，不仅可以消除采样中的混叠噪声和量化噪声，也避免了模拟滤波元器件本身引入噪声，从而提高了光谱仪的信噪比。
下面结合附图和具体实施例对本发明进行详细说明。
傅里叶光谱仪为了得到高精度的光谱，必须进行等光程差采样。通常采用He-Ne激光作为参考光源产生激光干涉信号，去除直流分量后激光干涉信号过零点触发采样得到均匀抽样干涉图。如图1所示，为本发明一种提高大气探测傅里叶光谱仪信噪比的方法的实现框图。首先利用高频时钟对傅里叶光谱仪激光干涉信号的过零点脉冲进行采样，并把每一个周期过零点脉冲细分成若干等时间间隔的子脉冲，利用细分后的子脉冲作为采样触发信号，实现对红外干涉信号的过采样，将采集到的红外干涉信号采样点的时间坐标信息变换成光程差坐标信息，并对过零点之间的采样数据进行等光程差插值重采样计算，对重采样后数据进行降采样数字滤波，实现提高红外干涉数据的信噪比。
1、过零点脉冲周期细分
如图2所示，以周期为T0的高频参考时钟Clk1作为测量参考时钟，在激光干涉信号的起始时刻(对应参考时钟的P点)和结束时刻(对应参考时钟的Q点)之间的时间间隔Ti可表示为：
Ti＝niT0+ta-tb (1)
式(1)中ni表示在第i个激光干涉信号过零点脉冲时间间隔Ti内参考时钟信号的整周期个数，在处理器中可由一个计数器来实现对ni的测量，即当激光干涉信号的过零点脉冲下降沿到来时打开计数器Counter1，在参考时钟Clk1上升沿触发Counter1计数加1，当下一个过零点脉冲下降沿到来时Counter1计数结束，保存Counter1计数结果并对其清零，同时计数器Counter2开始计数，计数器Counter1和Counter2二者交替计数实现对每一个周期过零点脉冲的测量。过零点脉冲下降沿到来的时刻是随机的，与参考时钟上升沿相差时间间隔分别为ta、tb，它的极限值是±T0，即Ti＝(ni±1)T0。把当前激光干涉信号过零点脉冲时间间隔计时定义为T2＝(n2±1)T0，前一次激光干涉信号过零点脉冲时间间隔计时定义为T1＝(n1±1)T0，因为光谱仪动镜的摆臂在连续两个过零点脉冲周期内速度不会突变，因此可以将前一次过零点脉冲的参考时钟计数值n1作为当前周期过零点脉冲被细分的时间间隔，而当前周期过零点脉冲的参考时钟计数值n2在计时完成后，作为下一周期过零点脉冲被细分的时间，依此类推。对每个过零点脉冲周期τ0等分，则新的子脉冲周期为int为取整。以周期为T0/2的高频时钟Clk2作为子脉冲产生时钟，由计数器Counter3来实现占空比为50％的子脉冲，如图3所示，实现方法如下(定义过零点脉冲为FPulse)：
(1)初始化计数器Counter3计数长度
(2)用Clk2时钟对激光干涉信号的过零点脉冲采样，检测到过零点脉冲下降沿时计数器Counter3清零，标志位Flag清零，子脉冲信号SPulse为高；
(3)在参考时钟Clk2上升沿触发Counter3计数加1；判断Counter3的计数值，当Counter3<A时，SPulse<＝SPulse AND(Not FPulse)，否则SPulse<＝(Not SPulse)AND(NotFPulse)，且Counter3计数值清零，标志位Flag加1；
(4)判断标志位Flag值，如果Flag小于2τ0则返回(3)，否则Counter3计数值保持不变，返回(2)；
2、红外干涉数据的采样点时间坐标信息变换成光程差坐标信息
由于大气探测光谱仪的双镜运动摆臂型式，加上端镜作用，使得整个光路的最大光程与机械行程之间形成了8：1的关系，当光谱仪摆臂转动θ角度时，则光程差为：
OPD＝8R sinθ(t) (2)
其中，OPD为光程差，R为光谱仪摆臂长。考虑到摆臂的摆动角度范围很小，sinθ(t)近似等于θ(t)，则可近似有如下关系：
OPD≈8Rθ(t)＝8Rωt (3)
其中，ω为角速度，t为时间。根据公式(3)可以把红外干涉数据采样点的时间信息转换成光程差信息。
3、多项式插值重采样
在利用激光干涉信号的过零点脉冲细分后的子脉冲对红外干涉信号采样时，是以时间刻度为基准的采样，由于光谱仪摆臂速度的抖动和过零点脉冲细分时误差的存在，会导致采样点的间隔并不是等光程差，因此需要将过零点周期内的采样数据按照等光程差重新采样，用多项式去逼近原始的红外干涉信号变化规律，再利用这一规律，计算插值区间内重采样的值，重新采样的方法如下。
把每一个激光干涉信号过零点采样周期作为一个插值区间，则插值区间内共有τ0个红外干涉信号的采样数据，加上插值区间下边界的一个采样数据，共τ0+1个数据作为本区间的插值节点进行多项式插值，则插值区间内的离散数据可表示为其中xi表示第i节点的光程差，yi表示第i节点的采样值即红外干涉数据幅值。则n阶拉格朗日插值多项式Ln是基函数l0,l1,...,ln的线性组合：
式中：
因本发明中τ0＝4，选择拉格朗日插值多项式的阶数n＝4。根据公式(4)和公式(5)，对于插值区间内任意光程差的位置，都可计算出其干涉数据的幅值。本发明中重采样的采样率不变，仍为每个插值区间有τ0个采样数据，但采样位置和原始数据发生变化，按照在插值区间内等光程差的间隔重新采样。
4、降采样数字滤波
对光谱仪采样系统而言，根据公式(4)进行的重采样是一种过采样处理，系统的随机噪声频带将被展开，需要用低通滤波器将信号的带外噪声滤除，才能提高系统的信噪比。同时，为减小数据量，低通滤波后的信号需要被降采样，降采样后数据采样频率下降，为使信号不发生混叠，需要用低通滤波器将二分之一采样频率外的噪声滤除，理论上，可将这两个低通滤波器合并成一个降采样低通滤波器。降采样数字滤波器的作用有两个：一是滤除过采样后的随机噪声，减小量化噪声的功率；二是相对于最终的采样率它必须起到抗混叠滤波器的作用。
低通滤波器的选择，在满足设计要求的前提下，需要从计算效率、实现难易程度以及系数存储量来考虑。算数平均滤波器隶属于FIR滤波器，是一种简单高效的低通滤波器，本发明采算数平均滤波器作为降采样滤波器。
由上述可知，根据公式(4)进行重采样后，每个插值区间内为τ0个采样数据，对插值区间内的数据进行数学平均作为其信号输出。其滤波公式定义如下：
其中yn为输出信号，xn+i为输入信号，N为滤波的点数，本发明中N＝τ0。公式(6)的滤波器实现的下抽取率与滤波器长度相等为N，即每个激光干涉信号的过零点采样周期只输出一个数据。本发明中N＝4，将4个连续的采样值累加，在完成累加后，将累加和右移2位即为本次滤波输出值，本发明中的滤波器并不需要乘法器来完成滤波运算，并且不需要存储器来保存滤波系数。
利用FPGA实现本发明的方法，并在某傅里叶光谱仪上进行测试。为了简单、客观的对本发明方法进行评价，这里采用实施本发明方法前后的系统噪声测试曲线、红外干涉信号的光谱曲线、红外干涉信号的光谱信噪比(SNR)三个指标进行比较。
图4给出傅里叶光谱仪在全黑条件下噪声，图4(a)为光谱仪采用传统方法设计采集到的系统噪声，噪声幅值主要集中在-8DN～8DN之间。图4(b)为光谱仪采用本发明方法采集到的系统的噪声，噪声幅值主要集中在-5DN～5DN之间，从图4(b)中可以看出，采用本发明方法系统的噪声得到明显的抑制。
图5给出了傅里叶光谱仪采集到的红外干涉信号经傅里叶变换后的光谱曲线，图5(a)为光谱仪采用传统方法设计，傅里叶变换后的光谱曲线，图(b)为光谱仪采用本发明方法采集红外干涉数据经傅里叶变换得到的光谱曲线，从两图对比可看出，虽然二者都能很好分辨出不同波数的吸收光谱，但图(b)的光谱噪声要低于图(a)(光谱平坦的部位幅值较窄)，由此可以看出系统的性能得到了改善。
图6为在相同辐照强度下，采用本发明方法前后红外干涉信号的光谱信噪比曲线，图6(a)为光谱仪采用传统方法设计，得到的红外干涉信号光谱信噪比曲线，图6(b)为光谱仪采用本发明方法设计，得到的红外干涉信号光谱信噪比曲线，由图6可以看出，在相同波数情况下，图6(b)中的光谱信噪比要较图6(a)中的光谱信噪比得到提高，说明采用本发明方法傅里叶光谱仪的信噪比得到提高。
本发明说明书中未作详细描述的内容属本领域专业技术人员的公知技术。
Claims (5)
1.一种提高大气探测傅里叶光谱仪信噪比的方法，其特征在于步骤如下：
(1)将周期为T0的高频时钟作为参考时钟信号，对傅里叶光谱仪激光干涉信号的过零点脉冲进行采样，记录过零点脉冲当前周期内参考时钟的整周期个数计数值为n2，过零点脉冲当前周期的前一周期内参考时钟的整周期个数计数值为n1，所述过零点脉冲周期即为傅里叶光谱仪激光干涉信号的过零点脉冲周期；
(2)根据步骤(1)的n1，将当前周期的傅里叶光谱仪激光干涉信号过零点脉冲细分成若干等时间间隔的子脉冲；
(3)对进入傅里叶光谱仪的红外干涉信号进行采样，红外干涉信号的采样触发信号为步骤(2)中细分后的子脉冲，采样后得到红外干涉数据；
(4)将步骤(3)中采集到的红外干涉数据的采样点时间坐标信息变换成光程差坐标信息；
(5)对步骤(4)中的红外干涉数据以每相邻的过零点之间为一数据段，进行多项式插值计算，并按等光程差重采样，得到重采样数据；
(6)对步骤(5)中重采样数据进行降采样数字滤波，得到去噪声的红外干涉数据。
2.根据权利要求1所述的一种提高大气探测傅里叶光谱仪信噪比的方法，其特征在于：所述步骤(2)中将当前周期的傅里叶光谱仪激光干涉信号过零点脉冲细分成若干等时间间隔的子脉冲的方法为：将前一周期内过零点脉冲的参考时钟计数值n1当作当前过零点脉冲周期的时间间隔，则对当前过零点脉冲周期τ0等分(本专利中τ0＝4)，新的子脉冲周期为int为取整。
3.根据权利要求1所述的一种提高大气探测傅里叶光谱仪信噪比的方法，其特征在于：所述步骤(4)中红外干涉数据的采样点时间坐标信息变换成光程差坐标信息的方法为：转换公式为OPD＝8Rωt，其中，OPD为光程差坐标信息，R为发出激光干涉信号的傅里叶光谱仪的摆臂长，ω为发出激光干涉信号的傅里叶光谱仪的摆臂转动角速度，t为时间坐标信息。
4.根据权利要求1所述的一种提高大气探测傅里叶光谱仪信噪比的方法，其特征在于：所述步骤(5)中对步骤(4)中的红外干涉数据以每相邻的过零点之间为一数据段，进行多项式插值计算，并按等光程差重采样，得到重采样数据方法为：把每一个周期过零点脉冲作为一个插值区间，对当前周期过零点脉冲τ0等分，则插值区间内共有τ0个红外干涉采样数据，加上插值区间下边界的一个红外干涉采样数据，共τ0+1个红外干涉采样数据作为本区间的插值点进行多项式插值，计算插值函数L(x)，并按照等光程差对插值函数L(x)进行重采样，得到重采样数据，重采样点数为τ1。
5.根据权利要求1所述的一种提高大气探测傅里叶光谱仪信噪比的方法，其特征在于：所述的步骤(6)中对步骤(5)中重采样数据进行降采样数字滤波，得到去噪声的红外干涉数据方法为：将每一个周期过零点脉冲内按照等光程差重采样的τ1个数据进行求均值，此均值作为当前周期过零点脉冲的输出值，则把采样数据由τ1个降采样为1个。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611045931.9A CN106841086B (zh) | 2016-11-22 | 2016-11-22 | 一种提高大气探测傅里叶光谱仪信噪比的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611045931.9A CN106841086B (zh) | 2016-11-22 | 2016-11-22 | 一种提高大气探测傅里叶光谱仪信噪比的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106841086A true CN106841086A (zh) | 2017-06-13 |
| CN106841086B CN106841086B (zh) | 2019-05-24 |
Family
ID=59145886
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611045931.9A Active CN106841086B (zh) | 2016-11-22 | 2016-11-22 | 一种提高大气探测傅里叶光谱仪信噪比的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106841086B (zh) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107870593A (zh) * | 2017-11-24 | 2018-04-03 | 中国航空工业集团公司洛阳电光设备研究所 | 一种高速大动态范围红外模拟信号采集电路 |
| CN107990983A (zh) * | 2017-11-21 | 2018-05-04 | 北京空间机电研究所 | 一种窄谱段干涉信号的数字欠采样方法 |
| CN108111759A (zh) * | 2017-12-23 | 2018-06-01 | 航天恒星科技有限公司 | 面向面阵ccd光电转换的仿真设计方法 |
| CN108449135A (zh) * | 2018-05-11 | 2018-08-24 | 华中科技大学 | 一种基于空间光干涉的光信噪比监测装置及监测方法 |
| CN110864808A (zh) * | 2019-11-12 | 2020-03-06 | 天津大学 | 一种基于高速采样的傅里叶变换光谱检测方法 |
| CN111398981A (zh) * | 2020-04-27 | 2020-07-10 | 武汉海达数云技术有限公司 | 一种圆光栅角度测量装置及测量方法、激光扫描仪 |
| CN112051442A (zh) * | 2020-08-05 | 2020-12-08 | 中电科仪器仪表有限公司 | 一种微波峰值功率测量中提高时间参数测量速度的方法 |
| CN113075152A (zh) * | 2021-03-26 | 2021-07-06 | 云南电网有限责任公司电力科学研究院 | 一种检测变压器油中溶解气体含量的红外光增强系统 |
| CN113899715A (zh) * | 2021-09-25 | 2022-01-07 | 杭州谱育科技发展有限公司 | 降低傅里叶变换红外光谱仪噪声的方法 |
| CN113899714A (zh) * | 2021-09-25 | 2022-01-07 | 杭州谱育科技发展有限公司 | 傅里叶变换红外光谱仪的数据处理方法 |
| CN114046727A (zh) * | 2021-11-05 | 2022-02-15 | 中国科学院合肥物质科学研究院 | 基于等效时钟法的干涉仪速度信息获取方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070263208A1 (en) * | 2006-01-10 | 2007-11-15 | The General Hospital Corporation | Systems and methods for generating data based on one or more spectrally-encoded endoscopy techniques |
| WO2008048087A1 (en) * | 2006-10-18 | 2008-04-24 | Erasmus University Medical Center | Method for identifying natural drug analogues |
| CN101527585A (zh) * | 2009-04-17 | 2009-09-09 | 哈尔滨工业大学深圳研究生院 | 一种实现电力系统数据自适应压缩的装置及方法 |
| CN102955949A (zh) * | 2011-08-25 | 2013-03-06 | 中国科学院西安光学精密机械研究所 | 基于干涉图匹配的目标识别方法 |
| CN203069274U (zh) * | 2013-01-22 | 2013-07-17 | 李剑平 | 一种激光干涉仪光程差定位系统 |
| CN104198424A (zh) * | 2014-08-29 | 2014-12-10 | 清华大学 | 一种提高傅里叶红外光谱仪偏低信噪比频段信噪比的方法 |
| CN104240242A (zh) * | 2014-09-05 | 2014-12-24 | 北京空间机电研究所 | 一种遥感图像最大信噪比测量方法 |
-
2016
- 2016-11-22 CN CN201611045931.9A patent/CN106841086B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070263208A1 (en) * | 2006-01-10 | 2007-11-15 | The General Hospital Corporation | Systems and methods for generating data based on one or more spectrally-encoded endoscopy techniques |
| WO2008048087A1 (en) * | 2006-10-18 | 2008-04-24 | Erasmus University Medical Center | Method for identifying natural drug analogues |
| CN101527585A (zh) * | 2009-04-17 | 2009-09-09 | 哈尔滨工业大学深圳研究生院 | 一种实现电力系统数据自适应压缩的装置及方法 |
| CN102955949A (zh) * | 2011-08-25 | 2013-03-06 | 中国科学院西安光学精密机械研究所 | 基于干涉图匹配的目标识别方法 |
| CN203069274U (zh) * | 2013-01-22 | 2013-07-17 | 李剑平 | 一种激光干涉仪光程差定位系统 |
| CN104198424A (zh) * | 2014-08-29 | 2014-12-10 | 清华大学 | 一种提高傅里叶红外光谱仪偏低信噪比频段信噪比的方法 |
| CN104240242A (zh) * | 2014-09-05 | 2014-12-24 | 北京空间机电研究所 | 一种遥感图像最大信噪比测量方法 |
Non-Patent Citations (1)
| Title |
|---|
| 邢廷等: "大气探测傅里叶变换光谱仪", 《遥感技术与应用》 * |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107990983A (zh) * | 2017-11-21 | 2018-05-04 | 北京空间机电研究所 | 一种窄谱段干涉信号的数字欠采样方法 |
| CN107990983B (zh) * | 2017-11-21 | 2020-05-12 | 北京空间机电研究所 | 一种窄谱段干涉信号的数字欠采样方法 |
| CN107870593A (zh) * | 2017-11-24 | 2018-04-03 | 中国航空工业集团公司洛阳电光设备研究所 | 一种高速大动态范围红外模拟信号采集电路 |
| CN108111759A (zh) * | 2017-12-23 | 2018-06-01 | 航天恒星科技有限公司 | 面向面阵ccd光电转换的仿真设计方法 |
| CN108449135A (zh) * | 2018-05-11 | 2018-08-24 | 华中科技大学 | 一种基于空间光干涉的光信噪比监测装置及监测方法 |
| CN110864808A (zh) * | 2019-11-12 | 2020-03-06 | 天津大学 | 一种基于高速采样的傅里叶变换光谱检测方法 |
| CN110864808B (zh) * | 2019-11-12 | 2021-07-27 | 天津大学 | 一种基于高速采样的傅里叶变换光谱检测方法 |
| CN111398981A (zh) * | 2020-04-27 | 2020-07-10 | 武汉海达数云技术有限公司 | 一种圆光栅角度测量装置及测量方法、激光扫描仪 |
| CN112051442B (zh) * | 2020-08-05 | 2023-08-25 | 中电科思仪科技股份有限公司 | 一种微波峰值功率测量中提高时间参数测量速度的方法 |
| CN112051442A (zh) * | 2020-08-05 | 2020-12-08 | 中电科仪器仪表有限公司 | 一种微波峰值功率测量中提高时间参数测量速度的方法 |
| CN113075152A (zh) * | 2021-03-26 | 2021-07-06 | 云南电网有限责任公司电力科学研究院 | 一种检测变压器油中溶解气体含量的红外光增强系统 |
| CN113899714A (zh) * | 2021-09-25 | 2022-01-07 | 杭州谱育科技发展有限公司 | 傅里叶变换红外光谱仪的数据处理方法 |
| CN113899715A (zh) * | 2021-09-25 | 2022-01-07 | 杭州谱育科技发展有限公司 | 降低傅里叶变换红外光谱仪噪声的方法 |
| CN113899715B (zh) * | 2021-09-25 | 2024-04-09 | 杭州谱育科技发展有限公司 | 降低傅里叶变换红外光谱仪噪声的方法 |
| CN113899714B (zh) * | 2021-09-25 | 2024-04-09 | 杭州谱育科技发展有限公司 | 傅里叶变换红外光谱仪的数据处理方法 |
| CN114046727A (zh) * | 2021-11-05 | 2022-02-15 | 中国科学院合肥物质科学研究院 | 基于等效时钟法的干涉仪速度信息获取方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106841086B (zh) | 2019-05-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106841086B (zh) | 一种提高大气探测傅里叶光谱仪信噪比的方法 | |
| CN104181577B (zh) | 一种基于全数字化技术的束流位置和相位测量系统及方法 | |
| CN105044456B (zh) | 一种基于正交子带的电网瞬时频率测量与跟踪方法 | |
| CN103196547A (zh) | 一种实现旋转机械振动信号同步阶比跟踪分析方法 | |
| CN107085144A (zh) | 一种快速测量电力系统谐波的方法 | |
| CN102109555B (zh) | 具有数字中频信号处理系统的频谱分析仪及其实现方法 | |
| CN106767745A (zh) | 一种光电传感器测角系统的信号处理方法 | |
| CN106772545B (zh) | 一种采用尖顶成形算法的数字脉冲幅度分析器 | |
| CN108181482A (zh) | 基于虚拟正弦波的实时低速检测装置 | |
| CN105180959B (zh) | 一种适用于腕式计步器的抗干扰计步方法 | |
| CN105067048B (zh) | 一种抗振动干扰的数字化涡街流量计的测量方法 | |
| CN102035554A (zh) | 电能计量和分析系统、方法和模数转换电路 | |
| CN102128981A (zh) | 脉冲信号周期测量方法 | |
| CN109214371A (zh) | 一种基于快速傅氏变换的正弦信号滤波系统 | |
| CN105021904B (zh) | 一种基于dds移相技术的快速相位噪声测量系统及测量方法 | |
| CN112781709A (zh) | 变速工况下设备振动信号早期故障分析和特征提取方法 | |
| CN109407501A (zh) | 一种基于相关信号处理的时间间隔测量方法 | |
| CN110837000B (zh) | 基于fpga的频率测量系统 | |
| CN110208601A (zh) | 基于fpga的瞬时测频方法及其数字接收机 | |
| KR101738822B1 (ko) | 임펄스 레이더를 이용한 타겟의 생체 정보 결정 장치 및 방법 | |
| CN203132688U (zh) | 一种实现旋转机械振动信号同步阶比跟踪分析的装置 | |
| US10924130B1 (en) | High resolution digital trigger detector | |
| CN112964931B (zh) | 基于双通道欠采样的非理想多阻尼谐波信号参数测量方法 | |
| CN102130682B (zh) | 一种鉴相电路 | |
| CN105759255B (zh) | 一种cic多相内插滤波超声相控阵波束延时方法 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |