Infrared Image Segmentation Based on Improved Spotted Hyena Optimizer
-
摘要: 针对斑点鬣狗优化算法（spotted hyena optimizer，SHO）容易陷入局部最优解、求解质量低等缺点，本文提出使用Lévy飞行和单纯形搜索算法改进SHO（spotted hyena optimizer based on simplex method and Lévy flight, Lévy_SM_SHO）。将Lévy_SM_SHO与Lévy飞行斑点鬣狗优化算法（spotted hyena optimizer based on Lévy flight, Lévy_SHO）、单纯形搜索斑点鬣狗优化算法（spotted hyena optimizer based on simplex method, SM_SHO）和SHO在测试函数上结果进行对比，实验证明改进算法能够取得较好的优化结果。并将Lévy_SM_SHO算法用于红外图像阈值分割问题，通过与粒子群算法（particle swarm optimization, PSO）分割结果对比，证明Lévy_SM_SHO算法能够取得较好的阈值分割结果。
-
关键词:
- Lévy飞行 /
- 单纯形搜索算法 /
- Lévy_SM_SHO /
- Lévy_SHO /
- SM_SHO
Abstract: Based on the shortcomings of the spotted hyena optimizer (SHO), falling into a local optimal solution or a low-quality solution is easy. In this study, the Lévy flight and simplex method are proposed to improve the SHO(Lévy_SM_SHO). Comparing Lévy_SM_SHO to Lévy flight spotted hyena optimizer (Lévy_SHO), simplex method spotted hyena optimizer (SM_SHO), and spotted hyena optimizer (SHO) on the test function, the experiment proves that the improved algorithm can achieve better optimization results. Finally, the Lévy_SM_SHO algorithm is applied to the infrared image threshold segmentation problem. By crosschecking the segmentation results with the particle swarm optimization algorithm (PSO), we proved that the Lévy_SM_SHO algorithm can achieve better threshold segmentation results.-
Keywords:
- lévy flight /
- simplex search algorithm /
- Lévy_SM_SHO /
- Lévy_SHO /
- SM_SHO
-
0. 引言
人体摔倒检测可以有效地检测视频中的摔倒行为，降低人在摔倒后无法及时呼救的风险[1-3]。现有的基于计算机视觉的人体摔倒检测方法大多使用可见光图像，这类方法受光照影响在昏暗环境下效果不佳。红外图像能够保护个人隐私，而且不受光照和恶劣天气的影响，适用于全天候的人体摔倒检测，本文将对红外图像摔倒检测算法进行研究。现有的人体摔倒检测方法大致可分为以下3类：1）基于Freeman链码的检测方法，文献[4]利用Freeman链码记录轮廓信息，计算实际Freeman链码与摔倒模板链码的欧式空间距离结合时间判定条件判断是否发生摔倒行为，但在大规模遮挡时容易出现误差；2）基于关键点的摔倒检测方法，文献[5]利用建立的红外图像摔倒数据集进行CenterNet网络训练，识别红外图像中的人体姿态，但没有突出摔倒动作的时序性，容易对缓慢坐地等类摔倒动作造成误检；3）基于宽高比和质心变化率的检测方法，文献[6]利用人体宽高比这一特征来判断是否有人摔倒，并使用质心变化率和有效面积比对容易造成误判的情况进行修正，从几何特征的角度判断人体是否摔倒，但对特定角度的摔倒检测效果不佳。
针对上述算法的不足，本文在现有的摔倒检测方法的基础上提出了一种基于改进Alphapose的红外图像摔倒检测算法，该算法使用改进后的Alphapose检测红外图像中的骨架关键点和人体姿态，利用得到的人体骨架关键点信息和姿态类别进行摔倒判定，兼顾了摔倒动作的时序性和前后动作的关联性。
1. 改进的Alphapose算法
1.1 YOLO v5s
人体检测的精度直接影响后续关键点检测的效果，是整个算法精度的保证。YOLO系列算法[7-9]将候选框的搜索和目标的识别定位合二为一，相对于RCNN[10-12]系列算法先搜索再识别的两阶段处理模式，YOLO算法的处理速度更快，适用于端到端的训练与预测，目前已发展至YOLO v5版本。YOLO v5的目标检测网络有4种版本，主要区别在于网络深度不同，如表 1所示。
表 1 YOLO v5各版本参数对比Table 1. Comparison of parameters between different versions of YOLO v5Model Size
/pixelsSpeedV100/ms Params(M) YOLO v5s 640 2.0 7.3 YOLO v5m 640 2.7 21.4 YOLO v5l 640 3.8 47.0 YOLO v5x 640 6.1 87.7 从表 1可以看到，YOLO v5s的模型更小、速度更快，方便部署在嵌入式设备上，本文算法使用YOLO v5s的网络结构。
YOLO v5s网络首先使用Focus模块对输入图像进行切片操作，如图 1所示。如果网络输入的图片大小为640×640×3，则将图像隔行隔列取值，生成320×320×12的子图，再通过32个卷积核卷积生成320×320×32的特征图。相对于普通下采样，Focus模块保留了原始图片的完整信息。
YOLO v5s的Backbone部分使用了CSP[13]模块，将输入的特征映射划分为两部分，其中一部分经过多个残差组件和卷积，然后与另一部分拼接，在减少计算量的同时可以保证准确率。
Neck部分采用了PANet[14]结构，主要用于生成特征金字塔。特征金字塔会增强模型对于不同缩放尺度对象的检测，从而能够识别不同大小和尺度的同一个物体。
1.2 Alphapose算法
Alphapose[15]算法是一种人体姿态估计算法，能够检测出图像中的人体骨架，其提出一种区域多人姿态估计框架（regional multi-person pose estimation，RMPE），主要由对称空间变换网络（symmetric spatial transformer network，SSTN）、参数化姿态非极大抑制（parametric pose no-maximum-suppression，PP-NMS）、姿态引导建议发生器（pose-guided proposals generator，PGPG）和并行单人姿态估计器（parallel single person pose estimation，Parallel SPPE）组成，其中PGPG用于生成大量训练样本，Parallel SPPE则充当额外的正则项，以避免局部最小值，而SSTN又由空间变换网络（STN）、单人姿态估计器（SPPE）和空间反变换网络（spatial de-transformer network，SDTN）组成。Alphapose首先使用目标检测算法检测单张图片，得到单人人体目标框，作为RMPE的输入，进入SSTN模块。检测到的人体目标框可能目标偏离中心或者目标框中人体占比过小，影响后续姿态估计的效果，所以需要先使用STN提取出高质量的单人区域，然后输入SPPE得到估计的姿态，再由SDTN将估计的姿态反变换到原始人体区域框中。冗余的人体区域框输入姿态估计网络，可能会检测到冗余的姿态。Alphapose提出PP-NMS用来消除多余的姿态。其定义了一种姿态距离来度量姿态间的相似度，并在此基础上建立了冗余姿态的消除标准。PP-NMS首先选取最大置信度的姿态作为参考，并且根据消除标准将靠近该参考的区域框进行消除，多次重复这个过程直到消除所有的冗余识别框。
1.3 改进的Alphapose算法
原Alphapose算法主要使用的人体检测器是YOLO v3，与之相比，YOLO v5s的模型更小，速度更快。本文算法使用YOLO v5s作为Alphapose的人体检测器，改进后的Alphapose算法结构如图 2所示。
改进后的算法使用YOLO v5s网络完成对人体姿态的直接分类，同时从输入的红外图像中提取人体区域框，输入人体姿态估计网络检测人体骨架关键点，将关键点信息与预测的姿态类别一并输出。
2. 基于改进Alphapose的红外图像摔倒检测方法
本文提出基于改进Alphapose的红外图像摔倒检测算法。该算法首先使用YOLO v5s网络对图像中的人体姿态进行分类，再提取人体目标框检测出单人的骨架关键点，最后利用关键点和姿态类别信息，通过关键点分析和摔倒判定，判断是否出现摔倒。
2.1 关键点分析
Alphapose算法检测出的人体骨架关键点如图 3所示。
图 3(a)是站立时的人体骨架，图 3(b)是摔倒后的人体骨架。图中编号为0~17的点代表检测得到的18个人体骨架关键点，其中点11和点12所在位置代表髋关节。
通过实验发现，站立或行走的人体目标在发生摔倒时，最直观的表现是髋关节位置在垂直方向上的迅速下跌，所以通过监测髋关节在序列帧中的移动速度可以及时发现产生的摔倒动作。但随着运动中的人体目标与摄像头的距离越来越远，其髋关节在图像上的位移速度也越来越小，容易造成漏检。针对该问题，本文提出以前一帧中人体目标框的高度作为参考对象，计算髋关节关键点的相对移动速度。
在对红外视频的检测中，设第i帧与第i－1帧均为单人红外图像，2≤i≤N，N表示红外视频的总帧数，且都是非遮挡状态，其髋关节关键点均被检出。设图像左上角为原点，水平向右为X轴正方向，垂直向下为Y轴正方向，建立直角坐标系。记第i帧人体骨架中编号为n的关键点的纵坐标为yin，取点11和点12的中心点M代表髋关节，则点M在第i帧中的纵坐标yiM如式(1)所示。
$$ y_i^M = (y_i^{11} + y_i^{12})/2 $$ (1) 在摔倒过程中，当前一帧人体目标框的高度明显小于其宽度时，即使下降速度不变，也会因为前一帧目标框高度减小，导致计算得到的速度偏大，所以本算法引入前一帧人体目标框的宽高比作为限制条件，防止前一帧目标框宽高比过大导致点M的微小起伏被误识别为摔倒。记第i帧中人体目标框的高度为Hi，宽度为Wi，则第i帧中的人体宽高比Pi如公式(2)所示：
$$ {P_i} = {{{W_i}} \mathord{\left/ {\vphantom {{{W_i}} {{H_i}}}} \right. } {{H_i}}} $$ (2) 记第i帧中点M在垂直方向上的相对移动速度为vi，如公式(3)所示：
$$ {v_i} = \frac{{y_i^M - y_{i - 1}^M}}{{{H_{i - 1}}}},\quad 2 \leqslant i \leqslant N $$ (3) 设置大于0的阈值Th，检测可能出现的摔倒行为，如公式(4)所示。
$$ {F_i} = \left\{ {\begin{array}{*{20}{l}} {1,\;\;{v_i} > \text {Th and }{P_{i - 1}} < 1.0,\quad 2 \leqslant i \leqslant N} \\ {0,\;\;else} \end{array}} \right.\; $$ (4) 当Fi等于1时，判断可能发生了摔倒。
2.2 摔倒判定
仅仅分析关键点的下落速度有可能将快速下蹲、俯身拾取等类摔倒动作误识别为摔倒，为了减少误判，需要结合关键点分析结果和姿态分类做进一步判定。
经实验发现，摔倒后人体姿态会在短时间内稳定，直观表现是坐在地上、躺下或手部撑地，所以本文算法在经过关键点分析判断有可能发生摔倒后，继续检测后续多帧内的人体姿态。
在统计后续的人体姿态中，若第i帧图像中人体姿态类别为li，则对应的人体图像得分为si，两者关系可用公式(5)表示：
$$s_i=\left\{\begin{array}{l} 1, l_i \in[\text { "sit_flat","lie","push_up"] } \\ 0, \text { else } \end{array}\right.$$ (5) 式(5)中：当图像中的人体姿态类别li为“sit_flat”（坐地）、“lie”（躺下）、“push_up”（撑地）这3种姿态中的任意一种时，该人体图像的得分记为1，否则得分记为0。
当Fi＝1，即第i帧检测到可能发生摔倒时，继续检测后续20帧图像的人体姿态，若累计图像得分大于10，最终判定为摔倒，记为FFi＝1，如公式(6)所示：
$$ F{F_i} = \left\{ {\begin{array}{*{20}{l}} {1,\quad {F_i} = 1\;\;\text{and}\;\;\sum\limits_{k = 1}^{20} {{s_{i + k}} > 10} } \\ {0,\quad \text{else}} \end{array}} \right. $$ (6) 在实时检测过程中，以实时视频流作为输入，检测每一帧当前画面中的人体姿态和骨架关键点。从第二帧开始，如果前一帧与当前帧均检测出同一人的髋关节位置，即“无遮挡”状态，则进行关键点分析，计算髋关节关键点在这相邻两帧间的相对移动速率和方向。当满足可能摔倒的条件后，进行摔倒判定，开始统计后续20帧图像中的姿态类别，如果满足最终摔倒条件，则判断为一次摔倒。实时检测的情况下一次完整的摔倒检测流程图如图 4所示。
3. 数据集
3.1 采集和分类
本文使用的数据集分为3部分。第一部分是自建的室内红外数据集，摄像头采用艾睿光电红外热成像模组Xmodule T3，分辨率384×288，样本共6787张，场景包括走廊、办公室、大厅等室内场景；第二部分是运动物体热红外成像数据集（moving object thermal infrared imagery dataset，MOTIID），样本取自其中三段室外行人行走的视频，获取样本571张，均为室外场景；第三部分使用的是开放的免费FLIR热数据集FLIR_ADAS，获取样本926张。
本文对红外图像下的人体姿态定义了8种类别，分别是“直坐”、“躺下”、“下蹲”、“站立”、“弯腰”、“撑地”、“坐地”和“遮挡”，如图 5所示。
![]()
图 5 姿态分类：(a) 直坐；(b) 躺下；(c) 下蹲；(d) 站立；(e) 弯腰；(f) 撑地；(g) 坐地；(h) 遮挡Figure 5. Posture classification: (a) Sit straight; (b) Lie; (c) Squat;(d) Stand; (e) Bend; (f) Push up; (g) Sit flat; (h) Cover3.2 数据预处理
本文采用了Mosaic数据增强的方式对训练数据进行预处理，如图 6所示，Mosaic数据增强通过每次随机抽取数据集中的4张图片，采用随机裁剪、随机缩放、随机组合的方式生成一张新的训练图片。通过Mosaic数据增强，能够丰富场景内容，增强样本多样性，提高目标检测算法的抗干扰性。
4. 实验分析
本实验的实验环境为Win10_x64平台，CPU为Intel Core i7-9750H，内存16 G，GPU为NVIDIA GeForce RTX 2060，使用框架为torch1.6.0+opencv-python4.4.0，输入图片尺寸384×288，学习率为0.01，epochs为100。实验分为两部分：进行红外人体骨架检测实验，以测试本文算法定位精度和实时性；进行红外视频摔倒判定实验，以验证本文策略的可行性。
4.1 红外人体骨架检测实验
为了排除遮挡状态对关键点检测的干扰，需准确识别遮挡与非遮挡状态，如图 7所示。
图 7(a)与图 7(d)分别是遮挡状态的两种标注方式示意图，图 7(a)将未被遮挡的人体部分标注为遮挡，其预测结果如图 7(b)和图 7(c)所示，图 7(d)在标注遮挡类别的目标框时加入部分遮挡物，其预测结果如图 7(e)和图 7(f)所示。从图 7(b)中可以看到，全身被识别为站立姿态的同时，上身也被识别为遮挡，即出现了冗余检测，而如图 7(e)和图 7(f)所示，在遮挡状态标记框中加入部分遮挡物后，站立和遮挡得到了明显的区分。经实验，在标记遮挡状态时，当选取的遮挡物在标记框纵向占比为0.2或横向占比为0.3时，对遮挡和非遮挡状态的区分效果较好。
![]()
图 7 遮挡状态标注与预测图：(a) 标注时不含遮挡物；(b) 无遮挡物标注方式检测无遮挡状态；(c) 无遮挡物标注方式检测遮挡状态；(d) 标注时含部分遮挡物；(e) 含遮挡物标注方式检测无遮挡状态；(f) 含遮挡物标注方式检测遮挡状态Figure 7. Annotated and predicted images under occlusion: (a) No obstructions included when labeling; (b) Detection of unobstructed state using unobstructed object annotation method; (c) Detection of occlusion status using unobstructed object annotation method; (d) Include partial obstructions when labeling; (e) Detection of unobstructed state using labeling method with occluded objects; (f) Detection of occlusion status using labeling method with occluded objects本实验将本文算法与原Alphapose算法进行对比，取500张图片作为测试集，部分检测结果如图 8所示。
![]()
图 8 实验结果：(a) Alphapose检测自建红外数据集；(b) Ours检测自建红外数据集；(c) Ours检测MOTIID红外测试集；(d) Ours检测FLIR_ADAS数据集Figure 8. Experimental results: (a) Alphapose for detecting self built infrared datasets; (b) Ours for detecting self built infrared datasets; (c) Ours for detecting MOTIID dataset; (d) Ours for detecting FLIR_ADAS datasets本实验采用DIoU（Distance-IoU）衡量算法的定位精度，DIoU反映了真实与预测之间的相关度，计算方法如公式(7)所示：
$$ \text{DIoU} = \text{IoU} - \frac{{{\rho ^2}(b,{b^{\text{gt}}})}}{{{c^2}}} $$ (7) 式(7)中：b和bgt分别表示预测框和真实框的中心点位置坐标；ρ(b, bgt)是两个中心点之间的欧氏距离；c表示真实框和预测框的最小外接矩形的对角线距离；IoU表示真实框与预测框之间的交并比。DIoU值越接近于1，则预测位置越接近真实位置。
分别使用原Alphapose算法和本文算法计算DIoU值，其平均值DIoU_Average和平均处理时间如表 2所示。
表 2 实验结果Table 2. Experimental resultsMethod DIoU_average Time_average/s Alphapose 0.910 0.085 Ours 0.926 0.068 由图 8(a)和图 8(b)可以看到，因为在标记遮挡姿态时加入了部分遮挡物，所以识别出的遮挡姿态目标框也向被遮挡方向进行了部分延伸，符合遮挡的特征，这也验证了本文针对遮挡姿态修改标记方法的有效性。如图 8(a)~图 8(b)所示，原Alphapose算法与本文算法的检测结果基本相近，但也出现少数偏差。原Alphapose算法对“躺下”姿态的定位不够准确，导致后续的骨架检测出现误差，而本文算法得到的预测框更加逼近其真实目标框，得到的人体区域更加紧凑，检出的人体骨架也较准确。结合表 2中的计算结果可知，本文算法对人体区域的定位较准，速度更快。如图 8(c)~图 8(d)所示，本文算法也可检测室外人体目标，但对于远距离的人体可能产生漏检，如图 8(d)所示，远处路灯下的人体没有被检出。同时，在人群聚集和多人重叠场景下，容易发生漏检现象。
4.2 红外视频摔倒判定实验
本实验使用本文算法检测实时红外视频流，从中截取12段热红外视频，共计36个摔倒动作和51个类摔倒动作，其中类摔倒动作包括快速下蹲、俯身拾取和缓慢坐地。记摔倒动作被识别为摔倒的个数为TP，摔倒动作被识别为非摔倒的个数为FN，类摔倒动作被识别非摔倒的个数为TN，类摔倒动作被识别为摔倒的个数为FP。改变阈值Th，计算不同阈值下的准确率，计算方法如公式(8)所示：
$$ \text{Accuracy} = \frac{{\text{TP} + \text{TN}}}{{\text{TP} + \text{TN} + \text{FP} + \text{FN}}} $$ (8) 从表 3可以看出，当Th大于或等于0.04时，漏检的摔倒个数逐渐增多，因为设定的阈值大于部分真实摔倒动作的相对位移速度，容易将真实的摔倒动作漏检。而当Th小于或等于0.03时，真实的摔倒动作已全部被正确检出，但也出现误识别的情况，如将俯身拾取动作误识别为摔倒，如图 9所示。
表 3 不同Th下检测结果Table 3. Experimental results under different Th valuesTh TP FN TN FP Accuracy 0.01 36 0 41 10 88.5% 0.02 36 0 47 4 95.4% 0.03 36 0 50 1 98.9% 0.04 34 2 51 0 97.7% 0.05 32 4 51 0 95.4% 0.06 29 7 51 0 92.0% 0.07 25 11 51 0 87.4% ![]()
图 9 摔倒检测结果：(a) Th=0.03时，真实摔倒动作检测结果；(b) Th=0.03时，伪摔倒动作检测结果；(c)Th=0.01时，伪摔倒动作检测结果Figure 9. Fall detection results: (a) Detection result of real fall movements when Th=0.03; (b) Detection result of false fall action when Th=0.03; (c) Detection result of false fall action when Th=0.01图 9(a)与图 9(b)分别是Th取0.03时，本算法对真实摔倒动作和俯身拾取动作的检测结果，图 9(c)是Th取0.01时，本算法对俯身拾取动作的检测结果。从图 9(c)可以看出，当Th取值偏小时，因为俯身拾取动作在下蹲过程中可能其速度大于设定的阈值，所以有可能将俯身拾取识别为摔倒。而当Th取0.03时，可有效过滤此类动作，并且能够检测出真实的摔倒动作。因为在本实验环境下，阈值Th为0.03时，准确率最高，所以取0.03作为本实验环境下Th的最佳值。因为本文策略以关键点瞬时位移速度大于阈值为前提，辅以后续姿态加以判定，而快速下蹲动作不满足其后续姿态判定中对撑地、坐地或躺下的要求，所以实验中没有出现快速下蹲被误识别为摔倒的情况。从实验结果来看，本算法对快速下蹲、缓慢坐地、俯身拾取这样的类摔倒动作具有一定的抗干扰性。
本实验的目的是对真实摔倒动作与类摔倒动作进行区分，可以根据当前Th得到的准确率变化情况更新Th值。如果本轮Th得到的准确率高于上一轮的准确率，则表明Th的变化方向是有效的，反之则表明上一轮Th值更接近最优值，转而反方向寻值。在迭代过程中，当准确率变化趋势发生改变时，通过改变步长逼近最佳阈值。阈值Th的更新计算公式如式(9)所示：
$$ \text{T}{\text{h}}_{j+1}=\text{T}{\text{h}}_{j}+a\ast （-2）{}^{-m}\text{，}j\ge 0 $$ (9) 式中：Thj表示第j轮迭代中的阈值；a*(－2)－m表示阈值的变化步长。本实验中a取值0.04，m的初始值设为1，记阈值为Thj时的检测准确率为Accj。
更新过程步骤如下：
1）设置初始阈值Th0＝0.05，此时j＝0，按式(8)计算对应的检测准确率Acc0；
2）设m＝1，a＝0.04，按公式(9)计算Thj+1，并计算对应的Accj+1；
3）若Accj+1＞Accj，则j＝j+1，再按公式(9)计算下一代阈值Thj+1，并计算Accj+1；否则，m＝m+1，j＝j+1，再按公式(9)计算Accj+1；
4）重复执行步骤3）。
5. 总结
针对现有人体摔倒检测方法对类摔倒动作容易误检的问题，本文提出基于改进Alphapose的红外图像摔倒检测算法，使用改进的Alphapose算法检测红外图像中的人体关键点和姿态类别，结合人体关键点信息和姿态类别进行摔倒判定。本文在原算法的基础上，使用YOLO v5s作为Alphapose的人体检测器和姿态分类器，在提取人体目标框的同时直接完成对人体姿态的分类，在摔倒与类摔倒场景中，对人体摔倒检测的准确率达到98%。本算法也存在着不足，如本算法利用人体关键点的位置信息做关键点分析，适用于无遮挡状态下的人体摔倒检测。阈值的选取与算法在不同设备上的运行速度有关，对阈值的自适应取值是后续研究的重点。
-
![]()
图 4 基于PSO-Otsu和Lévy_SM_SHO-Otsu算法的断路器瓷套二阈值分割结果
Figure 4. Two threshold segmentation results of ceramic sleeve of circuit breaker based on PSO-Otsu and Lévy_SM_SHO-Otsu algorithms
![]()
图 5 基于PSO-Otsu和Lévy_SM_SHO-Otsu算法的断路器静触头二阈值分割结果
Figure 5. Two threshold segmentation results of circuit breaker static contacts based on PSO-Otsu and Lévy_SM_SHO-Otsu algorithms
表 1 测试函数
Table 1 Test functions
Function Expression Dimension Search range Minimum F5 ${f_{\rm{5}}}{\rm{(}}x{\rm{) = }}\sum\limits_{i{\rm{ = 1}}}^n {{\rm{[100(}}{x_{i{\rm{ + 1}}}} - {x_i}^{\rm{2}}{{\rm{)}}^{\rm{2}}}{\rm{ + (}}{x_i} - {\rm{1}}{{\rm{)}}^{\rm{2}}}{\rm{]}}} $ 30 [-30, 30] 0 F6 ${f_6}(x) = \sum\limits_{i = 1}^n {{{(|{x_i} + 0.5|)}^2}} $ 30 [-100, 100] 0 F13 $\begin{gathered} {f_{13}}(x) = 0.1\{ {\sin ^3}(3{\rm{ \mathsf{ π} }}{x_1}) + \sum\limits_{i = 1}^n {{{({x_i} - 1)}^2}} [1 + {\sin ^2}(3{\rm{ \mathsf{ π} }}{x_i})] \\ \;\;\;\;\;\;\;\;\;\;\; + {({x_n} - 1)^2}[1 + {\sin ^2}(2{\rm{ \mathsf{ π} }}{x_n})]\} + \sum\limits_{i = 1}^n {u({x_i}, 5, 100, 4)} \\ \end{gathered} $
$u({x_i}, a, k, m) = \left\{ {\begin{array}{*{20}{c}} {k{{({x_i} - a)}^m}, \;\;\;{x_i} > a\;\;\;\;\;} \\ {0, \;\;\;\;\;\;\;\; - a \leqslant {x_i} \leqslant a\;\;} \\ {k{{( - {x_i} - a)}^m}, \;\;{x_i} < a\;\;} \end{array}} \right.$30 [-50, 50] 0 F16 ${f_{16}}(x) = 4{x_1}^2 - 2.1{x_1}^4 + \frac{1}{3}{x_1}^6 + {x_1}{x_2} - 4{x_2}^2 + 4{x_2}^4$ 2 [-5, 5] -1.0316 F17 ${f_{17}}(x) = {({x_2} - \frac{{5.1}}{{4{{\rm{ \mathsf{ π} }}^2}}}{x_1}^2 + \frac{5}{{\rm{ \mathsf{ π} }}}{x_1} - 6)^2} + 10(1 - \frac{1}{{8{\rm{ \mathsf{ π} }}}})\cos {x_1} + 10\;\;\;\;\;\;\;\;\;\;\;$ 2 [-5, 5] 0.398 F20 ${f_{20}}(x) = - \sum\limits_{i = 1}^4 {{c_i}} \exp ( - \sum\limits_{j = 1}^3 {{a_i}_j{{({x_j} - {p_{ij}})}^2}} )$ 6 [0,1] -3.3 
下载: 导出CSV
表 2 4种算法在测试函数上的测试结果
Table 2 Test results of 4 algorithms on test functions
Function Indicators SHO Lévy_SHO SM_SHO Lévy_SM_SHO F5 Maximum 28.9838 29.6218 29 29.0713 Minimum 28.7005 28.7027 28.5028 28.4771 Average 28.8738 28.9371 28.7862 28.7253 Standard 0.1025 0.0799 0.0947 0.1582 F6 Maximum 7.06 7.02 7.50 6.19 Minimum 0.326 0.0421 0.0182 0.003723 Average 5.23 4.710 3.26 3.19 Standard 1.93 2.29 2.92 2.25 F13 Maximum 2.99 3.2709 3 3.0271 Minimum 2.86 0.0046 2.7835 0.0023 Average 2.95 1.8008 2.8493 1.6320 Standard 0.0302 1.1731 0.0559 1.1147 F16 Maximum -0.11 -1.0101 0 -1.0094 Minimum -1.03 -1.0316 -1.0316 -1.0316 Average -0.94 -1.0252 -0.9882 -1.0304 Standard 0.147 0.0057 0.1444 0.0031 F17 Maximum 3.597 0.7162 0.51318 0.4690 Minimum 0.398 0.3982 0.3979 0.3979 Average 0.642 0.4282 0.4060 0.4031 Standard 0.484 0.0528 0.0201 0.0108 F20 Maximum -1.6 -2.0591 -2.7087 -2.9141 Minimum -3.1 -3.0991 -3.3031 -3.3047 Average -2.6 -2.7685 -3.167 -3.1990 Standard 0.3 0.2243 0.145831 0.09829
下载: 导出CSV
表 3 基于PSO-Otsu和Lévy_SM_SHO-Otsu的最佳阈值
Table 3 The optimal thresholds based on PSO-Otsu and Lévy_SM_SHO-Otsu algorithms
Infrared image The number of threshold PSO-Otsu Lévy_SM_SHO-Otsu threshold Abnormal temperature distribution in porcelain sleeve of circuit breaker 1 70 70 2 44, 145 46, 143 3 73, 74, 255 70, 76, 255 4 46, 138, 173, 255 46, 73, 144, 255 5 69, 88, 181, 255, 255 46, 54, 132, 139, 255 6 58, 142, 167, 181, 255, 255 53, 65, 76, 144, 255, 255 The circuit breaker still touches the hair to heat 1 63 63 2 51, 147 57, 146 3 44, 72, 255 63, 71, 255 4 38, 80, 123, 255 37, 63, 121, 255 5 58, 70, 135, 150, 255 55, 62, 72, 255, 255 6 54, 71, 97, 186, 246, 255 58, 72, 82, 138, 245, 255
下载: 导出CSV
表 4 基于PSO-Otsu和Lévy_SM_SHO-Otsu算法的适应度函数
Table 4 Fitness function based on PSO-Otsu and Lévy_SM_SHO-Otsu algorithms
Infrared image The number of threshold PSO-Otsu Lévy_SM_SHO-Otsu The value of fitness functions Load switch 1 851.7254 851.7254 2 1096.1 1096.2 3 1631.5 1633.5 4 1746.8 1877 5 2096.3 2197.8 6 2255.8 2614.6 Load switch 1 1282.8 1282.8 2 1486.6 1487 3 2305.5 2362.6 4 2521.7 2565.1 5 2763.8 3437.9 6 3182.6 3635.9
下载: 导出CSV
表 5 基于PSO-Otsu和Lévy_SM_SHO-Otsu算法的PSNR和SSIM值
Table 5 PSNR and SSIM values based on PSO-Otsu and Lévy_SM_SHO-Otsu algorithms
Infrared image The number of threshold PSO-Otsu Lévy_SM_SHO-Otsu PSNR SSIM PSNR SSIM Load switch 1 18.6778 0.0981 18.6778 0.0981 2 21.0977 0.1753 21.1041 0.1794 3 18.7711 0.0969 18.9324 0.1011 4 21.7095 0.1741 22.3763 0.1806 5 21.3552 0.1266 21.5929 0.1760 6 21.3714 0.1370 21.9312 0.1580 Load switch 1 19.1335 0.1567 19.1335 0.1567 2 20.7128 0.1686 21.2433 0.1720 3 20, 0060 0.1738 20.1006 0.1741 4 22.5955 0.1882 22.6117 0.2747 5 22.6896 0.1892 22.7556 0.1908 6 24.0231 0.1900 24.1233 0.1964
下载: 导出CSV
-
[1] Maryam M, Majid H, Fatemeh N. Air condition's PID controller fine-tuning using artificial neural networks and genetic algorithms[J]. Computers, 2018, 7(2): 32. DOI: 10.3390/computers7020032
[2] Ayala H V H, Coelho L D S. Tuning of PID controller based on a multiobjective genetic algorithm applied to a robotic manipulator[J]. Expert Systems with Applications, 2012, 39(10): 8968-8974. DOI: 10.1016/j.eswa.2012.02.027
[3] Beyer H G, Schwefel H P. Evolution strategies——a comprehensive introduction[J]. Natural Computing, 2002, 1: 3-52. DOI: 10.1023/A:1015059928466
[4] ZHU W, DUAN H. Chaotic predator–prey biogeography-based optimization approach for UCAV path planning[J]. Aerospace Science & Technology, 2014, 32(1): 153-161. http://or.nsfc.gov.cn/bitstream/00001903-5/99973/1/1000009426400.pdf
[5] Ghosh A, Das S, Chowdhury A, et. al. An improved differential evolution algorithm with fitness-based adaptation of the control parameters[J]. Information Sciences, 2011, 181(18): 3749-3765. DOI: 10.1016/j.ins.2011.03.010
[6] 孙云霞, 刘兆刚, 董灵波. 基于模拟退火算法逆转搜索的森林空间经营规划[J]. 林业科学, 2019, 55(11): 52-62. DOI: 10.11707/j.1001-7488.20191107 SUN Yunxia, LIU Zhaogang, DONG Lingbo. Spatial forest management planning based on reversion search technique of simulated annealing algorithm[J]. Scientia Silvae Sinicae, 2019, 55(11): 52-62 DOI: 10.11707/j.1001-7488.20191107
[7] Rashedi E, Nezamabadi-Pour H, Saryazdi S. BGSA: Binary gravitational search algorithm[J]. Natural Computing, 2010, 9(3): 727-745. DOI: 10.1007/s11047-009-9175-3
[8] PAN Q K, WANG L, GAO L. A Chaotic Harmony Search Algorithm for the Flow Shop Scheduling Problem with Limited Buffers[M]. Elsevier Science Publishers B. V., 2011.
[9] Shaikh N F, Doye D D. An adaptive central force optimization (ACFO) and feed forward back propagation neural network (FFBNN) based iris recognition system[J]. Journal of Intelligent and Fuzzy Systems, 2016, 30(4): 2083-2094. DOI: 10.3233/IFS-151921
[10] 金旭旸. 基于莱维飞行的水波优化算法[J]. 科技创新与生产力, 2019(5): 66-68. DOI: 10.3969/j.issn.1674-9146.2019.05.066 JIN Xuyang. Water wave optimization algorithm based on Lévy flight[J]. Sci-tech Innovation and Productivity, 2019(5): 66-68. DOI: 10.3969/j.issn.1674-9146.2019.05.066
[11] HUANG H, YANG X, HAO Z, et al. A novel ACO algorithm with adaptive parameter[C]//International Conference on Intelligent Computing on Lecture Notes in Computer Science, 2006, 4115: 12-21.
[12] NIU B, ZHU Y, HU K, et al. A novel particle swarm optimizer using optimal foraging theory[C]//International Conference on Intelligent Computing on Computational Intelligence and Bioinformatics, 2006, 4115: 61-71.
[13] Akay B, Karaboga D. A modified Artificial Bee Colony algorithm for real-parameter optimization[J]. Information Sciences, 2012, 192(1): DOI: 10.1016/j.ins.2010.07.015.
[14] QIANG Z, LI H, LIU C, et al. A new extreme learning machine optimized by firefly algorithm[C/OL]//Proceedings of the 2013 Sixth International Symposium on Computational Intelligence and Design of IEEE, 2013: https://ieeexplore.ieee.org/search/searchresult.jsp?newsearch=true&queryText=A%20new%20extreme%20learning%20machine%20optimized%20by%20firefly%20algorithm.
[15] Mukherjee A, Mukherjee V. A solution to optimal power flow with DC link placement problem using chaotic krill herd algorithm[C]// International Conference on Emerging Technological Trends of IEEE, 2016: DOI: 10.1109/ICETT.2016.7873756
[16] JIANG T, ZHANG C. Application of grey wolf optimization for solving combinatorial problems: job shop and flexible job shop scheduling cases[J]. IEEE Access, 2018, 6: 26231- 26240. DOI: 10.1109/ACCESS.2018.2833552
[17] Dhiman G, Kumar V. Spotted hyena optimizer: a novel bio-inspired based metaheuristic technique for engineering applications[J]. Advances in Engineering Software, 2017, 114: 48-70. DOI: 10.1016/j.advengsoft.2017.05.014
[18] Dhiman G, Kaur A. Spotted hyena optimizer for solving engineering design problems[C]// 2017 International Conference on Machine Learning and Data Science (MLDS)of IEEE, 2017: 114-119.
[19] Dhiman G, Kumar V. Spotted Hyena Optimizer for Solving Complex and Non-linear Constrained Engineering Problems[M]//Harmony Search and Nature Inspired Optimization Algorithms, Advances in Intelligent Systems and Computing, Singapore: Springer Singapore, 2019, 741: 857-867.
[20] 贾鹤鸣, 姜子超, 彭晓旭, 等. 基于改进鬣狗优化算法的多阈值彩色图像分割[J]. 计算机应用与软件, 2020, 37(5): 261-267. DOI: 10.3969/j.issn.1000-386x.2020.05.045 JIA Heming, JIANG Zichao, PENG Xiaoxu, et al. Multi-threshold color image segmentation based on improved hyena optimization algorithm[J]. Computer Applications and Software, 2020, 37(5): 261-267. DOI: 10.3969/j.issn.1000-386x.2020.05.045
[21] 孙悦, 何同祥. 基于莱维飞行的改进蚁群算法的PlD参数优化[J]. 仪器仪表用户, 2019, 26(6): 83-85. DOI: 10.3969/j.issn.1671-1041.2019.06.024 SUN Yue, HE Tongxiang. Optimization of PID parameters based on improved ant colony algorithm for Lévy flight[J]. Electronic Instrumentation Customers, 2019, 26(6): 83-85. DOI: 10.3969/j.issn.1671-1041.2019.06.024
[22] 张烈平, 何佳洁, 于滟琳, 等. 基于蚁群算法优化的布谷鸟搜索算法[J]. 微电子学与计算机, 2018, 35(12): 27-32. https://www.cnki.com.cn/Article/CJFDTOTAL-WXYJ201812005.htm ZHANG Lieping, HE Jiajie, YU Yanlin, et al. A cuckoo search algorithm based on ant colony algorithm optimization[J]. Microelectronics & Computer, 2018, 35(12): 27-32. https://www.cnki.com.cn/Article/CJFDTOTAL-WXYJ201812005.htm
[23] 张新, 李珂, 严大虎, 等. 改进入侵杂草算法求解柔性作业车间调度问题[J]. 系统仿真学报, 2018, 30(11): 446-453. https://www.cnki.com.cn/Article/CJFDTOTAL-XTFZ201811051.htm ZHANG Xin, LI Ke, YAN Dahu et al. Improved intrusion weed algorithm for solving flexible job shop scheduling problem[J]. Journal of System Simulation, 2018, 30(11): 446-453. https://www.cnki.com.cn/Article/CJFDTOTAL-XTFZ201811051.htm
[24] 赵洪, 李伟鹏, 刘铁军. 基于改进莱维飞行的狼群算法及其在翼型气动优化设计中的应用[J]. 科学技术与工程, 2019, 19(18): 315-323. DOI: 10.3969/j.issn.1671-1815.2019.18.048 ZHAO Hong, LI Weipeng, LIU Tiejun. An improved Lévy flight based grey wolf optimization algorithm for aerodynamic design problem[J]. Science Technology and Engineering, 2019, 19(18): 315-323. DOI: 10.3969/j.issn.1671-1815.2019.18.048
[25] 肖石林. 基于Lévy飞行的树种优化算法及在图像分割中的应用[D]. 南宁: 广西民族大学, 2019. XIAO Shilin. Tree Seed Optimization Algorithm Based on LévyFlight and Its Application in Image Segmentation[D]. Nanning: Guangxi University for Nationalities, 2019.
[26] 莫愿斌, 郑巧燕, 马彦追. 单纯形法的布谷鸟搜索算法及其在约束优化问题中的应用[J]. 计算机与应用化学, 2015(2): 213-218. https://www.cnki.com.cn/Article/CJFDTOTAL-JSYH201502019.htm MO Yuanbin, ZHENG Qiaoyan, MA Yanzhui. Cuckoo search based on simplex method and its application on constrained optimization problems[J]. Computers and Applied Chemistry, 2015(2): 213-218. https://www.cnki.com.cn/Article/CJFDTOTAL-JSYH201502019.htm
[27] 莫愿斌, 马彦追, 郑巧燕, 等. 单纯形法的改进萤火虫算法及其在非线性方程组求解中的应用[J]. 智能系统学报, 2014(6): 747-755. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNXT201406020.htm MO Yuanbin, MA Yanzhui, ZHENG Qiaoyan, et al. Improved firefly algorithm based on simplex method and its application in solving non-linear equation groups[J]. CAAI Transactions on Intelligent Systems, 2014(6): 747-755. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNXT201406020.htm
[28] 张红霞, 罗毅, 师瑞峰. 基于单纯形法的改进型人工鱼群算法[J]. 计算机应用, 2011(5): 1321-1323. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJY201105048.htm ZHANG Hongxia, LUO Yi, SHI Ruifeng. Artificial fish swarm algorithm based on simplex method[J]. Journal of Computer Applications, 2011(5): 1321-1323. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJY201105048.htm
[29] 肖辉辉. 基于单纯形法和自适应步长的花朵授粉算法[J]. 计算机工程与科学, 2016, 38(10): 2126-2133. DOI: 10.3969/j.issn.1007-130X.2016.10.025 XIAO Huihui. A flower pollination algorithm based on simplex method and self-adaptive step[J]. Computer Engineering & Science, 2016, 38(10): 2126-2133. DOI: 10.3969/j.issn.1007-130X.2016.10.025
计量
- 文章访问数: 228
- HTML全文浏览量: 72
- PDF下载量: 24
