CN102426171A - 一种基于小波理论的印刷品墨斑测量方法 - Google Patents

Abstract

本发明公开了一种基于小波理论的印刷品墨斑的测量方法，该方法按照以下方法进行印刷品墨斑的测量，通过扫描仪进行印刷墨斑的图像采集，按照预定的尺寸进行图像分析取样，利用不同的小波基进行多级小波分解图像，将图像分解为高频部分和低频部分，并除去高频噪声部分；采用离散小波逆变换对墨斑图像进行还原，得到印刷墨斑平滑图像，并计算平滑墨斑图像的灰度值的变化，即墨斑大小的程度。本发明基于小波理论的印刷品墨斑的测量方法采用小波理论来测量印刷墨斑，利用小波多级分解、离散小波逆变换和图像的变化，最终实现印刷墨斑的测量，方法简单，方便、精度高，很好解决了印刷品墨斑难测量及评价的问题。

Description

一种基于小波理论的印刷品墨斑测量方法

技术领域

本发明属于印刷品质量检测技术领域，涉及一种印刷品墨斑测量的方法，具体涉及一种基于小波理论的印刷品墨斑测量方法。

背景技术

印刷品墨斑，一般定义为：在印刷品的所有方向上，密度的空间频率小于0.4周期/mm的无周期性的波动，尤其是在实地区域由于油墨分布不均匀而形成的色彩不正常，深浅不一致的斑点和条纹。印刷墨斑是影响印刷品质量的重要因素之一，也是多色印刷中最常见的印刷缺陷。

目前对印刷墨斑的测量方法主要有ISO印刷墨斑测量法，比参数墨斑算法和人工视觉判断法。但这些方法对印刷墨斑进行分析及评价，仍然存在着相关的不足：ISO印刷墨斑法仅是从几何计算的方法对墨斑进行简单计算，且算法比较粗造，精度比较差；比参数墨斑测量法通过测量亮斑和暗斑的周长及其密度对比来表达印刷墨斑的程度，但是此方法由于未考虑噪音的影响及在50％的阈值处周长与人眼视觉评价时的不均匀性差，使得与人眼的相关性差，同时也降低了墨斑测量的精度；而人眼视觉判断法，是靠人工去逐一观察，十分费时，费力，而且因人的视觉差异，判断结果不太准确。

发明内容

本发明的目的是提供一种基于小波理论的印刷品墨斑测量方法，解决了现有测量方法精度低，计算方法粗糙，费时及费力的问题，同时考虑墨斑图像的噪声及将墨斑图像时域和频域联系来，精确的评价了印刷品墨斑的程度。

本发明所采用的技术方案是，一种基于小波理论的印刷品墨斑测量方法，具体按照以下步骤实施：

步骤1：用扫描仪采集印刷图像；

步骤2：按照取样图像尺寸为4cm-6cm，选择印刷图像上需要测量的部位，即为印刷墨斑图像；

步骤3：对步骤2得到的印刷墨斑图像，采用不同的小波基进行不同级数的分解，将印刷墨斑图像分解成高频图像和低频图像两部分，去除高频噪声，并保留低频图像；

步骤4：对步骤3得到的去除噪声后的图像通过二维离散平稳小波逆变换，进行图像重构，得到印刷墨斑的平滑墨斑图像；

步骤5：计算步骤4得到的平滑墨斑图像灰度值的变化，即墨斑大小PM：

$\mathrm{PM}=\underset{x,y}{\mathrm{\Σ}}{(I(x,y)-\underset{x,y}{\mathrm{\Σ}}I(x,y)/L)}^2/{100\·L}^2\·\underset{x,y}{\mathrm{\Σ}}I(x,y),$

其中：(x，y)为平滑图像的纵横坐标，I(x，y)为平滑图像的灰度值，L为离散图像矩阵的大小，即得印刷品的墨斑的大小，根据以下标准判断印刷品上是否有墨斑：在印刷品的所有方向上，比较PM与0.4的大小，小于0.4则有墨斑，大于0.4则没有墨斑，其中图像尺寸的单位为mm。

本发明的特点还在于，

其中的步骤3中采用不同的小波基进行不同级数的分解，具体按照以下步骤实施：分解级数J选为1到4级，小波基函数选为HAAR，SYM4，SYM5，BIOR3.7，DB1，DB3，DB4，DB5，首先选定特定的分解级数和小波基函数，利用Matlab软件小波工具箱中的二维离散平稳小波函数Swt2进行分解，分解完后，得到印刷墨斑图像的高频部分和低频部分，然后利用Size函数计算高频部分矩阵的大小，利用Zeros函数将高频矩阵信息置0，从而去除墨斑图像的高频部分。

其中的步骤4中的二维离散平稳小波逆变换，具体按照以下步骤实施：利用软件小波工具箱中的二维离散平稳小波逆变换函数Iswt2进行图像重构，得到印刷墨斑平滑图像。

本发明的有益效果是，本发明方法利用小波理论，采用不同的小波基函数对印刷墨斑图像进行多级分解及小波逆变换，计算平滑去噪的印刷品墨斑图像灰度值的变化，最终简单、方便、高精度的测量出印刷品墨斑。

附图说明

图1是本发明基于小波理论的印刷品墨斑测量方法的计算流程图。

具体实施方式

下面结合附图和具体实施方式对本发明进行详细说明。

本发明基于小波理论的印刷品墨斑测量方法，利用计算机存储图像并通过墨斑算法计算印刷品的墨斑，对于一个印刷图像，测量图像有没有墨斑及墨斑指数的大小，在印刷行业中，也以此来评价印刷质量。如图1所示，具体按照以下步骤实施：

(1)设置扫描仪参数，包括亮度、来源、扫描模式、分辨率及对比等参数，其中分辨率最小应该为1200DPI。其他参数的设置，必须使得能够采集到足够清晰的图像。用扫描仪采集印刷图像。

(2)印刷墨斑图像取样，选择图像上需要测量部位。为了准确测量墨斑，取样图像大小应该与图像分辨率相匹配，取样图像尺寸一般为4cm-6cm。

(3)采用不同的小波基进行不同级数的分解，将墨斑图像I分解成高频图像和低频图像两部分，去除高频噪声，并保留低频图像。分解级数J可选为1到4级，小波基函数可选为HAAR，SYM4，SYM5，BIOR3.7，DB1，DB3，DB4，DB5。具体算法步骤为：首先选定特定的分解级数和小波基函数，利用Matlab软件小波工具箱(WaveletToolbox)中的二维离散平稳小波函数Swt2进行分解，分解完后，将会得到印刷墨斑图像的高频部分和低频部分。然后利用Size函数计算高频部分矩阵的大小，利用Zeros函数将高频矩阵信息置0，从而去除墨斑图像的高频部分。

(4)对去除噪声的图像进行二维离散平稳小波逆变换，进行图像重构，得到印刷墨斑的平滑墨斑图像。具体步骤为：利用软件小波工具箱(Wavelet Toolbox)中的二维离散平稳小波逆变换函数Iswt2进行图像重构，便可得到印刷墨斑平滑图像I。

(5)计算平滑墨斑图像的灰度值的变化，即墨斑大小PM：

$\mathrm{PM}=\underset{x,y}{\mathrm{\Σ}}{(I(x,y)-\underset{x,y}{\mathrm{\Σ}}I(x,y)/L)}^2/{100\·L}^2\·\underset{x,y}{\mathrm{\Σ}}I(x,y)---\left(1\right)$

其中：(x，y)为平滑图像的纵横坐标，I(x，y)为平滑图像的灰度值，L为离散图像矩阵的大小，即得印刷品的墨斑的大小。根据以下标准判断印刷品上是否有墨斑：在印刷品的所有方向上，比较PM与0.4的大小，小于0.4则有墨斑，大于0.4则没有墨斑。

实施例

现以新闻纸的印刷墨斑为例，具体来说明本发明测量墨斑的方法及过程。

(1)调整扫描仪参数，放置好新闻纸印刷墨斑样张进行图像的正确采集。扫描仪参数设置如下：亮度50，来源选为一般，扫描模式为黑白照片，分辨率为1200DPI。

(2)取样新闻纸印刷图像，图像大小为4.27cm*4.27cm，图像像素点离散为2048*2048，图像的灰度等级为256，图像数量为70副。

(3)对每一副图像进行多尺度离散小波变换，采用不同的小波基进行不同级数的分解，将印刷墨斑图像分解成高频图像和低频图像两部分，去除高频噪声，并保留低频图像，并得到每副图像的低频系数及水平方向、垂直方向和对角方向的高频系数。具体算法为：

选定特定的分解级数J和小波基函数

利用Matlab软件小波工具箱(Wavelet Toolbox)中的二维离散平稳小波函数Swt2进行分解，Swt2函数的具体格式为：

Swt2函数的返回值便是印刷图像的低频部分A，高频部分的水平方向系数H，垂直方向系数V及对角方向系数D。从而得到印刷墨斑图像的高频部分和低频部分。

(4)将每一级上的高频系数去除，去除噪声，利用二维离散平稳小波逆变换函数，将去噪图像的低频系数还原，重构新的新闻纸墨斑平滑图像。具体算法为：

利用Size函数计算高频部分矩阵的大小，Size(D)，通过Zeros函数将高频矩阵信息置0，H＝V＝D＝Zeros(Size(D))，从而去除墨斑图像的高频部分；

利用软件小波工具箱(Wavelet Toolbox)中的二维离散平稳小波逆变换函数Iswt2进行图像重构，其具体格式为：

便可得到印刷墨斑平滑图像I。

(5)按照

$\mathrm{PM}=\underset{x,y}{\mathrm{\Σ}}{(I(x,y)-\underset{x,y}{\mathrm{\Σ}}I(x,y)/L)}^2/L^2\·\underset{x,y}{\mathrm{\Σ}}I(x,y)---\left(2\right)$

计算新闻纸墨斑平滑图像的墨斑，即得到该印刷品的墨斑值的大小。

(6)为了检验该方法的正确性，是否符合人眼的视觉评价系统，将本发明的方法跟人眼的视觉系统做了相关性分析，结果如下：

表1新闻纸墨斑与人眼视觉评价系统相似度

表1中数据表示采用该方法分析的新闻纸墨斑程度与人眼视觉系统的相似度。从结果表明，采用HAAR，SYM4，SYM5，BIOR3.7，DB1，DB3，DB4，DB5小波基并进行1-3级的分解，最终测量墨斑结果跟人眼有很好的相似度，可以代替人眼进行测量。其中对于新闻纸墨斑采用1-3级和DB3HAAR，SYM4，SYM5，BIOR3.7，DB1，DB3，DB4，DB5小波基进行测量，效果比较好。

除了新闻纸墨斑的评价外，将测量方法应用到喷墨纸，双胶纸上及铜版纸上印刷墨斑的测量，最终测量结果跟人眼有很好的相似度，反映了印刷墨斑的实际情况，最佳的处理参数应为下表2：

表2最佳参数值

Claims (3)

1.一种基于小波理论的印刷品墨斑测量方法，其特征在于，具体按照以下步骤实施：

步骤1：用扫描仪采集印刷图像；

步骤2：按照取样图像尺寸为4cm-6cm，选择印刷图像上需要测量的部位，即为印刷墨斑图像；

步骤3：对步骤2得到的印刷墨斑图像，采用不同的小波基进行不同级数的分解，将印刷墨斑图像分解成高频图像和低频图像两部分，去除高频噪声，并保留低频图像；

步骤4：对步骤3得到的去除噪声后的图像通过二维离散平稳小波逆变换，进行图像重构，得到印刷墨斑的平滑墨斑图像；

步骤5：计算步骤4得到的平滑墨斑图像的灰度值的变化，即墨斑大小PM：

$\mathrm{PM}=\underset{x,y}{\mathrm{\Σ}}{(I(x,y)-\underset{x,y}{\mathrm{\Σ}}I(x,y)/L)}^2/L^2\·\underset{x,y}{\mathrm{\Σ}}I(x,y),$

其中：(x，y)为平滑图像的纵横坐标，I(x，y)为平滑图像的灰度值，L为离散图像矩阵的大小，即得印刷品的墨斑的大小，根据以下标准判断印刷品上是否有墨斑：在印刷品的所有方向上，比较PM与0.4的大小，小于0.4则有墨斑，大于0.4则没有墨斑，其中图像尺寸的单位为mm。

2.根据权利要求1所述的基于小波理论的印刷品墨斑测量方法，其特征在于，所述的步骤3中采用不同的小波基进行不同级数的分解，具体按照以下步骤实施：分解级数J选为1到4级，小波基函数选为HAAR，SYM4，SYM5，BIOR3.7，DB1，DB3，DB4，DB5，首先选定特定的分解级数和小波基函数，利用Matlab软件小波工具箱中的二维离散平稳小波函数Swt2进行分解，分解完后，得到印刷墨斑图像的高频部分和低频部分，然后利用Size函数计算高频部分矩阵的大小，利用Zeros函数将高频矩阵信息置0，从而去除墨斑图像的高频部分。

3.根据权利要求1所述的基于小波理论的印刷品墨斑测量方法，其特征在于，所述的步骤4中的二维离散平稳小波逆变换，具体按照以下步骤实施：利用软件小波工具箱中的二维离散平稳小波逆变换函数Iswt2进行图像重构，得到印刷墨斑平滑图像。

Images