CN100545632C - 光纤光谱仪波长标定方法 - Google Patents
光纤光谱仪波长标定方法 Download PDFInfo
- Publication number
- CN100545632C CN100545632C CNB2007101779412A CN200710177941A CN100545632C CN 100545632 C CN100545632 C CN 100545632C CN B2007101779412 A CNB2007101779412 A CN B2007101779412A CN 200710177941 A CN200710177941 A CN 200710177941A CN 100545632 C CN100545632 C CN 100545632C
- Authority
- CN
- China
- Prior art keywords
- wavelength
- standard
- spectrometer
- optical fiber
- fiber spectrometer
- 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
Images
Landscapes
- Spectrometry And Color Measurement (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
Description
技术领域
本发明涉及一种光学仪器的标定方法,特别涉及一种光纤光谱仪波长的标定方法.
背景技术
光纤光谱仪是一种可用于紫外-可见-近红外光谱区的光谱测量的仪器,其测量波长范围一般在200~1100nm,广泛应用于化学分析、工业检测以及航空航天遥感等领域。近年来,随着光纤、CCD以及计算机技术的引入,光纤光谱仪的技术得到迅速的发展,可实现在线实时测量,更兼有信息采集、处理、存储等优点,从而使光谱测量技术在工业、科研领域中得到了广泛的应用。由于被测样品吸收或发射的光强是光波波长的函数,因此在实际测量中,光纤光谱仪的波长标定结果必须精确-即光纤光谱仪中各CCD像元所对应的光波波长必须精确确定,否则光谱测量的精确度就会降低。
下面是有关背景技术中的对光纤光谱仪波长进行标定的一些具体的方法,其中光纤光谱仪是目前最常用的是Quadratic(二次多项式)方法:
Linear(线性法):λi=a1+a2ki (1)
Quadratic(二次多项式法):λi=a1+a2ki+a3ki 2 (2)
Cubic(三次多项式法):λi=a1+a2ki+a3ki 2+a4ki 3 (3)
Trigonmetric1(三角法1):
Trigonmetric2(三角法2):
上式中:λi为波长,单位为纳米(nm);ai为待定系数;ki为像元编号;np为给定常数,其与CCD像元个数及光谱范围有关,通常可取np=2048。
采用Linear方法标定时,仅需要选取两条标准线状谱线来进行标定,虽然标定过程最简单,但标定结果误差很大,该方法目前在光纤光谱仪波长标定时很少采用。
目前光纤光谱仪几乎都采用Quadratic方法来进行波长标定,其标定过程比较简单,只需要选取三条标准谱线来进行标定,但标定结果误差较大,且误差大小与标定点的选择有直接关系;采用Quadratic方法标定时,用于标定的三条标准谱线有比较严格的要求:谱线跨度较大、且分布均匀,例如虽然在404.7nm~808.0nm范围内,误差较小;但在这个波长范围外,如在980.0nm处,误差很大,达到2.69nm。(其具体标定情况可参见后述的具体实施方式部分)
Cubic方法是在Quadratic方法基础上再增加一项高次项,可以有效的减小波长标定误差,但在980.0nm处误差最大,为0.35nm。采用Cubic方法标定时,需要采取四条标准谱线来进行标定,这在实际的标定过程中并不可取。
Trigonmetric1方法中引入一个sin函数,sin函数展开后,缺少偶次项。和Quadratic方法一样,用于标定的三条标准谱线也有比较严格的要求(如跨度大、需分布均匀),例如虽然在404.7nm~808.0nm范围内,误差较小;但在这个波长范围外,如在980.0nm处,误差很大,达到2.44nm;另外,在Trigonmetric1公式中,np为给定的参数,当np取不同值时,标定结果存在差异。
Trigonmetric2方法中分别引入sin函数和cos函数,可以有效的减小波长标定误差,在980.0nm处误差最大,为0.52nm;但采用Trigonmetric2方法标定时,需要采取四条标准谱线来进行标定,并且和Trigonmetric1方法存在相同的问题,即np为给定的参数,这在实际的标定过程中并不可取。
上述这些方法都存在着波长标定误差较大的缺陷。它们在其本质上都是多项式展开,波长标定的精确度随着展开级数的提高而提高,但用于标定的标准谱线也相应增多,计算量也相应增大;因此,如果只是提高多项式的级数并不能从根本上解决当前光纤光谱仪波长标定误差较大的问题,虽然曾有在多项式中引入sin、cos函数的方法,但其本质上仍然是多项式展开,并且人为的引入了相关参数。
发明内容
为了克服现有技术中对光纤光谱仪标定误差大,且用于标定的三条标准谱线要求高,如需跨度大、分布均匀等,而且需要一些给定的参数,并当给定参数取不同值时,标定结果也出现差异的问题。本发明提出一种标定误差非常小的光纤光谱仪标定方法,只需任意三条标准普线,就可以精确标定的方法,其提高了标定的灵活性,使用时更加方便、简洁。
根据本发明的标定方法,其包括以下步骤:
求出a1、a2、a3的解,其中d为光谱仪中的光栅常数,i为自然数,i=1,2,3…;
C)根据步骤B)中得出的a1、a2、a3值,再根据下式求解出对应于每个像元编号ki的标定计算波长值λi:
本发明的Formula(解析式法)方法相比于其他方法的优点有:
①波长标定精度得到极大的提高,标定误差可以控制在0.05nm以内;
②用于标定的三条标准谱线可以任意选择,不同于Quadratic、Trigonmetric1方法对标准谱线有比较严格的要求。③标定误差不仅在404.7nm~808.0nm范围内非常小,在这个波长范围外也非常小,如在980.0nm处,仅为0.05nm。
总之,与传统标定方法的标定结果进行比较。采用我们提出的方法来进行波长标定,其标定误差最大不超过0.06nm,如果用于标定的标准线状谱线峰值所对应的CCD像元读数足够精确,标定误差可以达到0.01nm以下,远远优于当前的标定方法,极大的提高了光谱测量的精确性。
附图说明
图1是光纤光谱仪构造原理图;
图2是图1所示光纤光谱仪的等效光路图;
图3是标定装置标定过程的示意图;
图4是光谱仪采集的标准线状谱线图;
图5是用本发明的方法和背景技术的方法进行波长标定比较的误差图(用于标定的标准线状谱线跨度大、分布均匀);
图6是用本发明的方法和当前常用的方法进行波长标定比较的误差图。(用于标定的标准线状谱线跨度小、分布不均匀);
具体实施方式
下面参考附图,具体描述本发明的实施例。如图1所示,为普通光纤光谱仪的构造:由光纤1、FC耦合探头2、球面反射镜3、光栅4、CCD线阵5、计算机操作系统6构成;光路采用水平对称式光路,光栅中心位于水平对称式光路的对称轴上,FC耦合探头2、CCD线阵5位于球面镜的焦平面上。
根据几何光学原理,可把光纤光谱仪光路简化为光路图2,在试验中,采用的是1级衍射条纹,故由光栅方程可知:
d(sinβ-sinα)=λ (6)
式中:β-入射角,α-衍射角,λ-波长(nm)
由图2可知,假设波长为λ的入射光经光栅衍射后,汇聚于CCD线阵P点上,P点的像元编号为k,根据几何原理,可以得到方程(7)
f tanα=kΔ+B (7)
式中:Δ-单个像元线度(mm);f-球面反射镜的焦距(mm);k-像元编号;B-常数
将(6)式代入(7)式,得:
为了方便计算,令: a3=sinβ;故(8)式可化简为:
式中:a1、a2、a3——待定常数。
根据公式(9),本发明光栅光谱仪标定方法的步骤如下:
求出a1、a2、a3的解,其中d为光谱仪中的光栅常数,i为自然数,i=1,2,3…;
C)根据步骤B)中得出的a1、a2、a3值,再根据下式求解出对应于每个像元编号ki的标定计算波长值λi:
从上述步骤可以看出,光纤光谱仪一旦组装完毕,a1、a2、a3就均为定值,不同波长的光谱谱线汇聚在与其相对应的CCD像元上,在本发明的该实施例中,可采用低压汞灯(其波长为404.7nm、435.8nm、546.1nm、577.0nm、579.0nm,)和激光器(其波长为532.0nm、632.8nm、808.0nm、980.0nm)为标准线状谱线,可通过计算机系统准确读出各标准线状谱线峰值所对应的像元。由方程(10)可知,方程中有三个未知数,只需要选择标准线状谱线中的三条谱线,将波长像元ki分别代入方程(10)中,即可求出a1、a2、a3的值(a1、a2、a3的计算值与真值的误差取决于光纤光谱仪的分辨本领以及标准线状谱线波长像元ki读数的准确性);最后,将a1、a2、a3的计算值代入公式(11)中,可以求出各个像元ki所对应的波长λi,即为标定波长,本文将该方法定义为解析式法。
本实施例中,两个球面反射镜的焦距f约为100mm,CCD线阵为2048点线阵,单个像元线度Δ约为14um,响应波长为200nm~1100nm;光栅常数d=2.5um。
一般来讲,用于标定的实验装置如图3所示,主要包括标准谱线光源、光纤光谱仪和计算机,计算机用来分析光纤光谱仪的光谱信号,确定标准线状谱线峰值与像元的对应关系。
下面具体描述本发明实施例中标定的详细步骤和结果。首先使用标准谱线光源(低压汞灯和激光器),通过光纤光谱仪采集标准光源的光谱。确定标准光源光谱的实际波长对应的光谱仪接收元件线阵CCD的像元编号ki,得到的对应关系如表1所示:
表1
标准谱线波长(nm) | 404.7 | 435.8 | 532.0 | 546.1 | 632.8 | 808.0 | 980.0 |
对应像元编号(k<sub>i</sub>) | 128.0 | 229.0 | 538.0 | 583.0 | 858.0 | 1409. | 1950. |
选取上述波长中的任意三个实际波长以及对应的像元编号ki,求解常数a1、a2、a3。这里选取实际波长为404.7nm、632.8nm、808.0nm的三条标准谱线来进行标定,分别将其波长像元编号ki代入公式(10)中,得到:
a1=0.17377735 a2=-0.0001274271 a3=0.31743
根据上面求解出的a1、a2、a3值,再将其代入公式(11)中,就计算出每个像元编号ki对应的计算波长λi,表2显示了部分ki及其对应的计算波长λi值。
表2
k<sub>i</sub> λ<sub>i</sub> | k<sub>i</sub> λ<sub>i</sub> | k<sub>i</sub> λ<sub>i</sub> |
0.10 365.546598960.20 365.577065380.30 365.607532000.40 365.637998820.50 365.668465830.60 365.698933040.70 365.729400450.80 365.759868050.90 365.790335851.00 365.820803841.10 365.851272031.20 365.881740421.30 365.912209001.40 365.942677781.50 365.973146761.60 366.003615931.70 366.034085301.80 366.064554861.90 366.09502463|| | ||1950.70 980.008952471950.80 980.040543801950.90 980.072135041951.00 980.103726201951.10 980.135317261951.20 980.166908231951.30 980.198499111951.40 980.230089911951.50 980.261680611951.60 980.293271221951.70 980.324861741951.80 980.356452181951.90 980.388042521952.00 980.419632771952.10 980.451222941952.20 980.482813011952.30 980.514402991952.40 980.545992891952.50 980.57758269|| | ||2048.20 1010.765566952048.30 1010.797063672048.40 1010.828560292048.50 1010.860056802048.60 1010.891553212048.70 1010.923049522048.80 1010.954545722048.90 1010.986041822049.00 1011.017537812049.10 1011.049033702049.20 1011.080529482049.30 1011.112025162049.40 1011.143520742049.50 1011.175016212049.60 1011.206511572049.70 1011.238006842049.80 1011.269501992049.90 1011.300997052050.00 1011.33249200 |
在评价光纤光谱仪标定误差时,通常采用预测标准误差(SEE)作为评价标准,如下式所示:
为了对本发明的标定方法与背景技术的标定方法进行比较,表3给出了各种标定方法的标定波长以及标定误差,其误差图如图5所示。
表3
从上表可以看出,本发明的解析式法的预测标准误差(SEE)为0.05,而线性法、二次多项式法、三次多项式法、三角法1、三角法2分别为1.55、1.36、0.21、1.25、0.30;本发明的误差远远低于背景技术的误差。
为了充分说明本发明解析式法的优越性,尤其是相对于当前通用的二次多项式法方法,我们采取632.8nm、808.0nm、980.0nm(分布不均、跨度小)三条标准线状谱线来进行标定,其标定结果如表4所示,根据表4绘制出的误差图,如图6所示:
表4
由图6可知,如果用于标定的标准线状谱线跨度较小且分布不均时,对Formula法标定结果没有影响,而对Quadratic法标定结果影响很大,虽然在632.8nm~980.0nm范围内误差较小,但在这个范围之外,误差急剧增大,在404.7nm处达到了4.33nm。
Claims (4)
2、如权利要求1所述的光纤光谱仪波长标定方法,其特征在于:使用的标准光源是波长为404.7nm、435.8nm、546.1nm、577.0nm、579.0nm的低压汞灯,和/或波长为532.0nm、632.8nm、808.0nm、980.0nm的激光器。
3、如权利要求2所述的光纤光谱仪波长标定方法,其特征在于:使用标准光源的三个实际波长为404.7nm、632.8nm、808.0nm的谱线,计算参数a1、a2、a3来进行标定。
4、如权利要求2所述的光纤光谱仪波长标定方法,其特征在于:使用标准光源的三个实际波长为632.8nm、808.0nm、980.0nm的谱线,计算参数a1、a2、a3来进行标定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007101779412A CN100545632C (zh) | 2007-11-22 | 2007-11-22 | 光纤光谱仪波长标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007101779412A CN100545632C (zh) | 2007-11-22 | 2007-11-22 | 光纤光谱仪波长标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101158636A CN101158636A (zh) | 2008-04-09 |
CN100545632C true CN100545632C (zh) | 2009-09-30 |
Family
ID=39306797
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007101779412A Expired - Fee Related CN100545632C (zh) | 2007-11-22 | 2007-11-22 | 光纤光谱仪波长标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100545632C (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102141440A (zh) * | 2010-12-30 | 2011-08-03 | 中国科学院长春光学精密机械与物理研究所 | 高分辨、微小型、宽光谱反射式光谱仪光学系统 |
CN102818630B (zh) * | 2012-07-27 | 2014-10-08 | 中国资源卫星应用中心 | 一种干涉型成像光谱仪的光谱定标方法 |
CN103115880B (zh) * | 2013-01-18 | 2015-03-11 | 山东大学 | 一种激光诱发听觉神经的光谱分析方法 |
CN103471813B (zh) * | 2013-08-27 | 2015-09-09 | 武汉光迅科技股份有限公司 | 多通道可调激光器定标装置及方法以及测试装置及方法 |
CN103616074B (zh) * | 2013-11-21 | 2015-06-10 | 中国科学院长春光学精密机械与物理研究所 | 数字微镜光栅光谱仪的波长标定方法 |
CN104266758B (zh) * | 2014-10-21 | 2016-05-25 | 中国科学院光电研究院 | 一种大孔径空间外差干涉光谱成像仪基准波数定标方法 |
CN109886351B (zh) * | 2019-03-04 | 2021-10-15 | 北京麦飞科技有限公司 | 高光谱数据和高分辨率图像融合方法 |
CN112595440B (zh) * | 2020-12-30 | 2022-11-01 | 西安和其光电科技股份有限公司 | 一种用于光纤测温仪的光谱标定方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0207498A2 (en) * | 1985-07-02 | 1987-01-07 | Shimadzu Corporation | Spectrophotometer |
US4643572A (en) * | 1983-08-04 | 1987-02-17 | Anritsu Electric Company, Limited | Driving circuit of a grating |
CN1125321A (zh) * | 1994-12-20 | 1996-06-26 | 中国科学院化工冶金研究所 | 智能光栅单色仪的波长校正装置及其方法 |
WO1997032187A1 (en) * | 1996-02-27 | 1997-09-04 | Amersham Pharmacia Biotech Ab | Calibration method |
WO2006066581A1 (en) * | 2004-12-21 | 2006-06-29 | Foss Analytical A/S | A method for standardising a spectrometer |
JP2007101365A (ja) * | 2005-10-04 | 2007-04-19 | Univ Of Tsukuba | 光コヒーレンストモグラフィーの構成機器の較正方法 |
-
2007
- 2007-11-22 CN CNB2007101779412A patent/CN100545632C/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4643572A (en) * | 1983-08-04 | 1987-02-17 | Anritsu Electric Company, Limited | Driving circuit of a grating |
EP0207498A2 (en) * | 1985-07-02 | 1987-01-07 | Shimadzu Corporation | Spectrophotometer |
CN1125321A (zh) * | 1994-12-20 | 1996-06-26 | 中国科学院化工冶金研究所 | 智能光栅单色仪的波长校正装置及其方法 |
WO1997032187A1 (en) * | 1996-02-27 | 1997-09-04 | Amersham Pharmacia Biotech Ab | Calibration method |
WO2006066581A1 (en) * | 2004-12-21 | 2006-06-29 | Foss Analytical A/S | A method for standardising a spectrometer |
JP2007101365A (ja) * | 2005-10-04 | 2007-04-19 | Univ Of Tsukuba | 光コヒーレンストモグラフィーの構成機器の較正方法 |
Non-Patent Citations (1)
Title |
---|
对PC2000-PC/104 型光谱仪的波长定标. 徐中民,禹秉熙.光学精密工程,第12卷第1期. 2004 * |
Also Published As
Publication number | Publication date |
---|---|
CN101158636A (zh) | 2008-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100545632C (zh) | 光纤光谱仪波长标定方法 | |
EP0800066B1 (en) | Precise calibration of wavelengths in a spectrometer | |
US6029115A (en) | Analyzing spectrometric data | |
CN102914268A (zh) | 膜厚测量装置 | |
CN102798342A (zh) | 一种用于光学散射测量的基于拟合误差插值的库匹配方法 | |
CN105628676A (zh) | 一种拉曼光谱修正系统及方法 | |
JP2001525931A (ja) | 光スペクトル測定用装置とその方法 | |
CN101813519B (zh) | 光谱仪杂散光校正方法 | |
CN102564591B (zh) | 一种光谱分析仪和光谱分析方法 | |
CN102778202B (zh) | 一种膜厚测量装置及方法 | |
White et al. | Compact diffraction grating laser wavemeter with sub-picometer accuracy and picowatt sensitivity using a webcam imaging sensor | |
US3606547A (en) | Spectrophotometer | |
ITAN20070019A1 (it) | Rifrattometro spettrofotometrico | |
Workman Jr | Calibration transfer, Part V: The mathematics of wavelength standards used for spectroscopy | |
Danyang et al. | Design and implementation of high sensitivity micro spectrometer based on area array CCD | |
CN113494890B (zh) | 基于fpi干涉仪的光纤光栅应变传感器精度测量装置及方法 | |
US11543294B2 (en) | Optical technique for material characterization | |
CN102589706B (zh) | 一种光纤束发射光谱层析中光纤传感单元的标定方法 | |
Sacco et al. | Graphene edge method for three‐dimensional probing of Raman microscopes focal volumes | |
CN112700373A (zh) | 一种基于非线性回归数据处理的提高光谱图像分辨率的方法、系统、设备及存储介质 | |
Wang et al. | A multimodal spectrometer for Raman scattering and near-infrared absorption measurement | |
Small | Improved out-of-plane BRDF measurement and modeling | |
JP2004177147A (ja) | 発光測定装置 | |
CN105180801B (zh) | 一种对比式抗干扰阶梯型角反射镜激光干涉仪及标定方法和测量方法 | |
JP2005121574A (ja) | 近赤外分光装置 |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090930 Termination date: 20131122 |