氧气A吸收带偏最小二乘基线拟合方法 下载: 868次
1 引言
氧气A带的基线拟合是依据黑体辐射理论,利用红外目标辐射光谱经大气衰减后测得的吸收光谱拟合出无衰减时的光谱曲线,即利用吸收带外的光谱数据拟合基线[1-2]。测量基线的拟合精度直接影响到测距的精度,利用带外光谱拟合基线主要考虑其是否为最佳拟合[3]。目前,关于光谱的基线拟合方法有插值法[4]、多项式拟合[5]、分段线性拟合方法[6]等,针对氧分子A吸收带内光谱波动大的情况,带外数据的基线拟合都存在误差平方过大的问题,这将影响光谱测量的精度,进而影响到光谱透射率的精度。而回归分析法根据测量数据分析变量之间的相互关系,更适用于氧分子A带基线的拟合。目前常用的回归拟合方法有多元线性回归(MLR)、主成分回归(PCR)、最小二乘回归和偏最小二乘回归(PLSR)等方法。最常用的是最小二乘回归方法[7],但最小二乘回归法的不稳健性会导致基线拟合偏差较大。使用MLR的前提是自变量必须相互独立,否则会出现多重共线性的问题,使得回归系数不稳定,导致计算误差变大[8]。PCR方法运算速度比MLR慢,且PCR仅考虑了自变量的信息,没考虑因变量的解释作用,会造成信息的丢失[9]。而PLSR既考虑了自变量的信息,又考虑了因变量对自变量的解释作用,稳定性强[10]。在相同环境,光谱仪分辨率为8 cm-1时,多次测量分别采用上述方法进行基线拟合,利用上述不同方法获取的拟合基线标准偏差分别为:最小二乘回归为0.32%,MLR法为0.51%,PCR法为0.42%,而PLSR法为0.23%。综上所述,采用偏最小二乘曲线拟合方法拟合氧分子A吸收带的基线,使数据点均匀分布在基线两侧,不仅考虑所计算的带外数据与拟合基线方差为最大,还考虑了需拟合的基线与带外数据相关性最大,有效区分系统信息和噪声,提高系统建模的准确性。
2 目标光谱的带平均透射率计算模型
氧分子光谱区间Δ
式中
如
根据(2)式,目标光谱
3 带外数据的基线拟合方法
3.1 PLSR数学模型
根据比尔定律,光强
式中
将光强
式中
同时考虑
式中
(3)~(6)式构成PLSR的基本数学模型。(4)式和(5)式为外部模型,(3)式为内部模型。外部模型各自进行拟合误差分解,内部模型则把两个外部模型通过最大程度的线性相关紧密联系。
经过不断迭代计算,若回归方程达到精度要求,则算法中止。否则继续提取第二对光谱信息,直到满足精度要求为止。
氧气吸收带外数据利用PLSR的方法可以得到被测目标光谱的基线,但是目标光谱测量过程中不可避免会出现奇异点,它会给PLSR相对较高的权重。因此过多地强调拟合基线和这些奇异点的关系,会导致拟合出现偏差,据此做出的基线可能会和实际情况有相当大的差距。因此必须剔除带外数据奇异点,获得新的回归参数,以期达到最佳拟合。
3.2 奇异点的剔除
基于上述基线拟合方法,采用Kennard-Stone (K-S)选择法[12]进行奇异点剔除。
基于带外光谱特征计算每个数据点与相邻数据点之间的欧氏距离[13],在带外均匀选择样本,其中,欧氏距离为:
式中
设有
1) 首先计算所有数据两两之间的欧式距离
2) 计算剩余
3) 计算剩余
4) 重复上述过程,直至选中
利用上述方法剔除了奇异点前后的效果如
综上可见,奇异点剔除前后基线拟合的曲线发生了变化,模型中存在的奇异点明显使残差增大,所以剔除奇异点后能有效提高基线拟合的准确性。利用相应的评价标准来判断是否达到最佳基线拟合。
3.3 基线拟合评价标准
采用基线拟合的方式估算原始光谱信息,不可避免地会引入测量结果的不确定性及误差。引入带平均吸收率的标准偏差来说明基线拟合的不确定度。
表 1. 奇异点剔除前后的平均残差
Table 1. Mean residual before and after eliminating singular points
|
根据(2)式可得带平均吸收率为:
式中,分子为被测目标辐射光强算术平均值,而分母是带外数据拟合基线的算术平均值计算带平均吸收率时,带外数据拟合基线至关重要,而基线拟合方法有很多,其评价标准可用统计分析方法计算不确定度表示[14],多次重复测量得到被测目标强度的平均值为:
总的带平均吸收率的不确定度用
式中,
将(8)式代入上式可得:
式中,
式中,
基线拟合是带平均透射率计算的关键,以(13)式为依据来评价带平均透射率计算的准确度,为预测模型误差的计算提供依据。
4 实验结果
为了验证氧气A带基线拟合方法精度对带平均透射率的影响,在空气温度为19.8 ℃,标准大气压下,在0~130 m内任取9个距离点进行测量,以500 W石英卤素灯作为光源,输入端采用直径为165 mm望远镜,采用不同分辨率分别多次扫描取平均后获得目标的衰减光谱。
为了分析基线拟合不确定度与不同分辨率的关系,在46.4 m处取不同分辨率下光谱数据进行氧分子A带的基线拟合,根据(12)式计算不同分辨率下的拟合标准偏差,结果如
表 2. 氧气A带46.4 m处基线拟合参数
Table 2. Baseline fitting parameters of oxygen A-band at 46.4 m
|
可见,同一距离点处不同分辨率下的平均标准偏差为0.23%,随着分辨率的降低,基线拟合不确定度变小,信噪比增大;但并非分辨率越小越好,当分辨率降低到32 cm-1时,虽然信噪比最高,但是不确定度却增加了。
在9个不同距离点处,不同分辨率下获得拟合不确定度与距离的关系如
图 3. 带平均吸收率基线拟合标准偏差
Fig. 3. Baseline fitting standard deviation with band average absorptivity
随着被测距离的增大,拟合不确定度随分辨率的降低而减小,总的拟合不确定度不超过0.4%。理论和实验证明,带平均透射率的计算方法还可消除吸收带内外仪器和散射带来的误差。
依据吸收带外数据拟合带内基线,计算带平均透射率的方法不仅是可行的,而且可以不依赖目标先验知识而获取目标光谱透射率。由
5 结论
本文以红外目标在氧分子A带的吸收峰为研究对象,利用带外数据对带内光谱进行基线拟合。为了减小基线拟合精度对带平均透射率的影响,首先剔除光谱采集中产生的奇异点,采用偏最小二乘拟合方法进行基线拟合。搭建了实验系统,通过对不同距离点处的吸收光谱拟合基线,计算了不同点、不同分辨率下的拟合不确定度。结果证明,同一距离点处不同分辨率下的平均标准偏差为0.23%,随着分辨率的降低,基线拟合不确定度变小,信噪比增大;随着被测距离的增大,拟合不确定度随分辨率的降低而减小,总的拟合不确定度不超过0.4%。该结论为后续测距系统的误差分析提供了依据。
表 3. 不同距离、不同分辨率的基线拟合不确定度Table 3 Baseline fitting uncertainty with different distances and resolutions%
Table 3.
|
[1] Hawks MR, Perram GP.Passive ranging of emissive targets using atmospheric oxygen absorption lines[C]. SPIE, 2005, 5811: 112-122.
[2] 吕航, 廖宁放, 吴文敏, 等. 基于自适应差分滤波的干涉图基线校正方法[J]. 光学学报, 2015, 35(10): 1030002.
[3] 刘霞, 王运锋. 基于最小二乘法的自动分段多项式曲线拟合方法研究[J]. 科学技术与工程, 2014, 14(3): 55-58.
[4] 宗鹏飞, 王志斌, 张记龙, 等. 基于红外被动测距的基线拟合算法研究[J]. 激光技术, 2013, 37(2): 174-177.
[5] 张宇峰, 戴景民, 陆晓东, 等. 红外光谱发射率测量系统的温漂修正方法[J]. 光谱学与光谱分析, 2016, 36(9): 3003-3007.
[6] 覃赵军, 陶站华, 刘军贤, 等. 分段式线性拟合校正拉曼光谱基线漂移[J]. 光谱学与光谱分析, 2013, 33(2): 383-386.
[7] 刁伟峰, 刘继桥, 竹孝鹏, 等. 全光纤相干多普勒激光雷达非线性最小二乘风速反演方法及实验研究[J]. 中国激光, 2015, 42(9): 0914003.
[9] 杜一平, 潘铁英, 张玉兰. 化学计量学应用[M]. 北京: 化学工业出版社, 2008.
[10] 李晋华, 杨志良, 王召巴, 等. 近红外漫透射技术检测玉米成分[J]. 红外技术, 2013, 35(11): 732-736.
[11] 李晋华, 王志斌, 陈媛媛, 等. 近红外氧气A带大气透射率的计算[J]. 激光与红外, 2013, 43(10): 1142-1145.
[12] Danielsson P E. Euclidean distance mapping[J]. Computer Graphics and Image Processing, 1980, 14(3): 227-248.
[13] 李慎安. 《测量不确定度评定与表示》计量技术规范问与答(上)[J]. 计量技术, 2002(11): 46-48.
[14] 曹宏燕. 分析测试中测量不确定度及评定第五部分测量不确定度评定中要注意的一些问题[J]. 冶金分析, 2006, 26(4): 89-94.
李晋华, 张敏娟, 张德龙. 氧气A吸收带偏最小二乘基线拟合方法[J]. 中国激光, 2017, 44(8): 0811001. Li Jinhua, Zhang Minjuan, Zhang Delong. Baseline Fitting of Partial Least Squares for Oxygen A Absorption Band[J]. Chinese Journal of Lasers, 2017, 44(8): 0811001.