CN103595894B - 打印机图像半色调输出算法质量评价方法 - Google Patents

打印机图像半色调输出算法质量评价方法 Download PDF

Info

Publication number
CN103595894B
CN103595894B CN201310626343.4A CN201310626343A CN103595894B CN 103595894 B CN103595894 B CN 103595894B CN 201310626343 A CN201310626343 A CN 201310626343A CN 103595894 B CN103595894 B CN 103595894B
Authority
CN
China
Prior art keywords
halftoning
value
luminance factor
munsell
color lump
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
CN201310626343.4A
Other languages
English (en)
Other versions
CN103595894A (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 University WHU
Original Assignee
Wuhan University WHU
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 WHU filed Critical Wuhan University WHU
Priority to CN201310626343.4A priority Critical patent/CN103595894B/zh
Publication of CN103595894A publication Critical patent/CN103595894A/zh
Application granted granted Critical
Publication of CN103595894B publication Critical patent/CN103595894B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了打印机图像半色调输出算法质量评价方法，通过将输出的半色调图像测量值转换为孟塞尔明度值，再对孟塞尔明度值与其对应的灰度值进行线性拟合，并使用误差项平方和、均方根误差、回归系数和调整后的回归系数四个指标进行一元线性回归分析，从而评价半色调算法对激光打印机非线性特性的鲁棒性。本发明达到了评价该半色调算法对激光打印机非线性特性的鲁棒性的目的，为进一步修改半色调算法及建立打印机模型提供基础。

Description

打印机图像半色调输出算法质量评价方法
技术领域
本发明属于数字图像处理领域，涉及一种输出半色调结果检查方法，具体涉及一种打印机图像半色调输出算法质量评价方法。
背景技术
由于打印机的非线性特性，半色调算法处理得到的半色调结果输出时，输出的半色调结果一般都会存在阶调非线性再现的情况。因此，半色调算法应当对打印机的非线性特性具有较强的鲁棒性。为了评价该半色调算法对激光打印机的鲁棒性，需要对打印机输出的半色调结果的线性程度进行评价，从而为进一步修改半色调算法及建立打印机模型提供基础。
目前的半色调算法评价方法，通常只是在半色调结果打印输出前，使用均方误差、峰值信噪比、通用质量评价参数等数字图像评价方法对其评价。但半色调结果打印输出后，这些数字图像评价方法均无法直接对输出半色调结果线性程度进行评价。
发明内容
为了解决上述的技术问题，本发明提出一种打印机图像半色调输出算法质量评价方法，使用色度仪或密度仪等测量设备对输出半色调结果进行测量，将测量结果转换为孟塞尔明度级，并对转换得到的孟塞尔明度级与其对应的灰度级进行一元线性回归分析，以评价半色调算法对激光打印机的鲁棒性。
本发明所采用的技术方案是：一种打印机图像半色调输出算法质量评价方法，其特征在于，包括以下步骤：
步骤1：设计不止一个灰度色块，其中，针对所述的灰度色块，用0表述黑，用1表述白，且其灰度值以0.05为增长步长从0开始逐一增长；
步骤2：使用半色调算法对各个灰度色块的中性灰色块进行半色调处理；
步骤3：使用打印机对半色调处理结果打印输出；
步骤4：对打印输出结果进行测量；
步骤5：将测量结果转换为CIE1931XYZ色度系统中的亮度因数Y；
步骤6：将亮度因数Y转换为孟塞尔明度V，根据亮度因数Y与孟塞尔明度V的函数关系，计算得到孟塞尔明度V；
步骤7：通过最小二乘法拟合转换得到的孟塞尔明度V与其对应的原始灰度色块的灰度值的线性关系；
步骤8：使用误差项平方和、均方根误差、决定系数和调整后的决定系数四个指标对孟塞尔明度V与其对应的原始灰度色块的灰度值线性程度进行评价。
作为优选，步骤1所述的设计不止一个灰度色块，其灰度色块的数量为20个。
作为优选，步骤4所述的对打印输出结果进行测量，是使用色度计或密度计对打印输出结果进行测量。
作为优选，步骤5所述的将测量结果转换为CIE1931XYZ色度系统中的亮度因数Y，若使用色度计对打印输出结果进行测量，则测量得到打印输出结果的亮度因数Y；若使用密度计对打印输出结果进行测量，则根据光学反射密度定义公式推导出的密度D与亮度因数Y的函数关系，计算出亮度因数Y：其中，Y0为孟塞尔系统第10级明度的亮度因数，Y0=102.75，ρ0为孟塞尔系统第10级明度的反射率，ρ0=1，D为某区域的光学反射密度，Y为该区域的亮度因数。
作为优选，步骤6所述的根据亮度因数Y与孟塞尔明度V的函数关系，计算孟塞尔明度V，其计算公式为：V=2.217Y0.352-1.324。
本发明相对于现有技术的创新点在于，将测量数据转换为孟塞尔明度值，通过对孟塞尔明度与对应灰度的线性程度评价半色调输出结果的线性程度。从而达到评价该半色调算法对激光打印机非线性特性的鲁棒性的目的，为进一步修改半色调算法及建立打印机模型提供基础。
附图说明
图1：是本发明的打印机图像半色调输出算法质量评价方法的流程图。
图2：是本发明实施例的待处理的20个灰度色块。
图3：是本发明实施例的亮度因数Y与孟塞尔明度V的转换关系示意图。
图4：是本发明实施例的20个灰度色块使用迟滞系数为1.5的绿噪声半色调方法进行半色调处理的结果图像。
图5：是本发明实施例的20个灰度色块通过迟滞系数为3.0的绿噪声半色调方法进行半色调处理后的结果图像。
图6：是本发明实施例的迟滞系数为1.5的绿噪声半色调方法所做的线性拟合示意图。
图7：是本发明实施例的迟滞系数为3.0的绿噪声半色调方法所做的线性拟合示意图。
具体实施方式
以下结合附图和具体实施例对本发明做进一步的阐述。
下面以待评价半色调方法为迟滞系数分别为1.5和3.0的绿噪声半色调方法、打印机为Pantum2000激光打印机、测量设备为Eyeone色度测量仪为例对本发明作进一步地描述。
请见图1，本发明所采用的技术方案是：一种打印机图像半色调输出算法质量评价方法，其特征在于，包括以下步骤：
步骤1：请见图2，本实施例首先设计20个灰度色块，其中，针对所述的灰度色块，用0表述黑，用1表述白，且其灰度值以0.05为增长步长从0增长到0.95，纸白为1。
步骤2：分别使用迟滞系数为1.5和3.0的绿噪声半色调方法对各个灰度色块的中性灰色块进行半色调处理；
请见图4，为20个灰度色块使用迟滞系数为1.5的绿噪声半色调方法进行半色调处理的结果图像；
请见图5，为20个灰度色块通过迟滞系数为3.0的绿噪声半色调方法进行半色调处理后的结果图像。
步骤3：使用Pantum2000激光打印机对半色调处理结果打印输出。
步骤4：使用Eyeone色度计对打印输出结果进行测量。
步骤5：将测量结果转换为CIE1931XYZ色度系统中的亮度因数Y。
步骤6：将亮度因数Y转换为孟塞尔明度V，根据亮度因数Y与孟塞尔明度V的函数关系，计算孟塞尔明度V，请见图3，为亮度因数Y与孟塞尔明度V的转换关系示意图，其计算公式为：
V=2.217Y0.352-1.324。
步骤7：通过最小二乘法拟合转换得到的孟塞尔明度V与其对应的原始灰度色块的灰度值的线性关系；
请见图6，为迟滞系数为1.5的绿噪声半色调方法所做的线性拟合示意图，图中最后一个点为纸白，其横坐标为1，纵坐标为其对应的孟塞尔明度值，其他每个点的横坐标为该色块的灰度值，纵坐标为其对应的孟塞尔明度值；
请见图7，为迟滞系数为3.0的绿噪声半色调方法所做的线性拟合示意图，图中最后一个点为纸白，其横坐标为1，纵坐标为其对应的孟塞尔明度值，其他每个点的横坐标为该色块的灰度值，纵坐标为其对应的孟塞尔明度值。
步骤8：使用误差项平方和、均方根误差、决定系数和调整后的决定系数四个指标对孟塞尔明度V与其对应的原始灰度色块的灰度值线性程度进行评价。
计算得到迟滞系数为1.5的绿噪声半色调处理的输出结果的误差项平方和为7.518，均方根误差为0.6290，决定系数为0.9316，调整后的决定系数为0.9280；迟滞系数为3.0的绿噪声半色调处理的输出结果的误差项平方和为0.5667，均方根误差为0.1727，决定系数为0.9944，调整后的决定系数为0.9941。从评价结果来看，迟滞系数3.0的绿噪声半色调处理输出结果的阶调再现明显比迟滞系数为1.5的绿噪声半色调方法更为线性，其对打印机非线性特性的鲁棒性明显优于迟滞系数为1.5的绿噪声半色调方法，符合绿噪声半色调算法中迟滞系数的理论作用。
本实施例的也可以使用密度计对打印输出结果进行测量，若使用密度计对打印输出结果进行测量，则根据光学反射密度定义公式推导出的密度D与亮度因数Y的函数关系，计算出亮度因数Y：
Y = Y 0 10 D ρ 0 ,
其中，Y0为孟塞尔系统第10级明度的亮度因数，Y0=102.75，ρ0为孟塞尔系统第10级明度的反射率，ρ0=1，D为某区域的光学反射密度，Y为该区域的亮度因数。
以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围，因此，凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (5)

1.一种打印机图像半色调输出算法质量评价方法，其特征在于，包括以下步骤：
步骤1：设计不止一个灰度色块，其中，针对所述的灰度色块，用0表述黑，用1表述白，且其灰度值以0.05为增长步长从0开始逐一增长；
步骤2：使用半色调算法对各个灰度色块的中性灰色块进行半色调处理；
步骤3：使用打印机对半色调处理结果打印输出；
步骤4：对打印输出结果进行测量；
步骤5：将测量结果转换为CIE1931XYZ色度系统中的亮度因数Y；
步骤6：将亮度因数Y转换为孟塞尔明度V，根据亮度因数Y与孟塞尔明度V的函数关系，计算得到孟塞尔明度V；
步骤7：通过最小二乘法拟合转换得到的孟塞尔明度V与其对应的原始灰度色块的灰度值的线性关系；
步骤8：使用误差项平方和、均方根误差、决定系数和调整后的决定系数四个指标对孟塞尔明度V与其对应的原始灰度色块的灰度值线性程度进行评价。
2.根据权利要求1所述的打印机图像半色调输出算法质量评价方法，其特征在于：步骤1所述的设计不止一个灰度色块，其灰度色块的数量为20个。
3.根据权利要求1所述的打印机图像半色调输出算法质量评价方法，其特征在于：步骤4所述的对打印输出结果进行测量，是使用色度计或密度计对打印输出结果进行测量。
4.根据权利要求3所述的打印机图像半色调输出算法质量评价方法，其特征在于：步骤5所述的将测量结果转换为CIE1931XYZ色度系统中的亮度因数Y，
若使用色度计对打印输出结果进行测量，则测量得到打印输出结果的亮度因数Y；
若使用密度计对打印输出结果进行测量，则根据光学反射密度定义公式推导出的密度D与亮度因数Y的函数关系，计算出亮度因数Y：
Y = Y 0 10 D ρ 0 ,
其中，Y0为孟塞尔系统第10级明度的亮度因数，Y0=102.75，ρ0为孟塞尔系统第10级明度的反射率，ρ0=1，D为某区域的光学反射密度，Y为该区域的亮度因数。
5.根据权利要求1所述的打印机图像半色调输出算法质量评价方法，其特征在于：步骤6所述的根据亮度因数Y与孟塞尔明度V的函数关系，计算孟塞尔明度V，其计算公式为：
V=2.217Y0.352-1.324。
CN201310626343.4A 2013-11-28 2013-11-28 打印机图像半色调输出算法质量评价方法 Expired - Fee Related CN103595894B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310626343.4A CN103595894B (zh) 2013-11-28 2013-11-28 打印机图像半色调输出算法质量评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310626343.4A CN103595894B (zh) 2013-11-28 2013-11-28 打印机图像半色调输出算法质量评价方法

Publications (2)

Publication Number Publication Date
CN103595894A CN103595894A (zh) 2014-02-19
CN103595894B true CN103595894B (zh) 2015-10-21

Family

ID=50085865

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310626343.4A Expired - Fee Related CN103595894B (zh) 2013-11-28 2013-11-28 打印机图像半色调输出算法质量评价方法

Country Status (1)

Country Link
CN (1) CN103595894B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104010107B (zh) * 2014-05-16 2016-10-05 电子科技大学 线扫描相机图像采集系统的测试设备及方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5675253B2 (ja) * 2009-12-28 2015-02-25 キヤノン株式会社 画像処理装置、画像処理方法及びコンピュータプログラム
JP5439426B2 (ja) * 2011-04-05 2014-03-12 富士フイルム株式会社 処理条件選択装置、処理条件選択方法及びプログラム

Also Published As

Publication number Publication date
CN103595894A (zh) 2014-02-19

Similar Documents

Publication Publication Date Title
JP5924091B2 (ja) 画像処理装置、画像処理システム、プログラム
CN106541702B (zh) 一种印刷质量校正系统及其方法
US8724197B2 (en) Image processing apparatuses, methods, and computer program products for color calibration superposing different color patches in different orders of gradation values
CN102307273A (zh) 一种综合目测、icc标准和自定义标准的数码打样质量评估方法
CN103076334A (zh) 一种定量评价数字印刷线条与文本感知质量的方法
CN103347142A (zh) 一种基于数字设备呈色适性表征的优化色表实现方法
EP2899964A2 (en) Method for generating lookup table and color conversion device
CN101587589A (zh) 喷墨打印质量模糊评价方法
JP2020102807A (ja) 色予測モデル構築方法、色予測方法、および色予測モデル構築プログラム
CN103595894B (zh) 打印机图像半色调输出算法质量评价方法
CN103895377B (zh) 多维度影像印艺
Nussbaum et al. Print quality evaluation and applied colour management in coldset offset newspaper print
JP5886462B1 (ja) 印刷プロセスをダイナミックに較正する方法
CN103604754B (zh) 一种报纸印刷的文字墨色的测量方法
US20230213387A1 (en) Generating Spectra from Tristimulus Values
JP7353763B2 (ja) 色予測方法および色予測プログラム
CN113306289A (zh) 一种哑膜表面的印刷颜色校正方法
CN103618845B (zh) 一种基于最小显色误差激光打印机模型绿噪声半色调算法
JP3831074B2 (ja) 画像評価方法
JP6794802B2 (ja) 刷版カーブ作成装置、方法及びプログラム
JP6127585B2 (ja) 画像処理装置、画像処理システム、画像処理方法
Su et al. Optimal parameters based stochastic dot model for tone compensation of dither matrix
US11368606B1 (en) Image processing apparatus and non-transitory computer readable medium
JP7119845B2 (ja) 画像処理装置及びプログラム
US9213923B2 (en) Printer calibration using limited range reflection scanners as input sources

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20151021

Termination date: 20191128

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