2. 武汉大学 遥感信息工程学院, 武汉 430079;
3. 北方民族大学 计算机科学与工程学院, 银川 750021
2. School of Remote Sensing and Information Engineering, Wuhan University, Wuhan 430079, China;
3. School of Computer Science and Engineering, North Minzu University, Yinchuan 750021, China
维持公共区域安全是保障社会和谐稳定和人民群众生命财产安全的一项重要工作。毫米波图像方法被广泛应用于机场和地铁站等公共区域的安检威胁识别任务中[1-2],该方法应用于威胁识别过程中时, 要求精度好、效率高和误报率低[3]。然而, 在威胁识别过程中采集到的毫米波图像存在形状、大小和位置未知的隐藏物体, 导致图像处理方法存在较高的误报率。毫米波图像包括有源毫米波图像和无源毫米波图像2种类型。有源毫米波图像通过区分目标物体周围辐射强度的差异来生成扫描图像[4-6], 能够以相对较快的速度生成图像, 但该方法生成图像的分辨率较低且识别精度较差。无源毫米波图像通过解析目标的毫米波能量来生成扫描图像, 可生成高分辨率的图像, 适用于公共场所的高精度威胁识别。文献[7]提出一种改进的双基址结构, 有效减少了扫描时间, 同时提高生成图像的质量。目前,关于无源毫米波图像的研究常采用超分辨率方法来提高图像分辨率[8-10], 以提高识别准确率。文献[11]提出一种基于高斯混合模型的方法, 并将期望最大化算法用于分类阶段,虽然其结果准确性有所提高, 但该方法会产生无关联区域。文献[12]提出一种双重阶段算法, 该算法的准确率高, 但是过程复杂且耗时长, 难以投入到实际应用中。文献[13]提出一种多重阶段算法, 通过优化图像处理速度, 有效减少算法耗时, 且可达到较高的准确率。
近年来,卷积神经网络(Convolutional Neural Network,CNN)已被广泛应用于识别、分类和分割等任务[14-16],且CNN具有通过非线性映射来区分数据独立特征的性能[17-19],因此在图像识别领域的应用较为广泛。文献[20]使用CNN-DL分类器检测毫米波图像中对象的可行性,并指出提高毫米波图像识别精度的2个重要因素:模型的深度和图像的完整性。由于训练样本较少、样本间类别不平衡问题显著,目前深度学习方法在威胁识别领域中的应用仍然较少。针对扫描对象的多角度视图, 多视图间的特征信息联系显得尤为重要。文献[21]并未考虑多视图间的直接联系,导致存在较高的漏检率。因此,在威胁识别任务中存在隐藏物体形状、大小和位置未知的问题,且不是每个角度都能检测到危险物品,多视图角度下的区域重叠问题将会导致模型的误报率增大。与此同时,在实际安检过程中,由于多数检测对象都为无威胁样本,仅有少数样本为危险样本,从而引起数据不均衡问题。
针对以上问题和挑战, 本文构建一种基于多视图的深度神经网络模型PSS-Net用于毫米波图像的威胁识别。该模型针对扫描图像多视角问题进行多视图建模, 从而增强各区域在不同角度之间的关联程度, 增强扫描对象的整体表征能力,并使用稠密性连接对单一视图的表征能力进行增强。与此同时, 该模型运用注意力机制以及焦点损失函数解决扫描对象样本中存在的数据不均衡问题。
1 面向威胁识别的多视图深度神经网络 1.1 简述针对输入的人体扫描对象生成16个视图, 将每个视图划分为17个区域, 随后各视图数据分别发送至残差连接卷积神经网络进行特征提取, 且将通过注意力机制提取的多视图卷积特征作为时序特征, 通过基于稠密连接的长短期记忆模型强化威胁区域的特征表达。基于焦点损失函数优化整个网络, 构成一个端到端的可训练框架, 从而对单个扫描对象的各个区域进行危险预测。本文方法的应用场景如图 1所示。
|
Download:
|
| 图 1 威胁识别人物场景 Fig. 1 Threat identification character scene | |
图像序列建模[22-23]问题使得基于毫米波图像的威胁识别任务更加困难。图像序列建模[24]基于各扫描对象通常会生成多角度视图。由于图像部分区域重叠, 使得识别难度增大, 在这种情况下, 多视图[23]之间的关联性显得尤为重要。与此同时, 长短期记忆(Long Short Term Memory, LSTM)网络[25]在处理序列问题上具有广泛应用, 可以较好地建立序列之间的关系。根据以上问题与解决方案, 构建一个多视图卷积神经网络模型架构, 如图 2所示。其中, 每个视图被发送到参数相同的预训练ResNet-50[26]网络中,以获取像素之间的关系并提取空间位置上的特征。此外, 在最终特征图上, 该模型采用一种基于稠密连接的长短期记忆注意力机制,特征映射由具有不同步幅和内核大小的卷积层处理。考虑到多视图信息序列的处理问题, 通过LSTM对扫描对象的多个视图进行空间邻域上的扫描,并对扫描信息建模, 从而对各个视图在不同扫描角度下的特征进行关联性建模, 提高区域的表征能力。PSS-Net模型架构中的参数设置如表 1所示。
|
Download:
|
| 图 2 PSS-Net模型架构 Fig. 2 Architecture of PSS-Net model | |
|
下载CSV 表 1 PSS-Net模型架构中的参数设置 Table 1 Parameter setting in PSS-Net model architecture |
利用模型对一个样本中的17个区域分别进行危险概率的预测, 即对于N个样本, 模型将输出17×N个概率。通常情况下, 模型考虑使用交叉熵损失函数作为目标函数对模型进行优化, 其表示方法如式(1)所示:
| $\begin{aligned} \text { Loss }=&-\frac{1}{17 N} \sum\limits_{n=1}^{N} \sum\limits_{i=1}^{17}\left[y_{i}^{(n)} \lg\;\hat{y}_{i}^{(n)}+\right.\\ &\left.\left(1-y_{i}^{(n)}\right) \lg \left(1-\hat{y}_{i}^{(n)}\right)\right] \end{aligned} $ | (1) |
其中, 当样本存在威胁时, yi(n)为1, 否则为0, ŷi(n)代表第n个扫描样本的第i个区域的预测结果。
由于在毫米波图像威胁识别任务中存在危险的样本数量较少, 从而带来数据不均衡问题, 导致模型的识别召回率较低。因此对损失函数进行加权, 以缓解样本不均衡带来的影响, 具体如式(2)所示:
| $\begin{aligned} \text { Loss }=&-\frac{1}{17 N} \sum\limits_{n=1}^{N} \sum\limits_{i=1}^{17}\left[\alpha y_{i}^{(n)} \lg \hat{y}_{i}^{(n)}+\right.\\ &\left.(1-\alpha)\left(1-y_{i}^{(n)}\right) \lg \left(1-\hat{y}_{i}^{(n)}\right)\right] \end{aligned} $ | (2) |
其中, α为样本调节因子。
上述问题仅解决了样本不均衡问题, 并未考虑难分类样本在任务中的重要性。文献[27]提出一种焦点损失函数, 可有效解决数据不均衡问题。在收敛速度和性能方面, 焦点损失函数优于先前的损失函数, 具体如式(3)所示:
| $\begin{aligned} \text { Loss }=&-\frac{1}{17 N} \sum\limits_{n=1}^{N} \sum\limits_{i=1}^{17}\left[\alpha y_{i}^{(n)}(1-)^{\gamma} \lg \hat{y}_{i}^{(n)}+\right.\\ &\left.(1-\alpha)\left(1-y_{i}^{(n)}\right)^{\gamma} \lg \left(1-\hat{y}_{i}^{(n)}\right)\right] \end{aligned} $ | (3) |
其中, γ与α为可调节的超参数(预设γ=2.0, α=0.750)。
PSS-Net模型选取焦点损失函数为最终损失函数, 并对模型进行优化。
1.4 注意力机制注意力机制可增强图像全局信息的表达, 降低局部信息对结果的影响, 使得连接特征映射更加紧凑。CNN注意力机制通过全连接层进行特征映射变换, 考虑了单个视图下的图像特征。LSTM注意力机制通过增加目标的权重并强调重要的区域, 模拟人类观察多角度物体的过程, 从而有效减少信息损失。LSTM注意力机制根据M(本文设置为16)个图像的全局信息重要性对每个扫描视图进行加权, 并输入序列x=(x1, x2, …, xM), LSTM注意力机制计算输出序列zi=(z1, z2, …, zM), 具体如式(4)所示:
| $z_{i}=\sum\limits_{j=1}^{M} \beta_{i j} x_{i} $ | (4) |
采用Softmax函数计算每个权重系数βij, 具体如式(5)所示:
| $\beta_{i j}=\frac{\exp \left(W x_{i j}\right)}{\sum\limits_{k=1}^{M} \exp \left(W x_{i k}\right)} $ | (5) |
其中, βij表示第j个视图的注意力权重, Wxij是线性变换层, 输入维度为768, 输出维度为17。
2 实验结果与分析 2.1 数据集本文方法旨在预测指定身体区域(来自17个全身区域)存在威胁的概率, 每次扫描包含17个标签, 对应于17个身体区域中是否存在威胁。
HD-AIT毫米波人体威胁扫描数据集中的每个人类扫描对象包含16个视图, 且每个视图被分成17个身体区域, 具体如图 3所示。该数据集包含由新一代毫米波扫描仪获得的11 470个扫描对象。不同对象穿着不同服装类型(含夏季服装和冬季服装), 具有不同的体重指数, 含有不同的威胁数量和不同类型的威胁。
|
Download:
|
| 图 3 人体扫描图像身体区域分布 Fig. 3 Body area distribution in human body scanning images | |
实验选择60%数据作为训练集, 20%数据作为验证集, 20%数据作为测试集。对于焦点损失函数, 预设γ=2.0, α=0.750, 并使用余弦学习率衰减策略对模型进行训练, 初始学习率设置为0.1。基于Ubuntu 16.04 LTS系统, CPU采用36核2.6 GHz Intel Core i9-9870XE, 内存大小为94 GB。使用4块GeForce GTX 1080Ti型号GPU用于加速训练, 总共训练150个Epoch, 训练共计120 h。PSS-Net模型在PyTorch框架上训练, 实验过程中模型在训练集和验证集上的损失可视化如图 4所示。从图 4可以看出, 当训练到80个~100个Epoch时, 损失函数已经降低到相对较低水平且基本收敛。
|
Download:
|
| 图 4 训练集与验证集的损失函数 Fig. 4 Loss function of training set and verification set | |
本文以准确率(Accuracy)、精确度(Precision)、召回率(Recall)、特异性(Specificity)和假阳性率(False Positive Rate, FPR)作为评估指标, 计算方法分别如式(6)~式(10)所示:
| ${\rm{Accuracy}} = ({\rm{TP}} + {\rm{TN}})/({\rm{TP}} + {\rm{TN}} + {\rm{FN}} + {\rm{FP}}) $ | (6) |
| ${\rm{Precision}} = {\rm{TP}}/({\rm{TP}} + {\rm{FP}}) $ | (7) |
| ${\rm{Recall}} = {\rm{TP}}/({\rm{TP}} + {\rm{FP}}) $ | (8) |
| ${\rm{Specificity}} = {\rm{TN}}/({\rm{TN}} + {\rm{FN}}) $ | (9) |
| ${\rm{FPR}} = 1.00 - {\rm{Specificity}} $ | (10) |
其中, TP表示非威胁样本中检测为无威胁样本的个数, FP表示非威胁区域中误检为威胁样本的个数, FN表示威胁样本中未能成功检测出威胁的个数, TN表示威胁样本中正确检测为威胁样本的个数。
2.3 测试结果由于损失函数中的超参数对实验结果有不同的影响, 实验研究了不同参数的α对焦点损失函数的性能影响, 结果如表 2所示。从表 2可以看出, 随着α的增大, 焦点损失函数的准确率和召回率均呈现先增大后降低的趋势, 当α为0.750时达到最优。实验进一步研究了不同的γ与α组合对焦点损失函数的性能影响, 结果如表 3所示, 在实验过程中, 当测试γ时, 对应的α均取实验中的最优值。从表 3可以看出, 当γ=2.0时性能达到最优, 因此实验选择γ=2.0, α=0.750作为损失函数的超参数。
|
下载CSV 表 2 超参数α对焦点损失函数的性能影响(γ=0) Table 2 Effect of super parameter α on the performance of focal loss function(γ=0) |
|
下载CSV 表 3 超参数α与γ对焦点损失函数的性能影响 Table 3 Effect of super parameter α and γ on the performance of focal loss function |
为了有效对比各种模型在威胁识别方面的性能, 实验对文献[26-28]中的基线模型与本文模型进行性能对比, 结果如表 4所示。从表 4可以看出, 本文提出的PSS-Net模型在测试集上的表现均高于其他基线模型, 且该任务主要识别危险图像。同时,PSS-Net模型在召回率上相比其他方法提升更多。
|
下载CSV 表 4 4种模型的测试结果对比 Table 4 Comparison of test results of four models |
每个区域的测试结果如表 5所示。由表 5可知, 误判集中在区域15(右脚踝)和区域16(左脚踝)中, 误判主要源于人体生理结构的差异, 考虑到数据集中多数样本不存在局部病变的情况, 因此当目标身体存在畸变时将干扰判断, 尤其是过于严重的生理形状病变可能导致误判。
|
下载CSV 表 5 各区域的测试结果 Table 5 Test results of each region |
为比较模型中不同模块对于性能的贡献, 本文设计了消融实验进行研究, 分别比较了LSTM、多视图卷积神经网络、注意力机制和焦点损失函数模块的测试结果, 如表 6所示。其中, 实验1、实验2分别仅使用以ResNet-50和LSTM为基础的模型, 实验3使用以ResNet-50为基础的多视图架构, 实验4在实验3的基础上加入LSTM以用作序列处理。在多视图架构和LSTM的模型上, 实验5加入注意力机制, 实验6同时加入了注意力机制和焦点损失函数。从表 6可以看出, 对于不包含多识图架构的模型(LSTM和ResNet-50), 从准确率角度而言, 其结果均超过0.9, 但召回率方面存在明显不足,加入多视图架构后, 测试结果显示召回率至少提升了0.28, 将注意力机制加入模型后, 准确率和召回率均达到了0.97, 最终PSS-Net模型在测试集上的准确率为0.997, 召回率为0.999。
|
下载CSV 表 6 消融研究测试结果 Table 6 Test results of ablation study |
本文提出一种基于多视图的神经网络模型PSS-Net用于毫米波图像威胁识别。多视图架构解决了扫描对象的多方位问题, 注意力机制显著提升召回率, 在实际应用中可有效降低误判率。该模型在HD-AIT毫米波人体威胁扫描数据集中取得较高的识别精度, 验证了其在威胁识别领域的有效性, 有助于危险物品的精细化识别, 对于提升公共区域安全具有重要意义。后续将针对不同对象的生理畸变进行研究, 以提高本文模型在几何变形严重场景下的识别准确率。
| [1] |
WANG Huijun.Research on millimeter wave passive imaging and image processing[D].Nanjing: Nanjing University of Science and Technology, 2007.(in Chinese) 王慧君.毫米波被动成像及图像处理的研究[D].南京: 南京理工大学, 2007. |
| [2] |
ALEXANDER N E, CARLOS C A, GONZALO R. Multispectral mm-wave imaging:materials and images[J]. Proceedings of SPIE-the International Society for Optical Engineering, 2008, 2: 39-50. |
| [3] |
OKA S, TOGO H, KUKUTSU N, et al. Latest trends in millimeter-wave imaging technology[J]. Progress in Electromagnetics Research Letters, 2008, 1: 197-204. DOI:10.2528/PIERL07120604 |
| [4] |
YUJIRI L, SHOUCRI M, MOFFA P. Passive millimeter-wave imaging[J]. IEEE Microwave Magazine, 2003, 4(3): 39-50. DOI:10.1109/MMW.2003.1237476 |
| [5] |
FETTERMAN M R, GRATA J, JUBIC G, et al. Simulation, acquisition and analysis of passive millimeter-wave images in remote sensing applications[J]. Optics Express, 2008, 16(25): 203-215. |
| [6] |
CHEN Yangkan.Research on passive millimeter wave image processing and feature analysis[D].Nanjing: Nanjing University of Science and Technology, 2010.(in Chinese) 陈阳侃.无源毫米波图像处理与特征分析研究[D].南京: 南京理工大学, 2010. |
| [7] |
FARZAD M K, SEYED M S T, ZAHRA K V, et al. Improved two-dimensional millimeter-wave imaging for concealed weapon detection through partial Fourier sampling[J]. Journal of Infrared, Millimeter, and Terahertz Waves, 2016, 37(3): 267-280. DOI:10.1007/s10762-015-0223-z |
| [8] |
AMIZIC B, SPINOULAS L, MOLINA R, et al.Compressive sampling with unknown blurring function: application to passive millimeter-wave imaging[C]//Proceedings of 19th IEEE International Conference on Image Processing.Washington D.C., USA: IEEE Press, 2012: 1-10.
|
| [9] |
SUN Qiansong.Research on super-resolution reconstruction and enhancement of passive millimeter wave image[D].Chengdu: University of Electronic Science and Technology of China, 2011.(in Chinese) 孙乾松.无源毫米波图像超分辨重建与增强技术研究[D].成都: 电子科技大学, 2011. |
| [10] |
ZHANG Renlin, LI Yan. Application of threshold shrinking iterative algorithm in passive millimeter wave image processing[J]. Journal of Jiujiang University (Natural Science Edition), 2013, 28(1): 24-26. (in Chinese) 张仁霖, 李炎. 阈值收缩迭代算法在被动毫米波图像处理中的应用[J]. 九江学院学报(自然科学版), 2013, 28(1): 24-26. |
| [11] |
HAWORTH C D, PETILLOT Y R, TRUCCO E. Image processing techniques for metallic object detection with millimetre-wave images[J]. Pattern Recognition Letters, 2006, 27(15): 1843-1851. DOI:10.1016/j.patrec.2006.02.003 |
| [12] |
MAQUEDA I G, DE LA BLANCA N P, MOLINA R, et al.Fast millimeter wave threat detection algorithm[C]//Proceedings of the 23rd European Signal Processing Conference.Washington D.C., USA: IEEE Press, 2015: 228-236.
|
| [13] |
MARTINEZ O, FERRAZ L, BINEFA X, et al.Concealed object detection and segmentation over millimetric waves images[C]//Proceedings of 2010 IEEE Computer Society Conference on Computer Vision and Pattern Recognition-Workshops.Washington D.C., USA: IEEE Press, 2010: 1-8.
|
| [14] |
LECUN Y, BENGIO Y, HINTON G. Deep learning[J]. Nature, 2015, 521(7553): 436-444. DOI:10.1038/nature14539 |
| [15] |
PENG Yong, LU Baoliang. Discriminative extreme learning machine with supervised sparsity preserving for image classification[J]. Neurocomputing, 2017, 261: 242-252. DOI:10.1016/j.neucom.2016.05.113 |
| [16] |
GU J X, WANG Z H, KUEN J, et al. Recent advances in convolutional neural networks[J]. Pattern Recognition, 2018, 77: 354-377. DOI:10.1016/j.patcog.2017.10.013 |
| [17] |
LEE C Y, GALLAGHER P W, TU Z W.Generalizing pooling functions in convolutional neural networks: mixed, gated, and tree[EB/OL].[2019-10-02].https://arxiv.org/abs/1509.08985.
|
| [18] |
CLEVERT D, UNTERTHINER T, HOCHREITER S.Fast and accurate deep network learning by Exponential Linear Units(ELUs)[EB/OL].[2019-10-02].https://arxiv.org/abs/1511.07289.
|
| [19] |
SHELHAMER E, LONG J, DARRELL T. Fully convolutional networks for semantic segmentation[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2017, 39(4): 640-651. DOI:10.1109/TPAMI.2016.2572683 |
| [20] |
LOPEZ-TAPIA S, MOLINA R, DE LA BLANCA N P. Deep CNNs for object detection using passive millimeter sensors[J]. IEEE Transactions on Circuits and Systems for Video Technology, 2019, 29(9): 2580-2589. DOI:10.1109/TCSVT.2017.2774927 |
| [21] |
GUO Lin, QIN Shiyin. High-performance detection of concealed forbidden objects on human body with deep neural networks based on passive millimeter wave and visible imagery[J]. Journal of Infrared, Millimeter, and Terahertz Waves, 2019, 40(3): 314-347. DOI:10.1007/s10762-018-0558-3 |
| [22] |
LI Yuanju. Survey of classification with imbalanced data[J]. Modern Computer, 2016(3): 30-33, 50. (in Chinese) 李元菊. 数据不平衡分类研究综述[J]. 现代计算机, 2016(3): 30-33, 50. |
| [23] |
DU Feng, SHI Wenkang, DENG Yong, et al. Fast infrared image segmentation method[J]. Journal of Infrared and Millimeter Waves, 2005, 24(5): 370-373. (in Chinese) 杜峰, 施文康, 邓勇, 等. 一种快速红外图像分割方法[J]. 红外与毫米波学报, 2005, 24(5): 370-373. |
| [24] |
SAK H, SENIOR A, BEAUFAYS F.Long short-term memory based recurrent neural network architectures for large vocabulary speech recognition[EB/OL].[2019-10-02].https://arxiv.org/abs/1402.1128.
|
| [25] |
YU Haiping, HE Fazhi, PAN Yiteng, et al.An efficient similarity-based level set model for medical image segmentation[EB/OL].[2019-10-02].https://www.jstage.jst.go.jp/article/jamdsm/10/8/10_2016jamdsm0100/_pdf.
|
| [26] |
HE Kaiming, ZHANG Xiangyu, REN Shaoqing, et al.Deep residual learning for image recognition[C]//Proceedings of 2016 IEEE Conference on Computer Vision and Pattern Recognition.Washington D.C., USA: IEEE Press, 2016: 1-10.
|
| [27] |
GUIMARAES A, TOFIGHI G.Detecting zones and threat on 3D body for security in airports using deep machine learning[EB/OL].[2019-10-02].https://arxiv.org/abs/1802.00565.
|
| [28] |
YAN Zhichong, FENG Shuai, LI Fangqi, et al. A multi-view deep learning approach for detecting threats on 3D human body[M]. Berlin, Germany: Springer, 2019: 286-294.
|
2020, Vol. 46

,