CN107356335B - 一种基于数字光阑技术的高通量傅里叶变换光谱检测方法 - Google Patents
一种基于数字光阑技术的高通量傅里叶变换光谱检测方法 Download PDFInfo
- Publication number
- CN107356335B CN107356335B CN201710508320.1A CN201710508320A CN107356335B CN 107356335 B CN107356335 B CN 107356335B CN 201710508320 A CN201710508320 A CN 201710508320A CN 107356335 B CN107356335 B CN 107356335B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- throughput
- fourier
- trasform spectroscopy
- jstop
- 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.)
- Active
Links
- 238000004611 spectroscopical analysis Methods 0.000 title claims abstract description 45
- 238000005516 engineering process Methods 0.000 title claims abstract description 28
- 238000001514 detection method Methods 0.000 title claims abstract description 16
- 238000001228 spectrum Methods 0.000 claims abstract description 69
- 230000004907 flux Effects 0.000 claims abstract description 37
- 230000003595 spectral effect Effects 0.000 claims abstract description 29
- 230000008439 repair process Effects 0.000 claims abstract description 18
- 238000011084 recovery Methods 0.000 claims abstract description 10
- 238000000034 method Methods 0.000 claims abstract description 6
- 230000005012 migration Effects 0.000 claims description 8
- 238000013508 migration Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000007476 Maximum Likelihood Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 2
- 241000406668 Loxodonta cyclotis Species 0.000 claims 1
- 230000003287 optical effect Effects 0.000 description 6
- 230000009466 transformation Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000000862 absorption spectrum Methods 0.000 description 1
- 238000001831 conversion spectrum Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009615 fourier-transform spectroscopy Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Spectrometry And Color Measurement (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于数字光阑技术的高通量傅里叶变换光谱检测方法,所述方法包括以下步骤:以大口径JStop条件下的高通量傅里叶变换光谱为修复对象,通过数字光阑,对高通量傅里叶变换光谱进行数字形式的修复处理,以此改善光谱的加宽和偏移状态,获得具有较高光通量和光谱分辨率的复原光谱;以小口径JStop条件下的高分辨傅里叶变换光谱为参考对象,通过复原光谱与参考对象的比对,验证数字光阑技术的可靠性和准确性。本发明突破常规傅里叶变换光谱仪中由JStop决定的光通量和光谱分辨率之间的相互限制,大幅度提高傅里叶变换光谱的光通量,从而进一步提高信噪比。
Description
技术领域
本发明涉及傅里叶变换光谱检测领域,尤其涉及一种基于数字光阑技术的高通量傅里叶变换光谱检测方法。
背景技术
傅里叶变换光谱仪是对物质进行定性分析和定量检测的重要光谱检测工具,目前,被广泛地应用在众多领域,包括:食品业、制药业、化工产业、安全领域以及环保领域等,尤其适用于微弱信号的测量。相比于色散型光谱仪,傅里叶变换光谱仪具有高通量的优势。与色散型光谱仪不同,傅里叶变换光谱仪中没有设置狭缝,而是设置有Jacquinot光阑(JStop[1][2]),在光谱分辨能力相同的情况下,傅里叶变换光谱仪的光通量是色散型光谱仪的约60倍。因此,傅里叶变换光谱往往具有较高的信噪比。
然而,JStop的设置会不可避免地引起光谱的加宽和偏移。不同于理想点光源形成的完全准直光束,JStop具有一定的孔径,穿过JStop的光源光束经过准直镜后,具有一定的发散角。由于发散角的存在,导致光束中的不同光线在迈克尔逊干涉仪中的传播路径不同,形成的光程差不一致。正是因为光程差的不一致,产生了光谱的加宽和偏移现象。虽然JStop孔径越大,光通量越大,但是光谱的加宽和偏移也会越明显。因此,常规的傅里叶变换光谱仪中均存在由JStop决定的光通量和光谱分辨率相互限制的问题。目前,还没有有效的优化方法去突破光通量和光谱分辨率的相互权衡,通常的选择则是牺牲光通量去获得较高的光谱分辨率。因此,急需有效的技术手段在保持较高光谱分辨率的同时实现光通量的大幅度提高。
参考文献
[1]李志刚,齐文宗,李福田.软X射线-真空紫外傅里叶变换光谱仪及其光谱分辨率研究[J].光学精密工程,1998,6(5):84-90.
[2]Saptari V.Fourier transform spectroscopy instrumentationengineering[M].SPIE press,2004.
发明内容
本发明提供了一种基于数字光阑技术的高通量傅里叶变换光谱检测方法,本发明突破常规傅里叶变换光谱仪中由Jacquinot光阑(JStop)决定的光通量和光谱分辨率之间的相互限制,大幅度提高傅里叶变换光谱的光通量,从而进一步提高信噪比,详见下文描述:
一种基于数字光阑技术的高通量傅里叶变换光谱检测方法,所述方法包括以下步骤:
以大口径JStop条件下的高通量傅里叶变换光谱为修复对象,通过数字光阑,对高通量傅里叶变换光谱进行数字形式的修复处理,以此改善光谱的加宽和偏移状态,获得具有较高光通量和光谱分辨率的复原光谱;
以小口径JStop条件下的高分辨傅里叶变换光谱为参考对象,通过复原光谱与参考对象的比对,验证数字光阑技术的可靠性和准确性。
所述通过数字光阑,对高通量傅里叶变换光谱进行数字形式的修复处理具体为:
通过反演计算,对加宽、偏移明显的高通量傅里叶变换光谱进行处理,获得具有较高光通量和光谱分辨率的傅里叶变换光谱。
所述通过反演计算,对高通量傅里叶变换光谱进行处理具体为:
获取复原光谱p0:
建立一个投影算子H,表示光谱加宽和偏移的规则,该投影算子H是光强分布和波数偏移步长的函数,表示为H(a,Δν):
pr=Hp0+e
使用迭代算法,寻找最大似然参数,直到收敛:
其中,β表示调整系数,Pr为实际采集的光谱;n表示波数偏移的次数;ai表示JStop处的光强分布;v为光谱波数坐标;Δν表示波数偏移步长;e为噪声;k为迭代次数。
本发明提供的技术方案的有益效果是:
1、本发明在维持常规傅里叶变换光谱仪硬件组成和结构不变的基础上,无需增加额外的元器件,无需考虑傅里叶变换光谱仪的具体光学参数,引进了数字光阑技术;
2、本发明在大口径JStop形成的高通量光谱采集条件下,对采集的傅里叶变换光谱进行数字修复处理,改善光谱的加宽和偏移问题,获得具有较高光通量和光谱分辨率的复原光谱,从而突破光通量和光谱分辨率的相互权衡;
3、本发明在保持较高光谱分辨率的同时,实现光通量的大幅度提高。
附图说明
图1是本发明提供的一种基于数字光阑技术的高通量傅里叶变换光谱检测方法的流程图;
图2是本发明提供的一种基于数字光阑技术的高通量傅里叶变换光谱检测方法的原理依据,即JStop对光谱加宽和偏移的影响;
图中,常规傅里叶变换光谱仪主要由光源S、光阑JStop、准直透镜L1、分束器BS、固定反射镜M1、可调反射镜M2、会聚透镜L2和探测器PD组成。
为了直观地显示JStop带来的影响,M1’是M1通过分束器BS所成的像。在傅里叶变换光谱仪中,探测器PD探测到的是干涉图,只有经过快速傅里叶变换(FFT),才能得到光谱图。
图3是本发明提供的一种基于数字光阑技术的高通量傅里叶变换光谱检测方法的实际效果图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面对本发明实施方式作进一步地详细描述。
实施例1
一种基于数字光阑技术的高通量傅里叶变换光谱检测方法,参见图1和图2,该方法包括以下步骤:
101:以大口径JStop条件下的高通量傅里叶变换光谱为修复对象,通过数字光阑,对高通量傅里叶变换光谱进行数字形式的修复处理,以此改善光谱的加宽和偏移状态,获得具有较高光通量和光谱分辨率的复原光谱;
102:以小口径JStop条件下的高分辨傅里叶变换光谱为参考对象,通过复原光谱与参考对象的比对,验证数字光阑技术的可靠性和准确性。
其中,步骤101中的通过数字光阑,对高通量傅里叶变换光谱进行数字形式的修复处理具体为:
通过反演计算,对加宽、偏移明显的高通量傅里叶变换光谱进行处理,获得具有较高光通量和光谱分辨率的傅里叶变换光谱。
综上所述,本发明实施例通过上述步骤101-步骤102突破常规傅里叶变换光谱仪中由JStop决定的光通量和光谱分辨率之间的相互限制,大幅度提高傅里叶变换光谱的光通量,从而进一步提高信噪比。
实施例2
下面结合图1和图2对实施例1中的方案进行进一步地介绍,详见下文描述:
本发明实施例是基于常规的傅里叶变换光谱仪进行的,无需改变常规傅里叶变换光谱仪的硬件组成和结构,无需增加额外的元器件,并且,无需考虑傅里叶变换光谱仪的具体光学参数。只需要利用数字修复处理的手段去改善大口径JStop引起的光谱加宽和偏移,从而获得同时具有较高光通量和光谱分辨率的复原光谱。
当傅里叶变换光谱仪设置有大口径JStop的时候,可以采集到高通量的傅里叶变换光谱,但是,光谱的加宽和偏移现象明显,导致光谱分辨率的显著恶化。
当傅里叶变换光谱仪设置有小口径JStop的时候,可以采集到高分辨的傅里叶变换光谱,但是,光谱的光通量较小,导致光谱的信噪比较差。
为了解决由JStop口径大小决定的光通量和光谱分辨率相互限制的问题,以大口径JStop条件下的高通量傅里叶变换光谱为修复对象,利用数字光阑技术,对光谱进行数字形式的修复处理,改善光谱的加宽和偏移状态。
具体而言,上述利用数字光阑技术,对光谱进行数字形式的修复处理,即主要是指通过反演计算,获得同时具有较高光通量和光谱分辨率的复原光谱。
为了验证上述反演计算所得的复原光谱是可靠的、准确的,以小口径JStop条件下的高分辨傅里叶变换光谱为参考,通过复原光谱与参考光谱的比对,确定数字光阑技术的可靠性和准确性。
一种基于数字光阑技术的高通量傅里叶变换光谱检测方法的实施是依据JStop引起光谱加宽和偏移的现象以及数量关系进行的,如图2所示。
在常规的傅里叶变换光谱仪中,当光源S是理想的点光源时,光源光束可以被完全准直,不存在发散角,最终经过快速傅里叶变换(FFT)可以得到理想的光谱图P0,中心波数为ν0。然而,由于光源S总是有一定的尺寸,为了限制S的大小,设置了一定口径的JStop,此时,准直光束存在最大发散角α,满足关系:α≈tanα=d/2f1。
其中,d是JStop的口径,f1是准直透镜的焦距。由于发散角α的存在,进入干涉仪的光束中的各条光线的传播路径不一样,从而导致形成的光程差不一致。对应于发散角α的光线经过FFT后得到的光谱不再是P0,而是变为Pα。两者的形态完全一样,只是Pα的中心波数να变小,为
由于JStop的设置,实际采集的光谱Pr是理想光谱P0和一系列对应于不同发散角的偏移光谱的叠加,相对于理想光谱P0,实际光谱Pr出现了加宽和偏移的现象。当JStop口径越大,光通量越大,然而,光谱的加宽和偏移现象也越明显。
为了突破常规傅里叶变换光谱仪中由JStop决定的光通量和光谱分辨率之间的相互限制,引进了数字光阑技术。数字光阑技术是一种通过数值反演计算修复高通量傅里叶变换光谱中的加宽和偏移状态的数字修复手段。
在数字光阑技术中,以p0表示同时具有较高光通量和光谱分辨率的复原光谱,以pr表示大口径JStop条件下的高通量傅里叶变换光谱,此时,可以通过求解公式(1),获取复原光谱p0:
公式(1)中,ai表示JStop处的光强分布,Δν表示波数偏移步长,n表示波数偏移的次数,v为光谱波数坐标,e为噪声。JStop处的光强分布ai可以事先通过阵列式探测器测量到,波数偏移步长Δν通常等于傅里叶变换光谱仪的采样波数间隔,偏移次数n由光束的最大发散角决定。
在此基础上,建立一个投影算子H,表示光谱加宽和偏移的规则。该投影算子H仅仅是光强分布ai和波数偏移步长Δν的函数,可以表示为H(a,Δν),与傅里叶变换光谱仪的具体光学参数无关。此时,可以将公式(1)转换为:
pr=Hp0+e (2)
为了求解公式(2),使用迭代算法,根据公式(3),去寻找最大似然参数,直到收敛:
公式(3)中,β表示调整系数,k为迭代次数。
综上所述,本发明实施例提供的一种基于数字光阑技术的高通量傅里叶变换光谱检测方法,具有显著的优势,不仅可以突破光通量和光谱分辨率之间的相互限制,大幅度提高光通量,而且基于现有的傅里叶变换光谱仪,无需增加额外的元器件,无需考虑傅里叶变换光谱仪的具体光学参数,具有较高的可操作性和较好的通用性。
实施例3
下面结合图3、以及实验数据对实施例1和2中的方案进行可行性验证,详见下文描述:
本发明实施例在维持常规傅里叶变换光谱仪硬件组成和结构不变的基础上,利用数字光阑技术,对光谱进行数字形式的修复处理,改善光谱的加宽和偏移状态,从而突破光通量和光谱分辨率的相互权衡。
图3示出了本发明实施例提供的一种基于数字光阑技术的高通量傅里叶变换光谱检测方法的实际效果图,详述如下:
利用傅里叶变换光谱仪采集水蒸气的中红外吸收光谱,其中,水蒸气的压强约为5kPa,傅里叶变换光谱仪的极限分辨率为0.08cm-1。为了提高原始光谱的信噪比,减小噪声对数字光阑技术的影响,采样时进行了100次平均。
分别设置JStop口径为1mm和4mm。其中,以4mm JStop口径时采集的光谱为修复对象,其具有的光通量是1mm JStop口径时的约12倍,但是其光谱加宽和偏移比较明显。以1mmJStop口径时采集的光谱为参考光谱,其具有较高的光谱分辨率,光谱加宽和偏移可以忽略,但是其光通量小,信噪比差。
为了简化数字修复处理,以高斯分布表示JStop处的光强分布,采用Landweber迭代算法,寻找最大似然参数,从而得到同时具有较高光通量和光谱分辨率的复原光谱。
其中,上述迭代算法,以及寻找最大似然参数的过程为本领域技术人员所公知,本发明实施例对此不做赘述。
复原光谱具有较大的光通量,其信噪比明显好于参考光谱,约为参考光谱的3倍。同时,复原光谱与参考光谱中的主要谱峰可以很好地匹配,说明复原光谱是可靠的、准确的。
本发明实施例对各器件的型号除做特殊说明的以外,其他器件的型号不做限制,只要能完成上述功能的器件均可。
本领域技术人员可以理解附图只是一个优选实施例的示意图,上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种基于数字光阑技术的高通量傅里叶变换光谱检测方法,其特征在于,所述方法包括以下步骤:
以大口径JStop条件下的高通量傅里叶变换光谱为修复对象,通过数字光阑,对高通量傅里叶变换光谱进行数字形式的修复处理,以此改善光谱的加宽和偏移状态,获得具有较高光通量和光谱分辨率的复原光谱;
以小口径JStop条件下的高分辨傅里叶变换光谱为参考对象,通过复原光谱与参考对象的比对,验证数字光阑技术的可靠性和准确性;
其中,所述通过数字光阑,对高通量傅里叶变换光谱进行数字形式的修复处理具体为:
通过反演计算,对加宽、偏移明显的高通量傅里叶变换光谱进行处理,获得同时具有较高光通量和光谱分辨率的傅里叶变换光谱;
进一步地,所述通过反演计算,对加宽、偏移明显的高通量傅里叶变换光谱进行处理具体为:
获取复原光谱p0:
建立一个投影算子H,表示光谱加宽和偏移的规则,该投影算子H是光强分布和波数偏移步长的函数,表示为H(a,Δν):
pr=Hp0+e
使用迭代算法,寻找最大似然参数,直到收敛:
其中,β表示调整系数,Pr为实际采集的光谱;n表示波数偏移的次数;ai表示JStop处的光强分布;v为光谱波数坐标;Δν表示波数偏移步长;e为噪声;k为迭代次数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710508320.1A CN107356335B (zh) | 2017-06-28 | 2017-06-28 | 一种基于数字光阑技术的高通量傅里叶变换光谱检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710508320.1A CN107356335B (zh) | 2017-06-28 | 2017-06-28 | 一种基于数字光阑技术的高通量傅里叶变换光谱检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107356335A CN107356335A (zh) | 2017-11-17 |
CN107356335B true CN107356335B (zh) | 2018-09-11 |
Family
ID=60272507
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710508320.1A Active CN107356335B (zh) | 2017-06-28 | 2017-06-28 | 一种基于数字光阑技术的高通量傅里叶变换光谱检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107356335B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019014522A1 (en) * | 2017-07-14 | 2019-01-17 | Thermo Electron Scientific Instruments Llc | ADVANCED INFRARED SPECTROSCOPY REFERENCE DETECTOR |
CN111076816B (zh) * | 2019-12-27 | 2021-01-15 | 中国科学院西安光学精密机械研究所 | 大孔径静态干涉成像光谱仪全视场光谱定标误差校正方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6940599B1 (en) * | 2002-02-08 | 2005-09-06 | Southwest Sciences Incorporated | Envelope functions for modulation spectroscopy |
CN1789933A (zh) * | 2005-11-10 | 2006-06-21 | 中国科学院上海技术物理研究所 | 红外静止型高光通量傅立叶变换成像光谱仪 |
JP5345785B2 (ja) * | 2008-01-08 | 2013-11-20 | Sumco Techxiv株式会社 | 分光吸収測定方法及び分光吸収測定装置 |
CN101782431A (zh) * | 2009-01-15 | 2010-07-21 | 中国科学院西安光学精密机械研究所 | 高分辨率反射式高速转镜干涉光谱仪 |
CN102095498B (zh) * | 2010-11-12 | 2012-10-31 | 北京工业大学 | 一种扫描式高精度傅立叶变换测量光谱的方法 |
-
2017
- 2017-06-28 CN CN201710508320.1A patent/CN107356335B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107356335A (zh) | 2017-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102297721B (zh) | 斜入射宽带偏振光谱仪和光学测量系统 | |
CN103954360B (zh) | 一种基于偏振阵列的光谱偏振装置及探测方法 | |
CN104634280B (zh) | 通用水平转台绝对角度和旋转角度的测量方法 | |
CN110709204B (zh) | 小样品曲率半径和热膨胀实时测量系统和方法 | |
CN107356335B (zh) | 一种基于数字光阑技术的高通量傅里叶变换光谱检测方法 | |
US10473525B2 (en) | Spatially resolved optical emission spectroscopy (OES) in plasma processing | |
CN105066910A (zh) | 电光晶体z轴偏离角测量装置及测量方法 | |
KR20160055673A (ko) | 나노구조의 임계 치수에 대한 광학 측정 시스템 및 측정 방법 | |
CN110494722A (zh) | 从空间输出中恢复光谱形状 | |
CN110864808B (zh) | 一种基于高速采样的傅里叶变换光谱检测方法 | |
WO2018146456A1 (en) | Compact interferometer | |
CN110567883B (zh) | 实时测量Mueller矩阵光谱的方法 | |
CN107462525A (zh) | 一种线性渐变滤光片的光谱特征参数测试方法 | |
CN104502067A (zh) | 一种傅里叶变换平面光栅相对衍射效率测试仪 | |
CN104977084B (zh) | 一种提高aotf成像空间分辨率和光谱分辨率的方法 | |
US8725447B2 (en) | Method for correcting a wave front analyser and analyser implementing said method | |
CN110824865B (zh) | 对准误差测量方法及装置 | |
CN110231098B (zh) | 超短激光脉冲多路延时同步测试方法 | |
JP4543213B2 (ja) | 分子配向温度計 | |
KR102600883B1 (ko) | 플라즈마 프로세싱에서의 공간적으로 분해된 광학적 방출 스펙트로스코피(oes) | |
CN115667841A (zh) | 用于测量距离的方法和装置 | |
WO2021110949A1 (en) | Apparatus for carrying out polarization resolved raman spectroscopy | |
US20150355023A1 (en) | System and Method of Alignment for Balanced Detection in a Spectral Domain Optical Coherence Tomography | |
Næss et al. | A fast map-making preconditioner for regular scanning patterns | |
CN108106732B (zh) | 拉曼光谱仪波数分辨率的标定方法及装置 |
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 |