CN104050472B - 一种灰度图像二值化的自适应全局阈值方法 - Google Patents

一种灰度图像二值化的自适应全局阈值方法 Download PDF

Info

Publication number
CN104050472B
CN104050472B CN201410258640.2A CN201410258640A CN104050472B CN 104050472 B CN104050472 B CN 104050472B CN 201410258640 A CN201410258640 A CN 201410258640A CN 104050472 B CN104050472 B CN 104050472B
Authority
CN
China
Prior art keywords
image
histogram
binary
cluster
curve
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.)
Expired - Fee Related
Application number
CN201410258640.2A
Other languages
English (en)
Other versions
CN104050472A (zh
Inventor
陆成刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang Aerospace Changfeng Science & Technology Development Co Ltd
Original Assignee
Zhejiang University of Technology ZJUT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201410258640.2A priority Critical patent/CN104050472B/zh
Publication of CN104050472A publication Critical patent/CN104050472A/zh
Application granted granted Critical
Publication of CN104050472B publication Critical patent/CN104050472B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Analysis (AREA)
  • Collating Specific Patterns (AREA)

Abstract

一种灰度图像二值化的自适应全局阈值方法,首先对图像进行高斯低通滤波,以便消除高频噪声干扰;计算图像的直方图;归一化直方图,得到归一化的直方图曲线,曲线上点的数目为256;对这256个点进行二元均值聚类;聚类迭代收敛后,输出两个聚类中心的横坐标之间的曲线部分中纵坐标数值最低的点的横坐标值乘以255倍;根据得到的阈值再做二值化处理。本发明从图像灰度直方图的分布出发,模拟人工肉眼判别阈值的经验,提供了一个基于归一化直方图的二维点集聚类的全局阈值选取方法,能够有效地进行二值化处理、计算量较小。

Description

一种灰度图像二值化的自适应全局阈值方法
技术领域
本发明属于图像处理技术领域,具体地说是涉及灰度图像的二值化处理方法。
背景技术
灰度图像的二值化处理旨在生成黑白图,有效地提取前景目标、剔除背景干扰,为进一步提取目标的边缘轮廓、乃至进行有效地为目标识别匹配提供必要的基础。灰度图的二值化处理在基于图像的模式识别应用里无处不在,车牌识别、光学字符识别(OCR)、人脸识别、指纹识别、和掌纹识别等等都需要用到图像二值化处理技术。
灰度图的二值化处理关键在于全局阈值的选取,自适应全局阈值的选取方法就是对于各种不同像素值分布的图像都能够自动找出合适的阈值进行二值化。获得阈值之后的图像二值化处理过程如(其中x,y是图像坐标;g(x,y)是图像灰度值; Tr是二值化阈值)。自适应二值化全局阈值方法一直是学术和技术领域研究的难点、热点。
在图像处理的经典理论里,二值化自适应阈值法有双峰法、P参数法、大津法(OTSU)、最大熵阈值法、和迭代法。其中双峰法对图像灰度呈准多峰分布时确定波谷有难度、P参数法也需要依赖目标占比的先验知识、大津法对背景和目标比较接近的图像处理失效、最大熵阈值法和迭代法的缺陷是计算复杂度略大。
发明内容
为解决现有技术中计算复杂、效果差的问题,提供一种灰度图像二值化的自适应全局阈值方法,这种方法能有效地进行二值化处理、计算简便。
为实现上述目的,本发明一种灰度图像二值化的自适应全局阈值方法,采用以下步骤:
第一步、首先对图像进行高斯低通滤波,以便消除高频噪声干扰;
第二步、计算图像的直方图,其中i是灰度级,h(i)是灰度级为i的图像像素的数目;
第三步、归一化直方图,得到,归一化的直方图曲线由这些二维点集组成,曲线上点的数目为256;
第四步、对上述256个点进行二元均值聚类,聚类的两个点集合P0和P1的初始中心各为,其中,即取从曲线两端开始的第一个纵坐标非零的点做初始聚类中心;
第五步、聚类迭代收敛后,输出两个聚类中心的横坐标之间的曲线部分中纵坐标数值最低的点的横坐标值乘以255倍;
第六步、根据步骤5得到的阈值再做二值化处理。
步骤三需要对直方图做归一化处理的理由是消除直方图纵坐标、横坐标由于物理意义不一样引发的差异性使得直方图曲线在聚类时依照欧式距离来计算。
步骤四采用固定机制的初始中心选取法,避免常用的聚类法使用随机初值中心造成结果的不稳定。
步骤五由于总的点数为256,经过较少的计算,聚类迭代即可收敛。
本发明的计算量少于迭代法、最大熵阈值法等在全图像像素数目上的迭代计算。
本发明从图像灰度直方图的分布出发,模拟人工肉眼判别阈值的经验,提出了一个基于归一化直方图的二维点集聚类的全局阈值选取方法,能够有效地进行二值化处理、计算复杂程度较小。
附图说明
图1为原始图像。
图2为经过高斯低通滤波后的结果图像。
图3为经本发明方法计算处理后的二值化结果。
图4为是原始直方图。
图5为归一化直方图。
具体实施方式
实施例一
以掌纹图预处理的二值化过程为例。
在学术文献里掌纹图的二值化都是根据经验设定一个固定阈值的方式处理的,在实际应用中会带来适应性的问题。
通过本发明步骤对掌纹图的二值化处理:
第一步、首先对图像进行高斯低通滤波,以便消除高频噪声干扰;
第二步、计算图像的直方图
其中i是灰度级,h(i)是灰度级为i的图像像素的数目;
本例中i为,i范围是0<=i<=255,
i是一个变量,h(i)表征的是一个函数
第三步、归一化直方图,得到,归一化的直方图曲线由这些二维点集组成,曲线上点的数目为256个;
第四步、对上述256个点进行二元均值聚类,聚类的两个点集合P0和P1的初始中心各为
其中
均值聚类算法收敛后的聚类中心坐标为(0.256863, 0.274890)和(0.758824,0.416831)
第五步、聚类迭代收敛后,输出两个聚类中心的横坐标之间的曲线部分中纵坐标数值最低的点的横坐标值乘以255倍;
本例中,经计算得到阈值为66
第六步、根据步骤5得到的阈值再做二值化处理。
通过本发明的步骤进行计算,对掌纹可以清晰地获取其边缘,二值化效果明显。

Claims (1)

1.一种灰度图像二值化的自适应全局阈值方法,其特征在于包括以下步骤:
第一步、首先对图像进行高斯低通滤波,用于消除高频噪声干扰;
第二步、计算图像的直方图
其中i是灰度级,h(i)是灰度级为i的图像像素的数目;
第三步、归一化直方图,得到归一化的直方图曲线由这些二维点集组成曲线上点的数目为256;
第四步、对上述256个点进行二元均值聚类,聚类的两个点集合P0和P1的初始中心各为(i0/255,h(i0)/Max0<=i<=255(h(i))),和(i1/255,h(i1)/Max0<=i<=255(h(i))),
其中
第五步、聚类迭代收敛后,输出两个聚类中心的横坐标之间的曲线部分中纵坐标数值最低的点的横坐标值乘以255倍;
第六步、根据第五步得到的阈值再做二值化处理。
CN201410258640.2A 2014-06-12 2014-06-12 一种灰度图像二值化的自适应全局阈值方法 Expired - Fee Related CN104050472B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410258640.2A CN104050472B (zh) 2014-06-12 2014-06-12 一种灰度图像二值化的自适应全局阈值方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410258640.2A CN104050472B (zh) 2014-06-12 2014-06-12 一种灰度图像二值化的自适应全局阈值方法

Publications (2)

Publication Number Publication Date
CN104050472A CN104050472A (zh) 2014-09-17
CN104050472B true CN104050472B (zh) 2018-02-27

Family

ID=51503284

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410258640.2A Expired - Fee Related CN104050472B (zh) 2014-06-12 2014-06-12 一种灰度图像二值化的自适应全局阈值方法

Country Status (1)

Country Link
CN (1) CN104050472B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104573692B (zh) * 2014-12-18 2020-04-10 厦门翼歌软件科技有限公司 一种基于模糊退化模型的车牌二值化方法
CN107561962A (zh) * 2016-06-30 2018-01-09 浙江大丰实业股份有限公司 智能化剧院设备控制平台
CN106529543B (zh) * 2016-11-02 2018-09-11 徐庆 一种动态计算多色级二值化自适应阈值的方法及其系统
CN109035281A (zh) * 2018-07-12 2018-12-18 南京工程学院 一种基于直方图分布的图像三值化方法
CN115410249A (zh) * 2022-07-29 2022-11-29 成都云天励飞技术有限公司 人脸识别模型训练方法、识别方法、装置、设备及介质
CN116092015B (zh) * 2023-04-06 2023-06-27 安徽乾劲企业管理有限公司 一种道路施工状态监测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101238987A (zh) * 2007-09-06 2008-08-13 深圳先进技术研究院 一种关于ct脑出血图像的处理方法
CN101447027A (zh) * 2008-12-25 2009-06-03 东莞市微模式软件有限公司 磁码字符区域的二值化方法及其应用
CN101923645A (zh) * 2009-06-09 2010-12-22 黑龙江大学 适应复杂应用场景中低质量虹膜图像的虹膜分割方法
CN103646392A (zh) * 2013-11-21 2014-03-19 华为技术有限公司 逆光检测方法及设备

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8781187B2 (en) * 2011-07-13 2014-07-15 Mckesson Financial Holdings Methods, apparatuses, and computer program products for identifying a region of interest within a mammogram image

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101238987A (zh) * 2007-09-06 2008-08-13 深圳先进技术研究院 一种关于ct脑出血图像的处理方法
CN101447027A (zh) * 2008-12-25 2009-06-03 东莞市微模式软件有限公司 磁码字符区域的二值化方法及其应用
CN101923645A (zh) * 2009-06-09 2010-12-22 黑龙江大学 适应复杂应用场景中低质量虹膜图像的虹膜分割方法
CN103646392A (zh) * 2013-11-21 2014-03-19 华为技术有限公司 逆光检测方法及设备

Also Published As

Publication number Publication date
CN104050472A (zh) 2014-09-17

Similar Documents

Publication Publication Date Title
CN104050472B (zh) 一种灰度图像二值化的自适应全局阈值方法
US9251614B1 (en) Background removal for document images
TW201740316A (zh) 圖像文字的識別方法和裝置
CN102306384B (zh) 一种基于单幅图像的颜色恒常性处理方法
CN106709964B (zh) 基于梯度修正和多向纹理提取的素描画生成方法及装置
CN108564549A (zh) 一种基于多尺度稠密连接网络的图像去雾方法
CN104463814B (zh) 基于局部纹理方向性的图像增强方法
Yang et al. Research and analysis of Image edge detection algorithm Based on the MATLAB
CN109472792A (zh) 结合局部熵的局部能量泛函与非凸正则项的图像分割方法
CN106127193B (zh) 一种人脸图像识别方法
CN102930538B (zh) 高斯势与空间直方图融合的自适应图像分割方法
CN106600610B (zh) 一种fcm图像分割方法及装置
CN109033987B (zh) 一种人脸图像阴阳脸的处理方法及系统
Alenezi Image Dehazing Based on Pixel Guided CNN with PAM via Graph Cut.
Lee et al. Multiscale morphology based illumination normalization with enhanced local textures for face recognition
JP2015204023A (ja) 被写体検出装置、被写体検出方法及びプログラム
CN110874825A (zh) 一种复合绝缘子表面水迹二值图像提取方法
Liu An improved image enhancement algorithm based on fuzzy set
Das et al. Histogram equalization techniques for contrast enhancement: a review
CN104063874B (zh) 基于权值模型的灰度图像二值化算法定量评估方法
El-Zaart et al. SAR images thresholding for oil spill detection
CN107784626B (zh) 一种基于单幅图像的三维数字凹浮雕生成方法
Wang et al. A novel finger vein pattern extraction method using oriented filtering technology
Wang et al. A new method of speeded up robust features image registration based on image preprocessing
CN108961290B (zh) 一种基于Otsu的Ratio算子自适应图像边缘检测方法

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
TR01 Transfer of patent right

Effective date of registration: 20180427

Address after: 312000 Shaoxing, Zhejiang Province, 66 -9, three to 1 level.

Patentee after: ZHEJIANG AEROSPACE CHANGFENG SCIENCE & TECHNOLOGY DEVELOPMENT CO., LTD.

Address before: 310014 Zhejiang University of Technology, 18, Chao Wang Road, Xiacheng District, Hangzhou, Zhejiang

Patentee before: Zhejiang University of Technology

TR01 Transfer of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180227

Termination date: 20210612

CF01 Termination of patent right due to non-payment of annual fee