姚建刚等

摘要：为解决目前绝缘子低零值检测方法漏判率高、操作繁琐的问题，提出了一种新的绝缘子串红外图像中绝缘子盘面和铁帽区域自动提取方法和状态识别模型.首先将现场拍摄的绝缘子红外图像进行灰度化处理、去噪、二值化；然后从二值图像中提取反映绝缘子的特征点集合；通过特征点对二值图进行角度校正；最后通过区域提取中的特定算法提取出绝缘子的盘面和铁帽区域.通过提取该区域内的绝对温度、纹理和相对温差率作为绝缘子状态识别的特征集.将用电压分布法测得的绝缘子状态信息作为输出向量，通过训练得到优化的识别模型，用于绝缘子状态识别.该方法经过了220 kV试验验证，证实了模型的有效性和实用性.

关键词：绝缘子；红外图像；二值形态学；Hough变换；自动提取；铁帽；纹理特性

中图分类号：TP301.6；TM85 文献标识码：A

绝缘子被广泛应用于输电线路中，是输电网络的重要组成部分.因此绝缘子的状态检测将是一个繁重的工作.相较于传统的人工登杆巡线检测方式，利用红外热像仪进行巡线拍摄[1]，并使用图像处理技术自动检测故障的巡线方式更为高效[2]，是目前智能巡线检查的主要发展方向.

同串绝缘子污湿状态相似，环境因素也一致，因此各绝缘子承载电压基本一致、发热情况也基本相同，同串才有可比性.正常绝缘子串以铁帽为主要发热中心，其温度分布和电压分布规律相似，也呈现不对称的马鞍型，但相邻绝缘子温差很小.低值绝缘子发热比正常绝缘子温度高，零值绝缘子发热比正常绝缘子温度要低，而表面污秽绝缘子其热像特征以瓷盘面为发热区，表现为瓷盘温升偏高而铁帽温升正常[3].因此，获取绝缘子盘面和铁帽区域的信息至关重要，而红外图像中绝缘子的自动识别提取则是其中的基础步骤.

目前有基于航拍图像的绝缘子提取方法，但航拍图像和红外热像图像有本质的不同，航拍图像保留了绝缘子本来的外貌特性和性质，如颜色.而红外热像图像是基于温度的信息.对于绝缘子的零低值检测温度才是最关键的信息.

1模型的构建

劣化绝缘子的红外识别主要由图像处理、特征集提取以及模式识别这三部分组成.

图像处理是在兼顾运行速度的前提下，采用合适的算法将绝缘子串有效地分割出来；特征参数的选取必须能反映出串中正常绝缘子与劣化绝缘子之间的差异；模式识别主要采用人工智能算法，通过挖掘特征参数所含信息来判断是否含有劣化绝缘子.图1为识别模型流程.

2绝缘子自动提取的基本原理

瓷质绝缘子的红外图像具有以下4个特征[7]：

1）瓷质绝缘子为轴对称结构，由铁帽、钢脚和瓷盘面3部分组成；

2）红外图像中的单片绝缘子盘面常呈现为椭圆形状，铁帽呈等腰梯形；

3）绝缘子一般成串出现，其数量随着输电线路的电压等级而不同，一般110 kV输电线路为7片，220 kV输电线路为14片；

4）同串绝缘子一般由同一型号、物理外形相同的绝缘子组成，在红外图像中表现为尺寸基本一致，且等间距排列.

基于绝缘子串的上述特征，本文提出绝缘子红外图像目标区域自动提取的方法，其流程图如图1中右半图所示，实现了自动提取绝缘子盘面和铁帽区域的功能，整个提取过程包括红外图像预处理、特征点提取、角度校正和区域提取4个部分.

2.1灰度化和去噪

在YUV颜色模型中，Y分量表示的是亮度；将RGB颜色转换成YUV颜色，只取其中的Y分量，即可表示为灰度图像.如下式[8]：

[Y U V]=[R G B]0.299-0.1480.6150.587-0.289-0.5150.1140.437-0.100

f=Y=0.299R+0.587G+0.114B

红外热像图像具有大噪声、低对比度的特点，普通的滤波方法虽能滤去部分噪声，但同时也模糊了图像本身，丢失了图像的部分边缘和细节.因此，本文采用自适应平滑滤波算法对其滤波去噪[9].

2.2OTSU分割

OTSU算法也称最大类间差法[10]，是图像分割中阈值选取的最佳算法，计算简单，不受图像亮度和对比度的影响.本文将采用该算法确定阈值T，然后对其进行二值化处理：

2.3二值形态学

数学形态学中二值图像的形态变换是一种针对集合的处理过程.将二值图像看成是集合，并用结构元素去度量和提取图像中的对应形状，去除不相干的结构，以达到图像分析和识别的目的.它的最基本运算有：腐蚀和膨胀[11].

骨架，可以理解为图像的中轴.集合A的骨架化：反复移除图像A的边界像素，但不允许原本连接的目标图像断裂，且操作保持欧拉数不变.所谓细化，就是从原来的图中去掉一些边界点，但仍要保持原来的形状.

先膨胀后腐蚀的过程称为闭运算，它具有填充细小空洞、连接邻近物体和平滑边界的作用.先腐蚀后膨胀的过程称为开运算，它具有消除细小物体、在纤细处分离物体和平滑较大物体边界的作用，如消除导线干扰[12].

在不同的应用场合，结构元素的选择及其相应的处理算法是不一样的，结构元素的大小、形状选择合适与否，将直接影响图像的形态运算结果.为了保证图像处理的最佳效果，本文中不同地方的膨胀或腐蚀所用到的结构元素都有所不同.

2.4绝缘子图像的倾斜角度校正

当红外图像中绝缘子倾斜时，利用绝缘子图像的几何特征，对绝缘子倾斜角度进行估计校正.校正的基本思路为：1）对二值化图像细化，获取绝缘子骨架；2）从骨架图中计算出交点，作为Hough变换的特征点集S；3）采用Hough变换对S进行直线拟合，计算出最长的拟合直线L和其倾斜角θ，即用该倾斜角对绝缘子图像进行倾斜校正[13]，能将倾斜校正误差控制在1°左右.

最简单的Hough变换是在图像中检测直线.在平面直角坐标系（x-y）中，用方程y=kx+b表示一条直线.该直线上任意一点（x，y）变换到（k-b）参数空间将只变成一个“点”，即点（k，b）.（k-b）参数空间中的一个局部峰值点就很有可能对应着原图像空间中的一条直线.对图像上所有的点进行Hough变换，最终所要检测的直线对应的一定是参数平面中出现频率最多的那个点.这样就在图像中检测出了直线.由于直线的斜率可能为无穷小，或者无穷大，则在（k-b）参数空间不便于对直线进行刻画和描述.所以，文献[14-15]提出了采用极坐标参数空间（ρ-θ）进行直线检测：

ρ=xcosθ+ysinθ

在Hough变换中，检测倾斜角的基本思路是由图像空间中的特征数据点去计算参数空间中的参数点的可能轨迹，并在一个累加器J（ρ，θ）中统计参考点，检测Hough域中曲线最频繁的交点（即J（ρ，θ）的最大值），如图2（b）所示，该交点对应着直角坐标系中最长直线，求得θ=2.04°.

2.5盘面和铁帽区域提取算法

2.5.1投影统计法

根据灰度图像的投影统计计算过程，对于绝缘子和导线像素点所构成的二值图像f′（x，y），其中m为图像f′（x，y）的高度，n为f′（x，y）的宽度，则垂直方向的投影信号表示为：

2.5.2区域定位提取算法

输入：倾斜校正后的二值图像f′.输出：盘面和铁帽区域的二值图像.算法步骤：

1）对图像f′的起始列进行逐行扫描，像素点间距d定义为一列0，1的数组中连续1的最大长度，记录在二维数组D中，终点（i，j）的长度记为D（i，j）=d.统计d出现的频数并记录在数组P中，即进行P（d）++操作；

2）对图像的下一列，同样按1）的方法进行操作，直至遍历整个图像；

3）滤去导线的干扰，导线的像素间距一般较小，且出现频数较高，设定阈值φ，令P（i≤φ）=0；

4）去除边缘干扰，规定在间距值10%以内的都认为是等间距，则P1（d）=∑1.1di=0.9dP（i）；

5）求出P1中频数最大的两极大值所对应的i，j，有P1（i）=max1，P2（j）=max2，（i

6）根据i，j分别反向计算求出铁帽和盘面的二值图像；

7）形态学处理，角度恢复；

8）结束.

关于上述算法的一些说明：对于步骤1），如有二值图像矩阵：

步骤3）中，对P中的结果再过滤，以剔除某些符合特征的电线等其他干扰.因为绝缘子的盘面和铁帽有一定的宽度，而其频率通常是整幅图像中最高的，根据这一特点将绝缘子盘面和铁帽区域提取出来.实际间距又不完全是个定值，所以规定间距值的10%以内都认为是等间距的.

在无步骤4）的情况下极值有时不明显，经过该操作处理后极值较容易得到，如图2（d）所示，求得两极大值在i=20和j=46处.

形态学处理能达到将不相关的细小区域滤除，而将相关的邻近区域连接并平滑边界的作用.

3状态识别

绝缘子一旦出现某种内部或外部故障，则故障发热通过热传导或其他形式热交换，改变绝缘子相应表面部位的温升或温度分布，从而表现出其红外图谱的差异.本文模型主要针对实用性设计，应尽量排除变量因素的干扰，对于某一变电站，某一时间段内（一般1～2 h内可以测试完），其外界环境对其绝缘子的影响基本不变（如环境温度，湿度，日照等）；本文模型具体到同串的每个绝缘子，对于同串来说，其绝缘子的外界环境也是基本一致；且不同变电站，不同环境变量因素太多，可比性不强，反而增大模型识别误差.为抽取故障特征集，分析异常绝缘子和正常绝缘子在温度和灰度特征上的差异，取热像图中绝缘子铁帽区域的平均温度作为该绝缘子的温度，并定义以下3个参数：绝对温度T、纹理H和相对温差率δt.

3.1绝对温度判断法

绝缘子表面绝对温度判断法是根据检测得到绝缘子表面绝对温度值，结合串中不同位置绝缘子的温度和温升极限的有关规定，分析判断绝缘子温度过热部位（铁帽和盘面）状态的正常与否[3].

这种方法的优点是分析简单直观，但由于受环境、距离和辐射率等方面的影响，绝对温度值不十分准确，但也能在一定程度上反映绝缘子的状态，其结果可对其他方法起一定的参考作用.

3.2图像纹理特征判断法

根据绝缘子串的红外热像图谱来判断设备是否正常.将红外图像灰度化处理后，利用其纹理特征来表征绝缘子之间的差异，图像纹理表现为二维空间的灰度变化模式，是一种区域特征，反映图中各像素之间空间分布的特性[16]，如方差H1、熵H2和平均梯度H3.

3.3相对温差率判断法

绝缘子在串中的位置不同所承载的电压也不同，环境温度也会影响红外诊断的结果，当环境温度低，尤其是承载电压小时，设备的温度值虽未超过相关标准，但在电压增长或环境温度上升后，会引发设备故障.可使用相对温差法解决上述问题.该方法是指同串中的某绝缘子与环境温度的温差（即温升）与相邻绝缘子的温升之比的百分数，相对温差率的计算如式：

3.4辅助判别法

尽管红外诊断是一种先进的诊断方法，上述各种判别方法也可以从不同的角度提高诊断的准确性，但是有些绝缘子的内部故障用红外成像很难诊断，因此需要一种可靠性高的检测方法（如电压分布法等）来对上述方法的准确性进行测试.由于该方法需要登杆作业，费时费力，只能作为一种对上述检测方法进行实时校正的辅助方法，从而将校验后的方便、准确模型应用于现场绝缘子的检测中.

4结果分析

上述分割方法抗干扰性强，为了验证本文方法的有效性，具体分析了一个试验算例，提取结果如

图2所示.并另外随机选取了江西省某地区2个220 kV变电

站（白沙站和珠珊站）的各50张现场拍摄的复杂背景的红外绝缘子图像和试验图像50张，利用本文的提取方法批量进行绝缘子盘面和铁帽区域的提取，其正确提取率都在90%以上.表1为本文方法从试验算例中提取出的铁帽区域特征集和识别结果.

根据随机选取的30张现场红外图像和10张试验红外图像分析可得以下结论：

1）良好绝缘子串的发热温度分布于其电压分布规律基本相同，呈不对称的马鞍形.当绝缘子串中含有低值绝缘子时，其绝缘子的承载电压将降低，且由于自身阻值的大大减小，低值绝缘子的发热远大于良好绝缘子的发热.如图3所示.