一种基于局部区域聚类的电力设备故障区域提取方法技术

本发明专利技术提出了一种基于局部区域聚类的电力设备故障区域提取方法。本发明专利技术构建电力故障图像的Mediodshift区域聚类模型；根据Mediodshift区域聚类特性调整邻域像素的灰度；对调整后邻域像素再进行局部区域聚类；阈值自高向低变化直至故障区域聚类完成；提高对整幅图像中故障区域像素聚类效率。本发明专利技术优点为局部区域Mediodshift聚类以及调整邻域相似灰度；自高向低的聚类阈值分割机制提升故障区域提取效率。

A method of power equipment fault region extraction based on local region clustering

【技术实现步骤摘要】
一种基于局部区域聚类的电力设备故障区域提取方法
本专利技术属于电力故障检测领域，尤其涉及一种基于局部区域聚类的电力设备故障区域提取方法。
技术介绍
红外热成像仪是一种在不停电的状态下对配电设备进行现状分析和故障诊断的工具，能有效的发现系统故障并快速成像记录，具有检测速度快、准确率高、费用低、灵敏性好和安全性高等优点。然而，传统的电力设备故障红外诊断时，通常需要运维人员对电力设备故障红外图像的认知以及人工方式采集、录入、上传等操作，使得整个检测过程效率低下，特别对于含有众多电力设备的变电站，人工诊断将变得异常繁琐。近年来，随着人工智能的兴起，一些智能化的电力设备运检装置得到了国内外学者的广泛研究，例如智能巡检机器人、巡检无人机等，其通过携带红外热像仪等装置，在线捕获电力设备红外成像信息，并传至控制室再进行人工分析。尽管这种方式减少了人工采集图片的工作量，但是仍然需要依赖人工诊断与处理。随着图像处理技术的发展以及智能化、自动化诊断的迫切需求，研究者开始着手采用计算机视觉、模式识别等技术解决电力设备故障智能预防、检测等问题。一般而言，智能化诊断包含了红外图像感兴趣区域的提取、分类、识别等主要过程，其中感兴趣区域(故障区域)的提取是整个系统中首要解决的关键内容，它在一定程度上决定了后续故障辨识的成败，也影响电力设备状态信息的分析。为此，研究针对电力设备故障的区域的高效自适应的红外自动检测方法非常重要。
技术实现思路
本专利技术主要解决了Mediodshift算法自适应迭代获取相似特性的局部区域。通过对区域聚类中心计算以及对灰度值介于两个均值之间的像素进行灰度调整，激励具有相似灰度值的像素进行局部聚类，惩罚具有相似灰度值的像素聚类，提升算法局部聚类性能，可不需要考虑迭代次数。本专利技术还解决了整幅红外图像故障区域提取效率的问题。通过引入阈值自高向低的分割机制，可直接对已经分割得到的区域进行邻域像素聚类，可不需要重新对整幅图像进行全局Mediodshift聚类。本专利技术提出的上述技术问题主要是通过下述技术方案得以解决的：本专利技术的目的是提供一种基于局部区域聚类的电力设备故障区域提取方法，通过mediodshift算法对局部区域像素的聚类，并提供了邻域像素的灰度调节机制，激励具有相似灰度值的像素进行聚类，提升了算法对故障区域灰度不均匀等情况下的自适应聚类能力；引入阈值自高向低的分割机制，解决了红外图像中用聚类方法提取故障区域速度慢的问题。本专利技术的方法为一种基于局部区域聚类的电力设备故障区域提取算法，具体包括以下步骤：步骤1：构建电力故障图像的Mediodshift区域聚类模型；步骤2：根据Mediodshift区域聚类特性调整邻域像素的灰度；步骤3：对调整后邻域像素进行局部区域聚类；步骤4：阈值自高向低变化，循环步骤2-4，直至故障区域聚类完成。作为优选，步骤1中所述Mediodshift区域聚类模型为：假设给定采样点xi∈Rd，i＝1，…，n，则点x处的概率密度估计为：其中，Φ(·)为高斯核函数，h为固定带宽，d表示维数，则聚类中心点可以通过最小化下式：迭代得到，其中下标k表示迭代次数；yk表示当前中心点的位置；yk+1为迭代到下一步的中心点，且选自样本集{xi}中满足上式的一个点；φ(x)＝-Φ′(x)，对等上式右边一阶导估计yk+1，可得由于在迭代过程中，yk是来自样本集{xi}中点。必然地，yk+1属于样本集；因此，在不断的迭代情况下，每一次得到的样本中心点都只需要计算在样本集中{xi}的结果，最终当前的样本中心会收敛到了概率密度极大的中心位置；由于二维图像空间平面和灰度空间是两个独立的空间，因此本文将每个空间的核进行相乘，组合得到新的核函数：其中Φ(·)为高斯核函数，hs为图像空间带宽，hr是灰度带宽，C为归一化常数，x为采样点，包含了图像空间特征xs以及灰度空间特征xr；作为优选，步骤2中所述根据Mediodshift区域聚类特性调整邻域像素的灰度，具体方法为：令m2(n)和m1(n)分别为目标区域Ω2(n)和背景区域Ω1(n)的均值，n为迭代次数，整个图像区间满足Ω＝Ω1∪Ω2；任取一个像素x∈N，如果该点的灰度值大于m2(n)，则必然会通过Medioshift归类为目标；而对于灰度值小于m1(n)，则迭代时Mediodshift会将其归类为背景；然而当该点的灰度值介于(m1(n)，m2(n))之间，则该像素的聚类会间接影响最终完整区域提取性能的优劣；因此，为了能够获取目标邻域相似灰度的像素，建立了如下的调整方法：其中，I(x)表示像素x的灰度值；ci,i＝1,2,3分别为系数，其累加和为1；sign()为符号项；由上式可知，右边由三部分构成，分别是激励项，自身项以及惩罚项。其中，自身项系数c2通常占据一定的比例；而其他两项由其区域特性决定。激励项c1，主要由该像素x的邻域中是否含有目标像素的多少决定，设置为：其中k为其邻域像素中属于目标的个数，当该像素邻域目标像素较多时，其本身为目标的概率也会因此而增加。而惩罚项因子c3主要由该像素灰度值是否处于目标和背景中间而进行惩罚，特别是低于中间值时，其得到的灰度值会降低，使得其不再被Mediodshift聚类；作为优选，步骤3中所述对调整后邻域像素进行局部区域聚类，具体为：在Mediodshift迭代聚类分割中，通常规定背景区域Ω1标记为0，目标区域Ω2标记为1，分别表示为Ω1＝{z|Rz(n)＝0,z∈Ω}和Ω2＝{z|Rz(n)＝1,z∈Ω}，其中R为整个图像Ω的聚类结果，n代表迭代时刻；其中，规定局部区域为在每次迭代过程中所产生的目标区域到背景区域的过渡区域；由此，可得到局部区域N为：N(n)＝{x|Mx(n)＝1}∩{x|Rx(n-1)＝0}式中为了获得邻域集合N，引入了数学形态学中的膨胀处理，由当前迭代得到的目标区域Ω2进行膨胀，得到M为：其中，表示膨胀符号，B为一个r×r的结构算子。通常，这个结构算子包含圆形、方形和棱形等；为了获得目标区域Ω2的邻域空间像素，在构造结构算子时采用方形，r为3的全1矩阵，以获取8-邻域空间像素；针对邻域像素x∈N，采用Medoidshift对局部区域N的像素进行分类判别；在规定标记区域为目标的情况下，当迭代得到的收敛点落入目标区域，则认为该候选像素为目标像素，即归类为Ω2；规定标记为目标区域像素不再重新分类，使得在迭代过程中满足背景区域则满足相对应地，会得到一组均值{mi(n)}i＝1,2；作为优选，步骤4中所述阈值自高向低变化直至局部区域聚类完成，具体为：引入了阈值分割机制，即采用当前的目标聚类均值m2(n)对整个图像进行分割，然后针对新的区域采用对邻域进行聚类，从而实现整幅图像中故障区本文档来自技高网...

【技术保护点】
1.一种基于局部区域聚类的电力设备故障区域提取算法，其特征在于，包括以下步骤：/n步骤1：构建电力故障图像的Mediodshift区域聚类模型；/n步骤2：根据Mediodshift区域聚类特性调整邻域像素的灰度；/n步骤3：对调整后邻域像素进行局部区域聚类；/n步骤4：阈值自高向低变化，循环步骤2-4，直至故障区域聚类完成。/n

【技术特征摘要】
1.一种基于局部区域聚类的电力设备故障区域提取算法，其特征在于，包括以下步骤：
步骤1：构建电力故障图像的Mediodshift区域聚类模型；
步骤2：根据Mediodshift区域聚类特性调整邻域像素的灰度；
步骤3：对调整后邻域像素进行局部区域聚类；
步骤4：阈值自高向低变化，循环步骤2-4，直至故障区域聚类完成。


2.根据权利要求1所述的基于局部区域聚类的电力设备故障区域提取算法，其特征在于：步骤1中所述Mediodshift区域聚类模型为：
假设给定采样点xi∈Rd，i＝1，…，n，则点x处的概率密度估计为：



其中，Φ(·)为高斯核函数，h为固定带宽，d表示维数，则聚类中心点可以通过最小化下式：



迭代得到，其中下标k表示迭代次数；yk表示当前中心点的位置；yk+1为迭代到下一步的中心点，且选自样本集{xi}中满足上式的一个点；φ(x)＝-Φ′(x)，对等上式右边一阶导估计yk+1，可得



由于在迭代过程中，yk是来自样本集{xi}中点，必然地，yk+1属于样本集；
因此，在不断的迭代情况下，每一次得到的样本中心点都只需要计算在样本集中{xi}的结果，最终当前的样本中心会收敛到了概率密度极大的中心位置；
由于二维图像空间平面和灰度空间是两个独立的空间，因此本文将每个空间的核进行相乘，组合得到新的核函数：



其中Φ(·)为高斯核函数，hs为图像空间带宽，hr是灰度带宽，C为归一化常数，x为采样点，包含了图像空间特征xs以及灰度空间特征xr。


3.根据权利要求1所述的基于局部区域聚类的电力设备故障区域提取算法，其特征在于：步骤2中所述根据Mediodshift区域聚类特性调整邻域像素的灰度，具体方法为：
令m2(n)和m1(n)分别为目标区域Ω2(n)和背景区域Ω1(n)的均值，n为迭代次数，整个图像区间满足Ω＝Ω1∪Ω2；
任取一个像素x∈N，如果该点的灰度值大于m2(n)，则必然会通过Medioshift归类为目标；而对于灰度值小于m1(n)，则迭代时Mediodshift会将其归类为背景；然而当该点的灰度值介于(m1(n)，m2(n))之间，则该像素的聚类会间接影响最终完整区域提取性能的优劣；
因此，为了能够获取目标邻域相似灰度的像素，建立了如下的调整方法：



其中，...

【专利技术属性】
技术研发人员：周正钦，罗传仙，龚浩，许晓路，江翼，吴念，周文，倪辉，朱诗沁，
申请(专利权)人：国网电力科学研究院武汉南瑞有限责任公司，
类型：发明
国别省市：湖北;42