1.一种数字印刷印品缺陷检测方法，其特征在于，包括如下步骤：

1）针对相应的数字印刷设备，求数字图像像素值到印后色值的映射关系；

2）在待印刷的彩色数字图像上抽样子块，求预期色值矩阵和实际的印刷色值矩阵；

3）比较预期色值矩阵和实际的印刷色值矩阵，并设置阈值，来判定印品质量合格与否。

2.根据权利要求1所述的一种数字印刷印品缺陷检测方法，其特征在于，所述步骤1）的具体步骤如下：

1.1）设计一幅总色阶行数为n的单色色阶图：设色阶图颜色深度为N，则 ，N取8或16；

1.2）依次读取单色色阶图中各色阶的色值，建立n×1阶矩阵S；设数字印刷设备的颜色通道数为m，建立元素全为1的1×m阶矩阵M；令矩阵 ；其中a1是n×m阶的矩阵，并且所有列向量都相同；

1.3）将单色色阶图在数字印刷设备的每个独立颜色通道上依次分别印刷，得到相应的m张色阶印品；

1.4）使用扫描仪依次扫描m张色阶印品，得到m个色阶印品扫描数字图像；

1.5）依次读取m个色阶印品扫描数字图像上各色阶的色值，得到n×m阶矩阵b1；b1中第x行第y列处数值表示第x个色阶的色值在第y个颜色通道上印刷出后，经过扫描得到的像素平均值；

1.6）读取矩阵a1和b1的第0列 和 ，使用常规的多项式插值操作或曲线拟合方法，得到映射关系 , 其中 是印前数字图像分色到颜色通道0的像素值，是印品扫描数字图像上相应位置处的颜色通道0上的色值；重复该操作，依次逐列读取矩阵a1和b1每一列，得到m个颜色通道上、从数字图像像素值到印后色值的映射关系：。

3.根据权利要求1所述的一种数字印刷印品缺陷检测方法，其特征在于，所述步骤2）的具体步骤如下：

2.1）记数字印刷设备待印刷的彩色数字图像的长为L像素、宽为W像素，并在数字图像上任选一个长为l像素宽为W像素的子图，记子图距离数字图像左边缘为i像素；

2.2）沿宽度方向上将子图连续划分为宽度都为w像素的子块，子块总数k是W/w的整数部分；

2.3）记数字印刷设备待印刷的彩色数字图像的长为L像素、宽为W像素，待印刷数字图像的每个子块都是彩色的，将每个子块按数字印刷设备的要求进行m个通道的分色处理，得到k×m个单色子块，分别计算各单色子块全部像素值的平均值，建立k×m阶的矩阵 ，矩阵中第x行第y列处数值是第x个子块分色到第y个颜色通道上的像素值的平均值；

2.4）应用求得的数字印刷设备m个颜色通道上、从数字图像像素值到印后色值的映射关系 ，将矩阵 中第x行第y列处数值代入 ，求得该子块在该通道上的预期印后扫描色值；同理，求出全部k个子块在m个通道上的预期印后扫描色值，建立k×m阶预期色值矩阵a2，a2中第x行第y列处数值表示第x个子块分色到第y个颜色通道上的像素值预期的经过印刷后再扫描回的色值；

2.5）完整的数字图像经过分色后在数字印刷设备的m个颜色通道上依次分别印刷，使用扫描仪扫描各通道印刷出的印品，得到m个长为L像素、宽为W像素的单色扫描数字图像，或者使用扫描仪扫描已印刷完成的彩色印品，得到彩色扫描数字图像，再经与印前相同的分色处理，得到m个长为L像素、宽为W像素的单色扫描数字图像；

2.6）在单色扫描数字图像上定位到距左边缘i像素位置，选取一个长为l像素宽为W像素的扫描子图，沿宽度方向上将扫描子图连续划分为宽度都为w像素的扫描子块，扫描子块总数是W/w的整数部分，扫描子块与子块总数相同，位置一一对应；读出每个扫描子块上的全部像素值并计算平均值，建立k×m阶的实际印刷色值矩阵b2，b2的第 x行第y列处数值表示数字图像经过该数字印刷设备实际印刷所得的印品上，距左边缘i像素、第x个子块位置，第y个颜色通道上的扫描仪扫描所得像素值的平均值。

4.根据权利要求1所述的一种数字印刷印品缺陷检测方法，其特征在于，所述步骤3）的具体步骤如下：阈值是给定的某一个确定的常数值Y，将预期色值矩阵a2和实际的印刷色值矩阵b2作差，得到的矩阵： ，分析矩阵Q,根据以下规则判定印品质量合格与否：

若矩阵Q中有P个值>Y，则该印品质量不合格；反之，合格；

求矩阵Q的均值 ， ，若该均值 ，则该印品质量不合格；反之，合格；

求矩阵Q的方差 :

若该方差 ，则该印品质量不合格；反之，合格；

求矩阵Q的标准差s:

若该标准差s>Y，则该印品的质量不合格；反之，合格；

若矩阵Q中的任意一个值|x|>Y，则该印品不合格；反之，合格；

用极差判定，极差是用来表示统计资料中的变异量数，其最大值与最小值之间的差距；

即最大值减最小值后所得之数据， ；若矩阵Q中的x>Y， 则该印品不合格；反之，合格。

5.根据权利要求1所述的一种数字印刷印品缺陷检测方法，其特征在于，所述步骤3）的具体步骤如下：阈值是给定的某一个确定的常数值Y，对于预期色值矩阵a2和实际的印刷色值矩阵b2，根据以下规则判定印品质量合格与否：求出矩阵a2的平均值 和矩阵b2的平均值 ,作差，若 ，则该印品不合格；反之，合格；

分别求出矩阵a2、b2各自的方差 、 ，将其作差，若 , 则该印品不合格；

反之，合格；

分别求出矩阵a2、b2各自的标准差 、 ，将其作差，若 , 则该印品不合格；反之，合格；

用极差判定，用矩阵a2中的最大值与矩阵b2中的最小值相减后，即 ，若该差值 , 则该印品不合格；反之，合格。

6.根据权利要求1所述的一种数字印刷印品缺陷检测方法，其特征在于，所述步骤3）的具体步骤如下：阈值是给定的某一个确定的常数值Y，根据预期色值矩阵a2和实际的印刷色值矩阵b2的均方误差MSE/均方根误差RMSE判定印品质量合格与否；

假设预期色值矩阵 ，

实际的印刷色值矩阵

则：

若 , 则该印品不合格；反之，合格。

7.根据权利要求1所述的一种数字印刷印品缺陷检测方法，其特征在于，所述步骤3）的具体步骤如下：阈值是给定的某一个确定的常数值Y，根据预期色值矩阵a2和实际的印刷色值矩阵b2的协方差判定印品质量合格与否；

在概率论和统计学中，协方差用于衡量两个变量的总体误差，而方差是协方差的一种特殊情况，即当两个变量是相同的情况；期望值分别为 与的两个实数随机变量X与Z之间的协方差定义为：

，它表示两组数据的相关性；

若 , 则该印品不合格；反之，合格。

8.根据权利要求1所述的一种数字印刷印品缺陷检测方法，其特征在于，所述步骤3）的具体步骤如下：阈值是给定的某一个确定的常数值Y，根据预期色值矩阵a2和实际的印刷色值矩阵b2的峰值信噪比判定印品质量合格与否：峰值信噪比PSNR表示信号最大可能功率和影响它的表示精度的破坏性噪声功率的比值；由于许多信号都有非常宽的动态范围，峰值信噪比常用对数分贝单位来表示；

其中，MAX表示图像颜色的最大数值，8bit图像最大取值为255；

若 ，则该印品不合格；反之，合格。

9.根据权利要求1所述的一种数字印刷印品缺陷检测方法，其特征在于，所述步骤3）的具体步骤如下：阈值是设定的某个常数矩阵C，根据矩阵之间的差值判定印品质量合格与否；

设阈值矩阵 ，矩阵 ，

将矩阵Q和矩阵C的对应位置一一对比，若 组数值对比中，有超过一半即 组数值大于阈值矩阵C中对应位置的数值，则该印品不合格；反之，合格。

10.根据权利要求1所述的一种数字印刷印品缺陷检测方法，其特征在于，所述步骤3）的具体步骤如下：阈值是设定的某个常数矩阵C，根据矩阵之间的平均值判定印品质量合格与否；

将预期色值矩阵a2和实际的印刷色值矩阵b2的对应位置上分别求均值，可得：, 阈值矩阵

, 将矩阵D和矩阵C的对应位置一一对比，若 组数值对比中，有超过一半即 组数值大于阈值矩阵C中对应位置的数值，则该印品不合格；反之，合格。