CN110017898A - 一种基于干涉成像光谱仪的高精度波长定标方法 - Google Patents
一种基于干涉成像光谱仪的高精度波长定标方法 Download PDFInfo
- Publication number
- CN110017898A CN110017898A CN201910248039.8A CN201910248039A CN110017898A CN 110017898 A CN110017898 A CN 110017898A CN 201910248039 A CN201910248039 A CN 201910248039A CN 110017898 A CN110017898 A CN 110017898A
- Authority
- CN
- China
- Prior art keywords
- signal
- interference
- imaging spectrometer
- calibrating method
- interference fringe
- 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 29
- 238000003384 imaging method Methods 0.000 title claims abstract description 19
- 238000001228 spectrum Methods 0.000 claims abstract description 22
- 230000003595 spectral effect Effects 0.000 claims abstract description 12
- 239000000284 extract Substances 0.000 claims abstract description 10
- 230000009466 transformation Effects 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011088 calibration curve Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000003716 rejuvenation Effects 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/2823—Imaging spectrometer
-
- 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)
Abstract
本发明公开了一种基于干涉成像光谱仪的高精度波长定标方法,包括以下内容:以某一波长的单色光作为光源,采集干涉条纹图;从干涉条纹图中提取某一行干涉条纹信息构成干涉信号I;对干涉信号I进行扩展获得信号长度为I信号长度M倍的信号Ie;获取信号Ie的频谱强度信息FIe;提取频谱强度信息FIe中峰值位置的横坐标Ke,并由Ke和M获取波数位置坐标K;以若干不同波长的单色光分别作为光源,重复前述过程;对所有波长与其对应的波数位置坐标的关系数据进行曲线拟合,由此获得波长定标结果。本发明能够在不增加太多计算复杂度的情况下显著提高波长定标中谱线位置的定位精度。
Description
技术领域
本发明属于干涉光谱信号处理领域,特别是一种基于干涉成像光谱仪的高精度波长定标方法。
背景技术
干涉光谱技术为:利用干涉信号与光谱信号之间的物理关系,通过对干涉信号进行傅立叶变换,即可复原目标的光谱强度信息。实际应用中,通过傅里叶变换获得的光谱曲线的横坐标是波数位置,需要将波数位置的值转换为波长的值,这样才能获得光谱强度与波长的关系曲线,才是完整的光谱测量的过程。所以需要对波数位置与波长的关系进行定标,也就是波长定标过程。
传统的波长定标过程为:首先采集单色光通过干涉光谱仪的干涉条纹图,再提取其中一行构成干涉信号,对该信号进行傅里叶变换,提取频谱中的峰值所在位置的横坐标,即为波数位置坐标。之后采用不同波长的单色光对干涉光谱仪进行照明,重复以上操作。最后对不同波长和谱线位置坐标的关系进行拟合,得到谱线位置定标曲线,完成波长定标。
传统的波长定标方法受常规傅里叶变换的限制,只能获得整数坐标下的频域波数坐标,定标精度低,对后续的辐射定标过程也有影响。
发明内容
本发明的目的在于提供一种解决传统波长定标方法精度低问题的基于干涉成像光谱仪的高精度波长定标方法。
实现本发明目的的技术解决方案为:一种基于干涉成像光谱仪的高精度波长定标方法,包括以下步骤:
步骤1、以单色光作为光源,采集其通过干涉成像光谱仪的干涉条纹图;
步骤2、从所述干涉条纹图中提取某一行干涉条纹信息构成干涉信号I;
步骤3、对干涉信号I进行扩展获得新信号Ie,所述新信号Ie的长度为干涉信号I长度的M倍;
步骤4、获取新信号Ie的频谱强度信息FIe;
步骤5、提取频谱强度信息FIe中峰值位置的横坐标Ke,并由Ke和M获取波数位置坐标K;
步骤6、以不同波长的单色光分别作为光源,重复上述步骤1~步骤5;
步骤7、对所有波长与其对应的波数位置坐标的关系数据进行曲线拟合,由此获得波长定标结果。
本发明与现有技术相比,其显著优点为:1)通过时域补零延拓实现频域插值的效果,相比于频域插值计算方法更简单,速度更快;2)根据系统的需要进行补零插值,可以获得更高精度的波数坐标值,以此进行曲线拟合可以实现亚纳米波长定标。
下面结合附图对本发明作进一步详细描述。
附图说明
图1为本发明基于干涉成像光谱仪的高精度波长定标方法的流程图。
图2为本发明实施例从干涉条纹图中提取的某一行干涉条纹信息构成干涉信号示意图。
图3为对图2放大获得的局部干涉信号示意图。
图4为本发明实施例中扩展后干涉信号的频谱强度信息示意图。
图5为本发明实施例中采用传统波长定标方法获得的原始干涉信号的频谱强度信息示意图。
具体实施方式
结合图1,本发明基于干涉成像光谱仪的高精度波长定标方法,包括以下步骤:
步骤1、以单色光作为光源,采集其通过干涉成像光谱仪的干涉条纹图;
步骤2、从所述干涉条纹图中提取某一行干涉条纹信息构成干涉信号I,I为:
式中,I(x)为采样点x处的干涉信号强度,x的取值为0至N-1,其中N为采样点数。
步骤3、对干涉信号I进行扩展获得新信号Ie,所述新信号Ie的长度为干涉信号I长度的M倍;
步骤4、获取新信号Ie的频谱强度信息FIe;
步骤5、提取频谱强度信息FIe中峰值位置的横坐标Ke,并由Ke和M获取波数位置坐标K;
步骤6、以不同波长的单色光分别作为光源,重复上述步骤1~步骤5;
步骤7、对所有波长与其对应的波数位置坐标的关系数据进行曲线拟合,由此获得波长定标结果。
进一步优选地,步骤3通过添加零值对干涉信号I进行扩展,获得长度为M*N的新信号Ie为:
进一步地,步骤4获取新信号Ie的频谱强度信息FIe,具体为:
对新信号Ie进行傅里叶变换并取模,即获取新信号Ie的频谱强度信息FIe,所用公式为:
FIe=abs(FT(Ie))。
进一步地,步骤5由Ke和M获取波数位置坐标K,所用公式为:
K=Ke/M。
进一步优选地,步骤7对所有波长与其对应的谱线位置坐标的关系数据进行曲线拟合,具体为:采用最小二乘法对所有波长与其对应的谱线位置坐标的关系数据进行曲线拟合。
下面结合实施例对本发明作进一步详细的描述。
实施例
本发明基于干涉成像光谱仪的高精度波长定标方法,包括以下内容:
1、以632.8nm波长的单色光作为光源,采集其通过干涉成像光谱仪的干涉条纹图。
2、从干涉条纹图中提取某一行干涉条纹信息构成干涉信号I如图2所示,由于周期较多,为了更清楚地说明问题,对干涉信号进行局部放大,如图3所示。
3、对干涉信号I进行扩展获得新信号Ie,新信号Ie的长度为干涉信号I长度的M倍,本实施例中M取值为10000。
4、获取新信号Ie的频谱强度信息FIe,如图4所示。
5、如图4所示,提取频谱强度信息FIe中的峰值Ke=1202613,并由Ke和M获取波数位置坐标K=Ke/M=120.2613。
而传统波长定标方法获得的频谱图如图5所示,可知,该方法得到的波数位置坐标值为121,与本发明对比可知,本发明的方法大大提高了波数位置定位精度。
由于采用本方法获得的波数位置坐标更精确,所以对其他波长的激光执行以上操作并进行曲线拟合,最终获得的波长定标结果也会更加精确。
在相同的干涉数据情况下,本发明能够在不增加太多计算复杂度的情况下显著提高波长定标中谱线位置的定位精度,为干涉光谱的精确测量提供了理论支持和技术保障,为干涉光谱技术向着超分辨的方向发展提供了一种可能。
Claims (6)
1.一种基于干涉成像光谱仪的高精度波长定标方法,其特征在于,包括以下步骤:
步骤1、以单色光作为光源,采集其通过干涉成像光谱仪的干涉条纹图;
步骤2、从所述干涉条纹图中提取某一行干涉条纹信息构成干涉信号I;
步骤3、对干涉信号I进行扩展获得新信号Ie,所述新信号Ie的长度为干涉信号I长度的M倍;
步骤4、获取新信号Ie的频谱强度信息FIe;
步骤5、提取频谱强度信息FIe中峰值位置的横坐标Ke,并由Ke和M获取波数位置坐标K;
步骤6、以不同波长的单色光分别作为光源,重复上述步骤1~步骤5;
步骤7、对所有波长与其对应的波数位置坐标的关系数据进行曲线拟合,由此获得波长定标结果。
2.根据权利要求1所述的基于干涉成像光谱仪的高精度波长定标方法,其特征在于,步骤2所述从干涉条纹图中提取某一行干涉条纹信息构成干涉信号I,具体为:
式中,I(x)为采样点x处的干涉信号强度,x的取值为0至N-1,其中N为采样点数。
3.根据权利要求1或2所述的基于干涉成像光谱仪的高精度波长定标方法,其特征在于,步骤3所述对干涉信号I进行扩展获得新信号Ie,具体为:通过添加零值对干涉信号I进行扩展,获得长度为M*N的新信号Ie为:
4.根据权利要求3所述的基于干涉成像光谱仪的高精度波长定标方法,其特征在于,步骤4所述获取新信号Ie的频谱强度信息FIe,具体为:
对新信号Ie进行傅里叶变换并取模,即获取新信号Ie的频谱强度信息FIe,所用公式为:
FIe=abs(FT(Ie))。
5.根据权利要求4所述的基于干涉成像光谱仪的高精度波长定标方法,其特征在于,步骤5所述由Ke和M获取波数位置坐标K,所用公式为:
K=Ke/M。
6.根据权利要求5所述的基于干涉成像光谱仪的高精度波长定标方法,其特征在于,步骤7所述对所有波长与其对应的谱线位置坐标的关系数据进行曲线拟合,具体为:采用最小二乘法对所有波长与其对应的谱线位置坐标的关系数据进行曲线拟合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910248039.8A CN110017898B (zh) | 2019-03-29 | 2019-03-29 | 一种基于干涉成像光谱仪的高精度波长定标方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910248039.8A CN110017898B (zh) | 2019-03-29 | 2019-03-29 | 一种基于干涉成像光谱仪的高精度波长定标方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110017898A true CN110017898A (zh) | 2019-07-16 |
CN110017898B CN110017898B (zh) | 2021-01-08 |
Family
ID=67190198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910248039.8A Active CN110017898B (zh) | 2019-03-29 | 2019-03-29 | 一种基于干涉成像光谱仪的高精度波长定标方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110017898B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114397262A (zh) * | 2021-12-24 | 2022-04-26 | 杭州春来科技有限公司 | 一种傅里叶变换红外光谱仪波数漂移的校正方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2619245Y (zh) * | 2003-05-12 | 2004-06-02 | 中国科学院西安光学精密机械研究所 | 一种用于空间调制干涉光谱成像仪中的星上定标系统 |
CN1635347A (zh) * | 2003-12-31 | 2005-07-06 | 中国科学院西安光学精密机械研究所 | 干涉型超光谱成像仪星上一次像面定标方法 |
CN2718553Y (zh) * | 2003-12-25 | 2005-08-17 | 中国科学院西安光学精密机械研究所 | Ccd像元光谱响应均匀性定标装置 |
CN101979971A (zh) * | 2010-09-10 | 2011-02-23 | 中国科学院安徽光学精密机械研究所 | 一种辐亮度定标方法及激光光源、黑体辐射源定标系统 |
CN102818630A (zh) * | 2012-07-27 | 2012-12-12 | 中国资源卫星应用中心 | 一种干涉型成像光谱仪的光谱定标方法 |
-
2019
- 2019-03-29 CN CN201910248039.8A patent/CN110017898B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2619245Y (zh) * | 2003-05-12 | 2004-06-02 | 中国科学院西安光学精密机械研究所 | 一种用于空间调制干涉光谱成像仪中的星上定标系统 |
CN2718553Y (zh) * | 2003-12-25 | 2005-08-17 | 中国科学院西安光学精密机械研究所 | Ccd像元光谱响应均匀性定标装置 |
CN1635347A (zh) * | 2003-12-31 | 2005-07-06 | 中国科学院西安光学精密机械研究所 | 干涉型超光谱成像仪星上一次像面定标方法 |
CN101979971A (zh) * | 2010-09-10 | 2011-02-23 | 中国科学院安徽光学精密机械研究所 | 一种辐亮度定标方法及激光光源、黑体辐射源定标系统 |
CN102818630A (zh) * | 2012-07-27 | 2012-12-12 | 中国资源卫星应用中心 | 一种干涉型成像光谱仪的光谱定标方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114397262A (zh) * | 2021-12-24 | 2022-04-26 | 杭州春来科技有限公司 | 一种傅里叶变换红外光谱仪波数漂移的校正方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110017898B (zh) | 2021-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109186767A (zh) | 一种软件触发的傅里叶变换光谱仪中干涉图采样方法 | |
JP2007516760A5 (zh) | ||
CN112051583B (zh) | Fmcw距离测量系统中拍频信号非线性校正方法 | |
CN105222892B (zh) | 一种紫外光谱仪的波长定标方法 | |
CN112415078A (zh) | 一种质谱数据谱图信号校准方法及装置 | |
CN110674835A (zh) | 一种太赫兹成像方法及系统和一种无损检测方法及系统 | |
CN110987181B (zh) | 窗扫型干涉高光谱成像系统的光程差在线定标方法 | |
CN110017898A (zh) | 一种基于干涉成像光谱仪的高精度波长定标方法 | |
Plankey et al. | Hadamard spectrometer for ultraviolet-visible spectrometry | |
WO2022233110A1 (zh) | 一种阶梯叠加式傅里叶变换微分方法 | |
CN112184785B (zh) | 基于mcd度量和vtm的多模态遥感图像配准方法 | |
CN104268896A (zh) | 基于光谱抽样直方图的超光谱降维匹配方法及系统 | |
CN102818630A (zh) | 一种干涉型成像光谱仪的光谱定标方法 | |
CN110864808A (zh) | 一种基于高速采样的傅里叶变换光谱检测方法 | |
Liu et al. | SCBSS signal de-noising method of integrating EEMD and ESMD for dynamic deflection of bridges using GBSAR | |
CN109655158B (zh) | 基于大气廓线与led的高光谱遥感器在轨光谱定标的方法 | |
Iwamoto et al. | Kinematic study of root elongation in Arabidopsis thaliana with a novel image-analysis program | |
CN110210081A (zh) | 一种SS-OCT系统k-clock延时校正算法 | |
CN109186761A (zh) | 一种编码光谱成像系统的定标方法及编码模板 | |
CN110836878B (zh) | 一种卷积插值耦合高斯混合模型快速三维荧光寻峰方法 | |
CN111126452B (zh) | 一种基于主成分分析的地物光谱曲线扩展方法及系统 | |
CN116701920B (zh) | 一种提取oct色散失配系数的方法 | |
CN110160450B (zh) | 基于白光干涉光谱的大台阶高度的快速测量方法 | |
WO2020228293A1 (zh) | 太赫兹光谱成像数据的处理方法及装置 | |
CN108414462B (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 |