ABSTRACT

In this paper, based on the traditional pedestrian detection framework of RGB channel template and convolution neural network, a pedestrian detection method combining multi-channel template and convolution neural network is proposed. The main contribution is: firstly, the motion detection method is used to detect the interest area, improve pedestrian detection efficiency and reduce the false alarm. Secondly, multi-channel template and maximum fusion strategy are adopted to reduce the possible pedestrian leakage detection based on RGB channel template.

Keywords:Channel Template, Convolutional Neural Network, Pedestrian Detection

结合多通道模板和卷积神经网络的行人检测方法研究

盘先跃，刘怡杉，谭科

国防科技大学电子科学与工程学院，湖南 长沙

收稿日期：2018年1月5日；录用日期：2018年1月23日；发布日期：2018年1月30日

摘 要

本文在传统的依靠RGB通道模板和卷积神经网络的行人检测框架下，提出了结合多通道模板和卷积神经网络的行人检测方法。主要工作是：一是采用运动侦测方法初步检测感兴趣区域，提高行人检测效率和降低虚警。二是采用多通道模板和最大值融合策略，降低单纯依靠RGB通道模板可能存在的行人漏检现象。

关键词 :通道模板，卷积神经网络，行人检测

Copyright © 2018 by authors and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

1. 引言

对于营区周界警戒应用而言，行人的入侵检测需要考虑复杂天候的影响，譬如，在夜晚照明不足或者采用红外补光照明的情况下，图像的色彩信息丢失，纹理特征显著性降低，此时基于视频分析方式可靠检测行人目标是非常困难的。现有的行人检测方法所使用的训练和测试数据集都是白天的彩色图像，导致现有行人检测方法在夜晚低照度场景、雾雪天气等低对比度场景以及其他类似的复杂天候场景的行人检测正确率偏低，尤其是漏检率偏高 [1] 。这可能导致采用现有方法研制的营区周界警戒系统的安全性大幅下降。为此，本文针对复杂天候条件，提出一种结合多通道模板和卷积神经网络的行人检测方法，主要工作是融合多通道模板的卷积神经网络计算结果，降低复杂天候条件下依靠单一通道模板可能存在的行人漏检现象。该方法的基本流程如图1所示，主要包括感兴趣区域提取、多通道模板提取、卷积神经网络计算、融合判决四个部分。详细描述如下。

2. 感兴趣区域提取

目前，基于视频分析的行人检测方法中检测正确率较高的方法基本都是采用深度学习算法实现的，这类算法的优点是行人检测正确率高，虚警率低。然而，深度学习算法往往比较耗时，运算效率偏低，这对于营区周界警戒系统来说往往是难以接受的。为了提高行人检测效率，本文首先采用运动侦测技术快速提取视频中的感兴趣区域，仅针对感兴趣区域采用深度学习算法检测行人目标。这样可以在不降低行人检测正确率的基础上大幅提高行人检测效率，对营区周界警戒系统而言很有意义。同时，由于运动侦测可以剔除许多没有运动目标的背景区域，从而也可以降低这些区域引起虚警的可能，进一步降低行人检测的虚警率。

目前，运动侦测方法主要分为三类：时间差分法、背景差分法和光流法。时间差分法是利用图像序列中连续帧(通常为两到三帧)之间的各像素差异来检测运动目标的方法，比较适用于动态环境。但对光照变化非常敏感，虚检现象较为严重。同时，对于低速运动目标容易出现孔洞效应。背景差分法是运动侦测最常用的方法之一，其基本原理是从参考背景图像中减去当前帧图像。而且背景图像在一段时间内实时更新，可以适应动态背景的变化。该方法通常要求摄像机静止，对于动摄像机的运动侦测效果不佳。光流法通过检测由于观察者和参考场景之间的相关运动而产生的光流来检测运动目标，对动摄像机的适应性较好，但运算效率偏低。针对营区周界警戒系统对运算效率和目标完整性要求较高的需求，这里选择背景差分法来进行感兴趣区域的提取。具体是采用高斯混合模型(Gaussian Mixture Model, GMM)来检测运动的感兴趣目标。

图1. 基于多通道模板和CNN的行人检测流程

假设图像中的像素可以用K个高斯分布的和来描述， $X_t$ 代表在时间t的像素值，其概率描述为

$P\left(X_t\right)={\displaystyle \sum _{k=0}^K{\omega }_{k,t}{\eta }_k\left(X_t,{\mu }_{k,t},{\Sigma }_{k,t}\right)}$ (1)

其中，K表示高斯分布的数目， ${\omega }_{k,t}$ 表示时间t的第k个高斯分布的权重，它反映了该分布可以作为背景模型的程度， ${\mu }_{k,t}$ 是表示时间第k个高斯分布的均值， ${\Sigma }_{k,t}$ 表示时间t第k个高斯分布的协方差矩阵，假设 $X_t$ 的分量是独立的，且方差相同，即 ${\Sigma }_{k,t}={\sigma }_{k,t}I$ ， ${\eta }_k$ 表示第k个高斯分布，描述如下

${\eta }_k\left(X_t,{\mu }_{k,t},{\Sigma }_{k,t}\right)=\frac{1}{\sqrt{{\left(2\text{π}\right)}^n\left|{\Sigma }_{k,t}\right|}}\exp \left(-\frac{1}{2}{\left(X_t-{\mu }_{k,t}\right)}^{\text{T}}{\Sigma }_{k,t}{}^{-1}\left(X_t-{\mu }_{k,t}\right)\right)$ (2)

在模型更新之后，K个高斯分布按照 $\frac{\omega }{\sigma }$ 由大到小的顺序重新排序，选取前B个最有可能代表背景的

高斯分布，公式为

$B=\underset{b}{\arg \min }\left({\displaystyle \sum _{k=1}^b{\omega }_{k,t}}>T_W\right),b=1,2,\cdots ,K$ (3)

其中， $T_W$ 是描述背景在场景中的最小先验概率的阈值。

由于光照等复杂场景的变化，背景模型也应是一个动态模型，需要根据场景的变化来自适应更新。模型更新是背景差分算法的重要环节，对于高斯背景模型而言，需要更新的是高斯分布的参数，公式为

${\omega }_{k,t}=\left(1-\alpha \right){\omega }_{k,t}+\alpha M_{k,t}$ (4)

其中， $\alpha$ 称为学习率，取值范围在0~1之间。 $M_{k,t}$ 是一个二元值，当观察像素与高斯分布匹配时， $M_{k,t}=1$ ，否则 $M_{k,t}=0$ 。观察像素与高斯分布匹配的条件为

$\left|X_t-{\mu }_{k,t-1}\right|\le 2.5{\sigma }_{k,t-1}$ (5)

如果没有发现与观察像素相匹配的高斯分布，则将观察像素值作为权重最小的高斯分布的平均值。

如果存在与观察像素相匹配的高斯分布，则相匹配的高斯分布的平均值和方差将被更新，公式为

$\{\begin{array}{l}{\mu }_{k,t}=\left(1-\rho \right){\mu }_{k,t-1}+\rho X_t\\ {\sigma }_{k,t}^2=\left(1-\rho \right){\sigma }_{k,t-1}^2+\rho {\left(X_t-{\mu }_{k,t-1}\right)}^2\end{array}$ (6)

其中， $\rho$ 为学习因子，计算公式为

$\rho =\alpha {\eta }_k\left(X_t|{\mu }_{k,t-1},{\Sigma }_{k,t-1}\right)$ (7)

其余不匹配的高斯分布保持不变。

在本文中，依据部分运动侦测实验结果以及高斯背景模型的经验参数，设置学习率 $\alpha =0.005$ 、阈值 $T_W=0.7$ 、高斯分布数量 $K=6$ 。

3. 多通道模板提取

对于上一小节获取到的每一个感兴趣区域，本文采用目前检测正确率高的深度学习框架来进行行人检测。相对于原始视频图像，感兴趣区域的尺寸很小，行人检测耗时也小。现有基于深度学习的行人检测方法主要是针对RGB格式的图像进行的。然而，营区周界警戒环境的天候复杂度高，单纯依靠RGB格式的图像信息容易造成漏检。为此，本文提出多通道模板的联合卷积神经网络检测框架，通过多通道模板增强行人特征的显著性，通过联合卷积神经网络计算各个通道模板的行人概率，最后融合各个通道模板的CNN计算结果形成最终判决，可以大幅降低复杂天候条件下的行人漏检问题，这对提高营区周界警戒系统的安全性具有极大的促进作用。

本文所述的多通道模板包括：RGB通道模板(记为RGB_CM)、亮度通道模板(记为Y_CM)、局部二元模式通道模板(记为LBP_CM)和梯度通道模板(记为G_CM)。其中，RGB_CM也即RGB颜色空间的图像模板，Y_CM也即灰度空间的图像模板，LBP_CM也即采用局部二元模式对灰度图像进行变换之后得到的图像模板，G_CM也即对灰度图像求取边缘得到的图像模板。选择这四种通道模板的依据是：首先，RGB_CM是目前行人检测常用的通道模板，该通道模板对于白天场景的行人检测效果很好，已经在行人检测的各个数据集上得到了验证。其次，选择Y_CM主要原因是摄像机在夜晚模式红外补光条件下图像彩色信息会丢失，此时采用亮度通道更能体现行人特征的显著性。另外，选择LBP_CM和G_CM的主要原因是考虑到营区周界警戒系统需要全天候运行，环境光照变化比较大，采用LBP通道和梯度通道可以消除光照变化的影响。另外，在可视条件差的雾雪天气，LBP通道和梯度通道更能凸显行人的显著特征。这里LBP通道主要从纹理的角度来凸显行人特征，梯度通道主要从边缘的角度凸显行人特征。通过这四个通道模板，可以增强单纯使用RGB通道模板的行人特征的显著性。图2展示了四种通道模板的示意图。

下面简要介绍如何依据RGB_CM构建Y_CM、LBP_CM和G_CM。

令 $f_{RGB}^{\left(R\right)}\left(x,y\right)$ 、 $f_{RGB}^{\left(G\right)}\left(x,y\right)$ 和 $f_{RGB}^{\left(B\right)}\left(x,y\right)$ 分别表示像素点 $\left(x,y\right)$ 在RGB_CM上的R通道值、G通道值和B通道值， $f_Y\left(x,y\right)$ 表示像素点 $\left(x,y\right)$ 在Y_CM上的像素值， $f_{LBP}\left(x,y\right)$ 表示像素点 $\left(x,y\right)$ 在LBP_CM上的像素值， $f_G\left(x,y\right)$ 表示像素点 $\left(x,y\right)$ 在G_CM上的像素值，那么有

$f_Y\left(x,y\right)=0.3f_{RGB}^{\left(R\right)}\left(x,y\right)+0.59f_{RGB}^{\left(G\right)}\left(x,y\right)+0.11f_{RGB}^{\left(B\right)}\left(x,y\right)$ (8)

图2. 四种通道模板示意图

$f_G\left(x,y\right)=\sqrt{{\left(f_Y\left(x+1,y\right)-f_Y\left(x-1,y\right)\right)}^2+{\left(f_Y\left(x,y+1\right)-f_Y\left(x,y-1\right)\right)}^2}$ (9)

$f_{LBP}\left(x,y\right)={\displaystyle \sum _{n=1}^8{f}_{LBP}^{(n)}\left(x,y\right)2^{n-1}}$ (10)

其中，

$f_{LBP}^{\left(n\right)}\left(x,y\right)=\{\begin{array}{l}1,f_Y^{\left(n\right)}\left(x,y\right)\ge f_Y\left(x,y\right)\\ 0,其他\end{array}$ (11)

这里， $f_Y^{\left(n\right)}\left(x,y\right)$ 表示像素点 $\left(x,y\right)$ 所在8邻域的第n个像素点的亮度值，本文将像素点 $\left(x,y\right)$ 右侧相邻像素点作为第1个邻域像素点，按照顺时针的方式依次为第2~第8个邻域像素点，如图3所示。

4. 卷积神经网络计算

卷积神经网络(Convolutional Neural Networks, CNN)是一个由多层网络组成的网络结构，每一层的输出由上一层的输出和当前层的参数决定 [2] ，可以表示为

$x^l=f^l\left(x^{l-1},{\theta }^l\right)$ (12)

其中， $x^{l-1}$ 和 $x^l$ 分别表示第l-1层和第l层的输出， ${\theta }^l$ 表示第l层的参数，可以表示为

${\theta }^l=\left[w^l,b^l\right]$ (13)

这里， $w^l$ 和 $b^l$ 分别表示第l层的权重和偏差参数。

$f^l$ 表示第l层的网络模型，也即该层网络所要执行的操作，主要分为四种类型：卷积、池化、全连接和回归。

CNN计算将输入的每一个图像映射到每一个类别的分类得分可以表示为

$p_c=\text{softmax}\left(x^{L-1}|{\theta }^L\right)$ (14)

其中，L为CNN的网络层数， $p_c\in \left[0,1\right]$ ， $c=\left\{0,1\right\}$ 。 $c=0$ 表示类别为非行人， $c=1$ 表示类别为行人。函数softmax的定义为

$\text{softmax}\left(x^{L-1}\right)=\frac{\exp \left(x_c^{L-1}\right)}{\exp \left(x_0^{L-1}\right)+\exp \left(x_1^{L-1}\right)}$ (15)

卷积的定义是：

$x_{m^{\prime },n^{\prime }}^l={\displaystyle \sum _m{\displaystyle \sum _n\left(w_{m,n,d^{\prime }}^l{x}_{m+m^{\prime },n+n^{\prime }}^{l-1}+b_{d^{\prime }}^l\right)}}$ (16)

图3. 邻域像素点序号示意图

其中m、n分别表示图像的高度和宽度， $d^{\prime }$ 表示滤波器的数量。

卷积之后紧跟激活函数，这里采用ReLU (Rectified Linear Unit)激活函数，定义为

$a\left(x_{m^{\prime },n^{\prime }}^l\right)=\max \left(0,x_{m^{\prime },n^{\prime }}^l\right)$ (17)

池化层采用最大值池化，定义为

$x_{m,n}^l=\underset{m^{\prime },n^{\prime }\in L\left(m,n\right)}{\max }\left(x_{m^{\prime },n^{\prime }}^{l-1}\right)$ (18)

这里， $L\left(m,n\right)$ 表示输入位置。

关于CNN架构的选择和初始化，本文采用常用的VGG网络 [3] 进行初始化(具体是将caffe平台下已训练好的prototxt和caffemodel文件加入数据层)，然后使用自己构建的数据进行微调。为了提高行人检测效率，将VGG网络的原始输入图像尺寸由224 × 234改为64 × 64。VGG网络共包含16层，即L = 16。卷积层的尺寸为3 × 3，池化层的尺寸为2 × 2。具体采用开源的Caffe架构下的VGG网络实现。

5. 融合判决

本文使用了四个多通道模板(RGB_CM、Y_CM、LBP_CM和G_CM)，每一个通道模板经过CNN计算之后都可以由公式(3)~(14)得到一个分类得分，分别记为 $p_c^{\left(RGB\_CM\right)}$ 、 $p_c^{\left(\text{Y}\_CM\right)}$ 、 $p_c^{\left(\text{LBP}\_CM\right)}$ 和 $p_c^{\left(\text{G}\_CM\right)}$ 。下面考虑如何融合这些分类得分形成最终的判决。

常用的融合方法有四种，包括平均、加权平均、多数投票和加权投票。但是，本文采用多通道模板的目标是降低恶劣天候下行人检测的漏检率。因此，本文采用极大值融合策略，选出各个通道模板对行人类别的分类得分的最大值，表示为

$p_1=\max \left(p_1^{\left(RGB\_CM\right)},p_1^{\left(\text{Y}\_CM\right)},p_1^{\left(\text{LBP}\_CM\right)},p_1^{\left(\text{G}\_CM\right)}\right)$ (19)

设定一个得分阈值 $T_s$ ，如果 $p_1>T_s$ ，则给出检测到行人目标的判决，否则给出未检测到行人目标的判决。在本文中，阈值 $T_s=0.5$ 。也即，四个多通道模板中只要有一个模板的行人类别分类得分大于设定的得分阈值，就判断检测到行人目标。这样可以大幅降低传统的RGB通道模板可能存在的漏检现象。

6. 实验与分析

行人检测是计算机视觉的研究热点之一，目前用于行人检测的公开测试数据集也较多，本节首先介绍行人检测领域的数据集以及针对营区周界警戒应用的自建数据集，然后介绍行人检测的常用评价指标，最后分析和对比本文方法与现有性能较好的行人检测方法的性能指标。

6.1. 行人检测数据集

目前，行人检测领域的公开测试数据集主要有6个，包括INRIA、ETH、TUD-Brussels、Daimler、KITTI和Caltech。这些数据集依据构建目标的不同都有不同的特点。其中，常用的是INRIA和Caltech两个数据集。简要介绍如下：

1) INRIA

INRIA是最早使用的行人检测数据集，该数据集包含的图像共有5264幅(其中包含行人的图像有3548幅)。尽管该数据集规模不大，但优势在于包含的行人图像是经过裁剪的，非常适用于行人检测器的训练。

2) Caltech

Caltech数据集是由穿越美国洛杉矶街头的车辆记录的2.5小时30 Hz视频构成的，视频分辨率为640 × 480。视频中注释了共计35万个包围约2300个独立行人的包围盒，也即人工标注的行人窗口数量约为350,000个。由于这些数据都是来自真实街区场景的自然行人，行人存在遮挡现象和姿态变化，因此检测难度较大。

对于营区周界警戒系统而言，还需要考虑复杂天候的影响，譬如夜晚可视程度差的情况。因此，在现有数据集的基础上，本文还自建了营区周界警戒数据集，采用视频录制和互联网收集的方式构建了共有50个视频片段，并标记视频片段中每帧图像中行人出现的位置，用于训练和测试。

6.2. 常用性能评价指标

行人检测的性能评价指标主要有两个：真正率(True Positive Rate, TPR)和假正率(False Positive Rate, FPR)，计算公式如下

$\text{TPR}=\frac{正确检测到的行人数量}{标记的行人数量}\times 100\%$ (20)

$\text{FPR}=\frac{错误检测到的行人数量}{错误检测到的行人数量+正确检测到的行人数量}\times 100\%$ (21)

这里，正确检测到行人的含义是：检测到的行人窗口与人工标记的行人窗口的重合度大于50%。错误检测到行人的含义是：检测到的行人窗口与人工标记的行人窗口的重合度不大于50%。重合度定义为

${\text{IoU}}^{\left(r,\tilde{r}\right)}=\frac{{\stackrel{⌢}{S}}^{\left(r,\tilde{r}\right)}}{S^{\left(r\right)}+S^{\left(\tilde{r}\right)}-{\stackrel{⌢}{S}}^{\left(r,\tilde{r}\right)}}\times 100\%$ (22)

其中， $r$ 和 $\tilde{r}$ 分别表示标记的行人窗口和检测到的行人窗口， $S^{\left(r\right)}$ 和 $S^{\left(\tilde{r}\right)}$ 分别表示 $r$ 和 $\tilde{r}$ 的面积， ${\stackrel{⌢}{S}}^{\left(r,\tilde{r}\right)}$ 是指 $r$ 与 $\tilde{r}$ 之间交叉区域的面积。

另外，行人检测效率也是行人检测性能的重要评价指标之一，本文采用平均检测耗时来描述行人检测效率，该值越小说明行人检测效率越高。平均检测耗时是指平均对每一帧视频图像进行行人检测所耗费的时间，不包含视频解码、显示以及分类器训练耗时。由于耗时测算与处理器和处理环境性能关联大，因此实验时需要在相同平台上进行耗时测试，本文使用的计算性平台参数如表1所示，后续不再赘述。

6.3. 不同方法的性能对比分析

为了评价本文提出的结合多通道模板和卷积神经网络的行人检测方法(简记为MC_CNN方法)的性能，将该方法与目前行人检测领域性能较好的方法进行对比实验。这里选取了三种方法进行对比实验。其中，文献 [4] 是一种基于特征的行人检测方法，采用HOG特征和SVM分类器实现行人的检测，简记为HOG_SVM方法。文献 [5] 是一种基于深度学习的行人检测方法，采用深度卷积神经网络(Deep Convolutional

表1. 计算机平台参数

Neural Networks)进行行人目标检测，简记为DCNN方法。文献 [6] 是一种结合运动检测和特征分类的行人检测方法，采用光流法检测运动区域，然后结合HOG特征和SVM分类器进行行人的检测，简记为OF_HOG_SVM方法。这三种方法可以从不同的侧面来印证本文方法的优越性。需要说明的是，四种方法都使用INRIA数据集、Caltech数据集的set00~set05子集以及自建数据集的前20个视频片段进行训练。测试时，首先在Caltech数据集的set06~set10子集上进行测试，然后再在剩余的自建数据集上进行测试，分别统计各个数据集上的行人检测性能指标。其中，图4展示了Caltech数据集上四种方法的性能指标。图5展示了自建数据集上四种方法的性能指标。表2给出了四种方法的平均检测耗时对比结果。下面结合实验结果进行分析。

1) FPR指标对比分析

由图4和图5可见，本文方法和文献 [5] 所述方法的FPR值明显低于其他两种方法，究其原因，主要是因为本文方法和文献 [5] 所述方法都是采用深度学习框架实现行人目标的检测，而文献 [4] 和文献 [6] 所述方法都是采用结合特征与分类器的方法实现行人目标的检测。这说明，深度学习框架在行人目标检测应用中的虚检率低，也即把非行人目标检测为行人目标的概率低。另外，本文方法的FPR值低于文献 [5] 所述方法，文献 [6] 所述方法的FPR值也低于文献 [4] 所述方法，这主要是由于本文方法和文献 [6] 所述方法采用运动侦测滤除了背景区域的干扰，降低了在这些区域将非行人目标检测为行人目标的概率，从而降低了FPR的值。对比图4和图5中各方法的FPR值可以发现，文献 [4] 所述方法在自建数据集下的FPR值增大明显，这是因为自建数据集下的干扰目标较多，容易造成虚检。而文献 [6] 所述方法通过运动侦测过滤了大量干扰目标，因此FPR值增大不显著。文献 [5] 所述方法采用深度学习框架进行行人目标的检测，虚检现象不显著，故FPR增加也不明显。本文方法通过感兴趣区域提取环节过滤了大量干扰目标，又结合卷积神经网络的深度学习框架实现了可靠的行人检测，因此FPR值也没有明显增加。

2) TPR指标对比分析

由图4和图5可见，本文方法的TPR值高于文献 [5] 所述方法，这是因为本文采用多通道模板增强了行人特征的显著性，在特征融合时采用最大值融合策略降低了行人目标的漏检率。本文方法的TPR值也高于文献 [4] 和文献 [6] 所述方法，主要原因是后两种方法在提取行人目标特征时受行人姿态和遮挡情况的影响较大，如果训练数据集中的行人姿态和遮挡情况不完备，那么在测试时难免存在行人漏检的现象。对比图4和图5中各方法的TPR值可以发现，文献 [4] 、文献 [5] 和文献 [6] 所述方法在自建数据集下的TPR值都有明显的下降，这是因为自建数据集中行人目标的光照、对比度等条件不理想，导致这三种方法存在不同程度的漏检。而本文方法融合了多通道模板，在RGB_CM、Y_CM、LBP_CM和G_CM四个通道模板上分别检测行人目标，避免某通道行人特征不显著造成的漏检现象，因此TPR值没有明显下降。

图4. Caltech数据集测试结果

图5. 自建数据集测试结果

3) 运算效率对比分析

由表2可见，文献 [5] 所述方法的平均检测耗时最多，其次是文献 [4] 所述方法，本文方法和文献 [6] 所述方法的平均检测耗时相当，明显小于其他两种方法，说明本文方法和文献 [6] 所述方法的运算效率高于其他两种方法。分析其原因，主要是因为本文方法和文献 [6] 所述方法都使用了运动检测方法过滤了大量背景区域，从而大幅降低了检测耗时。

表2. 两个数据集下四种方法的平均检测耗时对比

7. 结束语

本文在传统的依靠RGB通道模板和卷积神经网络的行人检测框架下，提出了结合多通道模板和卷积神经网络的行人检测方法。主要贡献在于：一是采用运动侦测方法初步检测感兴趣区域，提高行人检测效率和降低虚警。二是采用多通道模板和最大值融合策略，降低单纯依靠RGB通道模板可能存在的行人漏检现象。实验结果表明，本文方法的真正率高、假正率低，检测耗时少。后续研究重点是进一步提高检测效率，以期达到实时处理的效果。

文章引用

盘先跃,刘怡杉,谭 科. 结合多通道模板和卷积神经网络的行人检测方法研究
A Study of Pedestrian Detection Based on Multi-Channel Template and Convolution Neural Network[J]. 计算机科学与应用, 2018, 08(01): 97-106. http://dx.doi.org/10.12677/CSA.2018.81013

参考文献 (References)