Infrared Thermal Image Detection of Faulty Insulators in Distribution Lines Based on Multi-scale Template Matching
-
摘要: 瓷绝缘子在配电线路中应用广泛,受长期机电应力与户外恶劣环境影响,在运行中易发生劣化。红外热像法是一种重要的劣化绝缘子带电检测方法,具有检测方便、安全高效和非接触式的优点,已成为线路巡检的重要手段,但劣化绝缘子热像特征不明显,肉眼识别易出现误判。为此,本文首先对配电线路瓷绝缘子进行温度场仿真分析,然后提出了一种劣化绝缘子红外热像检测方法,采用多尺度模板匹配算法定位识别绝缘子,获取绝缘子红外图像中的坐标参数,并对其进行分割提取,通过最小二乘线性拟合提取绝缘子表面温度。结合相关标准与仿真分析结果,通过同类比较判断法对比多个绝缘子温度状态的差异,实现劣化绝缘子检测。Abstract: Porcelain insulators are widely used in power distribution lines, but they are susceptible to degradation during operation owing to long-term electromechanical stress and harsh outdoor environments. Infrared thermal imaging is an important live insulator degradation detection method. It has the advantages of convenient detection, safety, high efficiency, and non-contact operation. It has become an important method in power inspection. However, the thermal image characteristics of faulty insulators are not evident and cannot be recognized directly by the naked eye. Therefore, in this study, we first conduct a temperature field simulation analysis of porcelain insulators in distribution lines and then propose an infrared thermal image detection method for faulty insulators. A multi-scale template matching algorithm is used to locate and identify the insulators. The coordinate parameters of the insulator in the infrared image are obtained, the insulator is segmented and extracted by multi-scale template matching, and the temperature of the insulator is extracted by least-square linear fitting. Combined with the relevant standards and simulation analysis results, the differences in the temperature states among multiple insulators were compared using a similar comparison judgment method to detect faulty insulators.
-
0. 引言
高光谱图像目标探测分为有先验知识的目标探测和无先验知识的异常检测,其中异常目标是指相对于场景中的其他目标地物存在光谱和空间属性差异且分布较少的目标。异常检测算法通过建立背景信号模型,找到图像中与背景模型不同的像元,将该像元归为异常目标[1]。
基于统计的异常检测算法通常假设背景属于某种特定的概率分布,而不遵循这种概率分布的样本即为异常,这类算法包括RX(Reed-Xiaoli)算法[2]、UTD-RXD算法[3]、局部RX(Local RX)算法[4-5]、核RX算法(Kernel RX)[6]以及局部核RX检测器(Local Kernel RX Detector)[6]等。而基于稀疏表达的异常检测算法通过背景的字典表达,将稀疏回归重建误差作为衡量异常的度量[7-8];另外稀疏低秩分解[9-10]将具有稀疏特性的目标和具有低秩特性的背景分离开来;基于协同表示[11]的异常检测认为背景中的每个像元都可以由其空间邻域近似表示,而异常则不能。
视觉注意机制(visual attention mechanism)是人类视觉所特有的大脑信号处理机制[12],人类视觉通过扫描全局图像,经过大脑处理,会自动识别图像中的感兴趣部分,即显著性区域。视觉注意机制模型分为自底向上与自顶向下两种形式,其中自底向上的视觉注意机制模型通过计算底层特征,无监督地对图像进行显著性检测。
Christof Koch和Laurent Itti[13]提出了对后续研究影响深远的ITTI视觉显著性模型,该模型设计来源是早期灵长类动物的视觉神经系统,其模型所用特征为亮度、颜色和方向。Liang等[14]将高光谱影像转换为三色图像进行显著性目标检测,得到的结果优于RGB影像,表明了视觉注意机制对于高光谱影像分析与处理的可行性。Cao等[15]用光谱角距离特征替代了颜色特征,证明了光谱特征在检测方面的优势。刘嘉诚等[16]利用显著性结合RX算法对高光谱影像进行异常检测,但该方法在进行显著性计算时仅利用了光谱信息。Zhao等[17]将频率调谐显著性检测模型与光谱显著性相结合,对高光谱图像中的目标进行检测。Xiang等[18]将视觉显著性模型引入到高光谱中,提出了一种通过视觉显著性模型和具有自适应权重的背景减法来实现高光谱异常检测的方法,但由于该方法将高光谱数据降维至三波段来替代RGB三波段进行显著性的计算,因此极易造成一些关键信息的缺失。
鉴于RX算法[2]计算背景的均值协方差时将异常目标也包含在内,导致背景模型估计不准确,从而降低检测精度,这里将视觉注意机制与RX融合,提出融合视觉注意机制的高光谱RX异常检测算法。该方法从采样方式、选取波段数以及局部光谱特征3方面对视觉注意机制进行改进,获取高光谱图像的显著性结果,然后利用显著性结果对原数据进行加权,并重新估计背景模型参数,进而对图像进行RX异常检测。
1. 基本原理
本文提出的算法首先以上-下采样的模式建立高斯金字塔,将原图像作为中间层,分别进行上采样与下采样,有效避免了因高光谱数据尺寸过小在下采样过程中而造成的信息损失,也有效防止了因高光谱数据量大的特点引起在上采样过程中造成信息量过大的问题。其次,在选取波段时选取多个波段,将每3个波段作为一个波段子集进行非光谱特征显著性的计算,有效避免了因数据降维造成的信息损失。最后,计算局部光谱特征显著性,将非光谱特征显著性与局部光谱特征显著性进行线性相加,得到最终的显著性结果,利用显著性结果对原数据进行加权,重新估计背景均值协方差,对原图像进行RX异常检测。算法流程图如图1所示。
1.1 RX算法
RX算法[2]是最早的异常检测算法,该算法可认为是广义似然比检测衍生出的一种具有恒虚警率的异常点探测器。定义图像X为M×N×L大小的图像,xi表示第i个像元光谱,RX算法的检测算子为:
$$ δ_\text{RXD}(x_{i})=(x_{i}-μ_\text{b})^{T}C_\text{b}^{-1}(x_{i}-μ_\text{b}) $$ (1) $$ {\mu _{\text{b}}} = \frac{1}{{{\text{Num}}}}\sum\limits_{i = 1}^N {{x_i}} $$ (2) $$ {C_{\text{b}}} = \frac{1}{{{\text{Num}}}}\sum\limits_{i = 1}^N {({x_i} - {\mu _{\text{b}}}){{({x_i} - {\mu _{\text{b}}})}^{\text{T}}}} $$ (3) 式中:Num为像元点个数;μb表示均值;Cb表示协方差。δRXD(xi)值越大,表明该像元越可能是异常。
1.2 融合自底向上视觉注意机制的高光谱图像异常检测
1.2.1 高斯金字塔
高斯金字塔利用高斯核函数对图像进行采样,得到图像不同尺度的子图集合。将原图像作为基础图像,利用高斯核对其进行卷积,对卷积后的图像进行上采样得到图像G,将此图像作为输入,重复卷积和采样操作得到下一步的图像,反复迭代多次,得到上采样高斯金字塔,其次对原图像进行卷积后下采样,迭代次数与上采样高斯金字塔相同,得到下采样高斯金字塔,两者结合形成最终的高斯金字塔,其中原图像为中间层。二维高斯卷积函数[19]为:
$$ G(x,y) = \frac{1}{{2{\text{π }}{\sigma ^2}}}{{\text{e}}^{ - \frac{{{x^2} + {y^2}}}{{2{\sigma ^2}}}}} $$ (4) 式中:σ为高斯函数的标准差,最终得到M层高斯金字塔,记为G(M)。
1.2.2 局部光谱特征
设局部窗口内任意像元为xp=[xp1, xp2, ⋯, xpL]T,其中p∈(1,n),n表示局部窗口内的像元个数。定义ds(xi, xj)表示局部窗口中任意一个像元的光谱xi和局部窗口中心像元的光谱xj之间的光谱相似性。为使局部窗口内的像元遵循越靠近中心像元影响越大的原则,定义dp(xi,xj)为像元点xi和中心像元xj在局部窗口中的空间距离,以dp(xi,xj)计算邻域像元对中心像元的影响,光谱显著性dSaliency(xi,xj)计算公式如下:
$$ {d_{{\text{Saliency}}}}({x_i},{x_j}){\text{ = }}\frac{{{d_\text{s}}({x_i},{x_j})}}{{1 + c*{d_\text{p}}({x_i},{x_j})}} $$ (5) 式中:ds(xi,xj)和dp(xi,xj)采用欧氏距离[20]进行计算;c是用来调整显著性中位置距离和光谱距离之间比例的尺度参数,本文中设置为1。
为防止噪音或其他因素造成计算结果不准确,本研究以局部窗口的显著性平均值作为像元最终的显著值,定义dSaliency_eva为局部窗口的中心像元点与邻域像元光谱显著性的平均值,计算式如下:
$$ {d_{{\text{Saliency}\_\text{eva}}}} = \frac{{\sum\limits_{i = 1}^N {{d_{{\text{Saliency}}}}} }}{{{\text{sum}} - 1}} $$ (6) 式中:sum为局部窗口中除中心像元以外,显著性值非零的像元总数。在每一层金字塔进行计算,所得的光谱显著性结果记为D(M)。
1.2.3 非光谱特征
非光谱特征包括亮度特征、颜色特征和方向特征。在计算非光谱特征时,需要用3个波段进行计算,本研究选取K个波段子集,每个波段子集包括3个波段,共选取B=3K个波段,对每一个波段子集分别进行3种非光谱特征的计算。本文通过基于最优聚类框架(optimal clustering framework, OCF)的波段选取方法进行波段的选取[21]。OCF可以在合理的约束下获得特定形式目标函数的最优聚类结果,对每个聚类结果进行贡献度计算,相加作为整个聚类结果的度量。
亮度特征I提取公式如下:
$$ I\left( {M,K} \right) = {{\left( {{b_1}\left( {M,K} \right) + {b_2}\left( {M,K} \right) + {b_3}\left( {M,K} \right)} \right)} \mathord{\left/ {\vphantom {{\left( {{b_1}\left( {M,K} \right) + {b_2}\left( {M,K} \right) + {b_3}\left( {M,K} \right)} \right)} 3}} \right. } 3} $$ (7) 颜色特征提取时需要提取红R、绿G、蓝B、黄Y四种,计算公式如下:
$$ \begin{array}{l} R\left( {M,K} \right) = {{{b_1}\left( {M,K} \right) - \left( {{b_2}\left( {M,K} \right) + {b_3}\left( {M,K} \right)} \right)} \mathord{\left/ {\vphantom {{{b_1}\left( {M,K} \right) - \left( {{b_2}\left( {M,K} \right) + {b_3}\left( {M,K} \right)} \right)} 2}} \right. } 2} \hfill \\ G\left( {M,K} \right) = {{{b_2}\left( {M,K} \right) - \left( {{b_1}\left( {M,K} \right) + {b_3}\left( {M,K} \right)} \right)} \mathord{\left/ {\vphantom {{{b_2}\left( {M,K} \right) - \left( {{b_1}\left( {M,K} \right) + {b_3}\left( {M,K} \right)} \right)} 2}} \right. } 2} \hfill \\ B\left( {M,K} \right) = {{{b_3}\left( {M,K} \right) - \left( {{b_1}\left( {M,K} \right) + {b_2}\left( {M,K} \right)} \right)} \mathord{\left/ {\vphantom {{{b_3}\left( {M,K} \right) - \left( {{b_1}\left( {M,K} \right) + {b_2}\left( {M,K} \right)} \right)} 2}} \right. } 2} \hfill \\ Y\left( {M,K} \right) = {{\left( {{b_1}\left( {M,K} \right) + {b_2}\left( {M,K} \right)} \right)} \mathord{\left/ {\vphantom {{\left( {{b_1}\left( {M,K} \right) + {b_2}\left( {M,K} \right)} \right)} 2}} \right. } 2} - \hfill \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\; {{\left| {{b_1}\left( {M,K} \right) + {b_2}\left( {M,K} \right)} \right|} \mathord{\left/ {\vphantom {{\left| {{b_1}\left( {M,K} \right) + {b_2}\left( {M,K} \right)} \right|} 2}} \right. } 2} - {b_3}\left( {M,K} \right) \hfill \\ \end{array} $$ (8) 方向特征利用Gabor滤波器求取,二维Gabor滤波器[22]核函数为:
$$g(x, y ; \lambda, \theta, \psi, \sigma, \gamma)=\mathrm{e}^{-\frac{x^2+\gamma^2 y^2}{2 \sigma^2}} \mathrm{e}^{\text{i}\left(2 \pi \frac{x^{\prime}}{\lambda}+\psi\right)}$$ (9) 式中:λ表示波长;θ表示方向;ψ表示相位偏移;γ表示空间纵横比。对每层金字塔影像进行计算,得到的结果记为O(θ,M, K),其中θ=[0°, 45°, 90°, 135°]。
1.2.4 显著性图
得到各特征高斯金字塔后,利用中心-周围差分方法计算对应特征图[13],计算公式为:
$$ \begin{array}{l} I\left( {c,s,K} \right) = \left| {I\left( {c,K} \right)\Theta I\left( {s,K} \right)} \right| \hfill \\ RG\left( {c,s,K} \right) = \left| {\left( {R\left( {c,K} \right) - G\left( {c,K} \right)} \right)\Theta \left( {G\left( {s,K} \right) - R\left( {s,K} \right)} \right)} \right| \hfill \\ BY\left( {c,s,K} \right) = \left| {\left( {B\left( {c,K} \right) - Y\left( {c,K} \right)} \right)\Theta \left( {Y\left( {s,K} \right) - B\left( {s,K} \right)} \right)} \right| \hfill \\ O\left( {c,s,\theta ,K} \right){\text{ = }}\left| {O\left( {c,\theta ,K} \right)\Theta O\left( {s,\theta ,K} \right)} \right| \hfill \\ D\left( {c,s} \right) = \left| {D\left( c \right)\Theta D\left( {\text{s}} \right)} \right| \hfill \\ \end{array} $$ (10) 式中:Θ表示将图像调整尺寸后,进行相同矩阵大小的相减运算;c代表Center;s代表Surround;I表示亮度特征图;RG与BY表示颜色特征图;O表示方向特征图;D表示局部光谱特征图。
由于高斯特征金字塔的尺度不同,需要利用插值将图像统一到相同的大小,才能对各个特征进行跨尺度相加。因此本文将各尺度特征图统一至原始数据大小,对不同特征进行跨尺度相加后得到显著性图集$ {\overline I _k},{\overline C _k},{\overline O _k},\overline D $,其中$ {\overline C _k} $表示颜色显著性图集,然后将每个波段子集的非光谱特征图相加,得到最终的非光谱特征显著性图$ \overline I ,\overline C ,\overline O $。
对每个特征显著性图进行归一化后进行线性相加得到最终的显著性图,计算公式为:
$$ S{\text{ = }}{{\left( {\overline I + \overline C + \overline O + \overline D } \right)} \mathord{\left/ {\vphantom {{\left( {\overline I + \overline C + \overline O + \overline D } \right)} 4}} \right. } 4} $$ (11) 1.2.5 融合视觉注意机制的高光谱RX异常检测算法
对显著性图进行高斯加权[23]得到显著性权重,定义w(i,j)为权重;s(i,j)为像元显著值,公式为:
$$ w(i,j) = {{\text{e}}^{ - \frac{{{{\left( {s(i,j) - \mu } \right)}^2}}}{{2{\sigma _s}^2}}}} $$ (12) 式中:μ为图像显著性的均值;σs为图像显著性的方差。得到显著性权重图像W后,对原数据进行显著性加权,像元显著值越大,权重越小,公式为:
$$\boldsymbol{X}_{\mathrm{s}}=\boldsymbol{W} \odot \boldsymbol{X}$$ (13) 式中:⊙表示逐波段加权运算。加权后的数据背景被增强,异常被抑制,我们利用加权数据Xs进行背景均值μs与协方差Cs的重新估计,进而对图像X进行RX异常检测。
$$ {\delta _{{\text{HITTIRX}}}}({x_i}) = {({x_i} - {\mu _{\text{s}}})^{\text{T}}}{\boldsymbol{C}_s}^{ - 1}({x_i} - {\mu _{\text{s}}}) $$ (14) 得到最终的异常检测结果。
2. 实验与结果分析
2.1 实验数据
本文选取模拟实验数据LCVF[24]与airport-beach-urban(ABU)异常检测数据集[25]中的数据进行试验。ABU数据集中选取4个数据,分别为:airport4、beach4、urban1、urban3。模拟实验数据为在AVIRIS传感器所拍摄高光谱影像中加入模拟异常目标,植入的模拟异常目标具体介绍如表1所示。
表 1 LCVF模拟异常目标Table 1. LCVF simulates anomalous targetsColumn Number of targets/cells Proportion of target Proportion of background 1 5/20 0.8 0.2 2 5/20 0.6 0.4 3 5/5 0.4 0.6 4 6/6 0.2 0.8 表 2 异常检测数据介绍Table 2. Anomaly detection data introductionData Sensor Resolution/m Size Targert type Number of targets/cells LCVF AVIRIS 20 120×160×190 Simulated anomalous targets 21/51 airport4 AVIRIS 3.4 100×100×205 Airport 3/60 beach4 ROSIS-03 1.3 150×150×102 Car and bare soil 7/68 urban1 AVIRIS 17.2 100×100×204 Building 9/67 urban3 AVIRIS 3.5 100×100×191 Ship 11/52 2.2 参数设置
本文需要设置的参数包括:高斯金字塔层数、波段选取数、局部光谱特征计算窗口大小。首先是高斯金字塔层数,固定局部光谱特征计算窗口大小为5×5,波段选取数量为3进行实验,获取金字塔层数分别为3、5、7、9层的实验结果。其次,固定局部光谱特征计算窗口大小为5×5,金字塔层数为5层进行试验,获取波段数为B={b=3k|1≤k≤10}的实验结果。最终是确定局部光谱特征计算窗口的大小,在计算局部光谱特征时,需要对图像进行扩充,本文通过围绕边界进行镜像反射来补边,固定金字塔层数为5层,波段选取数量为3进行实验,获取窗口大小分别为3×3,5×5,7×7的实验结果。
本文采用接收者操作特征(receiver operating characteristic, ROC)曲线与ROC曲线下面积(the area under the ROC curve, AUC)作为检测结果的精度评定指标[26]。实验结果如表3所示。由实验得出参数应设置为:金字塔5层,波段数15个,计算局部光谱特征的窗口大小为3×3。此外,在进行计算显著性图中的Center-Surround操作时,选取s={1,2},c={2,3}。
表 3 参数实验结果Table 3. Parameters experimental resultsParameter Size AUC Parameter Size AUC Number of bands 3 0.9759 Number of pyramid layers 3 0.9081 6 0.9621 5 0.9759 9 0.9681 7 0.9709 12 0.9777 9 0.9624 15 0.9779 18 0.9772 Spectral feature window size 3×3 0.9769 21 0.9765 5×5 0.9759 24 0.9746 7×7 0.975 27 0.9726 30 0.9725 2.3 消融实验
为验证3个方面的改进是否均有效,本文进行了消融实验。表4展示了消融实验的统计结果,对勾(√)表示改进该模块,结果表明,本文所做改进有效。模型G(进行3方面改进)算法精度最高。相较于模型A(未进行改进)来说,仅改进采样方式(模型B)时精度提升0.68%;改进波段选取数量(模型C)时精度提高了1.69%;融入局部光谱特征(模型D)时提升了1.39%;在同时对采样方式和波段选取进行改进(模型E)后,精度提高2.30%;在改进采样方式的同时融入局部光谱特征(模型F)后,精度提升了2.14%;同时改进3个模块时,精度提高了2.46%。由于局部光谱特征是利用全波段进行计算,波段选取数量影响的是非光谱特征,故不进行同时改进波段选取数量与融入局部光谱特征的实验。
表 4 消融实验Table 4. Ablation experimentmodel Sampling mode Band selection Spectral feature AUC A 0.9625 B √ 0.9690 C √ 0.9788 D √ 0.9759 E √ √ 0.9846 F √ √ 0.9831 G √ √ √ 0.9862 2.4 结果与分析
实验将展示ITTI显著性图(ITTI saliency map, ITTIS)与高光谱ITTI显著性图(Hyperspectral ITTI saliency map, HITTIS)2种显著性图,以及RX算法、ITTI显著性模型RX算法(ITTI saliency model RX algorithm, ITTIRX)、高光谱ITTI显著性模型RX算法(Hyperspectral ITTI saliency model RX algorithm, HITTIRX)3种算法的检测结果,并对结果进行对比分析。
如图3(a)所示,在模拟数据LCVF的检测结果中,RX检测结果仅检测出第一列与第三列目标,在引入显著性后进行RX检测,虽并未检测出更多目标,但背景抑制效果较好,在引入HITTIS后,目标被全部检测,背景几乎被全部抑制,仅有右下角的一部分未被完全抑制。
在真实数据中,HITTIRX的检测结果也同样展现了较好的效果。如图3(b)所示,在airport4的检测结果中,RX算法检测出了3个飞机中较大的一架,在引入ITTIS后,背景较为纯净,但检测目标的效果并不好,HITTIRX的检测结果中,3架飞机被成功检测出,背景抑制较RX与ITTIRX也更好;如图3(c)所示,在beach4的检测结果中,RX算法检测出目标,但并未抑制桥梁、右上与下方背景,在引入ITTIS后,桥梁被抑制,在HITTIRX检测结果中,桥梁与右上背景被完全抑制,下方背景被大部分抑制;如图3(d)所示,在urban1的检测结果中,RX算法虽检测到了目标,但虚警率较高,在ITTIRX的检测结果中,虚警率有所下降,但图像中仍有一个面状背景未被抑制,在引入HITTIS后,虚警率最低,目标被较好地检测出来,ITTIRX中未被抑制的面状背景也被抑制;如图3(e)所示,在urban3的检测结果中,RX算法将目标全部检测,但同其他数据检测结果一样存在虚警率高的问题,ITTIRX检测中虽降低了虚警率,但目标检测效果不好,引入HITTIS后,目标全部被检测,背景抑制效果也最好。
综上所述,相较于RX算法,引入ITTIS后,背景抑制效果增强,在HITTIRX的检测结果中,可以明显看到相较于ITTIRX来说,背景信息抑制更彻底,检测的目标更全面准确。总而言之,与RX和ITTIRX相比,HITTIRX虚警率更低,检测准确度更高。
如图4所示,可以明显看出本文建议方法的ROC曲线相较于RX算法与ITTIRX算法来说更接近左上角,模拟数据中本算法无限接近于左上角,证明HITTIRX算法的可行性,真实数据的结果也表明其性能优于RX算法与ITTIRX算法。
表5统计了各算法检测结果的AUC值。从表中可以看出,在5个数据上,HITTIRX算法精度最高且均在0.98以上,直接引入视觉注意机制会降低部分数据的检测准确性,但在进行改进后,检测准确度上升,虚警率下降。HITTIRX相较于RX检测精度来说分别提高了19.99%、3.86%、2.70%、0.77%、3.67%,相较于BUS-RX来说分别提高了14.24%、12.84%、5.86%、1.03%、2.46%。
表 5 AUC统计结果Table 5. AUC statisticsRX ITTIRX HITTIRX LCVF 0.8333 0.8753 0.9999 Airport4 0.9526 0.8768 0.9894 Beach4 0.9588 0.9302 0.9847 Urban1 0.9907 0.9881 0.9983 Urban3 0.9513 0.9625 0.9862 3. 结论
不同于以往的异常检测方法,本文引入了视觉注意机制,将视觉注意机制模型进行改进,并将高光谱图像的局部光谱特征融入视觉注意机制进行显著性图的获取,解决了传统视觉注意机制模型不适用于高光谱图像的问题。对原图像加权后重新估计背景参数进行异常检测,提高了RX异常检测算法背景模型的精确度。通过实验与分析得出以下结论:①引入视觉注意机制对高光谱图像异常检测具有可行性;②改进金字塔采样方式与波段选取数量对提高基于视觉注意机制的高光谱图像异常目标检测精度是有效的;③在视觉注意机制中加入局部光谱特征能够提高高光谱图像异常检测的精度,具有较低虚警率。
-
表 1 材料属性
Table 1 Material parameters
Property Heat conductivity
/[W(m2·K)-1]Specific heat capacity
/[J(kg·℃)-1]Cap 1.22 816 Cement 2 80 Pin 4.5 465 Wire Aluminum 170 881 Steel 43 470 Insulation 0.4 800-1200 表 2 仿真结果
Table 2 The simulation results
Maximum temperature Global Insulator Normal insulator 33.4℃ 31.3℃ Low resistance insulator 33.4℃ 32.4℃ Zero resistance insulator 31.2℃ 29.7℃ 表 3 灰度值以及与其对应的温度
Table 3 Gray value and its corresponding temperature
Gray Value (x) 27 39 51 60 66 75 81 90 96 105 112 120 Temperature (y)/℃ 1 2 3 4 5 6 7 8 9 10 11 12 Gray Value (x) 132 141 150 159 165 177 186 195 207 219 225 235 Temperature (y)/℃ 13 14 15 16 17 18 19 20 21 22 23 24 表 4 P-10型绝缘子红外热像的温度提取结果
Table 4 Temperature extraction results of P-10 insulator infrared thermal image
Sequence number Original images Frequency maximum gray value Actual temperature /℃ 1 a.jpg 215 22.1556 2 b.jpg 216 22.2688 3 c.jpg 214 22.0424 4 d.jpg 215 22.1556 5 e.jpg 221 22.8347 6 f.jpg 206 21.1368 7 g.jpg 216 22.2688 8 h.jpg 221 22.8347 表 5 PS-15型绝缘子红外热像的温度提取结果
Table 5 Temperature extraction results of PS-15 insulator infrared thermal image
Sequence number Original images Frequency maximum gray value Actual temperature /℃ 1 i.jpg 217 22.3820 2 j.jpg 224 23.1744 3 k.jpg 215 22.1556 4 l.jpg 216 22.2688 -
[1] 王羽, 文习山, 蓝磊, 等. 提高架空配电线路耐雷水平的仿真分析[J]. 高电压技术, 2011, 37(10): 2471-2476. https://www.cnki.com.cn/Article/CJFDTOTAL-GDYJ201110022.htm WANG Yu, WEN Xishan, LAN Lei, et al. Simulation analysis on improving lightning withstanding level of overhead distribution lines[J]. High Voltage Engineering, 2011, 37(10): 2471-2476. https://www.cnki.com.cn/Article/CJFDTOTAL-GDYJ201110022.htm
[2] 邱志斌, 阮江军, 黄道春, 等. 输电线路悬式瓷绝缘子老化形式分析与试验研究[J]. 高电压技术, 2016, 42(4): 1259-1267. https://www.cnki.com.cn/Article/CJFDTOTAL-GDYJ201604027.htm QIU Zhibin, RUAN Jiangjun, HUANG Daochun, et al. Study on aging modes and test of transmission line porcelain suspension insulators[J]. High Voltage Engineering, 2016, 42(4): 1259-1267. https://www.cnki.com.cn/Article/CJFDTOTAL-GDYJ201604027.htm
[3] 汪小丁. 10 kV配电线路防雷分析与海宁地区防雷讨论[J]. 农村电气化, 2022(7): 24-26. https://www.cnki.com.cn/Article/CJFDTOTAL-NCDH202311002.htm WANG Xiaoding. Analysis and discussion on lightning protection of 10 kV distribution line[J]. Rural Electrification, 2022, (7): 24-26. https://www.cnki.com.cn/Article/CJFDTOTAL-NCDH202311002.htm
[4] 国家能源局. 输变电设备状态检修试验规程: DL/T 393-2010 [S]. 北京: 中国电力出版社, 2010. National Energy Administration. Regulations of Condition-based Maintenance & Test for Electric Equipment: DL/T 393-2010[S]. Beijing: China Electric Power Press, 2010.
[5] 国家能源局. 劣化悬式绝缘子检测规程: DL/T 626-2015 [S]. 北京: 中国电力出版社, 2015. National Energy Administration. Aging suspension insulators inspection rule: DL/T 626-2015[S]. Beijing: China Electric Power Press, 2015.
[6] 王力农, 简思亮, 宋斌, 等. 基于电场分布测量法的输电线路劣化绝缘子检测研究[J]. 电瓷避雷器, 2019(4): 199-205. https://www.cnki.com.cn/Article/CJFDTOTAL-DCPQ201904034.htm WANG Linong, JIAN Siliang, SONG Bin, et al. Research on the detection of faulty insulator in transmission line based on electric field distribution measurement[J]. Insulators and Surge Arresters, 2019(4): 199-205. https://www.cnki.com.cn/Article/CJFDTOTAL-DCPQ201904034.htm
[7] 卢航, 姚建刚, 付鹏. 基于总体最小二乘的Shearlet自适应零值绝缘子红外图像去噪[J]. 红外技术, 2015, 37(10): 842-846. http://hwjs.nvir.cn/article/id/hwjs201510007 LU Hang, YAO Jiangang, FU Peng, et al. Faulty insulators infrared thermal image adaptive denoising based on total least squares estimation and Shearlet transform[J]. Infrared Technology, 2015, 37(10): 842-846. http://hwjs.nvir.cn/article/id/hwjs201510007
[8] 程洋, 夏令志, 李志飞, 等. 基于红外成像法的零值瓷绝缘子检测[J]. 绝缘材料, 2019, 52(3): 74-79. https://www.cnki.com.cn/Article/CJFDTOTAL-JYCT201903014.htm CHENG Yang, XIA Lingzhi, LI Zhifei, et al. Detection of faulty porcelain insulator based on infrared imaging method[J]. Insulating Materials, 2019, 52(3): 74-79. https://www.cnki.com.cn/Article/CJFDTOTAL-JYCT201903014.htm
[9] 张美金, 屈秋帛. 基于GWO-SVM的红外热成像低零值绝缘子识别[J]. 红外技术, 2021, 43(4): 397-402. http://hwjs.nvir.cn/article/id/2e60cc20-25fd-4344-8e38-ff3768e03fa2 ZHANG Meijin, QU Qiubo. Infrared thermography low-zero insulator identification based on GWO-SVM[J]. Infrared Technology, 2021, 43(4): 397-402. http://hwjs.nvir.cn/article/id/2e60cc20-25fd-4344-8e38-ff3768e03fa2
[10] 王胜辉, 姜婷玥, 李伟, 等. 零值绝缘子串红外和紫外成像特性试验及耦合场仿真研究[J]. 高压电器, 2021, 57(12): 1-9. https://www.cnki.com.cn/Article/CJFDTOTAL-GYDQ202112001.htm WANG Shenghui, JIANG Tingyue, LI Wei, et al. Study on infrared and UV imaging characteristics test and coupling field simulation of zero-value insulator strings[J]. High Voltage Apparatus, 2021, 57(12): 1-9. https://www.cnki.com.cn/Article/CJFDTOTAL-GYDQ202112001.htm
[11] 王海涛, 郑雷, 李健, 等. 500 kV瓷绝缘子串的分体式低/零值检测机器人优化设计与实测[J]. 高电压技术, 2020, 46(9): 3035-3043. https://www.cnki.com.cn/Article/CJFDTOTAL-GDYJ202302015.htm WANG Haitao, ZHENG Lei, LI Jian, et al. Optimization design and measurement test of a split low/zero insulator detection robot for 500 kV porcelain insulator strings[J]. High Voltage Engineering, 2020, 46(9): 3035-3043. https://www.cnki.com.cn/Article/CJFDTOTAL-GDYJ202302015.htm
[12] 彭向阳, 梁福逊, 钱金菊, 等. 基于机载红外影像纹理特征的输电线路绝缘子自动定位[J]. 高电压技术, 2019, 45(3): 922-928. PENG Xiangyang, LIANG Fuxun, QIAN Jinju, et al. Automatic recognition of insulator from UAV infrared image based on periodic textural feature[J]. High Voltage Engineering, 2019, 45(3): 922-928.
[13] 张也, 彭子健, 付强, 等. 环境湿度对瓷质绝缘子串电压分布及红外热像检测的影响分析[J]. 电网技术, 2018, 42(4): 1342-1349. https://www.cnki.com.cn/Article/CJFDTOTAL-DWJS201804043.htm ZHANG Ye, PENG Zijian, FU Qiang, et al. Analysis of environment humidity influence on voltage distribution and infrared thermal image detection of porcelain insulator strings[J]. Power System Technology, 2018, 42(4): 1342-1349. https://www.cnki.com.cn/Article/CJFDTOTAL-DWJS201804043.htm
[14] 张晓春, 欧阳广泽, 何洪英, 等. 基于红外图像匹配的零值绝缘子检测[J]. 电测与仪表, 2019, 56(6): 100-105. https://www.cnki.com.cn/Article/CJFDTOTAL-DCYQ201906018.htm ZHANG Xiaochun, OUYANG Guangze, HE Hongying, et al. Zero -insulator detection based on infrared images matching[J]. Electrical Measurement & Instrumentation, 2019, 56(6): 100-105. https://www.cnki.com.cn/Article/CJFDTOTAL-DCYQ201906018.htm
[15] 李唐兵, 陈国锋, 饶斌斌, 等. 一种改进的劣化盘形悬式瓷绝缘子红外热像诊断方法[J]. 电瓷避雷器, 2014(6): 8-13. https://www.cnki.com.cn/Article/CJFDTOTAL-DCPQ201406002.htm LI Tangbing, CHEN Guofeng, RAO Binbin. A improved method of detecting deterioration insulator based on infrared thermal image[J]. Insulators and Surge Arresters, 2014(6): 8-13. https://www.cnki.com.cn/Article/CJFDTOTAL-DCPQ201406002.htm
[16] 国家能源局. 带电设备红外诊断应用规范: DL/T 664-2016 [S]. 北京: 中国电力出版社, 2016. National Energy Administration. Application rules of infrared diagnosis for live electric equipment: DL/T 664-2016 [S]. Beijing: China Electric Power Press, 2016.
[17] 刘鹏. 基于红外和紫外成像法的零值绝缘子发热与放电特性研究[D]. 北京: 华北电力大学, 2016. LIU Peng. Study on Zero Value Insulator Heating and Discharge Characteristic Based on the Infrared and Ultraviolet Imaging Method [D]. Beijing: North China Electric Power University, 2016.
[18] 莫运德, 马谦. 陶瓷材料导热率及其影响因素的探讨[J]. 中国陶瓷, 1983(3): 14-17. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGTC198304004.htm MO Yunde, MA Qian. Discussion on thermal conductivity of ceramic materials and its influencing factors[J]. China Ceramics, 1983(3): 14-17. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGTC198304004.htm
[19] 张俊双, 郑重, 钱文晓, 等. 计及多因素影响的直流电弧对瓷绝缘子热效应仿真分析[J]. 电瓷避雷器, 2016(2): 18-23. https://www.cnki.com.cn/Article/CJFDTOTAL-DCPQ201602004.htm ZHANG Junshuang, ZHENG Zhong, QIAN Wenxiao, et al. The simulation analysis of DC arc thermal effect on porcelain insulators considering multiple factors[J]. Insulators and Surge Arresters, 2016(2): 18-23. https://www.cnki.com.cn/Article/CJFDTOTAL-DCPQ201602004.htm
[20] 徐鑫, 方春华, 智李, 等. 连续激光作用下瓷质绝缘子温度和热应力分析[J]. 光电子激光, 2021, 32(1): 78-87. https://www.cnki.com.cn/Article/CJFDTOTAL-GDZJ202101012.htm XU Xin, FANG Chunhua, ZHI Li, et al. Analysis of temperature and thermal stress of porcelain insulator under continuous laser irradiation[J]. Journal of Optoelectronics·Laser, 2021, 32(1): 78-87. https://www.cnki.com.cn/Article/CJFDTOTAL-GDZJ202101012.htm
[21] 王畇浩. 基于可见光与红外航拍图像的输电线路绝缘子多故障检测研究[D]. 成都: 电子科技大学, 2019. WANG Yunhao. Research on Multi-fault Detection of Transmission Line Insulators Based on Visible and Infrared Aerial Images[D]. Chengdu: University of Electronic Science and Technology of China, 2019.
[22] HUANG S, TANG C, XU M, et al. BM3D-based total variation algorithm for speckle removal with structure-preserving in OCT ima ges[J]. Applied Optics, 2019, 58(23): 6233-6243.
[23] 刘飞飞, 马礼然. 快速模板匹配算法在口罩耳绳检测中的应用[J]. 传感器与微系统, 2022, 41(1): 157-160. https://www.cnki.com.cn/Article/CJFDTOTAL-CGQJ202201041.htm LIU Feifei, MA Liran. Application of fast template matching algorithm in mask ear band detection[J]. Transducer and Microsystem Technologies, 2022, 41(1): 157-160. https://www.cnki.com.cn/Article/CJFDTOTAL-CGQJ202201041.htm
[24] 卞艳, 宫雨生, 马国鹏, 等. 基于无人机遥感影像的水体提取方法[J]. 浙江大学学报: 工学版, 2022, 56(4): 764-774. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC202204016.htm BIAN Yan, GONG Yusheng, MA Guopeng. Water extraction from unmanned aerial vehicle remote sensing images[J]. Journal of Zhejiang University : Engineering Science, 2022, 56(4): 764-774. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC202204016.htm
[25] 张俊, 关胜晓. 基于改进的最大后验概率矢量量化和最小二乘支持向量机集成算法[J]. 计算机应用, 2015, 35(7): 2101-2104. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJY201507060.htm ZHANG Jun, GUAN Shengxiao. Integration algorithm of improved maximum a posteriori probability vector quantization and least squares support vector machine[J]. Journal of Computer Applications, 2015, 35(7): 2101-2104. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJY201507060.htm
[26] Otsu N. A threshold selection method from gray-level histograms[J]. IEEE Transactions on Systems on Man and Cybernetics, 1979, 9(1): 62-66.
[27] 聂茜茜, 肖斌, 毕秀丽, 等. 基于超像素级卷积神经网络的多聚焦图像融合算法[J]. 电子与信息学报, 2021, 43(4): 965-973. https://www.cnki.com.cn/Article/CJFDTOTAL-DZYX202104011.htm NIE Qianqian, XIAO Bin, BI Xiuli, et al. Multi-focus image fusion algorithm based on super pixel level convolutional neural network[J]. Journal of Electronics & Information Technology, 2021, 43(4): 965-973. https://www.cnki.com.cn/Article/CJFDTOTAL-DZYX202104011.htm
-
期刊类型引用(1)
1. 贾本康. 基于变压器噪声及故障分析的变电站轮式巡检机器人研究. 电力设备管理. 2024(17): 166-168 . 百度学术
其他类型引用(2)