基于直方图凹度分析的印刷网点图二值化算法研究.docVIP

基于直方图凹度分析的印刷网点图二值化算法研究 　　摘要： 在印刷工业上，为了用数字图像处理的方法测量印刷网点的参数，二值化是图像预处理过程中关键的一步，直接影响参数测量的精度。在分析已有数字图像二值化方法的基础上，结合印刷网点图的特征，提出一种基于直方图凹度分析的二值化方法。该算法根据印刷网点图像的灰度直方图特征，利用数学求导的思想求出灰度直方图中双峰之间的谷所对应的灰度值，将该值作为二值化的阈值。实验证明，该方法可以比较准确地求出印刷网点图的二值化阈值，和已有算法相比，扫描率低、速度快、计算准确，具有很好的应用前景。 　　关键词： 二值化； 阈值； 灰度直方图； 数学求导 　　引言在印刷工业中，印刷网点（见图1，以下简称网点）是印刷品的基本单元。它在印刷工业中起到决定印品颜色、层次和图像轮廓的作用。网点的主要特征参数有网点大小、网点间距和网点角度。能否准确选取合适的网点参数决定了印品质量的好坏。图1印刷网点放大图 　　Fig.1Enlarged print dot image传统的网点参数测量方法是采用光学密度计。随着数字图像处理技术的发展，图像处理算法理论上可以快速得到网点的参数[1]。因此，基于嵌入式平台的网点测量仪可以用图像处理的算法得到网点的参数，从而取代了传统网点测量仪器。网点图像的处理过程分为图像预处理和图像识别，而预处理的第一步是对采集到的网点图进行二值化处理。二值化是将图像的灰度值置成0或255的过程，即选择一定的灰度值作为二值化阈值，将图像中大于阈值的像素置为白色，小于阈值的像素置为黑色。二值化主要有全局阈值法、局部阈值自适应法和基于方向信息的二值化方法[2]。本文在分析上述算法的基础上，结合网点图的特点，提出一种基于直方图凹度分析的二值化阈值算法，即利用数学求导的思想计算二值化阈值。实验证明，该算法能比较准确地求出网点图像的二值化阈值。1网点测量仪的工作原理测量仪共由以下三部分组成：（1）网点图像的采集；（2）网点图像的数字图像处理和参数计算；（3）参数计算结果显示。测量仪的整体框图如图2所示。在系统中，图像的采集由CMOS图像传感器完成，网点图的处理及网点参数的计算由ARM核心板完成。当系统初始化完成后，可以通过CMOS图像传感器采集印品的显微数字图像，一帧图像数据采集完成后在屏幕上显示，图像数据则由MCU通过一定的算法处理，最后得到网点的参数：网点大小、网点间距和网点角度。网点图像处理的算法包括图像的灰度化、二值化、去噪、网点参数计算，并将结果显示在GUI上[3]。光学仪器第35卷 　　在图像处理部分，二值化是图像预处理的关键，二值化准确与否，将影响后期参数计算的精度，进而影响网点参数识别的精度。2网点图的灰度直方图特征灰度直方图是用于表达图像灰度分布情况的统计图表。其横坐标是灰度值r，纵坐标是出现这个灰度值的概率密度h（r）（对连续图像f（x，y）而言），或者是出现这个灰度的概率值h（r）（对数字图像f（m，n）而言）[4]。由图1可以看出，印刷网点图有以下特点：（1）由两种颜色组成：一种是背景色，一种是网点颜色；（2）由大小相近间距相同的点组成，并且点分布在互相垂直的两个方向；（3）不同的网点图的网点面积大小不同。运用MATLAB软件的imhist（）函数得到图3所示的灰度直方图，横坐标是灰度级别r，纵坐标是该灰度级出现的个数h（r）。图3印刷网点图的灰度直方图 　　Fig.3Gray histogram of print dot image灰度直方图能很好地反映图像的总体特征。由图3可以看出，印刷网点图的灰度直方图由明显的双峰组成，一个代表背景，一个代表目标点，即印刷网点。因此，二值化阈值如果选择在双峰之间的“谷”处，可以将网点图二值化成只有背景和目标点的黑白图像[5]。3基于直方图凹度分析求二值化阈值灰度直方图是一种离散的分布图，其包络线是一条连续的曲线，因此对直方图的“谷”的寻找便可以转化为求其包络线h（r）极小值的问题。根据高等数学函数极值的有关问题，曲线h（r）中的极小值点应同时满足下列条件：h（r）r=0（1） 　　2h（r）r2=0（2）满足上述条件的极小值点所对应的灰度级即可作为分割阈值[6]。但是，上述方法需要用到数值逼近等计算，算法复杂度高，不利于在嵌入式平台中的实现。为了在离散的直方图上寻找到“谷”所对应的灰度值，采用以下离散求导的方法来确定二值化阈值。 　　3.1确定阈值的灰度级范围一般灰度直方图起止处灰度级的像素个数很少，将像素个数h（r）小于一定值的灰度忽略（即h（r）置0）。这样做的目的在于减小这些小量像素灰度级对于整体阈值计算的影响，从而提高阈值的准确度。设经过上述处理后像素个数不为0的灰度级范围为m～n，灰度值r对应的像素个数为h（r）。首