基于计算机视觉的变压器复合绝缘子套管监测方法及系统的制作方法

基于计算机视觉的变压器复合绝缘子套管监测方法及系统的制作方法
【专利摘要】一种基于计算机视觉的变压器复合绝缘子套管监测方法及系统，在各变压器上空龙门架处分别安装摄像头，拍摄所得图像传输到服务器端，在服务器端进行图像采集、图像预处理、图像区域分割、套管污秽分级处理及套管工况缺陷识别。本发明用可见光图像识别与处理技术提高变压器运行可靠性，及时发现缺陷，对套管污秽情况进行分级；而且对套管工况进行监测判断，提高套管运行可靠性。本发明可广泛用于各种电压等级的变压器套管监测。
【专利说明】
基于计算机视觉的变压器复合绝缘子套管监测方法及系统
技术领域
[0001] 本发明属于电气设备在线监测领域，特别设及一种基于计算机视觉的变压器复合 绝缘子套管在线监测方法及系统。
【背景技术】
[0002] 套管是变压器的最重要的附件之一。它不仅把变压器内部的高、低压引线引到油 箱之外，实现与外部网络连接，同时还起着固定引线的作用。所W，套管必须满足电气强度 和机械强度的双重要求。套管又是载流元件之一，在变压器运行中将长期通过负载电流。近 年，随着电网规模的扩大和环境污染的加剧，变压器套管的事故频发，严重威胁着电力系统 的安全运行。因此对变压器套管进行实时的在线监测，非常有必要。
[0003] 对变压器套管污秽的在线监测目前研究较少，通常选取其介质损耗(简称介损)正 切值tanS、电容量C、末屏泄漏电流I等作为绝缘状况判断，由于变电站电磁干扰较大，运些 监测都需要提取弱信号作为依据，准确度难W保证。
[0004] 运行中变压器套管的工况监测更是运行部口关注的热点，目前尚无运方面的研究 与成果出现。

【发明内容】

[0005] 为了兼顾在线监测变压器套管的运行工况和污秽等级，本发明提出了一种基于计 算机视觉的变压器复合绝缘子套管在线监测方法及系统。
[0006] 本发明技术方案提供一种基于计算机视觉的变压器复合绝缘子套管监测方法，在 各变压器上空龙口架处分别安装摄像头，拍摄所得图像传输到服务器端，在服务器端进行 图像采集、图像预处理、图像区域分割、套管污秽分级处理及套管工况缺陷识别；
[0007] 所述图像采集，包括采集拍摄所得图像，包括绝缘子的憎水性图像和变压器套管 整体的工况图像；
[000引所述图像预处理采用基于同态増晰的局部直方图均衡化方式实现；
[0009] 所述图像区域分割，包括对预处理后的图像，采用改进的Canny算子检测图像边 缘，然后对边缘检测后的图像采用数学形态学运算进行修正，得到分割结果；
[0010] 所述采用改进的化nny算子检测图像边缘实现方式如下，
[0011] 利用高斯滤波器对预处理后的绝缘子憎水性图像进行平滑滤波处理；
[0012] 用邻域一阶偏导的有限差分来计算平滑后图像的梯度幅值和方向，包括使用3X3 邻域来计算梯度，包括计算水平和垂直方向、的梯度幅值，还计算45°和135°方向的梯度幅 值；
[0013] 对梯度幅值图像进行处理，抑制梯度幅值的非极大值；
[0014] 用OTSU算法求出最佳分割阔值，然后用该最佳阔值进行检测并连接边缘；
[0015] 所述套管污秽分级处理，包括根据憎水性图像的分割结果，提取四个与复合绝缘 子憎水性有关的量作为特征码，包括最大水珠面积比K，水珠覆盖率C，最大水珠偏屯、率e，最 大水珠形状因子f。，采用基于最小二乘支持向量机的方式得到检测样本所属的憎水等级；
[0016] 最大水珠面积比K计算如下，
[0017]
[0018] 式中，S表示最大水珠的面积，XY为整幅图像的面积；
[0019] 水珠覆盖率C计算如下，
[0020]
[0021] 式中，n表示识别出的水珠个数，Si为第i个水珠的面积；
[0022] 忌十水珠偏屯、率e计算如下，
[0023]
[0024] 式中，a表示最大水珠的长轴，b表示最大水珠的短轴；
[0025] 县+古曲巧状因子fc计算如下，
[0026]
[0027] 式中，L表示最大水珠的周长，S表示最大水珠的面积；
[0028] 所述套管工况缺陷识别，包括将工况图像作为待检测图像，根据工况图像的分割 结果，将待检测图像边缘矩阵减去标准样本图像边缘矩阵并取绝对值，进行标准样本图像 和待检测图像的对比；使用bwareaopen函数去除污秽和斑点，将标准样本图像和待检测图 像的边缘图像最大灰度差的绝对值和最大赌法求取的阔值相组合成为判断是否存在裂痕 的阔值，若未超过阔值则判定为无裂痕，若是则提出告警。
[0029] 而且，针对憎水性图像，对边缘检测后的图像采用数学形态学运算进行修正的实 现方式如下，
[0030] 采用扁平的直线型结构元素进行膨胀运算，使整个图像的边缘变得连续;然后紧 接着对图像进行区域填充，使得图像边界内部的孔桐得到填补;最后，利用菱形结构元素对 区域填充后的图形进行腐蚀运算，消除边缘附近的毛刺，使得图像边缘变得光滑，同时去掉 过小而无意义的水滴。
[0031] 而且，针对工况图像，对边缘检测后的图像采用数学形态学运算进行修正的实现 方式如下，
[0032] 对待检测图像和标准样本图像经预处理和图像边缘检测后的结果，采用方形结构 元素对两幅图像分别进行开运算，并用圆盘形结构元素实行膨胀运算得到图像矩阵A和B; 同时，采用方形结构元素对两幅图像分别进行闭运算，再用圆盘形结构元素进行腐蚀得到 图像矩阵Al和BI;
[0033] 分别用A减去AUB减去BI,得到标准样本图像和待检测图像的边缘；
[0034] 使用最大赌法求取阔值得到标准二值边缘图像，对两图像边缘进行中值滤波后还 使用bwareaopen函数进行后处理，去除图像中的多余杂点。
[0035] 而且，将工况图像作为待检测图像时，标准样本图像采用相同的光照环境下的拍 摄结果。
[0036] 而且，根据多幅不同光照环境下待检测图像和标准样本图像逐一进行对比，如果 都超过每次对比后的阔值，则判断存在裂痕，进行报警。
[0037] -种基于计算机视觉的变压器复合绝缘子套管监测系统，在各变压器上空龙口架 处分别安装摄像头，拍摄所得图像传输到服务器端，在服务器端设置图像采集模块、图像预 处理模块、图像区域分割模块、套管污秽分级处理模块及套管工况缺陷识别模块；
[0038] 所述图像采集模块，用于采集拍摄所得图像，包括绝缘子的憎水性图像和变压器 套管整体的工况图像；
[0039] 所述图像预处理模块，用于采用基于同态増晰的局部直方图均衡化方式实现预处 理；
[0040] 所述图像区域分割模块，用于对预处理后的图像，采用改进的化nny算子检测图像 边缘，然后对边缘检测后的图像采用数学形态学运算进行修正，得到分割结果；
[0041] 所述采用改进的化nny算子检测图像边缘实现方式如下，
[0042] 利用高斯滤波器对预处理后的绝缘子憎水性图像进行平滑滤波处理；
[0043] 用邻域一阶偏导的有限差分来计算平滑后图像的梯度幅值和方向，包括使用3 X 3 邻域来计算梯度，包括计算水平和垂直方向、的梯度幅值，还计算45°和135°方向的梯度幅 值；
[0044] 对梯度幅值图像进行处理，抑制梯度幅值的非极大值；
[0045] 用OTSU算法求出最佳分割阔值，然后用该最佳阔值进行检测并连接边缘；
[0046] 所述套管污秽分级处理模块，用于根据憎水性图像的分割结果，提取四个与复合 绝缘子憎水性有关的量作为特征码，包括最大水珠面积比K，水珠覆盖率C，最大水珠偏屯、率 e，最大水珠形状因子f。，采用基于最小二乘支持向量机的方式得到检测样本所属的憎水等 级；
[0047] 最大水珠面积比K计算如下，
[004引
[0049]式中，S表示最大水珠的面积，XY为整幅图像的面积；
[(K)加]水珠覆盖率C计算如下，
[0化1 ]
[0052] 式中，n表示识别出的水珠个数，Si为第i个水珠的面积；
[0053] 最大水珠偏屯、率e计算如下，
[0化4;
[0055] 式中，a表示最大水珠的长轴，b表示最大水珠的短轴；
[0化6] 县地巧状因子fc计算如下，
[0化7]
[005引式中，L表示最大水珠的周长，S表示最大水珠的面积；
[0059]所述套管工况缺陷识别模块，用于将工况图像作为待检测图像，根据工况图像的 分割结果，将待检测图像边缘矩阵减去标准样本图像边缘矩阵并取绝对值，进行标准样本 图像和待检测图像的对比；使用bwareaopen函数去除污秽和斑点，将标准样本图像和待检 测图像的边缘图像最大灰度差的绝对值和最大赌法求取的阔值相组合成为判断是否存在 裂痕的阔值，若未超过阔值则判定为无裂痕，若是则提出告警。
[0060] 而且，针对憎水性图像，对边缘检测后的图像采用数学形态学运算进行修正的实 现方式如下，
[0061] 采用扁平的直线型结构元素进行膨胀运算，使整个图像的边缘变得连续;然后紧 接着对图像进行区域填充，使得图像边界内部的孔桐得到填补;最后，利用菱形结构元素对 区域填充后的图形进行腐蚀运算，消除边缘附近的毛刺，使得图像边缘变得光滑，同时去掉 过小而无意义的水滴。
[0062] 而且，针对工况图像，对边缘检测后的图像采用数学形态学运算进行修正的实现 方式如下，
[0063] 对待检测图像和标准样本图像经预处理和图像边缘检测后的结果，采用方形结构 元素对两幅图像分别进行开运算，并用圆盘形结构元素实行膨胀运算得到图像矩阵A和B; 同时，采用方形结构元素对两幅图像分别进行闭运算，再用圆盘形结构元素进行腐蚀得到 图像矩阵Al和BI;
[0064] 分别用A减去AUB减去BI,得到标准样本图像和待检测图像的边缘；
[0065] 使用最大赌法求取阔值得到标准二值边缘图像，对两图像边缘进行中值滤波后还 使用bwareaopen函数进行后处理，去除图像中的多余杂点。
[0066] 而且，将工况图像作为待检测图像时，标准样本图像采用相同的光照环境下的拍 摄结果。
[0067] 而且，根据多幅不同光照环境下待检测图像和标准样本图像逐一进行对比，如果 都超过每次对比后的阔值，则判断存在裂痕，进行报警。
[0068] 本发明采用计算机视觉，在高压变电站的主变的龙口架上架设多台高清摄像设 备，24小时不间断监控现场变压器套管的运行工况，采用图像识别方式，判断设备缺陷。同 时系统能准确捕捉复合绝缘子套管水滴图像，通过复合绝缘子套管的憎水机理，分析套管 的套管憎水特性，通过智能方式对套管的污秽等级进行判断，为电力系统安全运行提供保 障。本系统可广泛用于各种电压等级的变压器套管监测。
【附图说明】
[0069] 图1为本发明实施例提供的一种基于计算机视觉的变压器复合绝缘子套管在线监 测系统线路结构示意图；
[0070] 图2为本发明实施例提供的一种基于计算机视觉的变压器复合绝缘子套管在线监 测方法流程图。
[0071 ]图3为本发明实施例的具体流程原理图。
【具体实施方式】
[0072]为了使本发明的目的、技术方案及优点更加清楚明白，W下结合附图及实施例，对 本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用W解释本发明，并不 用于限定本发明。
[0073] 本发明的技术方案为一种基于计算机视觉的变压器复合绝缘子套管在线监测系 统。系统由多个高清摄像头，光纤图像传输网络及服务器组成。
[0074] 参见图1，在各变压器上空龙口架处分别安装摄像头，为保证准确率，建议采用高 清摄像头。所述的高清摄像头经过A、B处控制箱中的光电转换装置通过光通信网络将变压 器套管的实时运行图像拍摄后传输到控制室中的主控柜服务器，可W采用后台软件方式处 理后，给出变压器套管的运行工况和污秽等级判断。
[0075] 实施例中，高清摄像设备采用200万像素 W上全景自动变焦进口摄像头，通过室外 双屏蔽双绞线与千兆光收发器连接转换，然后通过12忍单模光缆通过TCP/IP协议将实时图 像传输到变电站控制室中的服务器。
[0076] 具体实施时，建议摄像设备的安装位置既考虑套管的全方位工况监测，也考虑单 个盘面的污秽监测，可W在龙口架处选用多个摄像头分别对全方位工况和单个盘面的污秽 进行监测，建议设置一个摄像头拍用于工况监测的整体画面，得到工况图像，另设一个摄像 头拍盘面的画面，得到憎水性图像，通过设置摄像头的旋转的角度，能够自动拍摄针对盘面 的若干位置。
[0077] 在服务器端，可W采用可见光图像识别与处理技术提高变压器运行可靠性，及时 检测及分析套管盘面水滴情况，对套管污秽情况进行分级；实验证明本发明所提供识别处 理技术能够在60秒钟能给出判断。
[0078] 进一步地，可W对套管工况进行监测判断，给出缺陷报告;提高套管运行可靠性。
[0079] 具体实施时，可W在服务器端提供具备可视化的用户界面，W网站的形式运行，网 页形式展示，实现人机交互。服务器端同时支持电脑和移动设备等用户终端的远程登录访 问，方便现场作业。
[0080] 进一步地，通过服务器端和高清摄像设备之间的光通信网络，通过系统在线自动 控制多台网络摄像机进行定时、定点图像抓拍、保存和传输功能。
[0081] 具体实施时，系统具备对采集的图像进行读取、截取图片、预处理(去噪、增强、边 缘检测等）、特征量提取和智能识别等功能。
[0082] 具体实施时，服务器端可具备数据库管理功能，可同步读取、存入、查询、插入、删 除、更新后台数据。
[0083] 具体实施时，服务器端可具备检测结果报警功能。
[0084] 参见图2,本发明实施例提供的在线监测方法，依次包括W下步骤：
[0085] 步骤1:图像采集；
[00化]步骤2:图像预处理；
[0087] 步骤3:图像区域分割；
[0088] 步骤4:套管污秽分级处理，套管工况缺陷识别。
[0089] 参见图3,虽然基本步骤相同，但由于处理对象不同，具体实施时，可W设计流程 为，根据图像采集结果，
[0090] 如果是憎水性图像，对憎水性图像进行图像预处理，图像区域分割，然后进行套管 污秽分级处理；
[0091] 如果是工况图像，取相应标准样本图像，对工况图像和标准样本图像分别进行图 像预处理、图像区域分割，然后进行套管工况缺陷识别。
[0092] W下对实施例的具体步骤实现进行说明：
[0093] 图像采集，即接收传输到服务器端的图像，包括绝缘子的憎水性图像和变压器套 管整体的工况图像。一般是实时采集并传输到服务器端，不同光照条件下变压器套管整体 完好的套管图像可W预先采集并存储到服务器端，根据实时得到的工况图像，取用相同光 照条件下变压器套管整体完好的套管图像，作为标准样本图像即可。
[0094] 进一步地，变压器套管长期工作在户外，拍摄图像如果不进行预处理对后面的工 况识别和污秽识别精度影响较大，本发明采用基于同态増晰的局部直方图均衡化对拍摄的 图像进行预处理。同态増晰方法在频域范围内通过选取合适的传递函数可W分别对高频分 量与低频分量产生不同的影响，使得图像的边缘与细节更加明显并增强了图像的对比度， 同时对噪声也有一定的抑制，但是由于压缩了灰度动态范围使得图像灰度级范围变窄，图 像亮度因此受到了一定的影响。局部直方图均衡化对图像处理的效果相比传统的均衡化方 法更佳，但是如果直接对复杂的憎水性图像和工况图像进行处理，会使噪声放大，不利于后 面的图像分割。所W考虑到两种方法的优缺点，本发明提出采用基于同态増晰的局部直方 图均衡化方法，对图像进行增强处理。即先采用同态増晰的方法对图像进行处理，使图像的 边缘与细节更加明显，然后对处理后的图像进行局部直方图均衡化，使图像的灰度动态范 围能获得比较均匀的调整。通过实验比较，该方法不仅可W对图像进行滤波处理，还可W增 强图像对比度，调整图像的灰度动态范围。经过实验比较该方法的处理，图像细节和边缘更 加清晰，图像的边缘和背景也更易区分，为后面图像的准确分割创造了良好的条件。
[00M]进一步地，在套管图像预处理的基础上，通过分析本方明发现传统的化nny算子在 梯度幅值计算和阔值设定两方面存在不足，因此在计算梯度幅值时将像素对角线方向考虑 进来并将OTSU阔值算法引入到化nny算子。改进的化nny算子用于图像边缘检测分割时，明 显改善了边缘的质量，但仍然存在边缘断线问题。紧接着，还利用数学形态学来修正W解决 边缘断线问题。基于此，本发明提出了基于Canny算子的复合绝缘子憎水性图像的分割算 法，分割后的图像能满足后续边缘提取的要求。
[0096] OTSU算法即最大类间方差法，本发明针对套管图像具有随机性和模糊性的特点， 利用最大类间方差法寻找一个能使得方差最大的最优解。该算法基于最小二乘法原理并在 灰度直方图的基础上推导而来，具有统计意义上的最佳分割阔值。利用该最佳分割阔值把 灰度图像分割成两部分，此时两部分间的方差最大，即目标与背景的分离性最大。
[0097] OTSU算法:如一幅灰度图像可能具有的灰度级数为L，灰度级为i的像素个数为m， 其中i G {0，1，…，L-I }，那么图像总共含有的像素个数为：
[009引
[0099] 灰度值为i的像素出现的概率为：
[0100] pi = ru/n
[0101] 然后，用灰度值k作为阔值将图像分成Ao和Al两部分，其中区域Ao的灰度范围为0~ k-1，区域Al的灰度范围为k~kl。
[0102] 整幅图像的平均灰度级为：
[0103]
[0104] X域Al中像素出现的概率：
[0105]
[0106]
[0107] E：域Al中像素的平均灰度：
[010 引
[0109]
[0110] 所W两个区域的类间方差为：
[0111]
[0112]
[011引上式中参数wo、wi、ii0和m均与阔值k有关，故类间方差的大小仅由阔值k决定。方差 可用来度量灰度分布的均匀性，如果其值愈大，表明构成图像的不同部分的差别愈大。通过 计算不同k值对应的类间方差S2化），选取使得S2化)最大时的那个k值作为最佳分割阔值。
[0114] 传统的化nny算子在计算梯度幅值与方向时，使用的是2X2邻域一阶偏导的有限 差分。因为采用的点较少，虽然能够比较准确的定位边缘，但是对噪声十分敏感，所W可能 丢掉一些真实边缘却带来虚假的边缘。使用3X3邻域来计算梯度，不仅计算水平(X)和垂直 (y)方向的梯度幅值，还计算45°和135°方向的梯度幅值，即考虑到像素的对角线方向，因而 有效地降低了噪声带来的影响。
[0115] 经过改进的化nny算子用于图像边缘检测时的步骤如下：
[0116] (1)利用高斯滤波器对预处理后的绝缘子憎水性图像进行平滑滤波处理。
[0117] (2)用邻域一阶偏导的有限差分来计算平滑后图像的梯度幅值和方向。
[0118] (3)对梯度幅值图像进行处理，抑制梯度幅值的非极大值。
[0119] (4)用OTSU算法求出最佳分割阔值，然后用该最佳阔值进行检测并连接边缘。
[0120] 进一步地，套管图像的区域分割步骤还使用数学形态学边缘检测算法对W上边缘 图像进行增强处理。
[0121] 进一步地，为实现复合绝缘子套管的污秽分级，改进的化nny算子虽然能够大致提 取出憎水性图像的边缘轮廓，但是由于憎水性图像和套管破损、裂痕图像的复杂性，提取的 边缘因为灰度问题不免存在断裂的问题。运个容易对后面的识别造成误判，同时，图像分割 的最终目的是要计算出憎水性图像的特征码及判断是否有裂痕等。所W还要填充边缘轮廓 使得图像尽可能忠实原图像并消除掉小而无意义的水滴。基于上述的原因，本发明对边缘 检测后的图像进行膨胀、区域填充、腐蚀的数学形态学运算进行修正。
[0122] (1)膨胀运算
[0123] 膨胀运算是一种使图像中的目标"变粗"或"生长"的操作，它可W填补边缘的缝 隙，解决巧錄齡錶的问颗。,4 ? X亲元图像A被结构元素 S膨胀，定义为：
[0124]
[0125] 式中，Sxy表示将结构元素 S的原点移到点（x，y)时得到的图像；d)表示空集。图像A 被结构元素 S膨胀得到的图像可W认为是运样的点（x，y)构成的集合:如果结构元素 S的原 点移到点(x，y)时，S与A的交集不为空集。
[0126] (2)腐蚀运算
[0127] 腐蚀运算是一种"细化"或"收缩"图像中目标的操作，它能够消除边缘附近的毛 刺，使得图像边缘变得光滑，同时去掉小而无意义的水滴。^该及表示图像A被结构元素 S腐 蚀，定义为：
[012 引
[0129] 图像A被结构元素 S腐蚀得到的图像可W认为是运样的点（x，y)构成的集合:如果 结构元素 S的原点移到点(x，y)时，S完全包含于A中。
[0130] (3)区域填充
[0131] 通过使用数学形态学中的膨胀运算、交集运算、求补运算W及补集运算，填充边缘 内部的孔桐。从图像A边缘内部一点Xo开始，将其灰度记为1，按照下面的表达式填充整个区 域：
[0132] .....
[01削式中，S表示结构元素;Ae表示图像A的补集。结束条件为Xk = Xk-I,此时Xk U A包括图 像边缘及其内部填充的集合。
[0134] 本发明对经过W上改进化nny算子边缘检测后的图像，首先对其用扁平的直线型 结构元素进行膨胀运算，使整个图像的边缘变得连续;然后紧接着对图像进行区域填充，使 得图像边界内部的孔桐得到填补;最后，利用菱形结构元素对区域填充后的图形进行腐蚀 运算，消除边缘附近的毛刺，使得图像边缘变得光滑，同时去掉小而无意义的水滴。为后面 的污秽等级识别提供基础。
[0135] 进一步地，本发明的复合绝缘子套管的污秽分级考虑识别效果和工程实用采用基 于最小二乘支持向量机化east Squares Su卵ort Vector Machines,简称LS-SVM)的复合 绝缘子套管的憎水性等级识别处理，由于运行于现场的复合绝缘子套管的复杂性，其憎水 性图像包含的信息是十分丰富的。为了准确识别复合绝缘子套管的憎水性等级，找到能够 表征憎水性的特征码并将其提取出来是非常必要的。对于经过W上图像区域分割步骤所得 分割后的憎水性图像，提取如下四个与复合绝缘子憎水性有关的量作为特征码:最大水珠 积比K，水珠覆盖率C，最大水珠偏屯、率e，最大水珠形状因子f C。
[0136] 最大水珠面积比K
[0137]
[0138] 式中，S表示最大水珠的面积;XY为整幅图像的面积。
[0139] 水珠覆盖率C
[0140]
[0141] 式中，n表示识别出的水珠个数;Si为第i个水珠的面积。
[0142] 最大水珠偏屯、率e
[0143]
[0144] ，a表示最大水珠的长轴，即最大水珠边界上相距最远的两像素点之间的线段 长度;b表示最大水珠的短轴，即与长轴垂直且与之构成的外接矩形恰好能包围最大水珠边 界的线段长度。
[0145] 最大水珠形状因子f C
[0146]
[0147] 式中，L表示最大水珠的周长;S表示最大水珠的面积。
[0148] LS-SVM方法简化了计算的复杂性。另外，由于LS-SVM采用了最小二乘算法，因此运 算速度明显快于支持向量机的其它方法，实际工程应用效果良好。
[0149] LS-SVM估计算法基于如下思想:通过事先选择的非线性映射将输入向量映射到高 维特征空间，在运个空间中构造最优决策函数。在构造最优决策函数时，利用了结构风险最 小化原则。并巧妙的利用原空间的核函数取代了高维特征空间中的点积运算。
[0150] 设样本为n维向量，定义区域m的样本及其值表示为：
[0151] XI，......,XiGR"
[0152] R为实数。其中，Xi表示区域m的第i个样本的n维向量样本值，i = l，2，...l。
[0153] 首先用一非线性映射W ( X )把样本从原空间R n映射到特征空间
[0154] X是原样本值。
[01巧]在运个高维特征空间中构造最优决策函数.r(.Y) = VV，例.V) +
[0156] W是连接权值，
[0157] b是阔值。
[0158] 运样非线性估计函数转化为高维特征空间中的线性估计函数。利用结构风险最小 化原则，寻找W，b就是最小化。
[0159]
，
[0160] 其甲，W 3空制模型的复杂度，C为正规化参数。Remp为误差控制函数，实施例采 用e不敏感损失函数。常用的损失函数有线性e损失函数，二次e损失函数，huber损失函数。 具体实施时选取不同的损失函数，可构造不同形式的支持向量机。最小二乘支持向量机在 优化目标中的损失函数为误差Cl的二次项。故优化问题为：
[0161]
[0162]
[0163] 其中，yi为构造的最优决策函数。
[0164] 用拉格朗日法求解运个优化问题：
[01 化]
[0166] 其中ai，i = 1，. . .，1，是拉格朗日乘子，根据优化条件：
[0167]
[016 引
[0169] (4)
[0170] 定义核函数
I满足Mercer条件的对称函数。根据(4)尤化I、品挪掉化责龙偏純化古浩.
[0171]
[0172]
[0173]
[0174] 然后将检测样本化，C，e，fc)输入上述非线性方程，即xi = (化1，片，61，片1))，可得 到该检测样本所属的憎水等级。其中K表示最大水珠(或水迹)面积比，C表示水珠(或水迹） 覆盖率，e表示最大水珠(或水迹)偏屯、率，fc表示最大水珠(或水迹)形状因子。其中K(xi，xj) 为核函数，Si是拉格朗日乘子，i二1, . . ? , 1，j二1, . . ? , 1，b是阔值。
[0175] 实际应用效果见下表1。
[0176] 表1: LS-SVM 识别率
[0177]
[0178] 进一步地，为进行变压器套管工况识别，与憎水性图像的数学形态学处理不同，本 发明对完好的套管图像（即标准样本图像)和实时采集的套管图像（即工况图像)在经过前 面图像预处理和图像区域分割步骤中的边缘提取后，进行数学形态学处理时，本发明实施 例进行如下操作：
[0179] 首先，选取方形结构元素对两幅图像分别进行开运算，并用圆盘形结构元素对他 们实行膨胀运算得到图像矩阵A和B。同时，选取方形结构元素分别对他们进行闭运算，再用 圆盘形结构元素对他们进行腐蚀得到图像矩阵Al和B1。具体实施时，本领域技术人员可自 行预设结构元素，例如边长为8的方形结构元素。
[0180] 分别用A减去Al, B减去Bl可W得到标准样本图像和待检测图像的边缘。
[0181] 而后使用最大赌法求取阔值得到标准二值边缘图像。对两图像边缘进行中值滤波 后还使用bwareaopen函数进行后处理，去除图像中的多余杂点。
[0182] 进一步地，套管图像缺陷识别步骤还调整摄像头拍摄方向及固定时间点拍摄，在 尽量相同的光照环境下（如每天的对应时间），将待检测图像边缘矩阵减去标准样本图像边 缘矩阵并取绝对值，进行标准样本图像和待检测图像的对比，可W-定程度上排除绝缘子 上污秽，斑点，及光照投影对裂纹识别的影响，同时还可再次使用bwareaopen函数去除小面 积的污秽，斑点，仅保留可能存在的较狭长的裂缝图像。为减小同态增晰对光照差别造成的 放大崎变作用，比较标准样本图像和待检测图像的边缘图像最大灰度差的绝对值和预设的 阔值，成为判断是否存在裂痕的阔值的依据。若绝对值未超过阔值则判定为无裂痕，若超过 阔值则提出告警。进一步地，可根据多幅不同光照环境下正常图像和相应环境下待检测图 像逐一进行多次对比，进一步判断，最终如果都超过每次对比后的阔值，则判断存在裂痕， 进行报警。具体实施时，阔值可由本领域技术人员自行预设，建议采用最大赌法求取。本发 明试验效果非常好，套管如果破损、轻微裂痕都能很快识别。
[0183] 具体实施时，本发明所提供方法可基于软件技术实现自动运行流程，也可采用模 块化方式实现相应系统。
[0184] -种基于计算机视觉的变压器复合绝缘子套管监测系统，在各变压器上空龙口架 处分别安装摄像头，拍摄所得图像传输到服务器端，在服务器端设置图像采集模块、图像预 处理模块、图像区域分割模块、套管污秽分级处理模块及套管工况缺陷识别模块；
[0185] 所述图像采集模块，用于采集拍摄所得图像，包括绝缘子的憎水性图像和变压器 套管整体的工况图像；
[0186] 所述图像预处理模块，用于采用基于同态増晰的局部直方图均衡化方式实现预处 理；
[0187] 所述图像区域分割模块，用于对预处理后的图像，采用改进的化nny算子检测图像 边缘，然后对边缘检测后的图像采用数学形态学运算进行修正，得到分割结果；
[0188] 所述采用改进的化nny算子检测图像边缘实现方式如下，
[0189] 利用高斯滤波器对预处理后的绝缘子憎水性图像进行平滑滤波处理；
[0190] 用邻域一阶偏导的有限差分来计算平滑后图像的梯度幅值和方向，包括使用3X3 邻域来计算梯度，包括计算水平和垂直方向、的梯度幅值，还计算45°和135方向的梯度幅 值；
[0191] 对梯度幅值图像进行处理，抑制梯度幅值的非极大值；
[0192] 用OTSU算法求出最佳分割阔值，然后用该最佳阔值进行检测并连接边缘；
[0193] 所述套管污秽分级处理模块，用于根据憎水性图像的分割结果，提取四个与复合 绝缘子憎水性有关的量作为特征码，包括最大水珠面积比K，水珠覆盖率C，最大水珠偏屯、率 e，最大水珠形状因子f。，采用基于最小二乘支持向量机的方式得到检测样本所属的憎水等 级；
[0194] 島十水铁巧积比K计算如下，
[0195]
[0196] 式中，S表示最大水珠的面积，XY为整幅图像的面积；
[0197] 水珠覆盖率C计算如下，
[019 引
[0199] 式中，n表示识别出的水珠个数，Si为第i个水珠的面积；
[0200] 最大水珠偏屯、率e计算如下，
[0201]
[0202] 式中，a表示最大水珠的长轴，b表示最大水珠的短轴；
[0203] '?状因子fc计算如下，
[0204]
[0205] 式中，L表示最大水珠的周长，S表示最大水珠的面积；
[0206] 所述套管工况缺陷识别模块，用于将工况图像作为待检测图像，根据工况图像的 分割结果，将待检测图像边缘矩阵减去标准样本图像边缘矩阵并取绝对值，进行标准样本 图像和待检测图像的对比；使用bwareaopen函数去除污秽和斑点，将标准样本图像和待检 测图像的边缘图像最大灰度差的绝对值和最大赌法求取的阔值相组合成为判断是否存在 裂痕的阔值，若未超过阔值则判定为无裂痕，若是则提出告警。
[0207] 各模块具体实现可参见相应步骤，本发明不予寶述。
[0208] 本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领 域的技术人员可W对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替 代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
【主权项】
1. 一种基于计算机视觉的变压器复合绝缘子套管监测方法，其特征在于：在各变压器 上空龙门架处分别安装摄像头，拍摄所得图像传输到服务器端，在服务器端进行图像采集、 图像预处理、图像区域分割、套管污秽分级处理及套管工况缺陷识别； 所述图像采集，包括采集拍摄所得图像，包括绝缘子的憎水性图像和变压器套管整体 的工况图像； 所述图像预处理采用基于同态増晰的局部直方图均衡化方式实现； 所述图像区域分割，包括对预处理后的图像，采用改进的Canny算子检测图像边缘，然 后对边缘检测后的图像采用数学形态学运算进行修正，得到分割结果； 所述采用改进的Canny算子检测图像边缘实现方式如下， 利用高斯滤波器对预处理后的绝缘子憎水性图像进行平滑滤波处理； 用邻域一阶偏导的有限差分来计算平滑后图像的梯度幅值和方向，包括使用3X3邻域 来计算梯度，包括计算水平和垂直方向、的梯度幅值，还计算45°和135°方向的梯度幅值； 对梯度幅值图像进行处理，抑制梯度幅值的非极大值； 用OTSU算法求出最佳分割阈值，然后用该最佳阈值进行检测并连接边缘； 所述套管污秽分级处理，包括根据憎水性图像的分割结果，提取四个与复合绝缘子憎 水性有关的量作为特征码，包括最大水珠面积比K，水珠覆盖率C，最大水珠偏心率e，最大水 珠形状因子f。，采用基于最小二乘支持向量机的方式得到检测样本所属的憎水等级； 最大水珠面积比K计算如下，式中，S表示最大水珠的面积，XY为整幅图像的面积； 水珠覆盖率C计算如下，式中，η表示识别出的水珠个数，&为第i个水珠的面积； 最大水珠偏心率e计算如下，式中，a表;^最大水珠的长轴，b表;^最大水珠的短轴； 最大水珠形状因子f。计算如下，式中，L表示最大水珠的周长，S表示軍木的囬不只； 所述套管工况缺陷识别，包括将工况图像作为待检测图像，根据工况图像的分割结果， 将待检测图像边缘矩阵减去标准样本图像边缘矩阵并取绝对值，进行标准样本图像和待检 测图像的对比；使用bwareaopen函数去除污秽和斑点，将标准样本图像和待检测图像的边 缘图像最大灰度差的绝对值和最大熵法求取的阈值相组合成为判断是否存在裂痕的阈值， 若未超过阈值则判定为无裂痕，若是则提出告警。2. 根据权利要求1所述基于计算机视觉的变压器复合绝缘子套管监测方法，其特征在 于:针对憎水性图像，对边缘检测后的图像采用数学形态学运算进行修正的实现方式如下， 采用扁平的直线型结构元素进行膨胀运算，使整个图像的边缘变得连续;然后紧接着 对图像进行区域填充，使得图像边界内部的孔洞得到填补;最后，利用菱形结构元素对区域 填充后的图形进行腐蚀运算，消除边缘附近的毛刺，使得图像边缘变得光滑，同时去掉过小 而无意义的水滴。3. 根据权利要求1所述基于计算机视觉的变压器复合绝缘子套管监测方法，其特征在 于:针对工况图像，对边缘检测后的图像采用数学形态学运算进行修正的实现方式如下， 对待检测图像和标准样本图像经预处理和图像边缘检测后的结果，采用方形结构元素 对两幅图像分别进行开运算，并用圆盘形结构元素实行膨胀运算得到图像矩阵A和B;同时， 采用方形结构元素对两幅图像分别进行闭运算，再用圆盘形结构元素进行腐蚀得到图像矩 阵A1和B1; 分别用A减去A1、B减去B1，得到标准样本图像和待检测图像的边缘； 使用最大熵法求取阈值得到标准二值边缘图像，对两图像边缘进行中值滤波后还使用 bwareaopen函数进行后处理，去除图像中的多余杂点。4. 根据权利要求1或2或3所述基于计算机视觉的变压器复合绝缘子套管监测方法，其 特征在于:将工况图像作为待检测图像时，标准样本图像采用相同的光照环境下的拍摄结 果。5. 根据权利要求4所述基于计算机视觉的变压器复合绝缘子套管监测方法，其特征在 于:根据多幅不同光照环境下待检测图像和标准样本图像逐一进行对比，如果都超过每次 对比后的阈值，则判断存在裂痕，进行报警。6. -种基于计算机视觉的变压器复合绝缘子套管监测系统，其特征在于：在各变压器 上空龙门架处分别安装摄像头，拍摄所得图像传输到服务器端，在服务器端设置图像采集 模块、图像预处理模块、图像区域分割模块、套管污秽分级处理模块及套管工况缺陷识别模 块； 所述图像采集模块，用于采集拍摄所得图像，包括绝缘子的憎水性图像和变压器套管 整体的工况图像； 所述图像预处理模块，用于采用基于同态増晰的局部直方图均衡化方式实现预处理； 所述图像区域分割模块，用于对预处理后的图像，采用改进的Canny算子检测图像边 缘，然后对边缘检测后的图像采用数学形态学运算进行修正，得到分割结果； 所述采用改进的Canny算子检测图像边缘实现方式如下， 利用高斯滤波器对预处理后的绝缘子憎水性图像进行平滑滤波处理； 用邻域一阶偏导的有限差分来计算平滑后图像的梯度幅值和方向，包括使用3X3邻域 来计算梯度，包括计算水平和垂直方向、的梯度幅值，还计算45°和135°方向的梯度幅值； 对梯度幅值图像进行处理，抑制梯度幅值的非极大值； 用0TSU算法求出最佳分割阈值，然后用该最佳阈值进行检测并连接边缘； 所述套管污秽分级处理模块，包括根据憎水性图像的分割结果，提取四个与复合绝缘 子憎水性有关的量作为特征码，包括最大水珠面积比K，水珠覆盖率C，最大水珠偏心率e，最 大水珠形状因子f。，采用基于最小二乘支持向量机的方式得到检测样本所属的憎水等级； 最大水珠面积比K计算如下，式中，S表示最大水珠的面积，XY为整幅图像的面积； 水珠覆盖率C计算如下，式中，η表示识别出的水珠个数，Si为弟i个水秌的囬枳； 最大水珠偏心率e计算如下，式中，a表;^最大水珠的长轴，b表;^最大水珠的短轴；最大水珠形状因子f。计算如下， 式中，L表不最大水珠的周长，S表不取人7」Vi木tfj IMJ不只； 所述套管工况缺陷识别模块，用于将工况图像作为待检测图像，根据工况图像的分割 结果，将待检测图像边缘矩阵减去标准样本图像边缘矩阵并取绝对值，进行标准样本图像 和待检测图像的对比；使用bwareaopen函数去除污秽和斑点，将标准样本图像和待检测图 像的边缘图像最大灰度差的绝对值和最大熵法求取的阈值相组合成为判断是否存在裂痕 的阈值，若未超过阈值则判定为无裂痕，若是则提出告警。7. 根据权利要求6所述基于计算机视觉的变压器复合绝缘子套管监测系统，其特征在 于:针对憎水性图像，对边缘检测后的图像采用数学形态学运算进行修正的实现方式如下， 采用扁平的直线型结构元素进行膨胀运算，使整个图像的边缘变得连续;然后紧接着 对图像进行区域填充，使得图像边界内部的孔洞得到填补;最后，利用菱形结构元素对区域 填充后的图形进行腐蚀运算，消除边缘附近的毛刺，使得图像边缘变得光滑，同时去掉过小 而无意义的水滴。8. 根据权利要求6所述基于计算机视觉的变压器复合绝缘子套管监测系统，其特征在 于:针对工况图像，对边缘检测后的图像采用数学形态学运算进行修正的实现方式如下， 对待检测图像和标准样本图像经预处理和图像边缘检测后的结果，采用方形结构元素 对两幅图像分别进行开运算，并用圆盘形结构元素实行膨胀运算得到图像矩阵A和B;同时， 采用方形结构元素对两幅图像分别进行闭运算，再用圆盘形结构元素进行腐蚀得到图像矩 阵A1和B1; 分别用A减去A1、B减去B1，得到标准样本图像和待检测图像的边缘； 使用最大熵法求取阈值得到标准二值边缘图像，对两图像边缘进行中值滤波后还使用 bwareaopen函数进行后处理，去除图像中的多余杂点。9. 根据权利要求6或7或8所述基于计算机视觉的变压器复合绝缘子套管监测系统，其 特征在于:将工况图像作为待检测图像时，标准样本图像采用相同的光照环境下的拍摄结 果。10. 根据权利要求9所述基于计算机视觉的变压器复合绝缘子套管监测系统，其特征在 于:根据多幅不同光照环境下待检测图像和标准样本图像逐一进行对比，如果都超过每次 对比后的阈值，则判断存在裂痕，进行报警。
【文档编号】G06T7/00GK105957088SQ201610315778
【公开日】2016年9月21日
【申请日】2016年5月12日
【发明人】蔡日, 孙云莲, 王庆斌, 莫雄心, 黄嘉文, 黄雅鑫, 李和阳, 黎圣旻, 余军伟
【申请人】广东电网有限责任公司云浮供电局, 武汉大学