CN105205829A - 基于改进的二维Otsu的变电站红外图像分割方法 - Google Patents
基于改进的二维Otsu的变电站红外图像分割方法 Download PDFInfo
- Publication number
- CN105205829A CN105205829A CN201510698451.1A CN201510698451A CN105205829A CN 105205829 A CN105205829 A CN 105205829A CN 201510698451 A CN201510698451 A CN 201510698451A CN 105205829 A CN105205829 A CN 105205829A
- Authority
- CN
- China
- Prior art keywords
- infrared image
- dimentional otsu
- transformer station
- image segmentation
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10048—Infrared image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20036—Morphological image processing
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于改进的二维Otsu的变电站红外图像分割方法。解决了现有的二维Otsu算法中出现的误分割,目标分割不完整的问题。本发明是采用形态学的权重自适应算法对输入的图像进行去噪,最后使用改进的二维Otsu算法进行图像分割。本发明能够有效地将目标区域完整的分割出来,为后续目标区域的特征提取奠定良好的基础。
Description
技术领域
本发明属于电力系统变电站设备图像处理领域,具体涉及一种基于改进的二维Otsu的变电站红外图像分割方法。
背景技术
由于红外成像技术是对设备表面辐射的红外光线进行非接触、远距离热成像检测,不受电场干扰,因此具有直观、准确、灵敏度高、快速、安全、应用范围广等特点,已成为变电站设备健康状态监测和故障诊断的重要手段。
红外图像分割技术是变电站电气设备故障监测与诊断中智能软件模块的一个重要组成部分,也是系统完成故障监测与诊断的一个承上启下的重要环节,通过分割提取的电气设备特征可以为系统后期的电气设备目标提取与识别,智能判断和决策提供依据。
传统的二维Otsu阈值分割法虽然考虑了图像的局部空间信息,抑制了一部分噪声的干扰,分割效果也较为理想,但其没有考虑像素点与其周围领域像素灰度值的偏离情况,并且根据二维直方图粗略的将一小部分像素点看成边缘和噪声来处理,给图像的分割造成误差。
发明内容
本发明的目的在于提供一种基于改进的二维Otsu的变电站红外图像分割方法,能够有效地将目标区域完整的分割出来,为后续目标区域的特征提取奠定良好的基础。
为了实现上述目的,本发明采用的技术方案如下:
一种基于改进的二维Otsu的变电站红外图像分割方法,包括以下步骤:
包括以下步骤:
(1)输入变电站红外图像;
(2)采用形态学的权重自适应算法对输入的图像进行去噪:
(2a)根据输入的绝缘子串原始红外图像选择合适的形态学结构元素算子;
(2b)构建形态学滤波器并自适应确定其权值;
(2c)滤波处理,得到去噪后的图像;
(3)使用改进的二维Otsu算法进行图像分割:
(3a)求出被阈值对划分的两类图像像素对应的先验概率,进一步地得到相应的类内均值矢量;
(3b)求出二维直方图总均值矢量,假设二维直方图的副对角区域图像的概率可以忽略不计,于是可进一步得到一个关系式;
(3c)定义一个类间方差,结合步骤(3b)中的关系式,并以该类间方差的迹作为类间方差的测度,然后求得非改进的二维Otsu分割算法;
(3d)当步骤(3b)中的关系式并不成立时,须修正步骤(3c)中的类间方差测度关系式,进一步地得到改进后的二维Otsu分割算法,然后求得改进后的二维Otsu分割算法的最佳阈值。
具体地,所述步骤(2a)中选取的形态学结构元素算子为线性算子,每个线性算子具有不同的长度和角度。
具体地,所述步骤(2b)中选择构建串联滤波器来进行图像滤波,然后将上述不同形状的结构元素所构成的串联滤波器进行并联,结合自适应权值算法进一步构建串、并联复合滤波器。
具体地,将每次的串联处理结果与原始图像的差异值作为权值向量,分别记为β1,β2,…,βn,每次串联处理结果的权值分别记为α1,α2,…,αn,则自适应权值计算公式为:
具体地,每次串联滤波处理的结果为Ii,其中i=1,2,3,...,n,则输出去噪结果图像I为:
其中,αi为每次串联结果的权值,Ii为每次串联滤波的结果,式中i=1,2,3,...,n。
具体地,所述步骤(3a)中设图像灰度级从1到L,2维直方图任意一点的值定义为pij,pij表示二元组(i,j)发生的频率;假设图像被阈值对(s,t)划分为两类Co和Cb,对应的先验概率ωo(s,t)与ωb(s,t)分别为:
相应的类内均值矢量mo和mb分别为:
在公式(4)中,
具体地,所述步骤(3b)中二维直方图总均值矢量mT为:
根据假设得到的关系式为:
ωo+ωb≈1mT≈ωomo+ωbmb。(7)
具体地,所述步骤(3c)中类间方差σB的定义为:
σB=ωo[(mo-mT)(mo-mT)T]+ωb[(mb-mT)(mb-mT)T](8)
使用类间方差的迹作为类间方差的测度,有:
trσB=ωo[(moi-mTi)2+(moj-mTj)2]+ωb[(mbi-mTi)2+(mbj-mTj)2](9)将(7)带入(9)中,有:
此时二维Otsu分割算法的最佳阈值可以由(11)式确定:
具体地,所述步骤(3d)中对(10)式进行修正得到式(12):
此时改进的二维Otsu分割算法的最佳阈值可以由(13)式确
本发明与现有技术相比,具有以下优点及有益效果:
(1)本发明在图像分割前首先对图像进行去噪处理,缩小了像素点与图像领域像素灰度值的偏离程度,极大地减少了后续图像分割时的误差。
(2)本发明的分割阈值可以涵盖目标区域低灰度级部分,能够较好的将目标区域完整的分割出来,为后续目标区域的特征提取奠定良好的基础。
附图说明
图1为本发明仿真实验中待分割的断路器红外图像故障原图。
图2为本发明仿真实验中经过基于形态学的权重自适应去噪后的图像。
图3为本发明仿真实验中经过非改进的二维Otsu分割算法处理后的图像。
图4为本发明仿真实验中经过改进的二维Otsu分割算法处理后的图像。
具体实施方式
下面结合附图和实施例对本发明作进一步说明,本发明的实施方式包括但不限于下列实施例。
实施例
一种基于改进的二维Otsu的变电站红外图像分割方法,包括以下步骤:
(1)输入变电站红外图像;
(2)采用形态学的权重自适应算法对输入的图像进行去噪:
(2a)根据输入的绝缘子串原始红外图像选择合适的形态学结构元素算子:
(2b)构建形态学滤波器并自适应确定其权值;
(2c)滤波处理,得到去噪后的图像;
(3)使用改进的二维Otsu算法进行图像分割:
(3a)求出被阈值对划分的两类图像像素对应的先验概率,进一步得到相应的类内均值矢量;
(3b)求出二维直方图总均值矢量,假设二维直方图的副对角区域图像的概率可以忽略不计,于是可进一步得到一个关系式;
(3c)定义一个类间方差,结合步骤(3b)中的关系式,并以该类间方差的迹作为类间方差的测度,然后得到非改进的二维Otsu分割算法;
(3d)在很多情况下,步骤(3b)中的关系式并不成立,所以有必要修正步骤(3c)中的类间方差测度关系式,进一步得到改进后的二维Otsu分割算法,然后求得改进后的二维Otsu分割算法的最佳阈值。
所述步骤(2a)中选取的形态学结构元素算子为线性算子,每个线性算子具有不同的长度和角度。
所述步骤(2b)中形态学滤波器的选择方式为:先构建串联滤波器来进行图像滤波,然后将上述不同形状的结构元素所构成的串联滤波器进行并联,结合自适应权值算法进一步构建串、并联复合滤波器。
所述步骤(2b)中将每次的串联处理结果与原始图像的差异值作为权值向量,分别记为β1,β2,…,βn,每次串联处理结果的权值分别记为α1,α2,…,αn,则自适应权值计算公式为:
所述步骤(2c)中每次串联滤波的结果为Ii,其中i=1,2,3,...,n,则输出去噪结果图像I为:
所述步骤(3a)中设图像灰度级从1到L,2维直方图任意一点的值定义为pij,pij表示二元组(i,j)发生的频率。假设图像被阈值对(s,t)划分为两类Co和Cb,对应的先验概率ωo(s,t)与ωb(s,t)以及相应的类内均值矢量mo和mb分别为:
所述步骤(3b)中二维直方图总均值矢量mT为:
所述步骤(3b)中根据假设得到的关系式为:
ωo+ωb≈1mT≈ωomo+ωbmb。
所述步骤(3c)中类间方差σB的定义为:
σB=ωo[(mo-mT)(mo-mT)T]+ωb[(mb-mT)(mb-mT)T];
所述步骤(3c)中使用类间方差的迹作为类间方差的测度,有:
trσB=ωo[(moi-mTi)2+(moj-mTj)2]+ωb[(mbi-mTi)2+(mbj-mTj)2];
所述步骤(3c)中将步骤(3b)中根据假设得到的关系式带入类间方差的测度,有:
所述步骤(3d)中对非改进的二维Otsu分割算法公式进行修正得到改进的二维Otsu分割算法公式:
所述步骤(3d)中改进后的二维Otsu分割算法的最佳阈值计算公式为:
采用本发明的方法进行仿真实验。如图1所示,本发明仿真实验所要处理的红外图像故障原图,受到设备和拍摄的影响,该图像存在着许多噪声点。
如图2所示,对图像进行基于形态学的权重自适应去噪,去除这些噪声点的影响,这里选择的线性结构算子对图像进行形态学的腐蚀运算,线性结构算子长度分别为3、5、7,角度分别为-45、0、45.
如图3所示,采用非改进的二维Otus分割算法对图像进行分割,可见图像灰度级较低的部分没有分割出来,从而说明该算法不能将目标区域完整地分割出来。
如图4所示,采用改进的二维Otus分割算法对图像进行分割,可见该算法能够较好的将目标区域完整的分割出来,分割结果清晰并且没有错分现象。
按照上述实施例,便可很好地实现本发明。值得说明的是,基于上述结构设计的前提下,为解决同样的技术问题,即使在本发明上做出的一些无实质性的改动或润色,所采用的技术方案的实质仍然与本发明一样,故其也应当在本发明的保护范围内。
Claims (9)
1.一种基于改进的二维Otsu的变电站红外图像分割方法,其特征在于,包括以下步骤:
(1)输入变电站红外图像;
(2)采用形态学的权重自适应算法对输入的图像进行去噪:
(2a)根据输入的绝缘子串原始红外图像选择合适的形态学结构元素算子;
(2b)构建形态学滤波器并自适应确定其权值;
(2c)滤波处理,得到去噪后的图像;
(3)使用改进的二维Otsu算法进行图像分割:
(3a)求出被阈值对划分的两类图像像素对应的先验概率,进一步地得到相应的类内均值矢量;
(3b)求出二维直方图总均值矢量,假设二维直方图的副对角区域图像的概率可以忽略不计,于是可进一步得到一个关系式;
(3c)定义一个类间方差,结合步骤(3b)中的关系式,并以该类间方差的迹作为类间方差的测度,然后求得非改进的二维Otsu分割算法;
(3d)当步骤(3b)中的关系式并不成立时,须修正步骤(3c)中的类间方差测度关系式,进一步地得到改进后的二维Otsu分割算法,然后求得改进后的二维Otsu分割算法的最佳阈值。
2.根据权利要求1所述的基于改进的二维Otsu的变电站红外图像分割方法,其特征在于,所述步骤(2a)中选取的形态学结构元素算子为线性算子,每个线性算子具有不同的长度和角度。
3.根据权利要求1所述的基于改进的二维Otsu的变电站红外图像分割方法,其特征在于,所述步骤(2b)中选择构建串联滤波器来进行图像滤波,然后将上述不同形状的结构元素所构成的串联滤波器进行并联,结合自适应权值算法进一步构建串、并联复合滤波器。
4.根据权利要求3所述的基于改进的二维Otsu的变电站红外图像分割方法,其特征在于,将每次的串联处理结果与原始图像的差异值作为权值向量,分别记为β1,β2,…,βn,每次串联处理结果的权值分别记为α1,α2,…,αn,则自适应权值计算公式为:
。
5.根据权利要求4所述的基于改进的二维Otsu的变电站红外图像分割方法,其特征在于,每次串联滤波处理的结果为Ii,其中i=1,2,3,...,n,则输出去噪结果图像I为:
其中,αi为每次串联结果的权值,Ii为每次串联滤波的结果,式中i=1,2,3,...,n。
6.根据权利要求1所述的基于改进的二维Otsu的变电站红外图像分割方法,其特征在于,所述步骤(3a)中设图像灰度级从1到L,2维直方图任意一点的值定义为pij,pij表示二元组(i,j)发生的频率;假设图像被阈值对(s,t)划分为两类Co和Cb,对应的先验概率ωo(s,t)与ωb(s,t)分别为:
相应的类内均值矢量mo和mb分别为:
在公式(4)中,
7.根据权利要求1所述的基于改进的二维Otsu的变电站红外图像分割方法,其特征在于,所述步骤(3b)中二维直方图总均值矢量mT为:
根据假设得到的关系式为:
ωo+ωb≈1mT≈ωomo+ωbmb。(7)。
8.根据权利要求7所述的基于改进的二维Otsu的变电站红外图像分割方法,其特征在于,所述步骤(3c)中类间方差σB的定义为:
σB=ωo[(mo-mT)(mo-mT)T]+ωb[(mb-mT)(mb-mT)T](8)
使用类间方差的迹作为类间方差的测度,有:
trσB=ωo[(moi-mTi)2+(moj-mTj)2]+ωb[(mbi-mTi)2+(mbj-mTj)2](9)
将(7)带入(9)中,有:
此时二维Otsu分割算法的最佳阈值可以由(11)式确定:
9.根据权利要求8所述的基于改进的二维Otsu的变电站红外图像分割方法,其特征在于,所述步骤(3d)中对(10)式进行修正得到式(12):
此时改进的二维Otsu分割算法的最佳阈值可以由(13)式确
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510698451.1A CN105205829B (zh) | 2015-10-23 | 2015-10-23 | 基于改进的二维Otsu的变电站红外图像分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510698451.1A CN105205829B (zh) | 2015-10-23 | 2015-10-23 | 基于改进的二维Otsu的变电站红外图像分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105205829A true CN105205829A (zh) | 2015-12-30 |
CN105205829B CN105205829B (zh) | 2019-02-19 |
Family
ID=54953490
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510698451.1A Active CN105205829B (zh) | 2015-10-23 | 2015-10-23 | 基于改进的二维Otsu的变电站红外图像分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105205829B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107292902A (zh) * | 2017-07-07 | 2017-10-24 | 国家电网公司 | 一种结合果蝇优化算法的二维Otsu图像分割方法 |
CN108062508A (zh) * | 2017-10-13 | 2018-05-22 | 西安科技大学 | 变电站复杂背景红外图像中设备的提取方法 |
CN108319966A (zh) * | 2017-10-13 | 2018-07-24 | 西安科技大学 | 一种变电站复杂背景红外图像中设备的识别分类方法 |
CN109461148A (zh) * | 2018-10-30 | 2019-03-12 | 兰州交通大学 | 基于二维Otsu的钢轨缺陷分割自适应快速算法 |
-
2015
- 2015-10-23 CN CN201510698451.1A patent/CN105205829B/zh active Active
Non-Patent Citations (4)
Title |
---|
景晓军等: "一种基于二维最大类间方差的图像分割算法", 《通信学报》 * |
薛卫鹏: "基于Otsu算法与数学形态学的图像分割算法研究", 《中国优秀硕士学位论文全文数据库信息科技辑(月刊)》 * |
陈虎等: "基于数字形态学的图像去噪方法研究", 《工程图学学报》 * |
黄宵宁等: "直升机巡检航拍图像中绝缘子图像的提取算法", 《电网技术》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107292902A (zh) * | 2017-07-07 | 2017-10-24 | 国家电网公司 | 一种结合果蝇优化算法的二维Otsu图像分割方法 |
CN107292902B (zh) * | 2017-07-07 | 2020-05-19 | 国家电网公司 | 一种结合果蝇优化算法的二维Otsu图像分割方法 |
CN108062508A (zh) * | 2017-10-13 | 2018-05-22 | 西安科技大学 | 变电站复杂背景红外图像中设备的提取方法 |
CN108319966A (zh) * | 2017-10-13 | 2018-07-24 | 西安科技大学 | 一种变电站复杂背景红外图像中设备的识别分类方法 |
CN108319966B (zh) * | 2017-10-13 | 2019-08-09 | 西安科技大学 | 一种变电站复杂背景红外图像中设备的识别分类方法 |
CN109461148A (zh) * | 2018-10-30 | 2019-03-12 | 兰州交通大学 | 基于二维Otsu的钢轨缺陷分割自适应快速算法 |
Also Published As
Publication number | Publication date |
---|---|
CN105205829B (zh) | 2019-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108388896B (zh) | 一种基于动态时序卷积神经网络的车牌识别方法 | |
Yamaguchi et al. | Fast crack detection method for large-size concrete surface images using percolation-based image processing | |
CN102831582B (zh) | 一种微软体感装置深度图像增强方法 | |
CN102541954B (zh) | 一种商标检索方法及系统 | |
CN102682428B (zh) | 一种基于方向场的指纹图像计算机自动化修补方法 | |
CN109919002B (zh) | 黄色禁停线识别方法、装置、计算机设备及存储介质 | |
CN105205829A (zh) | 基于改进的二维Otsu的变电站红外图像分割方法 | |
CN104463870A (zh) | 一种图像显著区域检测方法 | |
CN110660065B (zh) | 一种红外故障检测识别算法 | |
CN104282027B (zh) | 一种基于Hough变换的圆检测方法 | |
CN112308872B (zh) | 基于多尺度Gabor一阶导数的图像边缘检测方法 | |
CN103955682A (zh) | 基于surf兴趣点的行为识别方法及装置 | |
CN103903238A (zh) | 图像特征的显著结构和相关结构融合方法 | |
CN104766319A (zh) | 提升夜间拍照图像配准精度的方法 | |
CN103150725B (zh) | 基于非局部均值的susan边缘检测方法及系统 | |
CN104077775A (zh) | 一种结合骨架特征点和形状上下文的形状匹配方法及装置 | |
CN104966054A (zh) | 无人机可见光图像中的弱小目标检测方法 | |
CN111696072A (zh) | 彩色图像的直线检测方法、装置、电子设备及存储介质 | |
CN103839248B (zh) | 图像分割方法和装置 | |
CN116311201A (zh) | 一种基于图像识别技术的变电站设备状态识别方法及系统 | |
CN113205494B (zh) | 基于自适应尺度图像块加权差测量的红外小目标检测方法和系统 | |
CN109636822A (zh) | 一种基于新构建隶属度函数的改进Canny自适应边缘提取方法 | |
CN104966296A (zh) | 滑窗N-Smoothlets图像边缘检测方法 | |
Zhao | Motion track enhancement method of sports video image based on otsu algorithm | |
CN108764110B (zh) | 基于hog特征行人检测器的递归错检校验方法、系统及设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |