发明名称 一种绝缘子图像缺陷检测方法
摘要 本发明公开了一种绝缘子图像缺陷检测方法。其特征在于包括如下步骤:首先利用彩色图像中绝缘子的颜色呈现聚类的特性,通过在Lab颜色空间中利用K均值聚类算法将绝缘子从复杂的背景图像中初步分割出来;然后利用绝缘子在图像形成的长条对称特性,采用基于主成分分析的连通区域判决方法定位出绝缘子的具体位置;在此基础上,通过统计分析绝缘子相邻盘片之间的距离识别出绝缘子缺陷区域。该方法受不同光线和复杂背景的影响较小,提高了绝缘子缺陷检测的准确性和鲁棒性,具有较高的工程实用价值。
申请公布号 CN103886610B 申请公布日期 2016.08.17
申请号 CN201410136596.8 申请日期 2014.04.05
申请人 东北电力大学 发明人 胡建平;谢琪;张杰;李鹏松
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 吉林市达利专利事务所 22102 代理人 陈传林
主权项 一种绝缘子图像缺陷检测方法,其特征是:包括以下步骤:1)基于Lab颜色空间的K均值聚类把直升机巡检获取的输电线路的彩色绝缘子图像从RGB空间转换到Lab空间,然后将表示颜色的色度分量a和b作为属性值向量进行K均值聚类分割,进而将其分割成多个类别,考虑到绝缘子的颜色呈现聚类特性,则其必然属于某一个类别;2)基于主成分分析的绝缘子定位经过基于K均值的聚类方法分割后存在几个类别,每个类别又以多个连通区域的形式出现,需要从众多的连通区域中找到绝缘子所在的连通区域;由于绝缘子整体呈现规则的对称长条特性,因此采用基于主成分分析的连通区域形状判决方法进行绝缘子精确定位,主要步骤如下:a)预处理:通过数学形态学操作对聚类分割后的图像进行腐蚀、膨胀去掉小的噪声,同时对图像中连通区域的孔洞进行填充,以保证分割后绝缘子的完整性;b)计算连通区域信息:利用主成分分析计算连通区域的轴向、宽、高信息;c)区域合并:通过对特定的连通区域进行合并,以保证由于缺陷导致绝缘子分开的部分连接在一起;d)连通区域形状判决:判定所有连通区域的形状得到绝缘子所在的精确位置;3)基于绝缘子盘片距离分析的缺陷检测将基于主成分分析的绝缘子定位所得的绝缘子二值图像进行旋转,使得绝缘子的主轴方向平行于竖直y轴方向,再对绝缘子的二值图像通过沿竖直y轴方向进行积分投影,得到积分投影曲线;通过绝缘子盘片完好处两相邻极大值之间的距离近似相等,而在绝缘子缺陷处两相邻极大值的距离明显变大的特性,精确定位出绝缘子的缺陷位置,并在图像中标记;所述步骤1)中把直升机巡检获取的输电线路的彩色绝缘子图像从RGB空间转换到Lab空间,具体计算公式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>X</mi><mo>=</mo><mn>0.4125</mn><msup><mi>R</mi><mo>&prime;</mo></msup><mo>+</mo><mn>0.3576</mn><msup><mi>G</mi><mo>&prime;</mo></msup><mo>+</mo><mn>0.1805</mn><msup><mi>B</mi><mo>&prime;</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><mi>Y</mi><mo>=</mo><mn>0.2126</mn><msup><mi>R</mi><mo>&prime;</mo></msup><mo>+</mo><mn>0.7152</mn><msup><mi>G</mi><mo>&prime;</mo></msup><mo>+</mo><mn>0.0722</mn><msup><mi>B</mi><mo>&prime;</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><mi>Z</mi><mo>=</mo><mn>0.0193</mn><msup><mi>R</mi><mo>&prime;</mo></msup><mo>+</mo><mn>0.1192</mn><msup><mi>G</mi><mo>&prime;</mo></msup><mo>+</mo><mn>0.9505</mn><msup><mi>B</mi><mo>&prime;</mo></msup></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000957423760000011.GIF" wi="966" he="198" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>L</mi><mo>=</mo><mn>116</mn><mi>f</mi><mrow><mo>(</mo><mi>Y</mi><mo>/</mo><msub><mi>Y</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mn>16</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>a</mi><mo>=</mo><mn>500</mn><mi>f</mi><mrow><mo>(</mo><mi>X</mi><mo>/</mo><msub><mi>X</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mi>f</mi><mrow><mo>(</mo><mi>Y</mi><mo>/</mo><msub><mi>Y</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>b</mi><mo>=</mo><mn>200</mn><mi>f</mi><mrow><mo>(</mo><mi>Y</mi><mo>/</mo><msub><mi>Y</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mi>f</mi><mrow><mo>(</mo><mi>Z</mi><mo>/</mo><msub><mi>Z</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000957423760000012.GIF" wi="934" he="198" /></maths>式中<img file="FDA0000957423760000013.GIF" wi="1270" he="135" />X<sub>0</sub>,Y<sub>0</sub>,Z<sub>0</sub>为标准光源D65的三刺激值,其值为X<sub>0</sub>=95.045,Y<sub>0</sub>=100,Z<sub>0</sub>=108.255;所述步骤2)中的步骤b)采用如下的方式计算连通区域的轴向、宽、高信息:对于连通区域A=(p<sub>1</sub>,…,p<sub>k</sub>,…,p<sub>N</sub>),p<sub>k</sub>=(x<sub>k</sub>,y<sub>k</sub>),k=1,…N表示该连通区域中的点,则其协方差矩阵为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>M</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msub><mi>p</mi><mi>k</mi></msub><mo>-</mo><mi>c</mi><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><msub><mi>p</mi><mi>k</mi></msub><mo>-</mo><mi>c</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000957423760000021.GIF" wi="1213" he="118" /></maths>式中c表示该连通区域的质心,(·)<sup>T</sup>表示向量的转置;通过计算协方差矩阵M的最大特征值对应的特征向量得到该连通区域的主轴,计算最小特征值对应的特征向量得到该连通区域的次轴,计算沿着两个轴线方向的最远距离得到连通区域的宽度w和高度h;所述的步骤2)的步骤c)采用如下的方式对特定的连通区域进行合并:记连通区域A<sub>i</sub>和A<sub>j</sub>的单位主轴方向分别为d<sub>i</sub>和d<sub>j</sub>,质心分别为c<sub>i</sub>和c<sub>j</sub>,宽度分别为w<sub>i</sub>和w<sub>j</sub>,如果下列关系式成立,则合并这两个连通区域:|w<sub>i</sub>‑w<sub>j</sub>|&lt;w<sub>T</sub>  (4)arccos(|(d<sub>i</sub>·d<sub>j</sub>)|)&lt;a<sub>T</sub>  (5)<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>arccos</mi><mrow><mo>(</mo><mo>|</mo><mo>(</mo><msub><mi>d</mi><mi>i</mi></msub><mo>&CenterDot;</mo><mfrac><mrow><msub><mi>c</mi><mi>i</mi></msub><mo>-</mo><msub><mi>c</mi><mi>j</mi></msub></mrow><mrow><mo>|</mo><mo>|</mo><msub><mi>c</mi><mi>i</mi></msub><mo>-</mo><msub><mi>c</mi><mi>j</mi></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>)</mo></mrow><mo>|</mo><mo>)</mo><mo>&lt;</mo><msub><mi>a</mi><mi>T</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000957423760000022.GIF" wi="1326" he="167" /></maths>式中(·)表示两个向量内积,w<sub>T</sub>和a<sub>T</sub>分别表示选定的宽度差异阈值和偏离角度阈值,设w<sub>T</sub>=3,a<sub>T</sub>=1,即合并的两个子串的宽度应该不超过3个像素,偏离角度不超过1弧度;所述的步骤2)中的步骤d)采用如下的方式判定连通区域的形状得到绝缘子所在的精确位置:对于连通区域A,设其宽和高分别为w和h,连通区域中点的个数为N,计算如下形状判决度量:高宽比k<sub>1</sub>=h/w;连通区域最小包围盒里物体像素所占比例k<sub>2</sub>=N/(w*h);连通区域的对称度k<sub>3</sub>,设连通区域中关于主轴的对称点落在连通区域中的点的数目为N<sub>1</sub>,则对称度为k<sub>3</sub>=N<sub>1</sub>/N;如果形状判决度量k<sub>1</sub>,k<sub>2</sub>,k<sub>3</sub>都能超过给定的阈值,设k<sub>1</sub>&gt;8,k<sub>2</sub>&gt;0.7,k<sub>3</sub>&gt;0.9,则判定该连通区域是绝缘子;所述的步骤3)中的积分投影方式计算如下:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>f</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>w</mi></munderover><mi>f</mi><mrow><mo>(</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mi>n</mi><mo>=</mo><mn>1</mn><mo>...</mo><mi>h</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000957423760000031.GIF" wi="1438" he="118" /></maths>式中w,h分别表示绝缘子所在连通区域的宽度和高度,f(j,n)表示位置(j,n)处的函数值,如果该点是绝缘子上的点,取值为1,如果是背景点,取值为0。
地址 132012 吉林省吉林市船营区长春路169号