基于FFWA的自适应Canny飞机蒙皮红外图像边缘检测

王 坤，刘沛伦，王 力

基于FFWA的自适应Canny飞机蒙皮红外图像边缘检测

王 坤，刘沛伦，王 力

（中国民航大学，电子信息与自动化学院，天津 300300）

针对传统自适应Canny算法阈值选取精度低、速度慢的不足，提出基于快速烟花算法（Fast Fireworks Algorithm，FFWA）的自适应Canny边缘检测算法，该算法采用最大类间方差法结合快速烟花算法对检测和连接边缘的高低阈值进行自动设定。快速烟花算法对传统烟花算法的爆炸半径，爆炸火花产生方式与选择策略进行改进。实验结果表明，快速烟花算法比传统烟花算法的计算时间节省了36%，在稳定性方面也有了可观的提升，基于快速烟花算法的自适应Canny边缘检测算法在精度保持不变的情况下，计算速度比改进前加快了49%，使飞机蒙皮损伤热像图的边缘检测效果更加理想。

边缘检测；烟花算法；Canny算法；最大类间方差法；自适应

0 引言

蒙皮是飞机运行中最容易出现损伤的零部件，因其直接与外部环境相接触，很容易产生裂纹、积水、腐蚀等故障。这些故障大都难以用肉眼察觉，却又直接影响着飞机的飞行安全[1]。因此，如何对飞机蒙皮的故障进行快速准确地边缘检测显得尤为重要。目前的检测手段主要有：射线检测、涡流检测、磁粉检测、渗透检测和超声检测等[2]。红外热成像检测有着范围广、速度快、稳定非接触等优点，已成为重要的检测手段[3]。但红外图像分辨率低且灰度较为集中，如何准确快速地对红外图像进行边缘提取尤为重要。

Canny算法[4]因其较优的信噪比和检测精度，在图像处理方面得到了广泛应用。但传统的Canny算法在双阈值的选择上缺乏自适应性。因此，众多自适应Canny算法被提出，李二森等[5]提出的自适应Canny边缘检测算法利用全局阈值分割法来寻找动态阈值，但此法精确度较低。杜磊等[6]采用图片信息熵最大时的阈值作为图像的自适应阈值。郭方方等[7]和宋人杰等[8]利用了最大类间方差法（Otsu）来确定Canny的分割阈值。这种方法虽然阈值精度较高，但计算量大，在处理大型图像时运算缓慢。考虑引入烟花算法来减少最大类间方差法（Otsu）的计算时间。

烟花算法（Fireworks Algorithm，FWA）[9-11]是Tan等在2010年提出的一种先进的、高效的群体智能算法，具有全局探索、局部暴发、分布并行等优点。但烟花算法也存在一定的不足，总体来说主要有以下4个方面：①爆炸半径机制缺乏自适应性，从而导致适应度值最好的烟花爆炸半径会被限制在极小的范围，致使具备最多爆炸火花数的最优烟花几乎丧失局部挖掘的能力，浪费计算资源。②当前适应度值最优的火花被直接保留为下一代的烟花，导致算法一旦陷入局部最优则很难跳出。③高斯变异火花在增加解的多样性方面作用有限，反而耗费了大量的计算时间。④每个烟花爆炸后依然作为下一代的火花，导致重复计算，降低了计算效率。

近些年来，科研工作者针对上述问题投入了大量精力，提出了众多改进型的烟花算法，并表现出良好的性能。王亮等[12]在算法的搜索过程中引入了柯西变异算子，并改进了选择策略，使新的算法具备了更多样的爆炸火花与更精确的收敛精度。Shaoqiu Zheng等[13]针对FWA的选择策略做出改进，给烟花之间加入了合作机制，使算法在全局探索和收敛速度方面有了进步。张水平等[14]的带有动态爆炸半径的增强型烟花算法（Enhanced Fireworks Algorithm with Dynamic Explosion Radius，EFWA-DER）改进了传统烟花算法的基本算子并加入爆炸半径动态调整的策略，提高了算法的优化性能。Li等[15]的导引烟花算法，在传统烟花算法中引入导引火花，提高了算法的搜索效率。陶小华等[16]提出的具有导向功能的改进烟花算法（improved FWA with directional function，IDFWA）充分利用了当代最优个体的位置信息以及历代的全局最优解，改进了传统烟花算法的搜索方式，消除了爆炸半径对搜索范围的限制。以此来使烟花具有更好的搜索能力。Li J.等[17]在改进型烟花算法（Improved Fireworks Algorithm，IFWA）中提出新的爆炸半径与火花产生方式，以此来提高算法的寻优能力。以上改进型烟花算法都从一定程度上使爆炸火花呈现出多样性，增强算法跳出局部最优的能力。但烟花算法陷入局部最优的主要原因在于算法在下一代的选择策略上都无条件地将当前适应度值最好的火花（并不一定是全局最优解范围内的火花）列入下一代烟花。而这类火花在下一代的爆炸半径又会非常小，所以算法会无法跳出局部最优，不利于全局探索。基于此，提出快速烟花算法（Fast Fireworks Algorithm，FFWA），该算法采用不同的爆炸半径产生方式，使核心烟花的爆炸半径更加合理，增强了局部挖掘能力；剔除了烟花爆炸序列中本烟花作为第一个爆炸火花的机制，减少了算法的运算量，提高了算法的收敛速度；并且在选择策略上增加了判断机制，使算法可以更加轻易地跳出局部最优。

本文主要致力于提升自适应Canny算法中分割阈值的计算速度，引入烟花算法来对图像的最大类间方差进行寻优，同时针对传统烟花算法会出现早熟和计算效率低的缺陷，对其爆炸半径和爆炸火花的生成方式与次代烟花的选择策略做出改变，提出快速烟花算法，最后通过实验分析了快速烟花算法和其他寻优算法的表现，对各类自适应Canny算法的计算速度与边缘提取效果进行了比较。

1 自适应Canny边缘检测

1.1 传统Canny边缘检测

传统的Canny算法主要分为4个部分。首先对图像的噪声进行滤波（通常为高斯滤波），然后对去噪后的图像计算其梯度幅值与方向，其次对图像的梯度幅值进行非极大值抑制，最后进行双阈值判定检测边缘。

Canny算法中，边缘的检测和连接需要根据双阈值（高阈值THh，低阈值THl）来判断，其中高阈值主要区分背景与前景，而低阈值主要平滑边缘，将不连续的部分连接起来。通常低阈值等于1/2～1/3的高阈值，但高阈值的选取则需要按照经验预设或通过多次实验不断改进。可利用最大类间方差法实现对THh的计算。

1.2 最大类间方差法（Otsu）

最大类间方差法由日本学者大津提出，是根据最小二乘原理得出的一种可以自动选取阈值的算法。在最大类间方差法中，用阈值将图像分为目标和背景两部分。目标和背景之间的方差值越大，表示组成图像的两部分区分度越高。当一部分背景被错分为目标或者是当目标被错分为背景，则会导致两部分的区分度降低。所以，类间方差最大，则表明错分的概率最小。假设要分割的图像的灰度范围是{0, 1, …,－1}。通过阈值图像像素将被按式(1)和(2)分为两个部分：

0＝{0, 1, …,} (1)

1＝{＋1,＋2, …,－1} (2)

0和1分别代表目标和背景。这类图像的误差在0和1之间：

()2＝0()*1()*[0()－1()]2(3)

式中：是阈值；0()是图像灰度值小于阈值的像素数量。1()是图像灰度值大于阈值的像素数量。0()是图像灰度值小于阈值的像素平均灰度值。1()是图像灰度值大于阈值的像素平均灰度值，使得()2最大值的max是最佳分割阈值。

max作为图像前景与背景的最佳分割阈值，在本质上与Canny算法中的高阈值THh作用是一样的，可以令THh＝max。通过多次实验对比得到，THh＝2.5THl时算法提取的边缘最准确与完整，令THl＝max/2.5，以此来连接边缘，获得边缘检测图像。但Otsu计算量较大，耗费时间较长，可通过结合快速烟花算法解决此问题。

2 快速烟花算法

2.1 传统烟花算法

燃放烟花爆竹是中国传统节日尤其是除夕的一项重要节日庆祝活动，点燃的烟花被发射到夜空中，爆炸产生火花继而照亮其临近的夜空，产生出一幅美丽的图案。在一个优化问题中，烟花被看作为最优化问题的解空间中一个可行解，那么烟花爆炸产生一定数量火花的过程即为其搜索邻域的过程[9]。烟花算法由爆炸算子、变异算子、映射规则和选择策略4部分组成。烟花的爆炸示意图如图1所示。

图1 烟花爆炸示意图

如图1所示，烟花爆炸产生火花，然后在所有的烟花与火花中按照选择策略选择下一代烟花进行爆炸操作，依次循环迭代，在解空间中寻找最优解。

2.1.1 爆炸算子

在可行解空间随机产生一些烟花，对其初始化后利用适应度函数计算出适应度值并进行评估，则第（＝1, 2, …,）个烟花爆炸产生的火花数目为：

式中：E是控制个烟花产生火花总数的参数；(x)为第个烟花的适应度值。max＝max(())表示烟花所对应适应度的最大值。是为避免除零操作引入的一个极小正数。为避免生成的S过多或过少，还需对S进行如下限制：

式中：，为爆炸数目限制因子。

根据烟花适应度的好坏，计算不同烟花的爆炸半径，适应度好的烟花在邻近区域生成较多的火花，适应度差的烟花在较大范围内产生较少的火花。计算公式如下：

式中：r为预设的最大爆炸半径；min＝min(())为烟花所对应适应度值的最小值；为一个极小的正数。

则x生成的第(＝1, 2, …,S)个火花的第(＝1, 2, …,)维度上坐标为：

式中：rand(－1,1)为－1～1之间的随机数。

2.1.2 变异算子

在传统烟花算法中产生第(＝1, 2,…,)个高斯火花的第(＝1, 2, …,)维度上的坐标为：

式中：Gaussian(1,1)代表平均值为1；标准方差为1的高斯分布。

2.1.3 映射规则

在产生爆炸火花或变异火花时，有可能更新后的坐标超出了取值范围，这时需要采用一定的映射规则来对越界火花进行修正：

式中：LB,k和UB,k分别表示该优化问题中烟花取值的下边界与上边界。

2.1.4 选择策略

在每一次迭代中，所有的个体（烟花和火花）中，适应度最好的烟花无条件保留到下一代，剩余的－1个个体采用轮盘赌的方式选取，每一个火花被选中的概率为：

式中：(X)为烟花X与其他烟花的距离之和；(X,X)为烟花X与烟花X之间的欧式距离；为所有烟花与火花数目之和。

烟花算法的具体执行步骤为：

步骤1 在解空间随机生成一定数量烟花，每个烟花都是一个可行解。

步骤2 根据每个烟花所在位置的适应度值，按照各自的爆炸数目（式(4)）爆炸半径（式(6)）进行爆炸。同时，随机选取一定数量的烟花进行高斯变异。

步骤3 计算所有的火花的适应度值并排序。

步骤4 判断是否达到迭代次数，如果是，则输出最优值，否则，选取下一代烟花。

步骤5 重复以上2, 3, 4步，直到达到迭代次数。

传统烟花算法的执行流程如图2所示。

图2 传统烟花算法流程图

如图2所示，在传统烟花算法中，火花分为爆炸火花和高斯火花，前者主要负责对解空间进行探索和挖掘，适应度好的烟花在小范围内生成较多的火花，适应度差的烟花在大范围内生成较少的火花。而高斯火花的作用主要是增加种群多样性。

2.2 快速烟花算法

烟花算法的选择策略为当前最优（适应度值最优）火花作为下一代烟花，其余－1个火花的选择方式有随机选取、轮盘赌的方式选取、按照适应度值排序选取等。这些选择策略都无一例外地直接保留了当前最优火花，而该火花（也是拥有爆炸火花数最多的烟花）在下一代的爆炸半径几乎等于0，这样不仅没有起到局部挖掘的作用，并且会在其他烟花没有找到更优点的情况下陷入局部值。因此，对烟花算法的选择策略增加一个判断机制，该机制如下：

每一代的当前最优火花必须满足以下条件之一时，才可以被选择为下一代烟花。

1）当前最优火花适应度优于上一代最优火花时。

2）当前最优火花适应度劣于上一代最优火花，且满足接受概率时。

为引入以上判断机制，需对烟花算法的爆炸半径计算方式与爆炸火花产生方式进行改变。

2.2.1 爆炸半径计算方式的改进

在烟花算法中，烟花主要分为两类，靠近极值点的烟花（当前适应度值较好的烟花）和远离极值点的烟花（当前适应度值较差的烟花），在算法的计算过程中，两类烟花分别发挥着不同的作用，前者主要在小范围内发挥其局部挖掘的能力，后者则在更大的范围内发挥探索的能力[18]。通过计算不难发现，传统烟花算法的爆炸半径计算公式会使拥有爆炸数量最多的最优烟花b的爆炸半径非常小甚至等于0，这会使b的爆炸火花出现大量重复计算，也是算法容易早熟的原因之一。因此，郑少秋[19]等在EFWA（enhanced fireworks algorithm）中增加了最小爆炸半径检测机制，使每一个烟花的爆炸半径不会小于min。增强了在算法早期烟花的探索能力。但这种限制最小爆炸半径的策略过分依赖于算法预先设定的最大迭代次数，而这是人为设定的。事实上，一个优秀的爆炸半径策略应该让两类烟花各自发挥其能力即核心烟花更好的发掘，非核心烟花更好的探索，并且也应该考虑到计算空间的整体信息和当前迭代的次数。基于此，提出一种新的爆炸半径计算公式：

式中：|X,k－X,k|为第维上当前烟花距离最优烟花的绝对距离，为当前迭代次数。

式(13)主要分为两部分，第一部分为当前烟花距离最优烟花的绝对距离与空间上下限的比值，全局信息的引入使爆炸半径的分配更具整体性。第二部分为当前迭代次数的倒数，这一部分主要是保证在算法前期各类烟花均具备良好的探索能力，并且随着算法深入，逐渐降低对爆炸半径的影响，使第一部分发挥优势。图3为FFWA在Sphere函数上自适应爆炸半径的取值。

图3 FFWA在Sphere函数上自适应爆炸半径的取值

由图3可知，爆炸半径的取值整体呈现非线性递减的趋势，在算法前期，核心烟花的爆炸半径足够大，可以快速锁定全局最优值的范围，后期爆炸半径在合理的范围内又足够小，充分发挥了其局部挖掘的能力。而非核心烟花根据与核心烟花的绝对距离与迭代次数的改变，自适应地选取爆炸半径，合理地在全局进行搜索，使计算过程快速收敛。

2.2.2 爆炸火花产生方式的改进

对于传统烟花算法（包括绝大多数改进型）为了保留代间信息，都选择将每一个烟花本身作为下一代爆炸火花的第一个火花。这导致每一代爆炸火花的适应度值评估都会因重复计算而耗费大量时间，当然并不能说这些重复计算是无意义的，因为在传统烟花算法对后代烟花的选择策略的前提下，为了保留当前最优与后期收敛，不得不这样做。但如果在选择策略中加入当前最优适应度值保存机制，那么就完全不必要进行这些重复计算。因此，新的爆炸方式中，每一个烟花的爆炸火花序列将不再存储该烟花本身，而是全部生成新的火花，改进后，每一代需要评估的火花数由原先的＝sum(S)＋个减少为¢=sum(S)个，从而使算法在没有降低全局探索能力的前提下，计算效率得到稳定的提升，且这种改进对各类适应度函数存在普遍适用性。

2.2.3 选择策略的改进

为了可以避免上一代的火花重复计算浪费时间，又可以使算法在每一代的迭代中保留上一代最优火花的信息，对每一代的最优火花进行一次判断选择机制，即每一代的最优火花与当代的最优烟花进行比对，如果适应度更好，则保留为下一代的烟花，反之，则按照一定概率（接受概率）接受。为使结果最终收敛，接受概率也随迭代次数逐渐收缩，最终接近于零。其余火花按照适应度排序择优选取。在理论上，这种新的选择策略既保证了最优值的保存，又使算法更加灵活，避免陷入局部最值，收敛稳定，此外，在改进后的选择策略中，其余烟花的选取不再依靠欧式距离的计算，又进一步降低了算法的运行时间。

快速烟花算法的具体步骤为：

步骤1 在解空间内随机设置一定数量的烟花，并设置接受概率与收缩系数。

步骤2 根据式(4)与式(12)规定的爆炸数目与爆炸半径进行爆炸。FFWA删除了烟花的高斯变异过程。并且每个烟花的爆炸序列中只储存爆炸火花的信息，不再对烟花本身进行储存。

步骤3 计算每个火花的适应度值并排序。

步骤4 判断是否达到迭代次数，如果是，则输出迭代结果，否则，判断当前最优火花是否优于上一代最优烟花，如果是，则将该火花选为下一代烟花，否则按照接受概率选取。剩余烟花按照适应度排序择优选取。

步骤5 对接受概率进行收缩。

步骤6 重复以上2，3，4步，直到达到迭代次数。

快速烟花算法流程如图4所示。

图4 快速烟花算法流程图

3 基于FFWA的自适应Canny边缘检测

本文用快速烟花算法对图像的类间方差进行寻优，并将求得的最大类间方差作为Canny算法中的高阈值THh，令低阈值THl等于THh/2.5，据此对图像进行边缘的检测与连接，获得图像的边缘。

算法的具体步骤为：

步骤1 对待检测的图像进行高斯去噪，降低噪声对后续计算的影响，防止误检测。

步骤2 对去噪后的图像计算梯度幅值与方向。

步骤3 对梯度幅值进行非极大值抑制，将局部最大值之外的所有梯度值抑制为0。

步骤4 初始化快速烟花算法，将Otsu设置为适应度函数，读入图像，循环迭代，输出最佳分割阈值m。

步骤5 将该阈值赋予Canny算法中的高阈值THh，低阈值THl＝THh/2.5。

步骤6 进行双阈值检测和连接边缘，得到边缘检测图像。

图5为基于FFWA的自适应Canny边缘检测算法流程图。

4 实验分析

实验主要分为两部分，首先对快速烟花算法，传统烟花算法，构造型烟花算法[20]，柯西烟花算法与动态搜索烟花算法[21]的计算性能做出对比，然后分别对基于上述算法的自适应Canny边缘检测结果进行对比。

实验系统环境为：Win7，Intel 2.20GHz CPU i5，8G RAM，编程环境为：Matlab 2016a。

4.1 快速烟花算法的性能对比

为了对改进算法的性能进行评估，实验选取了4个具有全局最优值的标准函数进行对比，如表1所示。此外，共设置了10组不同参数来评估参数对算法的影响，结果表明，不同组的参数设置对算法间的影响微乎其微，仅以表2中参数为例，进行实验。

表2中，为初始烟花数，为最大爆炸半径，为爆炸火花总数目，与为爆炸数目限制系数。4个测试函数如表1所示，每个测试函数分别独立运行20次，迭代次数设定为1000。由于初始烟花的位置对算法结果的影响很大，为了更加客观地进行比较，消除不确定因素，每次对比的各个算法都使用相同初始位置。并记录出现的最优值、最差值、平均值与平均运行时间4项数据对算法进行评估。测试对比结果见表3。

图5 基于FFWA的自适应Canny边缘检测算法流程图

由表3可知，FFWA在1,3,4三个函数中的计算精度均优于其他对照组，2中略差于CFWA。运行时间对比上，每一组的运行时间FFWA都是最低的。综合来看，FFWA的性能优于其他4个对照组。为了更加直观地对算法的收敛过程进行评估，给出4个标准函数的收敛曲线，如图6～图9所示。

表1 基准测试函数

表2 五种算法的参数设置

表3 测试函数测试结果对比

由图6～图9对比可知，因为爆炸半径的改进，使FFWA的核心烟花可以在算法早期快速定位在最优值附近。所以FFWA的收敛速度均优于其他4个对照组，而择优选取的选择策略又使得在算法后期有大量的火花在最优值附近挖掘，保证了算法的精度。值得指出的是，传统烟花算法因为高斯火花和映射关系的存在，使其在对最优值处于原点附近的函数计算时有着良好的表现，但红外图像普遍有着灰度集中，且最佳分割阈值不在零点附近的特点，这就需要算法具备更高的全局搜索能力和更高的计算精度。为验证FFWA在红外图像的性能表现，用红外热像仪对飞机损伤蒙皮实拍得到如图10～图13所示加热结束后冷却不同时刻飞机蒙皮损伤红外热像图。

图6 Sphere函数收敛曲线

图7 Ackley函数收敛曲线

图8 Griewank函数收敛曲线

图9 Six-hump函数收敛曲线

图10 积水损伤加热结束后不同时刻图像

图11 脱粘损伤加热结束后不同时刻图像

图12 裂痕损伤加热结束后不同时刻图像

分别用FFWA、FWA、IFWA、dynFWA和CFWA计算图10～图13中蒙皮损伤图像的边缘检测阈值。并记录在20次迭代过程中各个算法的平均耗时，平均迭代次数和总误差，总误差为20次试验中收敛阈值和遍历得到的真实值差的总和。实验结果在表4给出。

由表4可知，在不同损伤的蒙皮红外图像中，FWA的总误差较大，可见FWA的性能在一定程度上依赖于适应度函数的选取，当适应度函数的全局最优值离零点较远时，可靠性较差。而改进后的FFWA却可以在不同的适应度函数中有着较高的可靠性。计算时间上，FFWA的计算耗时明显低于其他对照组，其中平均比FWA节省了36%，比dynFWA节约了约20%。算法的可靠性方面FFWA，IFWA，CFWA在20次实验中均呈现出较高的可靠性，dynFWA与FWA出现了早熟的情况。

在大量实验中，不同图像中算法的收敛曲线差异性不大，由于篇幅限制，给出两个具有代表性的收敛曲线，分别是＝100s时的积水损伤阈值收敛曲线（图14）与＝100s时的脱粘损伤阈值收敛曲线（图15）。

由图14可知，FWA与dynFWA收敛于局部最优，FFWA，CFWA，IFWA都在计算初期收敛于最优值范围内，且收敛于最佳阈值附近。图15收敛曲线和图14类似，但值得指出的是dynFWA特殊的核心烟花爆炸半径机制使其依赖于初始烟花的位置，如果初始烟花在最优值周围，那么dynFWA无论在速度还是精度方面都有着很良好的表现，如果初始烟花离最优点位置较远，那么最优烟花的爆炸半径会越来越大，使得其火花落在最优值周围的概率很低，会出现收敛非常慢的现象。通过图14、图15可知，改进后的烟花算法可以在较低的代数内收敛于全局最优值。

图13 腐蚀损伤加热结束后不同时刻图像

表4 各个算法的实验结果

图14 积水损伤红外热像图阈值计算的收敛曲线

图15 脱粘损伤红外热像图阈值计算的收敛曲线

4.2 边缘检测结果的对比

针对大量冷却时间不同红外图像进行边缘检测，结果表明在＝100s时的损伤边缘特征最明显，边缘检测效果最优。为验证本文算法的有效性，分别用本文算法，基于全局阈值的自适应Canny算法（文献[5]算法），基于最大熵的自适应Canny算法（文献[6]算法），基于Otsu的自适应Canny算法（文献[7]算法）对＝100s时的积水、脱粘、裂痕与腐蚀图像进行边缘检测，检测结果如图16～图19所示。并以运行时间、分割阈值与Pratt品质因数（Pratt figure of merit，PFOM）[22]来作为算法的评价指标，品质因数用来评价边缘点的定位精度：

式中：A为检测到的边缘点；I为理想边缘点；为理想边缘点与检测到的边缘点连线的垂直距离。是对边缘错位时的惩罚系数（本次实验取0.1）Pratt值越大，则表明检测到的边缘越准确。理想边缘为人工标记的边缘图像。检测结果如表5所示。

由表5可知，本文算法和文献[7]算法具有同样的检测阈值和Pratt品质因数，但本文算法的计算耗时却比文献[7]算法减少了49%，比文献[6]算法减少了52%，比文献[5]算法减少了30%。从图16～图19中可以更加直观地进行对比。文献[5]算法只在腐蚀图中检测结果较好，其他3幅图像均有明显的误判，文献[6]算法基本可以检测出损伤边缘，本文算法与文献[7]算法一样，都对蒙皮损伤边缘做出了较为准确的检测。

5 结论

本文针对传统自适应Canny边缘检测算法检测阈值选取不准确和速度慢的不足，提出基于快速烟花算法的自适应Canny边缘检测，利用快速烟花算法对待检测图像的最大类间方差进行计算，使Canny算法分割阈值获取更加迅速，增强了传统自适应Canny算法的精确性与快速性，改善了Canny算法的检测效果。由于快速烟花算法对传统烟花算法的爆炸半径，产生爆炸火花方式与选择策略提出改进，使改进后的算法在速度与可靠性上又有了一定程度的改善。通过实验数据可知，所提方法具有可行性。

表5 4种算法检测对比结果

图16 飞机蒙皮积水损伤边缘检测结果

Fig.16 Edge detection results of hydrops

图17 飞机蒙皮脱粘损伤边缘检测结果

图18 飞机蒙皮裂痕损伤边缘检测结果

图19 飞机蒙皮腐蚀损伤边缘检测结果

[1] 王昊. 基于机器视觉的飞机蒙皮损伤检测与寿命分析方法[D]. 南京: 南京航空航天大学, 2013.

WANG Hao. Aircraft Skin Damage Detection and Life Cycle Analysis Method Based on Machine Vision[D]. Nanjing: Nanjing University College of Automation Engineering, 2013.

[2] CHENG Sheng, FENG Jianfei, DENG Faqing. A Summary of Nondestructive Testing Application in Chinese Civil Aircraft Industry[J]., 2017, 39(4): 76-79.

[3] Ammar K Al-Musawi, Fatih Anayi, Michael Packianather. Three-phase induction motor fault detection based on thermal image segmentation[J]., 2020, 104: 1350-4495.

[4] SHI H, Ward R. Canny edge based image expansion[C]//, 2002: 1-1. DOI: 10.1109/ ISCAS.2002.1009958.

[5] 李二森, 张保明, 周晓明, 等. 自适应Canny边缘检测算法研究[J]. 测绘科学, 2008, 66(6): 120-121.

LI Ersen, ZHANG Baoming, ZHOU Xiaoming, et al. Research on the algorithm of adaptive Canny edge detection[J]., 2008, 66(6): 120-121.

[6] 杜磊, 李立轻, 汪军, 等. 几种基于图像自适应阈值分割的织物疵点检测方法比较[J]. 纺织学报, 2014, 35(6): 56-60.

DU Lei, LI Liqing, WANG Jun, et al. Comparison of several fabric defect detection methods based on image adaptive threshold segmentation[J]., 2014, 35(6): 56-60.

[7] 郭方方, 严高师, 李旭东. 一种改进的基于Otsu算法的Canny红外边缘检测方法[J]. 红外, 2010, 31(7): 24-27.

GUO Fangfang, YAN Gaoshi, LI Xudong, et al. An Improved Canny Infrared Edge Detection Method Based on Otsu Algorithm[J]., 2010, 31(7): 24-27

[8] 宋人杰, 刘超, 王保军. 一种自适应的Canny边缘检测算法[J]. 南京邮电大学学报: 自然科学版, 2018, 176(3): 76-80.

SONG Renjie, LIU Chao, WANG Baojun. An adaptive Canny edge detection algorithm[J].: Natural Science Edition, 2018, 176(3): 76-80.

[9] TAN Y, ZHU Y. Fireworks Algorithm for Optimization[C]//, 2010: 12-15.

[10] 谭营. 烟花算法引论[M]. 北京: 科学出版社, 2015.

TAN Ying.[M]. Beijing: Science Press, 2015.

[11] TAN Ying, ZHENG Shaoqiu. Recent advances in fireworks algorithm[J]., 2014, 9(5) : 515-528．

[12] 王亮, 郭星. 基于柯西烟花算法的大规模服务组合优化[J]. 计算机工程与应用, 2018, 54(24): 34-40.

WANG Liang, GUO Xing. Large-scale service portfolio optimization based on Cauchy fireworks algorithm[J]., 2018, 54(24): 34-40.

[13] ZHENG S, LI J, Janecek A, et al. A Cooperative Framework for Fireworks Algorithm[J]., 2017, 14(1): 27-41.

[14] 张水平, 李殷俊, 高栋, 等. 带有动态爆炸半径的增强型烟花算法[J]. 计算机工程与应用, 2019, 961(18): 56-63.

ZHANG Shuiping, LI Yinjun, GAO Dong. Enhanced fireworks algorithm with dynamic explosion radius[J]., 2019, 961(18): 56-63.

[15] LI J, ZHENG S, TAN Y. The Effect of Information Utilization: Introducing a Novel Guiding Spark in the Fireworks Algorithm[J]., 2017, 21(1): 491-501.

[16] 陶小华, 陈基漓, 谢晓兰. 具有导向功能的改进烟花算法[J]. 计算机工程与设计, 2019, 40(12): 3479-3486.

TAO Xiaohua, CHEN Jili, XIE Xiaolan. Improved fireworks algorithm with guidance function[J]., 2019, 40(12): 3479-3486.

[17] YANG W, KE Liangjun. An improved fireworks algorithm for the capacitated vehicle routing problem[J]., 2019, 13(3): 552-564.

[18] LIU J, ZHENG S, TAN Y.[M]. Advances in Swarm Intelligence Berlin: Springer, 2013: 11-23.

[19] ZHENG S, Janecek A, Tan Y. Enhanced Fireworks Algorithm[C]//(CEC),., 2013: 2069-2077.

[20] RONG Cheng, BAI Yanping, ZHAO Yu, et al. Improved fireworks algorithm with information exchange for function optimization[J]., 2019, 163: 82-90.

[21] ZHENG S, Janecek A, LI J, et al. Dynamic search in fireworks algorithm[C]//2014(CEC), 2014: 3222-3229.

[22] 李姗姗, 陈莉, 张永新, 等. 结合四元数与最小核值相似区的边缘检测[J]. 中国图象图形学报, 2017, 22(7): 915-925.

LI Shanshan, CHEN Li, ZHANG Yongxin, et al. The Edge Detection algorithm Combining Dimension Univalue segment Depreciation and Quaternion[J]., 2017, 22(7): 915-925.

Infrared Image Adaptive Canny Edge-detection of Aircraft Skin Based on Fast Fireworks Algorithm

WANG Kun，LIU Peilun，WANG Li

(College of Electronic Information and Automation, Civil Aviation University of China, Tianjin 300300, China)

To address the low accuracy and slow speed of the traditional adaptive Canny algorithm in selecting the threshold value, an improved algorithm is herein proposed. The proposed algorithm utilizes the Otsu and fast fireworks algorithm (FFWA) to automatically set the high and low detection thresholds and the connecting edges. Consequently, the explosion radius, the production method, and selection strategy are improved as compared to those of the traditional fireworks algorithm, thereby increasing the speed and accuracy of the Otsu calculation. The experimental results show that the calculation speed of the fast fireworks algorithm increased by 36% as compared to that of the traditional fireworks algorithm. Moreover, the stability also increased considerably. The improved adaptive Canny algorithm not only maintains the same precision but also decreases the calculation time by 49%. This makes the result of edge detection of aircraft skin infrared images more ideal.

edge-detection, fireworks algorithm, Canny algorithm, Otsu, adaptive

TP274.52

A

1001-8891(2021)05-0443-12

2020-06-06；

2020-07-13.

王坤（1978-），女，博士，副教授，主要研究方向为图像处理、故障检测分析。E-mail：96288851@qq.com。

国家自然科学基金（U1733119）；中央高校基本科研业务费项目中国民航大学专项（3122018C001）。