[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN104063874A - 基于权值模型的灰度图像二值化算法定量评估方法 - Google Patents

基于权值模型的灰度图像二值化算法定量评估方法 Download PDF

Info

Publication number
CN104063874A
CN104063874A CN201410324856.4A CN201410324856A CN104063874A CN 104063874 A CN104063874 A CN 104063874A CN 201410324856 A CN201410324856 A CN 201410324856A CN 104063874 A CN104063874 A CN 104063874A
Authority
CN
China
Prior art keywords
gray level
image
level image
gray
matrix
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
Application number
CN201410324856.4A
Other languages
English (en)
Other versions
CN104063874B (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.)
WUHAN LANSHAN TECHNOLOGY CO.,LTD.
Original Assignee
Wuhan University of Science and Engineering WUSE
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 Wuhan University of Science and Engineering WUSE filed Critical Wuhan University of Science and Engineering WUSE
Priority to CN201410324856.4A priority Critical patent/CN104063874B/zh
Publication of CN104063874A publication Critical patent/CN104063874A/zh
Application granted granted Critical
Publication of CN104063874B publication Critical patent/CN104063874B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)
  • Facsimile Image Signal Circuits (AREA)
  • Image Analysis (AREA)

Abstract

本发明涉及一种基于权值模型的灰度图像二值化算法定量评估方法。其技术方案为:先对灰度图像A使用n(n为7~20的自然数)种灰度图像二值化算法进行二值化处理,再对生成的二值图像集进行处理,得到灰度图像A的前景点权值矩阵和背景点权值矩阵;然后用待评估灰度图像二值化算法对灰度图像A进行处理,得到灰度图像A的转换矩阵Ⅰ和Ⅱ;然后根据待评估的灰度图像二值化算法的加权正确率P的数学模型,得到待评估的灰度图像二值化算法的加权正确率P,定量评估待评估灰度图像二值化算法。本发明不仅适用性强,独立于参照集,且能对灰度图像二值化算法进行有效的定量评估。

Description

基于权值模型的灰度图像二值化算法定量评估方法
技术领域
本发明属于灰度图像二值化算法定量评估技术领域。具体涉及一种基于权值模型的灰度图像二值化算法定量评估方法。
背景技术
由于科技的发展和生活节奏的加快,现代人已经不满足仅通过文字获得知识,渴望用图片刺激眼球,激发求知欲,这种渴望加速了图像社会或视觉文化时代的来临,读图时代已经成为当今一种主导性的、全面覆盖性的文化景观。二值化处理作为图像处理中比较重要的预处理,其评估变得日趋重要。
目前,最常用的灰度图像二值化算法评估方法是通过定性比较同一幅灰度图像的不同二值化结果,由此来说明提出的算法的好坏,这样做存在很大的主观性。
虽然曾有人提出了一些灰度图像二值化算法定量评估方法(Hui, Zhang. Image segmentation evaluation: A survey of unsupervised methods [J]. Computer Vision and Image Understanding , 2008, (110): 260-280),但是这些方法中有些是针对图像的某一特征进行定量评估的,评估不够全面。如区域内一致性评估方法,虽然该灰度图像二值化算法定量评估方法能取得一定的评估效果,但是该方法只考虑了二值化处理需满足区域内一致性的特征,忽略了二值化处理需满足区域间不一致等特征,评估不够全面。
除此之外,有些灰度图像二值化算法定量评估方法(Y. Zhang, A survey on evaluation methods for image segmentation, Pattern Recognition 29 (8) (1996) 1335–1346.)依赖于参照集。如基于错误分类像素数目的灰度图像二值化算法定量评估方法,该方法定量评估时需要统计分类错误的像素点数目,统计过程依赖于参照集导致该评估方法依赖于参照集,但对于真实图像,只能人工获取参照集,不仅耗时还存在主观性,严重影响了灰度图像二值化算法的定量评估。
发明内容
本发明旨在克服现有技术的不足, 目的是提供一种适用性强和独立于参照集的基于权值模型的灰度图像二值化算法定量评估方法。
为实现上述目的,本发明所采用的技术方案的具体步骤是:
步骤一、选择n(n为7~ 20的自然数)种灰度图像二值化算法对任一幅灰度图像A分别进行二值化处理,得到各自对应的灰度图像二值化算法对灰度图像A处理后的二值图像,所述的各自对应的灰度图像二值化算法对灰度图像A处理后的二值图像构成二值图像集。
步骤二、灰度图像A的前景点权值矩阵和背景点权值矩阵的生成
分步一、先取灰度图像A中的任一像素点B,像素点B在灰度图像A中的坐标为(x,y)。再确定二值图像集的每副二值图像中与像素点B的位置相同的坐标为(x,y)的像素点,所述的二值图像集的每副二值图像中坐标为(x,y)的像素点构成像素点集C。然后在所述像素点集C中统计被认为是前景点的像素总数和被认为是背景点的像素总数。
所述的被认为是前景点的像素总数是指像素点B被认为是前景点的次数,所述的被认为是背景点的像素总数是指像素点B被认为是背景点的次数。
分步二、灰度图像A中其余像素点被认为是前景点的次数和背景点的次数的确定方法同分步一。
分步三、根据灰度图像A中的每个像素点被认为是前景点的次数,由前景点权值矩阵中的权值生成规则,生成灰度图像A中的每个像素点在前景点权值矩阵中对应的权值 。根据灰度图像A中的每个像素点被认为是背景点的次数,由背景点权值矩阵中的权值生成规则,生成灰度图像A中的每个像素点在背景点权值矩阵中对应的权值
分步四、灰度图像A中的所有像素点在前景点权值矩阵中对应的权值构成灰度图像A的前景点权值矩阵;灰度图像A中的所有像素点在背景点权值矩阵中对应的权值构成灰度图像A的背景点权值矩阵。
步骤三、用待评估的灰度图像二值化算法对灰度图像A进行二值化处理,得到待评估的灰度图像二值化算法对灰度图像A处理后的二值图像和待评估的灰度图像二值化算法对灰度图像A处理后的二值图像灰度值矩阵。用矩阵转换规则Ⅰ对所述二值图像灰度值矩阵进行处理,得到灰度图像A的转换矩阵Ⅰ。再用矩阵转换规则Ⅱ对所述二值图像灰度值矩阵进行处理,得到灰度图像A的转换矩阵Ⅱ。
步骤四、待评估的灰度图像二值化算法的加权正确率P的数学模型:
(1)
式(1)中:表示灰度图像A的灰度值矩阵的行数;
表示灰度图像A的灰度值矩阵的列数;
表示灰度图像A的前景点权值矩阵的第i行第j列的权值;
表示灰度图像A的背景点权值矩阵的第i行第j列的权值;
表示灰度图像A的转换矩阵Ⅰ的第i行第j列的值;
表示灰度图像A的转换矩阵Ⅱ的第i行第j列的值;
表示灰度图像A的前景点总数,即灰度图像A的前景点权值 矩阵中权值为正数的像素点的总数;
表示灰度图像A的背景点总数,即灰度图像A的背景点权值 矩阵中权值为正数的像素点的总数。
根据式(1)得到待评估的灰度图像二值化算法的加权正确率P
步骤五、根据待评估的灰度图像二值化算法的加权正确率P越高则待评估的灰度图像二值化算法越好的规则,对比待评估的灰度图像二值化算法的加权正确率P,得到评估结果。
所述前景点是指二值图像中灰度值为255的像素点。
所述背景点是指二值图像中灰度值为0的像素点。
所述前景点权值矩阵中的权值生成规则是:
①当灰度图像A中某一像素点被认为是前景点的次数等于所选灰度图像二值化算法数目的一半时,则前景点权值矩阵中与所述某一像素点的坐标对应位置的权值为0。
②当灰度图像A中某一像素点被认为是前景点的次数大于所选灰度图像二值化算法数目的一半时,则前景点权值矩阵中与所述某一像素点的坐标对应位置的权值为:
(2)
式(2)中:i,j为所述灰度图像A中某一像素点所对应的横坐标和纵坐标;
表示所述灰度图像A中某一像素点被认为是前景点的次数;
表示所选择的二值化算法的数目,为7~20的自然数;
表示所选择的二值化算法数目的修正值;
当所选择的二值化算法数目为奇数时,
当所选择的二值化算法数目为偶数时,
③当灰度图像A中某一像素点被认为是前景点的次数小于所选灰度图像二值化算法数目的一半时,则前景点权值矩阵中与所述某一像素点的坐标对应位置的权值为:
(3)
式(3)中:i,j为所述灰度图像A中某一像素点所对应的横坐标和纵坐标;
表示所述灰度图像A中某一像素点被认为是前景点的次数;
表示所选择的二值化算法数目为7~20的自然数)的修正值;
当所选择的二值化算法数目为奇数时,
当所选择的二值化算法数目为偶数时,
所述背景点权值矩阵中的权值生成规则是:
①当灰度图像A中某一像素点被认为是背景点的次数等于所选灰度图像二值化算法数目的一半时,则背景点权值矩阵中与所述某一像素点的坐标对应位置的权值为0。
②当灰度图像A中某一像素点被认为是背景点的次数大于所选灰度图像二值化算法数目的一半时,则背景点权值矩阵中与所述某一像素点的坐标对应位置的权值为:
(4)
式(4)中:i,j为所述灰度图像A中某一像素点所对应的横坐标和纵坐标;
表示所述灰度图像A中某一像素点被认为是背景点的次数;
表示所选择的二值化算法的数目,为7~20的自然数;
表示所选择的二值化算法数目的修正值;
当所选择的二值化算法数目为奇数时,
当所选择的二值化算法数目为偶数时,
③当灰度图像A中某一像素点被认为是背景点的次数小于所选灰度图像二值化算法数目的一半时,则背景点权值矩阵中与所述某一像素点的坐标对应位置的权值为:
(5)
式(5)中:i,j为所述灰度图像A中某一像素点所对应的横坐标和纵坐标;
表示所述灰度图像A中某一像素点被认为是背景点的次数;
表示所选择的二值化算法数目为7~20的自然数)的修正值;
当所选择的二值化算法数目为奇数时,
当所选择的二值化算法数目为偶数时,
所述矩阵转换规则Ⅰ是:
①确定二值图像灰度值矩阵中灰度值为255的位置,将转换矩阵Ⅰ中与二值图像灰度值矩阵中灰度值为255的相同位置的数值置1。
②确定二值图像灰度值矩阵中灰度值为0的位置,将转换矩阵Ⅰ中与二值图像灰度值矩阵中灰度值为0的相同位置的数值置0。
所述矩阵转换规则Ⅱ是:
①确定二值图像灰度值矩阵中灰度值为0的位置,将转换矩阵Ⅱ中与二值图像灰度值矩阵中灰度值为0的相同位置的数值置1。
②确定二值图像灰度值矩阵中灰度值为255的位置,将转换矩阵Ⅱ中与二值图像灰度值矩阵中灰度值为255的相同位置的数值置0。
由于采用上述技术方案,本发明与现有技术相比具有如下积极效果:
1、本发明能克服现有技术对灰度图像二值化算法进行定量评估的缺陷,现有技术中大多数的灰度图像二值化算法评估方法是通过定性比较同一幅灰度图像的不同二值化结果,由此来说明提出的算法的好坏,这样做存在很大的主观性。本发明通过确定灰度图像的前景点权值矩阵、背景点权值矩阵、转换矩阵Ⅰ和转换矩阵Ⅱ计算出待评估的灰度图像二值化算法的的加权正确率P,用待评估的灰度图像二值化算法的加权正确率P对灰度图像二值化算法进行定量评估,有效地克服了评估时存在的主观性缺陷。
2、本发明克服了对灰度图像二值化算法进行定量评估不全面的缺陷,目前有些灰度图像二值化算法定量评估方法只针对图像的某一特征进行定量评估,忽略了图像的其余特征的评估。本发明的评估对象不是二值图像的特征,而是通过确定二值图像集,在二值图像集中提取权值矩阵作为待评估的灰度图像二值化算法的加权正确率P的权值,有效地克服灰度图像二值化算法进行定量评估只针对某一特征的缺陷。
3、本发明克服了在没有参照集时对灰度图像二值化算法进行定量评估的缺陷,基于实验的灰度图像二值化算法定量评估方法往往依赖于参照集的选取,但对于真实图像,只能人工获取参照集,不仅耗时还存在主观性。本发明的定量评估方法独立于参照集,对于待评估的灰度图像二值化算法的评估不需要参照集作为评估标准,有效地克服了灰度图像二值化算法定量评估方法依赖于参照集的缺陷。
因此,本发明不仅适用性强和独立于参照集,且能进行有效的定量评估。
附图说明
图1是一种带钢缺陷灰度图像;
图2是对图1采用Otsu得到的二值图像;
图3是一种带钢缺陷灰度图像;
图4是对图3采用曲面拟合阈值法得到的二值图像;
图5是一种带钢缺陷灰度图像;
图6是对图5采用局部对比法得到的二值图像。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步的描述,并非对其保护范围的限制。
实施例1
一种基于权值模型的灰度图像二值化算法定量评估方法。本实施例所述的灰度图像二值化算法定量评估方法的具体步骤是:
步骤一、选择7种灰度图像二值化算法对任一幅灰度图像A分别进行二值化处理,得到各自对应的灰度图像二值化算法对灰度图像A处理后的二值图像,所述的各自对应的灰度图像二值化算法对灰度图像A处理后的二值图像构成二值图像集。
本实施例所述的灰度图像A是采用一幅如图1所示的带钢缺陷灰度图像。
本实施例所述的7种灰度图像二值化算法为:Otsu、迭代法、最大熵法、双峰法、Bernsen、最小误差法和平均灰度值法。
步骤二、灰度图像A的前景点权值矩阵和背景点权值矩阵的生成
分步一、先取灰度图像A中的任一像素点B,像素点B在灰度图像A中的坐标为(x,y)。再确定二值图像集的每副二值图像中与像素点B的位置相同的坐标为(x,y)的像素点,所述的二值图像集的每副二值图像中坐标为(x,y)的像素点构成像素点集C。然后在所述像素点集C中统计被认为是前景点的像素总数和被认为是背景点的像素总数。
所述的被认为是前景点的像素总数是指像素点B被认为是前景点的次数,所述的被认为是背景点的像素总数是指像素点B被认为是背景点的次数。
本实施例所述的前景点是指二值图像中灰度值为255的像素点,背景点是指二值图像中灰度值为0的像素点。取灰度图像A中坐标为(1,1)的像素点B,确定二值图像集的每幅二值图像中与像素点B的位置相同的坐标为(1,1)的像素点,得到二值图像集中与像素点B的坐标相同的像素点集C。 在所述像素点集C中统计得:在像素点集C中有0个像素点的灰度值为255,有7个像素点的灰度值为0。即像素点B被认为是前景点的次数为0,像素点B被认为是背景点的次数为7。
分步二、灰度图像A中其余像素点被认为是前景点的次数和背景点的次数的确定方法同分步一。
分步三、根据灰度图像A中的每个像素点被认为是前景点的次数,由前景点权值矩阵中的权值生成规则,生成灰度图像A中的每个像素点在前景点权值矩阵中对应的权值。根据灰度图像A中的每个像素点被认为是背景点的次数,由背景点权值矩阵中的权值生成规则,生成灰度图像A中的每个像素点在背景点权值矩阵中对应的权值
本实施例所述的像素点B在前景点权值矩阵中对应的权值为:
本实施例中灰度图像A中像素点B被认为是前景点的次数为0,所选择的灰度图像二值化算法数目为7,即灰度图像A中的坐标为(1,1)的像素点B被认为是前景点的次数小于所选灰度图像二值化算法数目的一半,则前景点权值矩阵中与像素点B的坐标对应位置的权值为:
式中:i,j为所述灰度图像A中某一像素点所对应的横坐标1和纵坐标1;
表示所述灰度图像A中某一像素点被认为是前景点的次数, =0;
表示所选择的二值化算法数目为7~20的自然数)的修正值;
本实施例所选择的二值化算法数目为7,
则灰度图像A中的像素点B在前景点权值矩阵中对应的权值为:
灰度图像A中其余像素点在前景点权值矩阵中对应的权值生成规则:
①当灰度图像A中某一像素点被认为是前景点的次数等于所选灰度图像二值化算法数目的一半时,则前景点权值矩阵中与所述某一像素点的坐标对应位置的权值为0。
②当灰度图像A中某一像素点被认为是前景点的次数大于所选灰度图像二值化算法数目的一半时,则前景点权值矩阵中与所述某一像素点的坐标对应位置的权值为:
(1)
式(1)中:i,j为所述灰度图像A中某一像素点所对应的横坐标和纵坐标;
表示所述灰度图像A中某一像素点被认为是前景点的次数;
表示所选择的二值化算法的数目,为7~20的自然数;
表示所选择的二值化算法数目的修正值;
当所选择的二值化算法数目为奇数时,
当所选择的二值化算法数目为偶数时,
③当灰度图像A中某一像素点被认为是前景点的次数小于所选灰度图像二值化算法数目的一半时,则前景点权值矩阵中与所述某一像素点的坐标对应位置的权值为:
(2)
式(2)中,i,j为所述灰度图像A中某一像素点所对应的横坐标和纵坐标;
表示所述灰度图像A中某一像素点被认为是前景点的次数;
表示所选择的二值化算法数目为7~20的自然数)的修正值;
当所选择的二值化算法数目为奇数时,
当所选择的二值化算法数目为偶数时,
本实施例所述的像素点B在背景点权值矩阵中对应的权值为:
本实施例中灰度图像A中像素点B被认为是背景点的次数为7,所选择的灰度图像二值化算法数目为7,即灰度图像A中的坐标为(1,1)的像素点B被认为是背景点的次数大于所选灰度图像二值化算法数目的一半,则背景点权值矩阵中与像素点B的坐标对应位置的权值为:
式中:i,j为所述灰度图像A中某一像素点所对应的横坐标1和纵坐标1;
表示所述灰度图像A中某一像素点被认为是背景点的次数, =7;
表示所选择的二值化算法的数目,为7;
表示所选择的二值化算法数目为7~20的自然数)的修正值,
本实施例所选择的二值化算法数目为7,
则灰度图像A中的像素点B在背景点权值矩阵中对应的权值为:
灰度图像A中其余像素点在背景点权值矩阵中对应的权值生成规则:
①当灰度图像A中某一像素点被认为是背景点的次数等于所选灰度图像二值化算法数目的一半时,则背景点权值矩阵中与所述某一像素点的坐标对应位置的权值为0。
②当灰度图像A中某一像素点被认为是背景点的次数大于所选灰度图像二值化算法数目的一半时,则背景点权值矩阵中与所述某一像素点的坐标对应位置的权值为:
(3)
式(3)中:i,j为所述灰度图像A中某一像素点所对应的横坐标和纵坐标;
表示所述灰度图像A中某一像素点被认为是背景点的次数;
表示所选择的二值化算法的数目,为7~20的自然数;
表示所选择的二值化算法数目的修正值;
当所选择的二值化算法数目为奇数时,
当所选择的二值化算法数目为偶数时,
③当灰度图像A中某一像素点被认为是背景点的次数小于所选灰度图像二值化算法数目的一半时,则背景点权值矩阵中与所述某一像素点的坐标对应位置的权值为:
(4)
式(4)中:i,j为所述灰度图像A中某一像素点所对应的横坐标和纵坐标;
表示所述灰度图像A中某一像素点被认为是背景点的次数;
表示所选择的二值化算法数目为7~20的自然数)的修正值;
当所选择的二值化算法数目为奇数时,
当所选择的二值化算法数目为偶数时,
分步四、灰度图像A中的所有像素点在前景点权值矩阵中对应的权值构成灰度图像A的前景点权值矩阵;灰度图像A中的所有像素点在背景点权值矩阵中对应的权值构成灰度图像A的背景点权值矩阵。
本实施例求得灰度图像A中每一个像素点在前景点权值矩阵中对应的权值,灰度图像A中的所有像素点在前景点权值矩阵中对应的权值构成灰度图像A的45行73列的前景点权值矩阵:
本实施例求得灰度图像A中每一个像素点在背景点权值矩阵中对应的权值,灰度图像A中的所有像素点在背景点权值矩阵中对应的权值构成灰度图像A的45行73列的背景点权值矩阵:
步骤三、用待评估的灰度图像二值化算法对灰度图像A进行二值化处理,得到待评估的灰度图像二值化算法对灰度图像A处理后的二值图像和待评估的灰度图像二值化算法对灰度图像A处理后的二值图像灰度值矩阵。用矩阵转换规则Ⅰ对所述二值图像灰度值矩阵进行处理,得到灰度图像A的转换矩阵Ⅰ。再用矩阵转换规则Ⅱ对所述二值图像灰度值矩阵进行处理,得到灰度图像A的转换矩阵Ⅱ。
本实施例的所述待评估的灰度图像二值化算法为Otsu。用Otsu对灰度图像A进行二值化处理,得到如图2所示的待评估的灰度图像二值化算法对灰度图像A处理后的二值图像和待评估的灰度图像二值化算法对灰度图像A处理后的二值图像灰度值矩阵,所述二值图像灰度值矩阵是一个大小为45行73列的一个矩阵:
本实施例所述矩阵转换规则Ⅰ是:
①确定二值图像灰度值矩阵中灰度值为255的位置,将转换矩阵Ⅰ中与二值图像灰度值矩阵中灰度值为255的相同位置的数值置1。如本实施例中二值图像灰度值矩阵的第44行第72列的灰度值为255,则根据所述规则将转换矩阵Ⅰ中第44行第72列的数值置1。
②确定二值图像灰度值矩阵中灰度值为0的位置,将转换矩阵Ⅰ中与二值图像灰度值矩阵中灰度值为0的相同位置的数值置0。如本实施例中二值图像灰度值矩阵的第1行第1列的灰度值为0,则根据所述规则将转换矩阵Ⅰ中第1行第1列的数值置0。
本实施例求得的灰度图像A的转换矩阵Ⅰ是一个大小为45行73列的一个矩阵:
本实施例所述矩阵转换规则Ⅱ是:
①确定二值图像灰度值矩阵中灰度值为0的位置,将转换矩阵Ⅱ中与二值图像灰度值矩阵中灰度值为0的相同位置的数值置1。如本实施例中二值图像灰度值矩阵的第1行第1列的灰度值为0,则根据所述规则将转换矩阵Ⅱ中第1行第1列的数值置1。
②确定二值图像灰度值矩阵中灰度值为255的位置,将转换矩阵Ⅱ中与二值图像灰度值矩阵中灰度值为255的相同位置的数值置0。如本实施例中二值图像灰度值矩阵的第44行第72列的灰度值为255,则根据所述规则将转换矩阵Ⅱ中第44行第72列的数值置0。
本实施例求得的灰度图像A的转换矩阵Ⅱ是一个大小为45行73列的一个矩阵:
步骤四、待评估的灰度图像二值化算法的加权正确率P的数学模型:
(5)
式(5)中:表示灰度图像A的灰度值矩阵的行数;
表示灰度图像A的灰度值矩阵的列数;
表示灰度图像A的前景点权值矩阵的第i行第j列的权值;
表示灰度图像A的背景点权值矩阵的第i行第j列的权值;
表示灰度图像A的转换矩阵Ⅰ的第i行第j列的值;
表示灰度图像A的转换矩阵Ⅱ的第i行第j列的值;
表示灰度图像A的前景点总数,即灰度图像A的前景点权值
矩阵中权值为正数的像素点的总数;
表示灰度图像A的背景点总数,即灰度图像A的背景点权值
矩阵中权值为正数的像素点的总数。
本实施例中灰度图像A的灰度值矩阵的行数为45,列数为73,灰度图像A的前景点总数为1642,背景点总数为1643。
根据式(5),本实施例得到的待评估的Otsu的加权正确率P为0.8089。
同理:本实施例得到的待评估的交叉熵阈值法的加权正确率P为0.6679;
本实施例得到的待评估的局部对比法的加权正确率P为0.7294。
步骤五、根据待评估的灰度值图像二值化算法的加权正确率P越高则待评估的灰度图像二值化算法越好的规则,对比待评估的灰度图像二值化算法的加权正确率P,则得到Otsu对灰度图像A即图1所示的带钢缺陷灰度图像的二值化效果最好。
实施例2
一种基于权值模型的灰度图像二值化算法定量评估方法。本实施例所述的灰度图像二值化算法定量评估方法的具体步骤是:
步骤一、选择14种灰度图像二值化算法对任一幅灰度图像A分别进行二值化处理,得到各自对应的灰度图像二值化算法对灰度图像A处理后的二值图像,所述的各自对应的灰度图像二值化算法对灰度图像A处理后的二值图像构成二值图像集。
本实施例所述的灰度图像A是采用一幅如图3所示的带钢缺陷灰度图像。
本实施例的所述14种灰度图像二值化算法为:Otsu、迭代法、最大熵法、双峰法、Bernsen、最小误差法、平均灰度值法、聚类阈值法、模糊聚类阈值法、模糊熵阈值法、交叉熵阈值法、边缘区域匹配阈值法、模糊相似度阈值法和共生阈值法。
步骤二、灰度图像A的前景点权值矩阵和背景点权值矩阵的生成
分步一、先取灰度图像A中的任一像素点B,像素点B在灰度图像A中的坐标为(x,y)。再确定二值图像集的每副二值图像中与像素点B的位置相同的坐标为(x,y)的像素点,所述的二值图像集的每副二值图像中坐标为(x,y)的像素点构成像素点集C。然后在所述像素点集C中统计被认为是前景点的像素总数和被认为是背景点的像素总数。
所述的被认为是前景点的像素总数是指像素点B被认为是前景点的次数,所述的被认为是背景点的像素总数是指像素点B被认为是背景点的次数。
本实施例所述的前景点是指二值图像中灰度值为255的像素点,背景点是指二值图像中灰度值为0的像素点。取灰度图像A中坐标为(1,2)的像素点B,确定二值图像集的每幅二值图像中与像素点B的位置相同的坐标为(1,2)的像素点,得到二值图像集中与像素点B的坐标相同的像素点集C。 在所述像素点集C中统计得:在像素点集C中有2个像素点的灰度值为255,有12个像素点的灰度值为0。即像素点B被认为是前景点的次数为2,像素点B被认为是背景点的次数为12。
分步二、灰度图像A中其余像素点被认为是前景点的次数和背景点的次数的确定方法同分步一。
分步三、根据灰度图像A中的每个像素点被认为是前景点的次数,由前景点权值矩阵中的权值生成规则,生成灰度图像A中的每个像素点在前景点权值矩阵中对应的权值。根据灰度图像A中的每个像素点被认为是背景点的次数,由背景点权值矩阵中的权值生成规则,生成灰度图像A中的每个像素点在背景点权值矩阵中对应的权值
本实施例所述的像素点B在前景点权值矩阵中对应的权值为:
本实施例中灰度图像A中像素点B被认为是前景点的次数为2,所选择的灰度图像二值化算法数目为14,即灰度图像A中的坐标为(1,2)的像素点B被认为是前景点的次数小于所选灰度图像二值化算法数目的一半,则前景点权值矩阵中与像素点B的坐标对应位置的权值为:
(1)
式(1)中:i,j为所述灰度图像A中某一像素点所对应的横坐标1和纵坐标2;
表示所述灰度图像A中某一像素点被认为是前景点的次数, =2;
表示所选择的二值化算法数目为7~20的自然数)的修正值;
本实施例所选择的二值化算法数目为14,
则灰度图像A中的像素点B在前景点权值矩阵中对应的权值为:
灰度图像A中其余像素点在前景点权值矩阵中对应的权值生成规则同实施例1的其余像素点在前景点权值矩阵中对应的权值生成规则。
本实施例所述的像素点B在背景点权值矩阵中对应的权值为:
本实施例中灰度图像A中像素点B被认为是背景点的次数为12,所选择的灰度图像二值化算法数目为14,即灰度图像A中的坐标为(1,2)的像素点B被认为是背景点的次数大于所选灰度图像二值化算法数目的一半,则背景点权值矩阵中与像素点B的坐标对应位置的权值为:
(2)
式(2)中:i,j为所述灰度图像A中某一像素点所对应的横坐标1和纵坐标2;
表示所述灰度图像A中某一像素点被认为是背景点的次数, =12;
表示所选择的二值化算法的数目,为14;
表示所选择的二值化算法数目为7~20的自然数)的修正值,
本实施例所选择的二值化算法数目为14,
则灰度图像A中的像素点B在背景点权值矩阵中对应的权值为:
灰度图像A中其余像素点在背景点权值矩阵中对应的权值生成规则同实施例1中的其余像素点在背景点权值矩阵中对应的权值生成规则。
分步四、灰度图像A中的所有像素点在前景点权值矩阵中对应的权值构成灰度图像A的前景点权值矩阵;灰度图像A中的所有像素点在背景点权值矩阵中对应的权值构成灰度图像A的背景点权值矩阵。
本实施例求得灰度图像A中每一个像素点在前景点权值矩阵中对应的权值,灰度图像A中的所有像素点在前景点权值矩阵中对应的权值构成灰度图像A的140行99列的前景点权值矩阵:
本实施例求得灰度图像A中每一个像素点在背景点权值矩阵中对应的权值,灰度图像A中的所有像素点在背景点权值矩阵中对应的权值构成灰度图像A的140行99列的背景点权值矩阵:
步骤三、用待评估的灰度图像二值化算法对灰度图像A进行二值化处理,得到待评估的灰度图像二值化算法对灰度图像A处理后的二值图像和待评估的灰度图像二值化算法对灰度图像A处理后的二值图像灰度值矩阵。用矩阵转换规则Ⅰ对所述二值图像灰度值矩阵进行处理,得到灰度图像A的转换矩阵Ⅰ。再用矩阵转换规则Ⅱ对所述二值图像灰度值矩阵进行处理,得到灰度图像A的转换矩阵Ⅱ。
本实施例的所述待评估的灰度图像二值化算法为曲面拟合阈值法。用曲面拟合阈值法对灰度图像A进行二值化处理,得到如图4所示的待评估的灰度图像二值化算法对灰度图像A处理后的二值图像和待评估的灰度图像二值化算法对灰度图像A处理后的二值图像灰度值矩阵,所述的二值图像灰度值矩阵是一个大小为140行99列的一个矩阵:
本实施例所述矩阵转换规则Ⅰ是:
①确定二值图像灰度值矩阵中灰度值为255的位置,将转换矩阵Ⅰ中与二值图像灰度值矩阵中灰度值为255的相同位置的数值置1。如本实施例中二值图像灰度值矩阵的第140行第85列为255,则根据所述规则将转换矩阵Ⅰ中第140行第85列置1。
②确定二值图像灰度值矩阵中灰度值为0的位置,将转换矩阵Ⅰ中与二值图像灰度值矩阵中灰度值为0的相同位置的数值置0。如本实施例中二值图像灰度值矩阵的第1行第1列为0,则根据所述规则将转换矩阵Ⅰ中第1行第1列置0。
本实施例求得的灰度图像A的转换矩阵Ⅰ是一个大小为140行99列的一个矩阵:
本实施例所述矩阵转换规则Ⅱ是:
①确定二值图像灰度值矩阵中灰度值为0的位置,将转换矩阵Ⅱ中与二值图像灰度值矩阵中灰度值为0的相同位置的数值置1。如本实施例中二值图像灰度值矩阵的第1行第1列为0,则根据所述规则将转换矩阵Ⅱ中第1行第1列置1。
②确定二值图像灰度值矩阵中灰度值为255的位置,将转换矩阵Ⅱ中与二值图像灰度值矩阵中灰度值为255的相同位置的数值置0。如本实施例中二值图像灰度值矩阵的第140行第85列为255,则根据所述规则将转换矩阵Ⅱ中第140行第85列置0。
本实施例求得的灰度图像A的转换矩阵Ⅱ是一个大小为140行99列的一个矩阵:
步骤四、待评估的灰度图像二值化算法的加权正确率P的数学模型:
(3)
式(3)中:表示灰度图像A的灰度值矩阵的行数;
表示灰度图像A的灰度值矩阵的列数;
表示灰度图像A的前景点权值矩阵的第i行第j列的权值;
表示灰度图像A的背景点权值矩阵的第i行第j列的权值;
表示灰度图像A的转换矩阵Ⅰ的第i行第j列的值;
表示灰度图像A的转换矩阵Ⅱ的第i行第j列的值;
表示灰度图像A的前景点总数,即灰度图像A的前景点权值矩阵中权值为正数的像素点的总数;
表示灰度图像A的背景点总数,即灰度图像A的背景点权值矩阵中权值为正数的像素点的总数。
本实施例中灰度图像A的灰度值矩阵的行数为140,列数为99,灰度图像A的前景点总数为7152,背景点总数为6318。
根据式(3),本实施例求得的待评估的的曲面拟合阈值法加权正确率P为0.7131。
同理:本实施例得到的待评估的共生阈值法的加权正确率P为0.8013;
本实施例得到的待评估的模糊熵阈值法的加权正确率P为0.5431。
步骤五、根据待评估的灰度值图像二值化算法的加权正确率P越高则待评估的灰度图像二值化算法越好的规则,对比待评估的灰度图像二值化算法的加权正确率P,则得到共生阈值法对灰度图像A即图3所示的带钢缺陷灰度图像的二值化效果最好。
实施例3
一种基于权值模型的灰度图像二值化算法定量评估方法。本实施例所述的灰度图像二值化算法定量评估方法的具体步骤是:
步骤一、选择20种灰度图像二值化算法对任一幅灰度图像A分别进行二值化处理,得到各自对应的灰度图像二值化算法对灰度图像A处理后的二值图像,所述的各自对应的灰度图像二值化算法对灰度图像A处理后的二值图像构成二值图像集。
本实施例所述的灰度图像A是采用一幅如图5所示的带钢缺陷灰度图像。
本实施例所述的20种灰度图像二值化算法为:Otsu、迭代法、最大熵法、双峰法、Bernsen、最小误差法、平均灰度值法、聚类阈值法、模糊聚类阈值法、模糊熵阈值法、交叉熵阈值法、边缘区域匹配阈值法、模糊相似度阈值法、共生阈值法、高阶熵阈值法、基于随机集的阈值法、局部对比法、曲面拟合阈值法、形状建模阈值法和凸包阈值法。
步骤二、灰度图像A的前景点权值矩阵和背景点权值矩阵的生成
分步一、先取灰度图像A中的任一像素点B,像素点B在灰度图像A中的坐标为(x,y)。再确定二值图像集的每副二值图像中与像素点B的位置相同的坐标为(x,y)的像素点,所述的二值图像集的每副二值图像中坐标为(x,y)的像素点构成像素点集C。然后在所述像素点集C中统计被认为是前景点的像素总数和被认为是背景点的像素总数。
所述的被认为是前景点的像素总数是指像素点B被认为是前景点的次数,所述的被认为是背景点的像素总数是指像素点B被认为是背景点的次数。
本实施例所述的前景点是指二值图像中灰度值为255的像素点,背景点是指二值图像中灰度值为0的像素点。取灰度图像A中坐标为(3,3)的像素点B,确定二值图像集的每幅二值图像中与像素点B的位置相同的坐标为(3,3)的像素点,得到二值图像集中与像素点B的坐标相同的像素点集C。 在所述像素点集C中统计得:在像素点集C中有17个像素点的灰度值为255,有3个像素点的灰度值为0。即像素点B被认为是前景点的次数为17,像素点B被认为是背景点的次数为3。
分步二、灰度图像A中其余像素点被认为是前景点的次数和背景点的次数的确定方法同分步一。
分步三、根据灰度图像A中的每个像素点被认为是前景点的次数,由前景点权值矩阵中的权值生成规则,生成灰度图像A中的每个像素点在前景点权值矩阵中对应的权值。根据灰度图像A中的每个像素点被认为是背景点的次数,由背景点权值矩阵中的权值生成规则,生成灰度图像A中的每个像素点在背景点权值矩阵中对应的权值
本实施例所述的像素点B在前景点权值矩阵中对应的权值为
本实施例中灰度图像A中像素点B被认为是前景点的次数为17,所选择的灰度图像二值化算法数目为20,即灰度图像A中的坐标为(3,3)的像素点B被认为是前景点的次数大于所选灰度图像二值化算法数目的一半,则前景点权值矩阵中与像素点B的坐标对应位置的权值为:
(1)
式(1)中:i,j为所述灰度图像A中某一像素点所对应的横坐标3和纵坐标3;
表示所述灰度图像A中某一像素点被认为是前景点的次数, =17;
表示所选择的二值化算法的数目,为20的自然数;
表示所选择的二值化算法数目的修正值;
本实施例所选择的二值化算法数目为20,
则灰度图像A中的像素点B在前景点权值矩阵中对应的权值为:
灰度图像A中其余像素点在前景点权值矩阵中对应的权值生成规则同实施例1的其余像素点在前景点权值矩阵中对应的权值生成规则。
本实施例所述的像素点B在背景点权值矩阵中对应的权值为:
本实施例中灰度图像A中像素点B被认为是背景点的次数为3,所选择的灰度图像二值化算法数目为20,即灰度图像A中的坐标为(3,3)的像素点B被认为是背景点的次数小于所选灰度图像二值化算法数目的一半,则背景点权值矩阵中与像素点B的坐标对应位置的权值为:
(2)
式(2)中:i,j为所述灰度图像A中某一像素点所对应的横坐标3和纵坐标3;
表示所述灰度图像A中某一像素点被认为是背景点的次数, =3;
表示所选择的二值化算法数目为7~20的自然数)的修正值;
本实施例所选择的二值化算法数目为20,
则灰度图像A中的像素点B在背景点权值矩阵中对应的权值为:
灰度图像A中其余像素点在背景点权值矩阵中对应的权值生成规则同实施例1的其余像素点在背景点权值矩阵中对应的权值生成规则。
分步四、灰度图像A中的所有像素点在前景点权值矩阵中对应的权值构成灰度图像A的前景点权值矩阵;灰度图像A中的所有像素点在背景点权值矩阵中对应的权值构成灰度图像A的背景点权值矩阵。
本实施例求得灰度图像A中每一个像素点在前景点权值矩阵中对应的权值,灰度图像A中的所有像素点在前景点权值矩阵中对应的权值构成灰度图像A的147行60列的前景点权值矩阵:
本实施例求得灰度图像A中每一个像素点在背景点权值矩阵中对应的权值,灰度图像A中的所有像素点在背景点权值矩阵中对应的权值构成灰度图像A的147行60列的背景点权值矩阵:
步骤三、用待评估的灰度图像二值化算法对灰度图像A进行二值化处理,得到待评估的灰度图像二值化算法对灰度图像A处理后的二值图像和待评估的灰度图像二值化算法对灰度图像A处理后的二值图像灰度值矩阵。用矩阵转换规则Ⅰ对所述二值图像灰度值矩阵进行处理,得到灰度图像A的转换矩阵Ⅰ。再用矩阵转换规则Ⅱ对所述二值图像灰度值矩阵进行处理,得到灰度图像A的转换矩阵Ⅱ。
本实施例的所述待评估的灰度图像二值化算法为局部对比法。用局部对比法对灰度图像A进行二值化处理,得到如图6所示的待评估的灰度图像二值化算法对灰度图像A处理后的二值图像和待评估的灰度图像二值化算法对灰度图像A处理后的二值图像灰度值矩阵,所述的二值图像灰度值矩阵是一个大小为147行60列的一个矩阵:
本实施例所述矩阵转换规则Ⅰ是:
①确定二值图像灰度值矩阵中灰度值为255的位置,将转换矩阵Ⅰ中与二值图像灰度值矩阵中灰度值为255的相同位置的数值置1。如本实施例中二值图像灰度值矩阵的第1行第1列为255,则根据所述规则将转换矩阵Ⅰ中第1行第1列置1。
②确定二值图像灰度值矩阵中灰度值为0的位置,将转换矩阵Ⅰ中与二值图像灰度值矩阵中灰度值为0的相同位置的数值置0。如本实施例中二值图像灰度值矩阵的第147行第52列为0,则根据所述规则将转换矩阵Ⅰ中第147行第52列置0。
本实施例求得的灰度图像A的转换矩阵Ⅰ是一个大小为147行60列的一个矩阵:
本实施例所述矩阵转换规则Ⅱ是:
①确定二值图像的灰度值矩阵中灰度值为0的位置,将转换矩阵Ⅱ中与二值图像灰度值矩阵中灰度值为0的相同位置的数值置1。如本实施例中二值图像灰度值矩阵的第147行第52列为0,则根据所述规则将转换矩阵Ⅱ中第147行第52列置1。
②确定二值图像的灰度值矩阵中灰度值为255的位置,将转换矩阵Ⅱ中与二值图像灰度值矩阵中灰度值为255的相同位置的数值置0。如本实施例中二值图像灰度值矩阵的第1行第1列为255,则根据所述规则将转换矩阵Ⅱ中第1行第1列置0。
本实施例求得的灰度图像A的转换矩阵Ⅱ是一个大小为147行60列的一个矩阵:
步骤四、待评估的灰度图像二值化算法的加权正确率P的数学模型:
(3)
式(3)中:表示灰度图像A的灰度值矩阵的行数;
表示灰度图像A的灰度值矩阵的列数;
表示灰度图像A的前景点权值矩阵的第i行第j列的权值;
表示灰度图像A的背景点权值矩阵的第i行第j列的权值;
表示灰度图像A的转换矩阵Ⅰ的第i行第j列的值;
表示灰度图像A的转换矩阵Ⅱ的第i行第j列的值;
表示灰度图像A的前景点总数,即灰度图像A的前景点权值
矩阵中权值为正数的像素点的总数;
表示灰度图像A的背景点总数,即灰度图像A的背景点权值
矩阵中权值为正数的像素点的总数。
本实施例中灰度图像A的灰度值矩阵的行数为147,列数为60,灰度图像A的前景点总数为4191,背景点总数为4629。
根据式(3),本实施例得到的待评估的局部对比法的加权正确率P为0.2637。
同理:本实施例得到的待评估的高阶熵阈值法的加权正确率P为0.5673;
本实施例得到的待评估的形状建模阈值法的加权正确率P为0.6328。
步骤五、根据待评估的灰度值图像二值化算法的加权正确率P越高则待评估的灰度图像二值化算法越好的规则,对比待评估的灰度图像二值化算法的加权正确率P,则得到形状建模阈值法对灰度图像A即图5所示的带钢缺陷灰度图像的二值化效果最好。
本具体实施方式与现有技术相比具有如下积极效果:
1、本具体实施方式能克服现有技术对灰度图像二值化算法进行定量评估的缺陷,现有技术中大多数的灰度图像二值化算法评估方法是通过定性比较同一幅灰度图像的不同二值化结果,由此来说明提出的算法的好坏,这样做存在很大的主观性。本具体实施方式通过确定灰度图像的前景点权值矩阵、背景点权值矩阵、转换矩阵Ⅰ和转换矩阵Ⅱ计算出待评估的灰度图像二值化算法的加权正确率P,用待评估的灰度图像二值化算法的加权正确率P对灰度图像二值化算法进行定量评估,有效地克服了评估时存在的主观性缺陷。
2、本具体实施方式克服了对灰度图像二值化算法进行定量评估不全面的缺陷,目前有些灰度图像二值化算法定量评估方法只针对图像的某一特征进行定量评估,忽略了图像的其余特征的评估。本具体实施方式的评估对象不是二值图像的特征,而是通过确定二值图像集,在二值图像集中提取权值矩阵作为待评估的灰度图像二值化算法的加权正确率P的权值,有效地克服灰度图像二值化算法进行定量评估只针对某一特征的缺陷。
3、本具体实施方式克服了在没有参照集时对灰度图像二值化算法进行定量评估的缺陷,基于实验的灰度图像二值化算法定量评估方法往往依赖于参照集的选取,但对于真实图像,只能人工获取参照集,不仅耗时还存在主观性。本具体实施方式的定量评估方法独立于参照集,对于待评估的灰度图像二值化算法的评估不需要参照集作为评估标准,有效地克服了灰度图像二值化算法定量评估方法依赖于参照集的缺陷。
因此,本具体实施方式不仅适用性强和独立于参照集,且能进行有效的定量评估。

Claims (7)

1.一种基于权值模型的灰度图像二值化算法定量评估方法,其特征在于所述方法是:
步骤一、选择n(n为7~ 20的自然数)种灰度图像二值化算法对任一幅灰度图像A分别进行二值化处理,得到各自对应的灰度图像二值化算法对灰度图像A处理后的二值图像,所述的各自对应的灰度图像二值化算法对灰度图像A处理后的二值图像构成二值图像集;
步骤二、灰度图像A的前景点权值矩阵和背景点权值矩阵的生成
分步一、先取灰度图像A中的任一像素点B,像素点B在灰度图像A中的坐标为(x,y);再确定二值图像集的每副二值图像中与像素点B的位置相同的坐标为(x,y)的像素点,所述的二值图像集的每副二值图像中坐标为(x,y)的像素点构成像素点集C;然后在所述像素点集C中统计被认为是前景点的像素总数和被认为是背景点的像素总数;
所述的被认为是前景点的像素总数是指像素点B被认为是前景点的次数,所述的被认为是背景点的像素总数是指像素点B被认为是背景点的次数;
分步二、灰度图像A中其余像素点被认为是前景点的次数和背景点的次数的确定方法同分步一;
分步三、根据灰度图像A中的每个像素点被认为是前景点的次数,由前景点权值矩阵中的权值生成规则,生成灰度图像A中的每个像素点在前景点权值矩阵中对应的权值 ;根据灰度图像A中的每个像素点被认为是背景点的次数,由背景点权值矩阵中的权值生成规则,生成灰度图像A中的每个像素点在背景点权值矩阵中对应的权值
分步四、灰度图像A中的所有像素点在前景点权值矩阵中对应的权值构成灰度图像A的前景点权值矩阵;灰度图像A中的所有像素点在背景点权值矩阵中对应的权值构成灰度图像A的背景点权值矩阵;
步骤三、用待评估的灰度图像二值化算法对灰度图像A进行二值化处理,得到待评估的灰度图像二值化算法对灰度图像A处理后的二值图像和待评估的灰度图像二值化算法对灰度图像A处理后的二值图像灰度值矩阵;用矩阵转换规则Ⅰ对所述二值图像灰度值矩阵进行处理,得到灰度图像A的转换矩阵Ⅰ;再用矩阵转换规则Ⅱ对所述二值图像灰度值矩阵进行处理,得到灰度图像A的转换矩阵Ⅱ;
步骤四、待评估的灰度图像二值化算法的加权正确率P的数学模型:
(1)
式(1)中:表示灰度图像A的灰度值矩阵的行数,
表示灰度图像A的灰度值矩阵的列数,
表示灰度图像A的前景点权值矩阵的第i行第j列的权值,
表示灰度图像A的背景点权值矩阵的第i行第j列的权值,
表示灰度图像A的转换矩阵Ⅰ的第i行第j列的值,
表示灰度图像A的转换矩阵Ⅱ的第i行第j列的值,
表示灰度图像A的前景点总数,即灰度图像A的前景点权值 矩阵中权值为正数的像素点的总数,
表示灰度图像A的背景点总数,即灰度图像A的背景点权值 矩阵中权值为正数的像素点的总数;
根据式(1)得到待评估的灰度图像二值化算法的加权正确率P
步骤五、根据待评估的灰度图像二值化算法的加权正确率P越高则待评估的灰度图像二值化算法越好的规则,对比待评估的灰度图像二值化算法的加权正确率P,得到评估结果。
2.根据权利要求1所述的基于权值模型的灰度图像二值化算法定量评估方法,其特征在于所述前景点是指二值图像中灰度值为255的像素点。
3.根据权利要求1所述的基于权值模型的灰度图像二值化算法定量评估方法,其特征在于所述背景点是指二值图像中灰度值为0的像素点。
4.根据权利要求1所述的基于权值模型的灰度图像二值化算法定量评估方法,其特征在于所述前景点权值矩阵中的权值生成规则是:
①当灰度图像A中某一像素点被认为是前景点的次数等于所选灰度图像二值化算法数目的一半时,则前景点权值矩阵中与所述某一像素点的坐标对应位置的权值为0;
②当灰度图像A中某一像素点被认为是前景点的次数大于所选灰度图像二值化算法数目的一半时,则前景点权值矩阵中与所述某一像素点的坐标对应位置的权值为:
(2)
式(2)中:i,j为所述灰度图像A中某一像素点所对应的横坐标和纵坐标,
表示所述灰度图像A中某一像素点被认为是前景点的次数,
表示所选择的二值化算法的数目,为7~20的自然数,
表示所选择的二值化算法数目的修正值,
当所选择的二值化算法数目为奇数时,
当所选择的二值化算法数目为偶数时,
③当灰度图像A中某一像素点被认为是前景点的次数小于所选灰度图像二值化算法数目的一半时,则前景点权值矩阵中与所述某一像素点的坐标对应位置的权值为:
(3)
式(3)中:i,j为所述灰度图像A中某一像素点所对应的横坐标和纵坐标,
表示所述灰度图像A中某一像素点被认为是前景点的次数,
表示所选择的二值化算法数目为7~20的自然数)的修正值,
当所选择的二值化算法数目为奇数时,
当所选择的二值化算法数目为偶数时,
5.根据权利要求1所述的基于权值模型的灰度图像二值化算法定量评估方法,其特征在于所述背景点权值矩阵中的权值生成规则是:
①当灰度图像A中某一像素点被认为是背景点的次数等于所选灰度图像二值化算法数目的一半时,则背景点权值矩阵中与所述某一像素点的坐标对应位置的权值为0;
②当灰度图像A中某一像素点被认为是背景点的次数大于所选灰度图像二值化算法数目的一半时,则背景点权值矩阵中与所述某一像素点的坐标对应位置的权值为:
(4)
式(4)中:i,j为所述灰度图像A中某一像素点所对应的横坐标和纵坐标,
表示所述灰度图像A中某一像素点被认为是背景点的次数,
表示所选择的二值化算法的数目,为7~20的自然数,
表示所选择的二值化算法数目的修正值,
当所选择的二值化算法数目为奇数时,
当所选择的二值化算法数目为偶数时,
③当灰度图像A中某一像素点被认为是背景点的次数小于所选灰度图像二值化算法数目的一半时,则背景点权值矩阵中与所述某一像素点的坐标对应位置的权值为:
(5)
式(5)中:i,j为所述灰度图像A中某一像素点所对应的横坐标和纵坐标,
表示所述灰度图像A中某一像素点被认为是背景点的次数,
表示所选择的二值化算法数目为7~20的自然数)的修正值,
当所选择的二值化算法数目为奇数时,
当所选择的二值化算法数目为偶数时,
6.根据权利要求1所述的基于权值模型的灰度图像二值化算法定量评估方法,其特征在于所述矩阵转换规则Ⅰ是:
①确定二值图像灰度值矩阵中灰度值为255的位置,将转换矩阵Ⅰ中与二值图像灰度值矩阵中灰度值为255的相同位置的数值置1;
②确定二值图像灰度值矩阵中灰度值为0的位置,将转换矩阵Ⅰ中与二值图像灰度值矩阵中灰度值为0的相同位置的数值置0。
7.根据权利要求1所述的基于权值模型的灰度图像二值化算法定量评估方法,其特征在于所述矩阵转换规则Ⅱ是:
①确定二值图像灰度值矩阵中灰度值为0的位置,将转换矩阵Ⅱ中与二值图像灰度值矩阵中灰度值为0的相同位置的数值置1;
②确定二值图像灰度值矩阵中灰度值为255的位置,将转换矩阵Ⅱ中与二值图像灰度值矩阵中灰度值为255的相同位置的数值置0。
CN201410324856.4A 2014-07-09 2014-07-09 基于权值模型的灰度图像二值化算法定量评估方法 Active CN104063874B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410324856.4A CN104063874B (zh) 2014-07-09 2014-07-09 基于权值模型的灰度图像二值化算法定量评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410324856.4A CN104063874B (zh) 2014-07-09 2014-07-09 基于权值模型的灰度图像二值化算法定量评估方法

Publications (2)

Publication Number Publication Date
CN104063874A true CN104063874A (zh) 2014-09-24
CN104063874B CN104063874B (zh) 2017-01-25

Family

ID=51551566

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410324856.4A Active CN104063874B (zh) 2014-07-09 2014-07-09 基于权值模型的灰度图像二值化算法定量评估方法

Country Status (1)

Country Link
CN (1) CN104063874B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105608708A (zh) * 2016-03-02 2016-05-25 浙江宇视科技有限公司 基于分块平面拟合的图像二值化方法及装置
CN109255277A (zh) * 2018-08-31 2019-01-22 阿里巴巴集团控股有限公司 一种二维码解析方法及装置
CN112415879A (zh) * 2020-11-25 2021-02-26 奥比中光科技集团股份有限公司 一种结构光场调控方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101577003B (zh) * 2009-06-05 2011-08-17 北京航空航天大学 一种基于改进交叉视觉皮质模型的图像分割方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105608708A (zh) * 2016-03-02 2016-05-25 浙江宇视科技有限公司 基于分块平面拟合的图像二值化方法及装置
CN109255277A (zh) * 2018-08-31 2019-01-22 阿里巴巴集团控股有限公司 一种二维码解析方法及装置
CN112415879A (zh) * 2020-11-25 2021-02-26 奥比中光科技集团股份有限公司 一种结构光场调控方法及系统
CN112415879B (zh) * 2020-11-25 2022-05-24 奥比中光科技集团股份有限公司 一种结构光场调控方法及系统

Also Published As

Publication number Publication date
CN104063874B (zh) 2017-01-25

Similar Documents

Publication Publication Date Title
CN105956582B (zh) 一种基于三维数据的人脸识别系统
Yousefi Image binarization using Otsu thresholding algorithm
CN109215034B (zh) 一种基于空间金字塔掩盖池化的弱监督图像语义分割方法
Poletti et al. A review of thresholding strategies applied to human chromosome segmentation
CN108022233A (zh) 一种基于改进型Canny算子的工件边缘提取方法
CN112017164A (zh) 基于深度阈值卷积模型的土石料级配检测方法
CN104636721B (zh) 一种基于轮廓与边缘纹理特征融合的掌纹识别方法
CN109461141A (zh) 一种工件缺胶检测方法
CN107767387B (zh) 基于可变感受野尺度全局调制的轮廓检测方法
CN105741279A (zh) 基于粗糙集的快速抑制模糊聚类的图像分割方法
CN109472792A (zh) 结合局部熵的局部能量泛函与非凸正则项的图像分割方法
CN110378924A (zh) 基于局部熵的水平集图像分割方法
CN111754538B (zh) Usb表面缺陷检测的阈值分割方法
CN102298773A (zh) 一种形状自适应的非局部均值去噪方法
CN105160355A (zh) 一种基于区域相关和视觉单词的遥感图像变化检测方法
CN107038432B (zh) 基于频率信息的指纹图像方向场提取方法
CN109859187B (zh) 一种爆堆矿岩颗粒图像分割方法
CN108090492B (zh) 基于尺度线索抑制的轮廓检测方法
CN106600610B (zh) 一种fcm图像分割方法及装置
CN103971362B (zh) 基于直方图和精英遗传聚类算法的sar图像变化检测
CN104050472A (zh) 一种灰度图像二值化的自适应全局阈值方法
CN104063874A (zh) 基于权值模型的灰度图像二值化算法定量评估方法
CN105205807A (zh) 基于稀疏自动编码机的遥感图像变化检测方法
CN103279960A (zh) 一种基于x射线背散射图像的人体隐藏物的图像分割方法
CN103093241B (zh) 基于同质化处理的光学遥感图像非匀质云层判别方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20220819

Address after: 2-518, Innovation Building, Huashi Science and Technology Park, East Lake High-tech Zone, Wuhan City, Hubei Province 430073

Patentee after: WUHAN LANSHAN TECHNOLOGY CO.,LTD.

Address before: 430081 construction of Qingshan District, Hubei, Wuhan

Patentee before: WUHAN University OF SCIENCE AND TECHNOLOGY

TR01 Transfer of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Quantitative evaluation method for grayscale image binarization algorithm based on weight model

Granted publication date: 20170125

Pledgee: Guanggu Branch of Wuhan Rural Commercial Bank Co.,Ltd.

Pledgor: WUHAN LANSHAN TECHNOLOGY CO.,LTD.

Registration number: Y2024980003781

PE01 Entry into force of the registration of the contract for pledge of patent right