基于统计窄谱带模型的油料火焰辐射光谱模拟

彭吴迪, 刘礼喜, 陈志莉, 唐瑾, 陈林, 胡天佑, 王皓文

彭吴迪, 刘礼喜, 陈志莉, 唐瑾, 陈林, 胡天佑, 王皓文. 基于统计窄谱带模型的油料火焰辐射光谱模拟[J]. 红外技术, 2022, 44(3): 217-224.
引用本文: 彭吴迪, 刘礼喜, 陈志莉, 唐瑾, 陈林, 胡天佑, 王皓文. 基于统计窄谱带模型的油料火焰辐射光谱模拟[J]. 红外技术, 2022, 44(3): 217-224.
PENG Wudi, LIU Lixi, CHEN Zhili, TANG Jin, CHEN Lin, HU Tianyou, WANG Haowen. Oil Fire Radiation Calculation Based on a Statistical Narrow-Band Model[J]. Infrared Technology , 2022, 44(3): 217-224.
Citation: PENG Wudi, LIU Lixi, CHEN Zhili, TANG Jin, CHEN Lin, HU Tianyou, WANG Haowen. Oil Fire Radiation Calculation Based on a Statistical Narrow-Band Model[J]. Infrared Technology , 2022, 44(3): 217-224.

基于统计窄谱带模型的油料火焰辐射光谱模拟

基金项目: 

国家自然科学基金 21976043

桂林理工大学科研启动基金 GUTQDJJ20172017075

详细信息
    作者简介:

    彭吴迪(1998-),男,硕士研究生,主要从事环境遥感监测方面的研究。E-mail:1981552937@qq.com

    通讯作者:

    陈志莉(1971-),女,教授,博士生导师,主要从事环境遥感监测方面的研究。E-mail:zhilichen518@foxmail.com

  • 中图分类号: O433

Oil Fire Radiation Calculation Based on a Statistical Narrow-Band Model

  • 摘要: 近年来油料火灾污染事故频发危害性极大,通过分析光谱特性来提取火灾信息已成为研究油料火灾事故的重要途径。目前国内外学者已建立了多种气体辐射和炭黑辐射的模型对燃料燃烧进行研究,但少有对火焰光谱进行建模深入分析燃烧污染产物光谱特征信息。本文搭建了油料火焰光谱测试平台,测量了单一尺度下酒精、92号汽油、95号汽油和0号柴油的火焰光谱,以及多尺度下0号柴油的火焰光谱。实验结果表明3种油料的火焰光谱相似,随着尺度的增大辐射亮度呈非线性增大。基于统计窄谱带法(statistical narrow band,SNB)构建了油料火焰光谱辐射模型,通过实验数据验证曲线拟合度达0.895。利用该光谱辐射模型计算出油料火焰大尺度下的平均辐射亮度与平均透过率、不同烟气浓度下的平均透过率,能为遥感探测火灾污染及反演污染物浓度提供帮助。
    Abstract: In recent years, oil fire pollution accidents have occurred frequently and caused significant harm. It has become important to study oil fire accidents by extracting fire information through the analysis of spectral characteristics. Many domestic and foreign scholars have established a variety of gas radiation and carbon black radiation models to study fuel combustion; however, few scholars have directly modeled flame spectra to analyze and extract spectral characteristic information on combustion pollution products. In this study, a test platform for the flame spectra of oil was constructed, and the flame spectra of alcohol, 92 gasoline, 95 gasoline, and 0 diesel were measured at a single scale, as was the flame spectra of 0 diesel at multiple scales. The experimental results demonstrate that the flame spectra of the three oils are similar and the radiance increases nonlinearly with an increase in scale. Based on the statistical narrow-band method, a spectral radiation model for the oil flame was established, and a curve fitting degree of 0.895 was obtained based on experimental data. The spectral radiation model can be used to calculate the average radiance and transmittance of oil flames on a large scale, as well as the average transmittance at different flue gas concentrations, which can be helpful for remote fire pollution detection and pollutant concentration inversion.
  • 根据国际划分标准,一个国家或地区65岁及以上人口占比达到7%,就意味着进入了老龄化社会[1]。研究表明,独居尤其高龄老人成为意外伤害的高危人群,跌倒是导致老年人残疾,甚至死亡的主要原因[2]。因此研究室内安全监测算法显得十分重要。

    目前,国内外已有很多针对室内跌倒检测算法的研究,使用的设备可以分为摄像装置和传感装置。基于摄像装置的检测算法通过对采集到的视频帧逐帧分析,完成动作识别[3]。但是摄像装置不仅会暴露使用者的隐私,还极易受到光照的影响,在黑暗环境中性能会大大降低。传感装置可分为需穿戴传感器和无需穿戴传感器。需穿戴传感器检测算法通过穿戴加速度传感器[4]等设备采集源数据,研究摔倒时突然的加速检测摔倒,但是需要老人时刻穿戴着传感器装置,十分不便。无需穿戴的传感器主要有压力传感器、声学传感器、二值传感器和红外传感器。基于压力传感器[5]和声学传感器[6]的检测算法是将传感器放置于地面,提取不同动作声音或者对地板产生的压力信号检测摔倒,但在日常生活中,其他外物时常产生声音干扰,对地板也可能产生压力从而造成干扰,检测准确率较低。基于二值传感器的检测算法[7]通过多个二值传感器同时采集数据,研究数据波形的时域相关性,检测跌倒,但是需要先采集固定地点和时间段的日常动作数据作为参照,否则需要再次采取源数据,鲁棒性极低。

    为了解决上述问题,有人提出使用被动红外传感器。被动红外传感器通过接收外界发出的红外辐射进行工作,采集到的数据为探测区域的温度数据。此类传感器不受光强弱变化的影响,无需人为穿戴,同时保护了用户的隐私。日本的Shota Mashiyama团队[8]采用分辨率为8×8的Grid-eye红外传感,对日常动作和静止等5个动作进行分类,摔倒、行走和静止这3种能明显区分的动作识别率很高,但是行走和站立这两种相似度较高的动作识别率相对较低。Akira Hayashida团队[9]则考虑将红外传感器装置应用在不同的家庭场景中检测摔倒。该团队对比了不同室温、家庭环境和光照条件下的摔倒检测正确率,验证鲁棒性。但该算法需设定较多的阈值,阈值的选取很大程度上决定了算法的性能表现,同时大量的阈值选取算法也增加了算法的计算量。杨任兵等人[10]将8×8的红外传感器置于天花板设计分类算法,将温度和特征进行了结合,提高了摔倒的检测率,但是所提取的特征存在在不同动作的反映上有着相似性,所以一些坐下的动作会被错判为摔倒,虚警率较大,且8×8传感器的广角较小,探测区无法完全覆盖面积较大的家庭场景。王召军等人[11]采用HTPA 32×32红外传感器,提取5种特征对3种动作进行分类,且根据传感器的安装角度对特征数据进行校准,最终识别率较高,但是实际生活中无法精准得到夹角的准确值,对结果会有一定的影响。

    本文采用HEIMANN型号为HTPA80x64dR 1L5.0/1.0的红外阵列传感器提出一种识别弯腰、坐下、站起、行走和摔倒5种日常动作的方法。该方法不仅对活动人体进行了检测,也进行了前景提取,通过背景减除法提取前景,根据奈曼-皮尔逊准则确定初始最佳判决阈值,完成了人体区域提取。之后从时间、温度、形变和轨迹4个方面提取6个特征,最后对比了不同分类器的识别结果,也将本文方法与其他方法的结果进行了对比。

    本文人体动作识别方法流程如图 1所示,分为活动人体检测、前景提取、特征提取和分类4个部分。计算每帧数据的最大温度方差,若该方差大于阈值,则在探测区存在人体活动。在接收到温度分布数据之后,为了方便进行特征提取,先通过前景提取算法提取人体轮廓。活动人体检测和前景提取完成后,提取特征,利用K近邻(K-Nearest Neighbor,KNN)分类器,与数据库的训练集数据进行对比,判断接收到的数据的动作类别。

    图  1  动作识别方法流程图
    Figure  1.  Flowchart of the proposed activity recognition approach

    红外探测区无人活动时,每个像素点的温度值在时域上波动较小,而当人在探测区中活动时,时域温度波动较大。根据这一特点,可以通过提取每帧的最大温度方差来判断探测区是否有人体活动,以便进行后续操作。

    由于制作工艺和测试环境等原因,噪声严重影响到温度数据的准确性。A. A. Trfimova等人在文献[8]的基础上采用卡尔曼滤波[12]的方式对数据先进行预处理,但是卡尔曼滤波需要大量环境数据,计算量较大,无法实现实时的去噪。经过分析,测试数据的噪声近似服从高斯分布,因此采用大小为3×3的高斯滤波模板对测试数据预处理。

    aijk为第k帧的第(i, j)个像素点的温度值。第k帧每个像素点的方差计算公式为:

    $$ {v_{ijk}} = \frac{1}{\tau }\sum\limits_{l = k - (\tau - 1)}^k {{{({a_i}_{jl} - \overline {{a_i}_{jk}} )}^2}} , {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} k \geqslant \tau $$ (1)

    式中:τ为滑动窗口大小,本方法通过计算窗口τ为15帧内的每个像素点的方差来进行人体活动检测。$\overline {{a_i}_{jk}} $为第k帧第(i, j)个像素点前连续τ帧的平均温度,其公式如下:

    $$ \overline {{a_i}_{jk}} = \frac{1}{\tau }\sum\limits_{l = k - (\tau - 1)}^k {{a_{ijl}}} , {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} k \geqslant \tau $$ (2)

    经过公式(1),可以得到每一帧的温度方差分布矩阵:

    $$ {\boldsymbol{v}_k} = \left[ {\begin{array}{*{20}{c}} {{v_{{a_{11k}}}}}&{{v_{{a_1}_{2k}}}}& \cdots &{{v_{{a_{1\;{\text{80}}k}}}}} \\ {{v_{{a_{21k}}}}}&{{v_{{a_{22k}}}}}& \cdots &{{v_{{a_{2\;{\text{80}}k}}}}} \\ \vdots & \vdots &{}& \vdots \\ {{v_{{a_{{\text{64}}\;1k}}}}}&{{v_{{a_{{\text{64}}\;2k}}}}}& \cdots &{{v_{{a_{{\text{64}}\;{\text{80}}k}}}}} \end{array}} \right] $$ (3)

    根据之前的分析,当有人体在探测区活动时,探测区的温度会发生较大波动,于是通过公式(4)提取每帧的最大温度方差:

    $$ {v_k}_{\max } = \max ({\boldsymbol{v}_k}) $$ (4)

    vkmax大于阈值vth则说明探测区有人体活动。阈值vth的选取决定着活动人体检测的准确率和特征提取的正确性。由于采集数据的时间很长,环境温度和器件自身的温度都会有微小改变,所以本文确定阈值vth的方法为在探测区无人的环境下,实验前、实验中和实验后3个时间段中,运行实验装置直至最大帧数,各重复3次,共9组数据,分别计算9组数据的最大温度方差,选取最大值作为阈值vth。为了研究环境温度和传感器安装空间对阈值的影响,分别在校内实验室和家里客厅的不同环境温度下采集数据,进行对比实验,结果如表 1所示。

    表  1  不同环境温度、安装空间下的阈值对比
    Table  1.  Comparison of thresholds under different ambient temperatures and installation spaces
    Ambient temperature Laboratory Living room
    22℃ 1.66 1.65
    26℃ 1.71 1.73
    30℃ 1.75 1.76
    下载: 导出CSV 
    | 显示表格

    表 1可以看出,相比较而言,传感器安装空间对阈值变化影响很小,环境温度变化对它的影响更大一些,所以根据实验时的环境温度,本文将最终阈值设为1.7。

    在提取出的前景基础上进一步进行特征提取,可以极大提高识别准确率。常见的前景提取算法有背景减除法[13],帧间差分法[14],单高斯模型法[15]和混合高斯模型法等,这些算法往往采用经验值[16]设置阈值提取前景,本文则采用奈曼-皮尔逊准则设置阈值。

    像素点为背景点的状态设为H0,为前景点的状态设为H1。奈曼-皮尔逊准则无需已知信源先验概率与代价因子,该准则是在P(H1/H0)=α的约束条件下,检测概率达到最大。

    利用拉格朗日乘子μ构建目标函数:

    $$ J=P\left(H_{0} / H_{1}\right)+\mu\left[P\left(H_{1} / H_{0}\right)-\alpha\right] $$ (5)

    根据要求,需求目标函数J的最小值,将公式(5)转化为积分运算,且由于:

    $$ \int\limits_{{Z_1}} {P(x/{H_{\text{0}}})} {\text{d}}x + \int\limits_{{Z_{\text{0}}}} {P(x/{H_{\text{0}}})} {\text{d}}x = {\text{1}} $$ (6)

    式中:Z0Z1分别为判决区域,公式(5)可写为以下形式:

    $$ J = \mu (1 - \alpha ) + \int\limits_{{Z_0}} {[P(x/{H_1})} - \mu P(x/{H_{\text{0}}})]{\text{d}}x $$ (7)

    由公式(7)可知,若要使J达到最小,被积函数部分应该取负值,可得以下关系:

    $$ \frac{{P(x/{H_1})}}{{P(x/{H_0})}}\mathop \gtrless \limits_{{H_0}}^{{H_1}} \mu $$ (8)

    其中,判决门限μ可由约束条件得到:

    $$ P({H_1}/{H_0}) = \int\limits_{{Z_1}} {P(x/{H_0}){\text{d}}x = \int_\mu ^{ + \infty } {P(x/{H_0}){\text{d}}x = \alpha } } $$ (9)

    通过每帧数据得到的判决门限μ即为每帧的最佳判决阈值xkth。本文采用虚警率α为0.01约束下获得每帧的最佳判决阈值xkth,前景提取公式如下:

    $$ f({x_k}) = \left\{ \begin{array}{l} {x_k}(i, j)\;\;, {\text{if}}\;\;{x_k}(i, j) \geqslant {x_{k{\text{th}}}} \hfill \\ 0\;\;\;\;\;\;\;\;\;\;, {\text{else}} \hfill \\ \end{array} \right. $$ (10)

    若第k帧的第(i, j)像素点的温度值xk(i, j)大于阈值xkth,则温度值不变,否则置为0。该算法结果如图 2(a)所示,该方法极干净地滤除了背景点,只剩下少数背景游离点,经过形态学处理后,如图 2(b)所示,前景图像只剩下了完整的人体轮廓。

    图  2  前景提取算法结果
    Figure  2.  Foreground extraction results

    本文从时间、温度、形变和轨迹4个方向考虑,提取每个数据集的运动时间kact、最大温度方差Vmax、运动前后宽高比比值rAR、运动前后重心到脚距离比值rgravity、运动速度的最大值velmax和运动前后垂直角度差值diffv这6种特征。具体方法如下:

    1)运动帧数kact:每帧最大温度分布方差vkmax大于等于阈值vth时,则表明运动开始,该帧计为kstart,否则视为运动结束,计为kend。运动帧数的定义如下:

    $$ k_{\text {act }}=k_{\text {end }}-k_{\text {start }}+1 $$ (11)

    行走动作耗费的时间最长,摔倒是一种突发情况,运动帧数最短,如图 3所示。

    图  3  各动作的最大温度方差和运动持续帧
    Figure  3.  Maximum temperature variance and motion duration frames for each action

    2)每个数据集的最大温度方差Vmax:该特征指的是在k帧内,每帧最大温度分布方差vkmax中的最大值,即:

    $$ V_{{\rm{max}} }={\rm{max}} \left(v_{1 {\rm{max}} }, v_{2 {\rm{max}} }, \cdots, v_{k {\rm{max}} }\right), k>k_{\text {act }} $$ (12)

    不同动作发生时探测区温度变化的剧烈程度也不同。摔倒动作瞬时发生,在温度数据上则表现在其Vmax相比于其他动作较高,不同动作的vkmaxVmax对比如图 3所示。

    3)运动前后宽高比比值rAR:该特征为人体外接矩形宽高比在kstartkend时的比率,表现了人体在运动前后的形变特征,其公式如下:

    $$ {r_{{\text{AR}}}} = \frac{{{r_{{\text{end}}}}}}{{{r_{{\text{start}}}}}} $$ (13)

    式中:rendrstart分别为运动结束帧和运动开始帧的宽高比。人体外接矩形如图 4所示。

    图  4  人体外接矩形
    Figure  4.  Human body circumscribed rectangle

    图中的人体外接矩形左上点akl=(ikl, jkl),人体外接矩形右下点akr=(ikr, jkr),l和r分别代表着第k帧外接矩形的左边界和右边界,于是人体外接矩形的宽高比可表示为:

    $$ r = \frac{{j_k^{\text{r}} - j_k^{\text{l}} + 1}}{{i_k^{\text{r}} - i_k^{\text{l}} + 1}} $$ (14)

    摔倒发生后人体的外接矩阵宽高比明显要远大于运动前的宽高比,因此摔倒动作的rAR是远大于1。弯腰动作后人体的宽度增加,高度降低,所以其rAR次之。

    4)运动前后重心到脚距离比值rgravity:表示运动前后人体重心到脚的距离的比率,与rAR一样表现了人体在运动前后的形变特征。人体区域的重心表示为xk=(ik, jk),其中:

    $$ {i_k} = \frac{{\sum\limits_{m = l}^r {\sum\limits_{n = b}^t {i{Q_k}(n, m)} } }}{{\sum\limits_{m = l}^r {\sum\limits_{n = b}^t {{Q_k}(n, m)} } }} $$ (15)
    $$ {j_k} = \frac{{\sum\limits_{m = l}^r {\sum\limits_{n = b}^t {j{Q_k}(n, m)} } }}{{\sum\limits_{m = l}^r {\sum\limits_{n = b}^t {{Q_k}(n, m)} } }} $$ (16)

    式中:tbl、和r分别为人体外接矩形的上下左右边界的坐标值,Qk(n, m)的定义如下:

    $$ {Q_k}(n, m) = \left\{ \begin{gathered} {\text{1}}\;\;, \;{\text{if}}\;\;{x_k} \ne 0 \hfill \\ 0\;\;, \;{\text{if}}\;\;{x_k} = 0 \hfill \\ \end{gathered} \right. $$ (17)

    那么重心到人脚的距离Ih为:

    $$ l_{\mathrm{h}}=i_{k}-i_{k}^{\mathrm{b}} $$ (18)

    式中:ikb为人体外接矩阵下边界的坐标值,运动前后重心到脚距离比值rgravity的定义为:

    $$ r_{\text {gravity }}=l_{\text {hend }} / l_{\text {hstart }} $$ (19)

    摔倒动作的lhend远小于lhstart,所以摔倒动作的rgravity远小于1,而行走运动前后人外接矩形并没有发生太大的变化,所以行走动作的rgravity是在1上下波动。

    5)运动速度的最大值velmax:将重心点xk=(ik, jk)作为追踪点,时间窗口τ等于传感器本身的帧率,计算重心点τ帧的欧氏距离,从而得到运动的速度,其公式如下:

    $$ {\text{ve}}{{\text{l}}_k} = \frac{{\sqrt {{{({i_k} - {i_{k - \tau + 1}})}^2} + {{({j_k} - {j_{k - \tau + 1}})}^2}} }}{\tau } $$ (20)

    式中:ikjk分别为第k帧重心点的坐标。通过公式(20)这样就可以得到动作的速度波动情况,将最大值作为特征:

    $$ {\text{ve}}{{\text{l}}_{\max }} = \max ({\text{ve}}{{\text{l}}_k}) $$ (21)

    摔倒动作是瞬时发生的,相比于其他动作有着极高的velmax,而行走动作相比于其他动作较为稳定,且速度稍高,所以行走动作的velmax仅次于摔倒动作。

    6)运动前后垂直角度差值diffv:该特征从不同动作的轨迹变化出发,将人体区域的重心点xk=(ik, jk)作为追踪点,由于传感器放置在墙壁之上,若将探测场的画面作为平面,在三维空间中,传感器则位于画面的下边界上。若人体在下边界重心点所在的垂直线上进行站立坐下动作,则垂直角度的变化均为0,无法判别动作,所以将整张图像左下角点,即坐标为(63, 0)的点作为顶点,与重心点形成角度,观察运动前后垂直角度的变化,如图 5所示。

    图  5  垂直角度示意图
    Figure  5.  Vertical angle diagram

    图 5可得垂直角度的定义如下:

    $$ \tan \theta = \frac{{{\text{63}} - {i_k}}}{{{j_k}}} $$ (22)

    则运动前后的垂直角度差值定义为:

    $$ {\text{dif}}{{\text{f}}_v} = {\theta _{{\text{end}}}} - {\theta _{{\text{start}}}} $$ (23)

    在摔倒动作发生前后,人的重心点会更加贴近于画面下边界,所以垂直角度会有极大的降低,摔倒动作的diffv较大且为负。

    本文将KNN算法、支持向量机(Support Vector Machine, SVM)算法、随机森林(Random Forest, RF)算法、决策树(Decision Tree, DT)算法和神经网络(Neural Network, NN)算法的结果进行对比。KNN算法是给定测试样本,基于距离度量找出特征空间中距离最靠近该样本的K个训练样本,选择这K个样本中出现最多的类别作为预测结果。距离度量和K值的选择均会影响到分类的准确性,距离度量会影响到“邻居”的选择,而K值如果过小容易产生过拟合的现象,如果过大则会大大增加算法的计算时间。因此本文采用交叉验证法[17]来确定K值和距离度量。

    本实验使用德国HEIMANN型号为HTPA80x64 dR1L5.0的热电堆阵列传感器,实验环境选择室内普通实验室,背景较为单一,室温为26℃,将传感器置于墙壁,距离地面2.6 m。测试人员共10名(7男3女),被测人员在探测区依次完成弯腰、站起、坐下、摔倒和行走5种动作,每种动作重复30次,最终每个动作获得300组数据,数据集一共1500组,按照动作将数据集等比例划分为训练集和测试集进行实验。

    本文采用交叉验证的方法确定KNN算法最优的K值和距离度量函数。首先将数据集参考动作类别按照一定比例划分为训练集和测试集。距离度量也在欧氏距离和曼哈顿距离两者中进行选择。K值以3为起始值,2为步长,分别在距离度量为欧氏距离和曼哈顿距离的条件下,依次选取不同的K值进行分类计算相对应的准确率。交叉验证结果如表 2表 3所示,选择曼哈顿距离作为距离度量,K取3进行实验。

    表  2  欧氏距离交叉验证结果
    Table  2.  Euclidean distance cross-validation results
    Test: Train K=3 K=5 K=7 K=9 K=11 K=13 K=15
    1:1 96.67% 96.8% 96.4% 96% 95.6% 95.2% 95.2%
    1:4 97.33% 97.67% 97.33% 97% 96.67% 96.33% 96%
    1:9 94.67% 95.33% 95.33% 95.33% 94% 93.33% 93.33%
    下载: 导出CSV 
    | 显示表格
    表  3  曼哈顿距离交叉验证结果
    Table  3.  Manhattan distance cross-validation results
    Test: Train K=3 K=5 K=7 K=9 K=11 K=13 K=15
    1:1 97.73% 97.73% 97.20% 96.53% 96.93% 96.53% 96.4%
    1:4 98% 97.67% 97.67% 97.67% 96.67% 96% 96%
    1:9 96.67% 96% 96% 96% 95.33% 95.33% 94.68%
    下载: 导出CSV 
    | 显示表格

    同样采用交叉验证对余下分类器的实验参数进行设定。经过验证,SVM分类器惩罚参数c最终选择为2,径向基函数的自带参数gamma,实验对比后取0.08,训练集780组,测试集720组。DT分类器采用信息增益作为分裂节点的评价指标,最优分裂作为拆分策略,剪枝时树的深度设为4,叶节点处所需最少样本数设为1,拆分最小样本数为3,样本权重相同,训练集900组,测试集600组。RF分类器中,最优结果下随机森林建立10棵决策树,深度为10,余下参数为默认参数,训练集975组,测试集525组。为了进行下一阶段研究,本文选择简单的神经网络结构构建分类器与以上几种分类器进行结果对比。NN分类器的结构为输入层,中间层以及输出层,输入为提取的特征,输出为动作类别,中间层为8层全连接层,使用sigmoid函数作为激活函数,网络训练总次数设定为500轮,依次改变学习率,最终最优学习率为0.01,训练集1125组,测试集375组。

    将KNN算法、SVM算法、RF算法、DT算法和NN算法的识别准确率进行对比,结果如表 4所示。

    表  4  5种分类算法结果对比
    Table  4.  Comparison of the results of five classification algorithms
    KNN SVM RF DT NN
    Bend 96.7% 92.8% 96.6% 97.9% 97.5%
    Stand 98.3% 97.3% 95.4% 96.3% 90.9%
    Sit 91.7% 75.8% 88.7% 91.3% 83.3%
    Fall 100% 99.1% 98.2% 96.3% 100%
    Walk 98.3% 95.6% 97.4% 96.7% 95.4%
    Accuracy 97% 92.12% 95.26% 95.7% 93.42%
    下载: 导出CSV 
    | 显示表格

    表 4可得KNN算法平均识别准确率高于其他算法,摔倒动作的识别准确率达到了100%,其他动作的识别准确率也在97%以上,最终选择KNN算法。KNN算法识别准确率的混淆矩阵图如图 6所示。

    图  6  准确率混淆矩阵图
    Figure  6.  Accuracy confusion matrix

    图 6可以发现,坐下和弯腰这两种有着下降趋势的动作易被误判,与其他动作相比,坐下仅在运动前后宽高比比值这一特征上与弯腰表现差异较大,所以坐下动作有8%的测试集数据错判为弯腰,弯腰动作有3%的测试集数据错判为坐下,误判率较高。为了进一步验证本文方法性能,将本文方法与文献[10]和文献[11]的方法的结果进行了对比,结果如表 5所示。

    表  5  3种方法结果对比
    Table  5.  Comparison of the results of three algorithms
    Ours Method in paper [10] Method in paper [11]
    Bend 96.7% 88.7% 94.7%
    Stand 98.3% 87.3% 93.5%
    Sit 91.7% 78.8% 90.5%
    Fall 100% 97.6% 96.7%
    Walk 98.3% 95.4% 100%
    Accuracy 97% 89.6% 95.1%
    下载: 导出CSV 
    | 显示表格

    表 5可以看出,本文方法准确率有了很大的提高,并且很好地区分了站起和坐下动作,这说明本文提取到的特征确实与人体动作相关联,能以较高的准确率区分日常的生活动作。

    本文提出了一种基于低分辨率红外传感器的人体动作识别方法。本文方法对活动人体进行了检测,也对前景提取进行了研究,根据得到的数据,采用背景减除法提取前景,基于奈曼-皮尔逊准则确定初始最佳判决阈值,使背景错判为前景的出错的概率达到最小。在以上工作基础上,从时间、温度、形变和轨迹4个方面提取6个相关特征,比较了不同分类器的结果,最终选用KNN分类器对弯腰、站起、坐下、摔倒和行走这5个日常动作进行分类,平均准确率达到97%。本文方法对摔倒动作有着100%的识别准确率,相比于其他方法,站起和坐下这两种易错分的动作的识别准确率也有提升,达到了初期构想的要求。未来,可采用多个红外传感器搭建系统,以提高细小差别动作的识别率。本文提出的方法有着较繁琐的计算机视觉方面的数据处理,避免不了阈值的设计,且计算量较大,今后可结合深度学习,搭建神经网络以进行更深入的研究。

  • 图  1   油料火焰光谱测试平台

    Figure  1.   Oil flame spectrum testing platform

    图  2   柴油与酒精的火焰光谱对比

    Figure  2.   Flame spectral comparison of diesel and alcohol

    图  3   6 cm尺度下柴油、92号汽油、95号汽油火焰光谱对比

    Figure  3.   Flame spectral comparison of diesel oil, No. 92 gasoline and No. 95 gasoline at the scale of 6cm

    图  4   6 cm、14 cm、18 cm和22 cm四种尺度下柴油火焰光谱对比

    Figure  4.   Spectral comparison of diesel flame at 6 cm, 14 cm, 18 cm and 22 cm scales

    图  5   柴油火焰光谱实测值与SNB模拟值对比

    Figure  5.   Comparison of measured values of diesel flame spectrum with simulated values of SNB

    图  6   实验数据、逐线法模拟值、SNB模拟值对比

    Figure  6.   Comparison of experimental data, line by line simulation value and SNB simulation value

    图  7   SNB模拟不同燃烧尺度、不同烟气浓度下火焰平均透过率

    Figure  7.   SNB simulates average flame transmittance at different combustion scales and flue gas concentrations

    图  8   SNB模拟不同尺度下火焰光谱辐射亮度

    Figure  8.   SNB simulates flame spectral luminance at different scales

    表  1   92号汽油,95号汽油,0号柴油,酒精燃烧60 s时的组分体积分数与温度

    Table  1   Volume fraction and temperature of gasoline 92, gasoline 95, diesel 0, alcohol in 60 s combustion

    Type 92#gasoline 95#gasoline 0#diesel Alcohol
    CO 0.0067% 0.0066% 0.0074% 0.0047%
    CO2 0.0821% 0.0892% 0.0854% 0.0823%
    T/K 1155 1165 1202 977
    下载: 导出CSV

    表  2   LBL, SNB和实验数据在2.5 μm, 2.7 μm, 4.3 μm, 5.5 μm,6.3 μm特征峰处辐射亮度比较

    Table  2   LBL, SNB and experimental data were compared at characteristic peaks of 2.5, 2.7, 4.3, 5.5, 6.3 μm

    Spectral radiance/(W·cm-2·μm-1·sr-1) 2.5 μm 2.7 μm 4.3 μm 5.5 μm 6.3 μm
    LBL 0.1236 0.2358 0.2936 0.0261 0.0142
    SNB 0.1194 0.2245 0.2956 0.0381 0.0206
    Experiment 0.1466 0.1976 0.3001 0.0299 0.0178
    下载: 导出CSV

    表  3   SNB, LBL拟合评价指标分析

    Table  3   Analysis of SNB and LBL fitting evaluation indexes

    Evaluation RMSE MAE R2
    LBL-experiment 0.0113 0.00769 0.9486
    SNB-experiment 0.0162 0.0109 0.8950
    SNB-LBL 0.0215 0.0129 0.7760
    下载: 导出CSV
  • [1] 周锋, 阮桢. 石油化工储罐火灾事故统计与分析[C]//2016中国消防协会科学技术年会论文集, 2016: 30-32.

    ZHOU Feng, RUAN Zhen. Statistics and analysis of fire accidents in petrochemical storage tanks[C]//Proceedings of 2016 China Fire Protection Association Science and Technology Annual Meeting, 2016: 30-32.

    [2] 刘洪涛, 陈志莉, 刘强, 等. 基于小波分解的油料火焰光谱特性分析研究[J]. 光学学报, 2016, 36(1): 320-326. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB201601040.htm

    LIU Hongtao, CHEN Zhili, LIU Qiang, et al. Analysis of oil flame spectrum characteristics based on wavelet decomposition[J]. Acta Optica Sinica, 2016, 36(1): 320-326. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB201601040.htm

    [3]

    Wakatsuki K, Fuss S P, Hamins A, et al. A technique for extrapolating absorption coefficient measurements to high temperatures[J]. Proceedings of the Combustion Institute, 2005, 30(1): 1565-1573. DOI: 10.1016/j.proci.2004.08.078

    [4]

    RAJ P K. LNG pool fire spectral data and calculation of emission power[J]. Journal of Hazardous Materials, 2007, 140(3): 720-729.

    [5]

    RAJ P K. Large hydrocarbon fuel pool fires: Physical characteristics and thermal emission variations with height[J]. Journal of Hazardous Materials, 2007, 140(1): 280-292. https://www.sciencedirect.com/science/article/abs/pii/S0304389406009885

    [6]

    Maoult Y L, Sentenac T, Orteu J J, et al. Fire Detection: A New Approach Based on a Low Cost CCD Camera in the Near Infrared[J]. Process Safety and Environmental Protection, 2007, 85(3): 193-206. DOI: 10.1205/psep06035

    [7] 徐建楠, 蒋新生, 张昌华, 等. 不同工况下汽油蒸气爆炸着火延迟与机理分析[J]. 化工学报, 2019, 70(1): 398-407. https://www.cnki.com.cn/Article/CJFDTOTAL-HGSZ201901045.htm

    XU Jiannan, JIANG Xinshen, ZHANG C H, et al. The ignition delay and mechanism analysis of gasoline vapor explosion under different working conditions[J]. CIESC Journal, 2019, 70(1): 398-407. https://www.cnki.com.cn/Article/CJFDTOTAL-HGSZ201901045.htm

    [8]

    Boudreau S, Levasseur S, Perilla C, et al. Chemical detection with hyperspectral lidar using dual frequency combs[J]. Optics Express, 2013, 21(6): 7411-8. DOI: 10.1364/OE.21.007411

    [9]

    Klassen M S, Sutula J A, Holton M M, et al. Transmission throughand breakage of single and multi-pane glazing due to radiant exposure: state of research[J]. Fire Technology, 2010, 46(4): 821-832. DOI: 10.1007/s10694-010-0150-4

    [10]

    Grosshandler W L. The effect of soot on pyrometric measurements of coal particle temperature[J]. Combustion & Flame, 1982, 55(1): 59-71. https://www.sciencedirect.com/science/article/pii/0010218084901494

    [11]

    Rothman L S, Gordon I E, Barber R J, et al. HITEMP, the high-temperature molecular spectroscopic database[J]. Journal of Quantitative Spectroscopy & Radiative Transfer, 2010, 111(15): 2139-2150. https://www.sciencedirect.com/science/article/pii/S002240731000169X

    [12]

    Malkmus W. Random lorentz band model with exponential-tailed S-1 line intensity distribution function[J]. Optical Society of America, 1967, 57(3): 323-329. DOI: 10.1364/JOSA.57.000323

    [13]

    CHU H Q, LIU F S, ZHOU H, Calculations of gas thermal radiation transfer in one dimensional planar enclosure using LBL and SNB models[J]. Heat Mass Transfer, 2011, 54(21): 4736-4745.

    [14]

    Hamins A, Anderson D T, Miller J H. Mechanistic studies of toluene destruction in diffusion flames[J]. Combustion Science & Technology, 1990, 71(4-6): 175-195.

    [15]

    Gordon I E, Rothman L S, Hill C, et al. The HITRAN2016 molecular spectroscopic database[J]. Journal of Quantitative Spectroscopy & Radiative Transfer, 2017, 130(11): 4-50. https://www.sciencedirect.com/science/article/pii/S0022407317301073

    [16]

    Elsasser W M. Heat Transfer by Infrared Radiation in the Atmosphere[M]. Harvard: Harvard University Press, 1943: 43-45.

    [17]

    Goody R M. A statistical model for water vapour absorption[J]. Quarterly Journal of the Royal Meteorological Society, 1952, 78(336): 165-169. DOI: 10.1002/qj.49707833604

    [18]

    Soufiani A, Taine J. High temperature gas radiative property parameters of statistical narrow band model for H2O, CO2 and CO and correlated k model for H2O and CO2[J]. Heat Mass Transfer, 1997, 40(4): 987–991. DOI: 10.1016/0017-9310(96)00129-9

  • 期刊类型引用(6)

    1. 王龙飞,胡雨旺,张泽光,刘悦,薛常喜. 非球面光学元件的快速制造技术. 光电工程. 2024(01): 3-34 . 百度学术
    2. 刘斌,王孝坤,程强,刘忠凯,李凌众,蔡梦雪,张洛嘉,李文涵,王金成,张学军. 复杂曲面光学元件高精度面形检测技术概述. 南通大学学报(自然科学版). 2024(01): 1-27 . 百度学术
    3. 席建普,周卓林,罗来华,张通通,聂子凯. 自由曲面轮廓仪机架多目标优化设计. 制造技术与机床. 2024(05): 122-127 . 百度学术
    4. 简恩霖,肖红军,张芷芊,杨皓桦,吴南寿,曾亚光. 谱域OCT技术在非球面镜面型测量中的应用. 现代信息科技. 2024(09): 58-61 . 百度学术
    5. 叶卉,王泽华,许世炜,姜晨. 基于便携式粗糙度仪的表面轮廓测量实验平台设计. 实验技术与管理. 2024(07): 140-147 . 百度学术
    6. 王兆明,栗孟娟,于秋跃,李春林,赵子跃,王经华,吕天斌,张兆健,于长锁. 两面共体非球面反射镜光轴一致性高精度测量方法研究(特邀). 红外与激光工程. 2023(09): 39-46 . 百度学术

    其他类型引用(2)

图(8)  /  表(3)
计量
  • 文章访问数:  158
  • HTML全文浏览量:  72
  • PDF下载量:  32
  • 被引次数: 8
出版历程
  • 收稿日期:  2021-03-07
  • 修回日期:  2021-04-24
  • 刊出日期:  2022-03-19

目录

/

返回文章
返回