CN105608457A - 直方图灰度矩阈值分割法 - Google Patents
直方图灰度矩阈值分割法 Download PDFInfo
- Publication number
- CN105608457A CN105608457A CN201510949163.9A CN201510949163A CN105608457A CN 105608457 A CN105608457 A CN 105608457A CN 201510949163 A CN201510949163 A CN 201510949163A CN 105608457 A CN105608457 A CN 105608457A
- Authority
- CN
- China
- Prior art keywords
- histogram
- gray
- image
- hist
- infrared 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
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/28—Quantising the image, e.g. histogram thresholding for discrimination between background and foreground patterns
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
- Image Input (AREA)
Abstract
为了对红外目标图像进行准确分割并实现对其识别与跟踪,本发明提出了一种直方图灰度矩阈值分割法。该方法借鉴力学杠杆原理,引入直方图灰度矩的概念,通过对直方图灰度区间进行遍历来查找使得直方图灰度矩平衡的阈值,以该阈值为最佳阈值实现对红外图像的二值化分割处理。该方法与经典的最大类间方差法(也称Otsu法)具有相同的算法时间复杂度,但对于红外图像,本发明所提方法具有较优的分割效果。
Description
技术领域
本发明涉及一种适用于红外目标识别与跟踪前预处理的直方图灰度矩阈值分割法,属于图像处理技术领域。
背景技术
图像分割常常作为图像特征分析、目标识别等技术的首要预处理操作,被广泛应用于科学研究和工程领域。在各类图像分割方法中,阈值分割法因其简单稳定、计算量小成为一种常用的分割方法。该方法依据图像区域相似性,一般将图像分割为前景区域和背景区域,其中前景区域常常是后续图像处理的感兴趣目标区域。
常用的阈值分割法有基于一维直方图的迭代选择阈值法、最小均方误差法、最大类间方差法以及基于二维直方图的二维最大熵法等。其中最大类间方差法(也称Otsu法)是由Otsu在文献一(OTSUN.Athresholdselectionmethodfromgraylevelhistogram[J].IEEETransactionsonSystems,ManandCybernetics,1979,9(1):62-66.)提出。Otsu法因其计算简单、性能稳定且分割效果较好成为了最受欢迎的阈值分割法之一,但是该方法对于图像直方图双峰不明显或者出现多峰情况时分割效果欠佳。因此之后许多学者都是以Otsu法为基础进行算法的改进。另外由Kapur等人在文献二(KapurJN,SahooPK,WongAKC.Anewmethodforgray-levelpicturethresholdingusingtheentropyofthehistogram[J].ComputerVision,GraphicsandImageProcessing,1985,29(3):273-285.)提出的最大熵法以及由Abutaleb在文献三(AbutalebAS.Automaticthresholdingofgray-levelpicturesusingtwo-dimensionentropy[J].ComputerVision,GraphicsandImageProcessing,1989,47(1):22-32.)将一维最大熵法推广至二维后,图像分割的效果得到了很大的改善,然而该方法计算量大。正如杜峰等人在文献四(《一种快速红外图像分割方法》文章编号:1001-9014(2005)05-0370-04)为了降低二维最大熵算法的计算量又引入了粒子群寻优策略,虽然该方法既对红外图像实现了较好的分割效果,又一定程度上降低了算法运算量,但是分割结果不稳定,而且在硬件平台上实现其计算量仍是一个需要解决的问题。
在红外目标识别与跟踪领域,红外图像的目标区域一般具有尺寸较小、灰度分布较窄、图像直方图双峰并不明显的特点。为了快速、稳定且最大可能地将目标区域从背景中分割出来,此时的Otsu法或最大熵法仍无法很好地满足该要求。
发明内容
针对上述问题,本发明的目的是提供一种直方图灰度矩阈值分割法,该方法将力学杠杆的平衡思想应用于红外图像的灰度直方图中,定义了直方图灰度矩的概念。通过对图像的灰度直方图进行逐级遍历来查找使得直方图杠杆平衡的最佳阈值,然后对图像进行二值化处理,从而实现了对红外目标图像的阈值分割。
本发明通过以下方案实现:
一种直方图灰度矩阈值分割法,其特征在于,包括如下步骤:
(1)图像输入:读入一幅原始的8bits红外图像I;
(2)统计:统计红外图像I的灰度直方图Hist(i),i表示红外图像I的灰度变化范围,i=0,1,2…255,Hist(i)表示灰度值为i的像素个数;
(3)计算红外图像I的直方图灰度矩L(T):
(a).将步骤(2)建立的灰度直方图Hist(i)模拟为一个杠杆,其中T表示直方图杠杆的支点,T为整数且取值范围为;
(b).定义红外图像I灰度值为i的直方图灰度矩,如下式:
;(1)
那么直方图灰度矩总和为:
(2)
其中在直方图支点的左半部分称为正矩,右半部分称为负矩,分别定义如下:
(3)
将公式(1)代入式(2)中,红外图像I的直方图灰度矩总和如下式:
(4)
其中。
对红外目标图像I的灰度直方图Hist(i)作对数变换后,直方图灰度矩最终公式定义如下:
(5)
(c).按照公式(5)依次计算不同支点T下的直方图灰度矩L(T);
(d).从T=0到T=255依序遍历L(T),直到发现L(T)大于或等于0为止,此时的T即为直方图灰度矩的平衡点;
(4)二值化分割:将直方图灰度矩L(T)的平衡点作为最佳分割阈值,对红外图像I进行二值化分割,从而完成红外图像I的阈值分割。
所述将步骤(2)建立的灰度直方图Hist(i)模拟为一个杠杆的具体做法是基于如下定义:
(1)定义任意一个灰度值为i的像素具有大小为i的重量,Hist(i)则表示该8bits红外图像I中重量为i的像素总数。
(2)构造矩形条:在直方图Hist(i)中,将灰度值为i的像素统计定义成一根具有重量大小为i×Hist(i)的矩形条,该矩形条的高度为Hist(i),宽度为1。
(3)构造直方图杠杆:对于8bits红外图像I而言,其直方图共有256根矩形条,将这些矩形条按图像灰度值从小到大排列在一根不计重量的水平轻质杠杆上,整个结构就组成了直方图杠杆。
分析公式(5)可知,直方图灰度矩与成正比例线性关系,对于一幅红外图像而言,理论上总存在一个使其直方图杠杆平衡的支点,使得直方图灰度矩,此时也是与关系曲线的过零点。而在实际应用中由于图像直方图的离散化,支点的取值只能为整数,此时直方图灰度矩在平衡点处并非严格等于0,即在处恰恰大于或等于0。对于小于时,大于时,总之在远离点两边直方图杠杆均会失衡。
对于如何评价图像分割的效果还没有一个统一的定量标准,一个广泛的评估原则就是看是否分割出感兴趣或期望的目标区域。因此评价图像分割的效果往往具有一定的主观判断因素。在红外目标识别与跟踪领域,如何能实现快速、稳定、自适应且最大可能地将目标区域从背景中分割出来才是解决后续图像分析的关键。利用公式(5)对红外目标图像的灰度直方图进行遍历搜索求取使得直方图灰度矩平衡的最佳阈值,然后对图像进行二值化处理,实现方法如下:
(6)
其中为原始红外图像,为阈值分割后的二值化图,i,j为图像像素索引下标,即为直方图灰度矩平衡点,也就是最佳分割阈值。
本发明的优点:相比现有技术而言,本发明通过对红外图像的灰度直方图进行逐级遍历来查找使得直方图杠杆平衡的最佳阈值,然后对图像进行二值化处理,从而实现了对红外目标图像的阈值分割,具有性能稳定、实现容易、计算量小且分割效果较Otsu法好的优势,能较好地提取红外图像的目标区域,既方便在MATLAB平台仿真,也容易实现嵌入式编程。
附图说明
图1为本发明所述灰度直方图的模拟杠杆示意图;
图2为第一幅原始红外图像;
图3为将本发明方法公式(5)应用于图2的L(T)与T的关系曲线图,其中横坐标标示出了使直方图灰度矩平衡的最佳阈值;
图4为现有Otsu法对图2阈值分割的二值化图,最佳阈值为77;
图5为本发明方法对图2阈值分割的二值化图,最佳阈值为101;
图6为第二幅原始红外图像;
图7为将本发明方法公式(5)应用于图6的L(T)与T的关系曲线图,其中横坐标标示出了使直方图灰度矩平衡的最佳阈值;
图8为现有Otsu法对图6阈值分割的二值化图,最佳阈值为79;
图9为本发明方法对图6阈值分割的二值化图,最佳阈值为103。
具体实施方式
现结合具体实施例及附图,来进一步阐述本发明。
实施例一
(1)图像输入:读入一幅分辨率为240×180的8bits的原始红外图像I,如图2或图6;
(2)统计:计算该红外图像I的行数row,列数col,并统计红外图像I的灰度直方图Hist(i),i表示红外图像I的灰度变化范围,i=0,1,2…255,Hist(i)表示灰度值为i的像素个数;
(3)计算红外图像I的直方图灰度矩:
(a).将步骤(2)建立的灰度直方图Hist(i)按图像灰度值从小到大排列模拟成一根直方图杠杆,如图1所示。T表示直方图杠杆的支点,T为整数且取值范围为;
(b).按照公式依次计算不同支点T下的直方图灰度矩L(T);
(c).从T=0到T=255依序遍历L(T),直到发现L(T)大于或等于0为止,此时的T即为直方图灰度矩的平衡点;
(4)二值化分割:将直方图灰度矩L(T)的平衡点作为最佳分割阈值,然后对红外图像I进行二值化分割,实现方法如下:
其中为原始红外图像,为阈值分割后的二值化图,i,j为图像像素索引下标,即为直方图灰度矩平衡点,也就是最佳分割阈值,从而完成红外图像I的阈值分割。
采用本发明的方法和现有Otsu法对图2和图6所示的两幅原始红外图像进行分割处理,两种方法处理的效果以二值化图所示,其中图4为现有Otsu法对图2进行阈值分割的二值化图(最佳阈值为77),图5为本发明方法对图2进行阈值分割的二值化图(最佳阈值为101);图8为现有Otsu法对图6进行阈值分割的二值化图(最佳阈值为79);图9为本发明方法对图6进行阈值分割的二值化图(最佳阈值为103)。由两种方法对两幅原始红外图像进行阈值分割的效果对比来看,本发明与经典的最大类间方差法(也称Otsu法)具有相同的算法时间复杂度,但对于红外图像,本发明所提方法具有较优的分割效果。
Claims (2)
1.直方图灰度矩阈值分割法,其特征在于,包括如下步骤:
(1)图像输入:读入一幅原始的8bits红外图像I;
(2)统计:统计红外图像I的灰度直方图Hist(i),i表示红外图像I的灰度变化范围,i=0,1,2…255,Hist(i)表示灰度值为i的像素个数;
(3)计算红外图像I的直方图灰度矩L(T):
(a).将步骤(2)建立的灰度直方图Hist(i)模拟为一个杠杆,其中T表示直方图杠杆的支点,T为整数且取值范围为;
(b).定义红外图像I灰度值为i的直方图灰度矩,如下式:
;
定义红外图像I的直方图灰度矩总和,如下式:
;
对红外目标图像I的灰度直方图Hist(i)作对数变换后,直方图灰度矩最终公式定义如下:
;
(c).按照公式依次计算不同支点T下的直方图灰度矩L(T);
(d).从T=0到T=255依序遍历L(T),直到发现L(T)大于或等于0为止,此时的T即为直方图灰度矩的平衡点;
(4)二值化分割:将直方图灰度矩L(T)的平衡点作为最佳分割阈值,对红外图像I进行二值化分割,从而完成红外图像I的阈值分割。
2.如权利要求1所述的直方图灰度矩阈值分割法,其特征在于,所述将步骤(2)建立的灰度直方图Hist(i)模拟为一个杠杆的具体做法是基于如下定义:
(1)定义任意一个灰度值为i的像素具有大小为i的重量,Hist(i)则表示该8bits红外图像I中重量为i的像素总数;
(2)构造矩形条:在直方图Hist(i)中,将灰度值为i的像素统计定义为一根具有重量大小为i×Hist(i)的矩形条,该矩形条的高度为Hist(i),宽度为1;
(3)构造直方图杠杆:对于8bits红外图像I而言,其直方图共有256根矩形条,将这些矩形条按图像灰度值从小到大排列在一根不计重量的水平轻质杠杆上,整个结构就组成了直方图杠杆。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510949163.9A CN105608457B (zh) | 2015-12-18 | 2015-12-18 | 直方图灰度矩阈值分割法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510949163.9A CN105608457B (zh) | 2015-12-18 | 2015-12-18 | 直方图灰度矩阈值分割法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105608457A true CN105608457A (zh) | 2016-05-25 |
CN105608457B CN105608457B (zh) | 2018-12-28 |
Family
ID=55988383
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510949163.9A Active CN105608457B (zh) | 2015-12-18 | 2015-12-18 | 直方图灰度矩阈值分割法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105608457B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106803238A (zh) * | 2016-12-26 | 2017-06-06 | 杭州施强教育科技有限公司 | 一种答题卡图像降噪处理方法 |
CN106815814A (zh) * | 2016-12-26 | 2017-06-09 | 杭州施强教育科技有限公司 | 一种应用于阅卷系统的图像污染处理方法 |
CN106845468A (zh) * | 2016-12-26 | 2017-06-13 | 杭州施强教育科技有限公司 | 一种提高答题卡图像识别正确率的处理方法 |
CN107194910A (zh) * | 2017-04-02 | 2017-09-22 | 南京汇川图像视觉技术有限公司 | 一种基于混合阈值的提环检测方法 |
CN109816627A (zh) * | 2018-12-13 | 2019-05-28 | 中国科学院西安光学精密机械研究所 | 平面玻璃元件油墨区弱小缺陷目标检测方法 |
CN109850518A (zh) * | 2018-11-12 | 2019-06-07 | 太原理工大学 | 一种基于红外图像的实时矿用胶带预警撕裂检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005072205A2 (en) * | 2004-01-15 | 2005-08-11 | Chemicon International, Inc. | Image analysis and assay system |
CN101059870A (zh) * | 2007-05-23 | 2007-10-24 | 华中科技大学 | 一种基于属性直方图的图像分割方法 |
CN101710425A (zh) * | 2009-12-25 | 2010-05-19 | 南京航空航天大学 | 基于图像灰度梯度和灰度统计直方图的自适应预分割方法 |
CN104732519A (zh) * | 2015-01-20 | 2015-06-24 | 中国科学院半导体研究所 | 鲁棒的全局阈值分割方法 |
-
2015
- 2015-12-18 CN CN201510949163.9A patent/CN105608457B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005072205A2 (en) * | 2004-01-15 | 2005-08-11 | Chemicon International, Inc. | Image analysis and assay system |
CN101059870A (zh) * | 2007-05-23 | 2007-10-24 | 华中科技大学 | 一种基于属性直方图的图像分割方法 |
CN101710425A (zh) * | 2009-12-25 | 2010-05-19 | 南京航空航天大学 | 基于图像灰度梯度和灰度统计直方图的自适应预分割方法 |
CN104732519A (zh) * | 2015-01-20 | 2015-06-24 | 中国科学院半导体研究所 | 鲁棒的全局阈值分割方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106803238A (zh) * | 2016-12-26 | 2017-06-06 | 杭州施强教育科技有限公司 | 一种答题卡图像降噪处理方法 |
CN106815814A (zh) * | 2016-12-26 | 2017-06-09 | 杭州施强教育科技有限公司 | 一种应用于阅卷系统的图像污染处理方法 |
CN106845468A (zh) * | 2016-12-26 | 2017-06-13 | 杭州施强教育科技有限公司 | 一种提高答题卡图像识别正确率的处理方法 |
CN107194910A (zh) * | 2017-04-02 | 2017-09-22 | 南京汇川图像视觉技术有限公司 | 一种基于混合阈值的提环检测方法 |
CN107194910B (zh) * | 2017-04-02 | 2020-06-09 | 南京汇川图像视觉技术有限公司 | 一种基于混合阈值的提环检测方法 |
CN109850518A (zh) * | 2018-11-12 | 2019-06-07 | 太原理工大学 | 一种基于红外图像的实时矿用胶带预警撕裂检测方法 |
CN109850518B (zh) * | 2018-11-12 | 2022-01-28 | 太原理工大学 | 一种基于红外图像的实时矿用胶带预警撕裂检测方法 |
CN109816627A (zh) * | 2018-12-13 | 2019-05-28 | 中国科学院西安光学精密机械研究所 | 平面玻璃元件油墨区弱小缺陷目标检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105608457B (zh) | 2018-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105608457A (zh) | 直方图灰度矩阈值分割法 | |
Yi et al. | An improved tiny-yolov3 pedestrian detection algorithm | |
CN106547880B (zh) | 一种融合地理区域知识的多维度地理场景识别方法 | |
Jiang et al. | Self-supervised relative depth learning for urban scene understanding | |
CN111126202A (zh) | 基于空洞特征金字塔网络的光学遥感图像目标检测方法 | |
CN108665485B (zh) | 一种基于相关滤波与孪生卷积网络融合的目标跟踪方法 | |
CN110349103A (zh) | 一种基于深度神经网络和跳跃连接的无干净标签图像去噪方法 | |
CN107103326A (zh) | 基于超像素聚类的协同显著性检测方法 | |
CN113487629B (zh) | 一种基于结构化场景和文本描述的图像属性编辑方法 | |
CN111080675A (zh) | 一种基于时空约束相关滤波的目标跟踪方法 | |
CN106055653A (zh) | 基于图像语义注释的视频浓缩对象检索方法 | |
Wang et al. | Hybrid attention network based on progressive embedding scale-context for crowd counting | |
CN111652240B (zh) | 一种基于cnn的图像局部特征检测与描述方法 | |
CN111582091B (zh) | 基于多分支卷积神经网络的行人识别方法 | |
Li et al. | Robust deep neural networks for road extraction from remote sensing images | |
CN114092517A (zh) | 基于传统和深度学习算法的多目标跟踪方法 | |
Cai et al. | Identification of grape leaf diseases based on VN-BWT and Siamese DWOAM-DRNet | |
CN106504219B (zh) | 有约束的路径形态学高分辨率遥感影像道路增强方法 | |
CN111008630A (zh) | 一种基于弱监督学习的目标定位方法 | |
Feng et al. | Local complexity difference matting based on weight map and alpha mattes | |
Zhong et al. | SCPNet: Self-constrained parallelism network for keypoint-based lightweight object detection | |
CN116433909A (zh) | 基于相似度加权多教师网络模型的半监督图像语义分割方法 | |
Wang et al. | Feature enhancement: predict more detailed and crisper edges | |
Tabib et al. | Categorization and selection of crowdsourced images towards 3d reconstruction of heritage sites | |
Zhao | Application of a clustering algorithm in sports video image extraction and processing |
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 |