基于PCA-ELM的红外多光谱辐射测温
沈阳航空航天大学 自动化学院, 沈阳 110136
Infrared Multispectral Radiation Temperature Measurement Based on PCA-ELM
School of Automation, Shenyang Aerospace University, Shenyang 110136, China
通讯作者: 傅 莉,女,教授,电话(Tel.): 15840355071;E-mail:ffulli@163.com.
责任编辑: 石易文
收稿日期: 2020-01-22
基金资助: |
|
Received: 2020-01-22
作者简介 About authors
席剑辉(1975-),女,辽宁省沈阳市人,副教授,主要研究方向为复杂系统模型辨识、故障检测与诊断、红外辐射测试与分析等 。
在目标发射率未知的情况下,建立一种基于主元分析(PCA)与极限学习机(ELM)相结合的红外多光谱测温方法.分析目标温度与辐射亮度谱的非线性数学模型,确定初始输入向量包含温度估计所需的充分信息;引入PCA方法从输入向量中提取相互独立的主元成分,降低神经网络输入维数;基于ELM网络对样本数据充分学习,最终建立PCA-ELM目标红外测温模型.利用黑体和未知发射率材料涂层目标作为测试目标源,验证该方法的有效性.
关键词:
In the case of unknown target emissivity, an infrared multispectral radiation temperature measurement method based on principal component analysis (PCA) and extreme learning machine (ELM) is established. The nonlinear mathematic model of target temperature and radiance spectrum is analyzed to find a set of initial input vectors, which include sufficient information to estimate temperature. The PCA method is used to extract the independent principle components in input vectors. This method can also reduce the input dimension for neural network. Based on ELM network, the sample data is sufficiently learned to build the target infrared temperature measurement model by PCA-ELM. The effectiveness of the proposed method is verified by using the blackbody and the coating material with unknown emissivity as test target sources.
Keywords:
本文引用格式
席剑辉, 姜瀚, 陈博, 傅莉.
XI Jianhui, JIANG Han, CHEN Bo, FU Li.
但在实际测量中,发射率与目标材料的成分、表面状况、所处的温度和考察的波长等因素均有关,无法建立完备的目标发射率数据库,因此在测温过程中,目标的发射率往往都是未知的.目前,针对目标发射率未知情况下的红外测温方法一般有比色测温法[6,7] 以及多光谱测温法[8,9].比色测温法是通过测量两个波段的辐射能量之比计算温度的.通常选取两个相近波段,认为两者的发射率近似相等可以互相消除.若所选波段的发射率发生较大变化时,其相应的测温误差也会较大.多光谱测温法通过测量多个目标光谱的辐射亮度建立测温模型.最初,Svet等[10]引入发射率非线性模型,建立一组与模型未知参数个数相等的方程组,通过求解方程组估计目标温度和光谱发射率.所建立的方程组数量与结构对估计结果的影响较大,线性方程难以表征目标温度与红外光谱信息之间的强非线性关系.建立非线性高阶方程组会导致方程组维数迅速上升,对样本数量和求解算法提出更高的要求. 测温过程中,大气、背景等复杂影响因子要求模型辨识方法应该具有良好的稳健性和适应性.近年来,自适应的智能模型辨识方法经常被用于多光谱测温[11,12].
1 目标辐射温度建模
当温度为T时,目标辐射亮度的典型曲线如图1所示.其中:λj为波长;Lj (λj, T) 为辐射亮度.当λj 分别取固定值为λ1,λ2,…,λp时,固定温度下对应的辐射亮度值分别为L1,L2,…,Lp.
图1
图1
当温度为T时的典型目标辐射亮度曲线
Fig.1
Radiance curve of a typical target at a temperature of T
由Planck定律和光谱发射率定义可得,当温度为T时,目标的辐射亮度为
式中:ε(λj, T)为光谱发射率;C1和C2分别为第一热辐射常数、第二热辐射常数.在同一温度T下,不同波长λ1,λ2,…,λp对应的亮度为
式中:ε(λ1, T),ε(λ2, T),…,ε(λp, T)为未知的光谱发射率,只与T和λj有关.因此,L1,L2,…,Lp可以看成p个关于温度T和λj的隐函数,通过选用固定波长下的多光谱亮度信息,就可以对温度进行估计,建立基于多光谱亮度信息的目标温度模型.
由式(2)可知:① 方程组具有强非线性,常规模型辨识方法很难获取亮度谱与温度之间的关系;② 方程组个数p不少于等式右边展开后未知参数的个数,p值过大时,式(2)求解困难.神经网络不需要具体知道系统动态特性,无需建立初始数学模型,尤其适用于难以用确定数学模型表达的非线性系统建模[13].所设计的测温模型结构如图2所示,其中z1,z2,…,zm (m<p)为主元成分.模型输入端为L1,L2,…,Lp,输出端为T.通过在已知温度下测量亮度谱曲线,获取输入输出样本对;再利用神经网络的强学习能力寻找输入输出间的潜在规律.为使输入包含足够的建模信息,输入亮度L1,L2,…,Lp往往取值密集,造成输入维数p数值高、建模困难.输入变量之间存在较强的相关性和冗余信息,导致网络学习精度下降,可泛化性弱.因此,在神经网络输入端引入PCA层提取输入变量中的主元成分z1,z2,…,zm (m<p),通过降低输入维数,以提高建模精度.
图2
2 PCA-ELM红外光谱测温方法
在不同温度T下测量光谱亮度,总共测量n个不同温度.设第k个温度下,对应p个波长的光谱辐射亮度输入为L(k)=[L1(k) L2(k) … Lp(k)]T (k=1,2,…,n).若输入变量间方差越大,相关性越小,则越具有好的区分能力.PCA 方法正是利用这一规律,以确定具有较大方差的主元变量为目标,对输入变量进行线性组合,使获得的输入向量具有更好的模式区分能力.
ELM是一种单隐层前馈网络,网络结构简单,输入层与隐含层连接权值以及隐含层神经元阈值在训练过程中无需调整,训练时间快,因此本文神经网络选用ELM神经网络.
2.1 PCA的原理和步骤
设有n个不同温度下的目标辐射亮度向量,每个向量包含p个不同波长下对应的亮度值,构成一个p×n阶的矩阵Λ,
令z(k)=[z1(k) z2(k) … zm(k)]T (m<p),为Λ中蕴含的主元变量,为便于表达,把变量k省略,则不同温度下,都满足下式:
式中:
(1) 将原始数据标准化.
(2) 建立变量的相关系数矩阵R=[rij]p×p.
式中:
(3) 求R的特征根λ1≥λ2≥…≥λp>0及其相应的单位特征向量.
a1=
a2=
ap=
(4) 计算前m个特征值的累计贡献率γm.
(5) 一般当γm>85%~95%时,可得m个主成分zi,令i=1,2,…,m
2.2 ELM的原理
图3
对第k个学习样本,有网络输入为z(k)=[z1(k) z2(k) … zm(k)]T,网络输出为对应该输入的第k个温度T(k).共测量n个温度,则有k=1,2,…,n;m为PCA确定的主元维数.ELM网络隐含层的输出为
式中:g(·)为激活函数,ELM对激活函数没有特别要求,常见的激活函数有sigmoid函数和tanh函数;Wl为输入权重;βl为隐含层第l个神经元到输出的权重;bl为隐含层第l个神经元的阈值.单隐层神经网络的训练目标是要使得拟合的误差最小,即网络输出T(k)与理想输出O(k)差的范数趋向于0,可表示为
寻找bl、Wl、βl使得
设Y为隐含层节点输出矩阵,通过矩阵初等变换得到矩阵的秩为h,则矩阵Y可分解为
式中:F为n×h满秩矩阵;G为h×N满秩矩阵.
最终找到
输入权重
式中:GH和FH为G和F的共轭转置矩阵,则有:
3 具体实例
实验采用加拿大ABB BOMEN 的MR170 型光谱辐射计,测量目标在不同温度下的辐射特性.辐射计的分辨率为32/cm,探测器的光谱范围为3~14μm波段.
3.1 黑体测温及验证
实验室现有3种不同型号的黑体,发射率出厂值以及测温范围略有不同.黑体A的发射率为0.98、有效测温范围为50~250℃;黑体B的发射率为0.96、有效测温范围为300~600℃;黑体D的发射率为0.99、有效测温范围为700~1200℃.分别将3种黑体作为测试目标源,利用光谱辐射计分别测量不同温度下的辐射特性,得到43组样本数据,其中黑体A有28组、黑体B有9组、黑体D有6组.3种黑体的发射率有差异、温度范围没有交叉,现将这3种黑体的数据同时做为训练样本,验证所提方法的测温精度.测试曲线采用文献[13]的方法进行校准补偿,最终得到测量波段光谱亮度样本,随机选取7组不同温度下的亮度曲线如图4所示.
图4
图4
黑体在不同温度下的光谱辐射亮度
Fig.4
Spectral radiance of blackbody at different temperatures
光谱亮度会受到杂散辐射、大气吸收散射等影响,导致某些固定波段内的值会产生较大震荡.选择测温训练样本时,可以避开这些波段.在同一温度下,选择有效辐射亮度值141个,样本可对测量波段全覆盖,使得构成的输入向量中包含温度估计需要的光谱信息,最终获取样本共43组,部分样本数据如表1所示.
表1 黑体光谱亮度样本
Tab.1
T/℃ | L1/(W·m-2· sr-1·μm-1) | L2/(W·m-2· sr-1·μm-1) | … | L141/(W·m-2· sr-1·μm-1) |
---|---|---|---|---|
80 | 0.629 | 0.671 | … | 17 |
85 | 0.755 | 0.781 | … | 17 |
︙ | ︙ | ︙ | … | ︙ |
1300 | 24158 | 23740 | … | 416 |
在不同温度段随机选取35组光谱亮度做为训练样本,可得式(3)的亮度矩阵Λ是一个35×141阶矩阵.根据PCA计算主元成分,获得的特征值和特征向量如表2所示.
表2 黑体光谱亮度PCA结果
Tab.2
特征值 | 特征向量 |
---|---|
λ1=138.7186 | a1=[ |
λ2=2.2762 | a2=[ |
λ3=0.0051 | a3=[ |
︙ | ︙ |
λ141=0 | a141=[ |
由式(6)可知,计算前两个特征值的贡献率γ1=98.3%,γ2=1.614%,累计贡献率达99.9%,确定主成分为2,得到降维后的网络输入为z1(k)=a1·L(k),z2(k)=a2·L(k)(k=1,2,…,n).
将z1和z2输入ELM网络,网络输出为T,隐含层节点数为8个.训练结束后,将8组测试样本输入网络得到的验证结果如图5所示.
图5
图5
黑体温度的ELM网络输出值与实测值比较
Fig.5
Comparison of blackbody temperature between estimated values obtained by ELM network and practical values
式中:Lbb(λj,T)为温度T时的标准黑体光谱辐射亮度.理想黑体光谱发射率不随波长变化,表现为常值,所测黑体为哈尔滨工业大学自动检测与过程控制系统研究所研制的面源黑体,发射率出厂值为0.98.根据表3得到的温度可获得理想黑体辐射亮度,根据式(16)计算的光谱发射率如图6所示,其中T(1),T(2),…,T(8)为8个温度.由图6可知,不同温度下光谱发射率在0.98上下浮动,且比较恒定.计算测量波段的平均发射率
表3 黑体的测温结果及其相对误差
Tab.3
T1/℃ | T2/℃ | |||
---|---|---|---|---|
95 | 97.4762 | 2.61 | 0.9318 | 4.8 |
115 | 119.1212 | 3.58 | 0.9673 | 1.3 |
135 | 131.8413 | 2.34 | 1.0110 | 3.2 |
155 | 150.7401 | 2.75 | 1.0217 | 4.3 |
180 | 176.3031 | 2.05 | 1.0151 | 3.6 |
450 | 455.9431 | 1.32 | 0.9410 | 4.0 |
520 | 517.7082 | 0.44 | 0.9688 | 1.1 |
1000 | 1025.3000 | 2.53 | 0.9523 | 2.8 |
图6
图6
黑体测试样本的光谱发射率估计
Fig.6
Estimation of spectral emissivity for testing samples of blackbody
通过黑体实验可以验证所提测温方法的合理性,且测温范围广,主要取决于所选网络训练样本的温度范围.
3.2 目标涂层的测温及验证
将所提方法应用于飞机机身表面涂层的温度估计,如图7所示.涂层温度通过外接温度控制器控制,利用光谱辐射计测量表面涂层的辐射特性,得到41组不同温度下的光谱亮度辐射数据,温度范围为80~180 ℃.
图7
与黑体实验相同,在每条温度曲线上,选取141个不同波长对应的辐射亮度值作为初始输入向量.部分样本数据如表4所示.
表4 目标涂层光谱亮度部分样本
Tab.4
T/℃ | L1/(W·m-2· sr-1·mm-1) | L2/(W·m-2· sr-1·mm-1) | … | L141/(W·m-2· sr-1·mm-1) |
---|---|---|---|---|
80 | 0.516 | 0.567 | … | 13.852 |
83 | 0.562 | 0.614 | … | 14.117 |
︙ | ︙ | ︙ | ︙ | ︙ |
180 | 2.417 | 2.233 | … | 17.000 |
随机选取34组训练样本,可得式(3)的亮度矩阵Λ是一个34×141阶矩阵.根据PCA计算主元成分,得到特征值和特征向量如表5所示.计算累计贡献率,确定主成分为1,得到降维后的网络输入为z1(k)=a1·L(k),k=1,2,…,n.得到第k个样本降维后的输入变量为
表5 目标涂层光谱亮度PCA结果
Tab.5
特征值 | 特征向量 |
---|---|
λ1=125.3228 | a1=[ |
λ2=12.0106 | a2=[ |
λ3=2.6217 | a3=[ |
︙ | ︙ |
λ141=0 | a141=[ |
z1(k)=0.0475L1(k)+0.0386L2(k)+…+ 0.0796L141(k)
将降维后的输入向量送入ELM网络,输出为温度,隐含层节点数为20,训练ELM网络.目标涂层7组测试样本的仿真结果如图8所示.
图8
图8
目标涂层温度的ELM网络输出与实测值比较
Fig.8
Comparison of coating target temperature between estimated values obtained by ELM network and practical values
表6 涂层的测温结果以及相对误差
Tab.6
T1/℃ | T2/ ℃ | ||
---|---|---|---|
88 | 88.8148 | 0.93 | 0.8197 |
98 | 98.7474 | 0.71 | 0.8026 |
113 | 112.3211 | 0.60 | 0.7718 |
130 | 129.3978 | 0.46 | 0.6684 |
143 | 143.6489 | 0.45 | 0.5941 |
158 | 157.9227 | 0.04 | 0.5364 |
168 | 167.5493 | 0.27 | 0.5143 |
4 结语
本文提出一种未知目标发射率的红外辐射测温方法.由红外辐射基本定律,目标温度与光谱辐射亮度之间存在着强非线性关系,很难应用常规辨识建模方法.建立一种基于主元分析与极限学习机相结合的测温模型.首先,建立能够充分包含温度估计信息的高维输入向量;利用PCA实现输入降维,确定线性无关的主元输入变量;利用ELM网络自适应学习多光谱亮度与目标温度之间的潜在规律,建立目标温度测温模型,最终构建未知发射率情况下的目标红外多光谱测温方法框架.用黑体和未知发射率的隐身材料涂层分别作为测试目标,获得的测温结果与实际温度的相对误差都很小,说明了这种测温方法的合理性.
参考文献
基于多光谱测温优化的材料光谱发射率测量
[J].
Material spectral emissivity measurement optimized by multi-spectral temperature measured
[J].
Automatic compensation of emissivity in three-wavelength pyrometers
[J].
基于傅里叶红外光谱仪的光谱发射率测量装置的研制
[J].
Development of spectral emissivity measurement system based on Fourier transform infrared spectrometer(ftir)
[J].
基于标校的双波段比色测温法
[J].
Correction-based dual-waveband color comparison thermometric method
[J].
辐射测温的广义有效亮度温度
[J].
Generalized effective radiance temperature in radiation thermometry
[J].
基于双温双波段法的郎伯体红外测温技术
[J].
Infrared temperature measurement technology on lambertian based on the dual temperature and dual-band method
[J].
基于彩色CCD的双色与三色比色测温法比较研究
[J].
Comparative study on two-color and three-color colorimetric temperature measurement based on colored CCD
[J].
Spectral emissivity and temperature measurements of selective bodies using multiband fiber-optic radiometry
[J].DOI:10.1063/1.1784551 URL [本文引用: 1]
Temperature measurements of high-temperature semi-transparent infrared material using multi-wavelength pyrometry
[J].
Radiative temperature-measurement in the presence of attenuation due to brown smoke
[J].
Variables selection methods in near-infrared spectroscopy
[J].DOI:10.1016/j.aca.2010.03.048 URL [本文引用: 1]
基于遗传神经网络的多光谱辐射测温法
[J].
Multi-spectral thermometry based on GA-BP algorithm
[J].Considering some defects of back-propagation neural network (BP), a new algorithm combining genetic algorithm (GA) with BP was described. The application of GA-BP to the data processing of multi-spectral thermometry was proposed. The simulation experiments, based on GA-BP algorithm and BP neural network respectively, show that the recognition precision of trained emissivity samples is +/-5 K and +/-10 K respectively, and that of untrained emissivity samples is +/-10 K and +/-20 K respectively. No matter GA-BP algorithm or BP neural network is used, in general, the recognition precision of trained emissivity samples is higher than that of untrained emissivity samples. The recognition precision of true temperature is lower near the edge of sample sets. The GA-BP algorithm was more efficient than the BP neural network in the true temperature measurement.
红外辐射亮度的RBF网络建模及其光谱发射率估计
[J].
Modeling infrared radiance and calculating spectral emissivity based on RBF network
[J].
Research on pipeline blocking state recognition algorithm based on mixed domain feature and KPCA-ELM
[J].DOI:10.1504/IJCSM.2018.095498 URL [本文引用: 2]
A model-based patient selection tool to identify who may be at risk of exceeding dose tolerances during pancreatic SBRT
[J].DOI:10.1016/j.radonc.2019.09.016 URL [本文引用: 3]
Optimal principal component analysis of STEM XEDS spectrum images
[J].DOI:10.1186/s40679-019-0065-1 URL [本文引用: 2]
基于在线序列超限学习机和主成分分析的蒸汽冷却型燃料电池系统快速故障诊断方法
[J].
Fast fault diagnosis method of evaporatively cooled fuel cell system based on online sequential extreme learning machine and principal component analysis
[J].
/
〈 |
|
〉 |
