CN102180003B - 一种打印机光谱特征化模型的修正方法 - Google Patents

一种打印机光谱特征化模型的修正方法 Download PDF

Info

Publication number
CN102180003B
CN102180003B CN 201110042430 CN201110042430A CN102180003B CN 102180003 B CN102180003 B CN 102180003B CN 201110042430 CN201110042430 CN 201110042430 CN 201110042430 A CN201110042430 A CN 201110042430A CN 102180003 B CN102180003 B CN 102180003B
Authority
CN
China
Prior art keywords
characterization model
printer
spectrum
spectrum characterization
correction
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
CN 201110042430
Other languages
English (en)
Other versions
CN102180003A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN 201110042430 priority Critical patent/CN102180003B/zh
Publication of CN102180003A publication Critical patent/CN102180003A/zh
Application granted granted Critical
Publication of CN102180003B publication Critical patent/CN102180003B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种打印机光谱特征化模型的修正方法。它可显著改善由不同墨水、纸张、温度以及空气湿度等打印环境变化所引起的光谱特征化模型精度降低的情况,适用于多种打印机光谱特征化模型。通过选取20-50个基本覆盖打印机色域的修正样本(墨水配方和光谱反射比),在新的打印环境下打印并测量,再经比较分析,求得相应的正向和反向波长修正系数,用于修正原光谱特征化模型的输入或输出参数,最终实现原正向和反向光谱特征化模型在新打印环境下的准确预测。本发明避免了重新建模所带来的浪费与不便,同时又克服了由打印环境变化造成的模型精度显著降低的问题,且易于实施。

Description

一种打印机光谱特征化模型的修正方法
技术领域
本发明涉及一种打印机光谱特征化模型的修正方法,具体是涉及修正由不同墨水、纸张、温度和空气湿度等打印环境变化引起的光谱特征化模型精度降低的方法。
背景技术
打印机光谱特征化模型通常包括正和反两个方向。正向光谱特征化模型是指建立从与设备相关的墨水配方(如,青、品红、黄、黑的数字化设备驱动值)到与设备无关的光谱反射比之间的关系;反之,则称为反向光谱特征化模型。由于人们常常希望通过打印可以获得某一种颜色,因此反向模型一般更为实用。打印机光谱特征化模型的建立往往需要打印并测量几百甚至几千个训练样本。然而,所建立的模型仅仅适用于建模时的打印环境,包括墨水、纸张、温度以及空气湿度等。一旦改变墨水(即使是同一型号的墨水)、纸张、打印时的温度和空气湿度,模型的精度会急剧地降低,需要再次建模。但是,建模需重新打印大量的训练样本,这又是对墨水、纸张等耗材的浪费,也不符合节能环保的理念。对于上述问题,目前学术界和工业界都还没有提出由于墨水、纸张、温度已经空气湿度等打印环境变化所导致光谱特征化模型精度降低的修正方法。
发明内容
本发明的目的是为了解决背景技术中所述的问题,提供一种打印机光谱特征化模型的修正方法。
打印机光谱特征化模型的修正方法,其特征在于包括正向打印机光谱特征化模型的修正和反向打印机光谱特征化模型的修正,其中:
正向打印机光谱特征化模型的修正的步骤如下:
1)选取20-50个基本覆盖打印机色域的不同于训练样本的修正样本,将其墨水配方作为输入墨水配方,经原正向光谱特征化模型预测相应的光谱反射比R pre
2)在新的打印环境下,通过打印机输出得到修正样本的硬拷贝,并测得其正向目标光谱反射比R target,正
3)采用最小二乘法,分别用截距为0的线性方程,拟合所有修正样本的R preR target,正在各波长的数据,并计算得出正向波长修正系数α λ; 
4)对于任一输入墨水配方,先利用原正向光谱特征化模型进行预测,然后用正向波长修正系数α λ修正原正向光谱特征化模型所预测的光谱反射比,得到基于新打印环境的预测光谱反射比;
反向打印机光谱特征化模型的修正的步骤如下:
5)选取20-50个基本覆盖打印机色域且不同于训练样本的修正样本,将其光谱反射比作为目标光谱反射比R target,反,通过原反向光谱特征化模型预测相应的墨水配方;
6)在新的打印环境下,通过打印机输出预测的墨水配方而得到修正样本的硬拷贝,并测得复现光谱反射比R rep
7)采用最小二乘法,分别用,其截距为0的线性方程拟合所有修正样本的R target,反R rep在各波长的对应关系,并计算得出反向波长修正系数β λ; 
8)对于任一输入目标光谱反射比,以反向波长修正系数β λ来修正目标光谱反射比R target,反,求得修正后的目标光谱反射比Rtarget,反,并将作为原反向特征化模型的输入目标光谱反射比,再经原反向光谱特征化模型预测在新打印环境下的墨水配方,实现原反向光谱特征化模型在新的打印环境下的准确复现。
本发明通过打印并测量20-50个修正样本以获得相应的波长修正系数,然后再对打印机光谱特征化模型的输入或输出参数进行修正,即可有效地补偿因打印环境变化所造成的打印机光谱特征化模型精度之降低,克服了重新建立特征化模型所造成的浪费与不便。
附图说明
图1是正向打印机光谱特征化模型的修正流程图;
图2是反向打印机光谱特征化模型的修正流程图;
图3是波长400nm的正向波长修正系数计算示意图。
具体实施方式
打印机光谱特征化模型的修正方法包括正向打印机光谱特征化模型的修正和反向打印机光谱特征化模型的修正、其中:
正向打印机光谱特征化模型的修正的步骤如下:
1)选取20-50个基本覆盖打印机色域的不同于训练样本的修正样本,将其墨水配方作为输入墨水配方,经原正向光谱特征化模型预测相应的光谱反射比R pre
2)在新的打印环境下,通过打印机输出得到修正样本的硬拷贝,并测得其正向目标光谱反射比R target,正
3)采用最小二乘法,分别用截距为0的线性方程,拟合所有修正样本的R preR target,正在各波长的数据,并计算得出正向波长修正系数α λ; 
4)对于任一输入墨水配方,先利用原正向光谱特征化模型进行预测,然后用正向波长修正系数α λ修正原正向光谱特征化模型所预测的光谱反射比,得到基于新打印环境的预测光谱反射比;
反向打印机光谱特征化模型的修正的步骤如下:
5)选取20-50个基本覆盖打印机色域且不同于训练样本的修正样本,将其光谱反射比作为目标光谱反射比R target,反,通过原反向光谱特征化模型预测相应的墨水配方;
6)在新的打印环境下,通过打印机输出预测的墨水配方而得到修正样本的硬拷贝,并测得复现光谱反射比R rep
7)采用最小二乘法,分别用,其截距为0的线性方程拟合所有修正样本的R target,反R rep在各波长的对应关系,并计算得出反向波长修正系数β λ; 
8)对于任一输入目标光谱反射比,以反向波长修正系数β λ来修正目标光谱反射比R target,反,求得修正后的目标光谱反射比Rtarget,反,并将作为原反向特征化模型的输入目标光谱反射比,再经原反向光谱特征化模型预测在新打印环境下的墨水配方,实现原反向光谱特征化模型在新的打印环境下的准确复现。
实施例
以一台CMYK(青、品、黄、黑)四墨打印机为例,对上述特征化模型修正方法的具体实施方法进行阐述。需要说明的是,本发明并不局限于四墨打印机,其可适用于各种多墨打印机基于光谱的特征化模型之修正。
基于某一打印环境的打印机光谱特征化模型的建立是进行模型修正的前提。打印机的正向光谱特征化模型可采用模型法和经验法两类。模型法是通过物理模型模拟设备的呈色机理,建立墨水配方与光谱反射比之间的对应关系,如Kubelka-Munk模型、Clapper-Yule模型、Neugebauer模型和Cellular Yule-Nielsen Neugebauer模型等,其反向模型往往是基于迭代计算的。经验法则通过数据拟合或插值等数学方法直接建立设备相关色与光谱反射比的对应关系,而不太注重设备的呈色机理,主要有多项式模型和三维查找表等,其反向模型也可通过多项式或查找插值的方法建立。在模型修正之前,首先假定基于某一打印环境的正向和反向打印机光谱特征化模型已经建立,其模型可采用模型法和经验法中的任一方法。下面分别对正向和反向光谱特征化模型的修正予以阐述。
在新的打印环境中,对于原正向光谱特征化模型的修正其主要步骤如下:
1)选取30个基本覆盖打印机色域的不同于训练样本的修正样本,其墨水配方将作为原正向光谱特征化模型的输入墨水配方以预测相应的光谱反射比R pre
2)通过打印机输出修正样本的硬拷贝,经分光光度计测得正向目标光谱反射比R target,正
3)分别用截距为0的线性方程拟合所有修正样本的R preR target,正在各波长的数据,并计算得出正向波长修正系数α λ,如以波长400nm为例,其对应的正向波长修正系数α 400R pre,400R target,400的关系可表示为
Figure DEST_PATH_IMAGE002
式中,R target,正,400为所有修正样本在波长400nm处的正向目标光谱反射比,R pre,400为所有修正样本在波长400nm处的预测光谱反射比,α 400为波长在400nm处的正向波长修正系数。利用最小二乘法求得最优的α400,如图3所示。
对于其他波长也做同样处理,以得到对应的正向波长修正系数α λ
4)对于任一墨水配方,利用原正向光谱特征化模型进行预测,然后将正向波长修正系数α λ用于修正原正向光谱特征化模型所预测的光谱反射比,即可求得基于新打印环境的预测光谱反射比。以波长400nm为例,可表示为
式中,R’ pre,400为修正后在波长400nm处的正向预测光谱反射比,即新打印环境下在400nm处的预测光谱反射比,R pre,400为在波长400nm处的原预测光谱反射比。
对于其他波长也做同样处理,可得到新打印环境下相应的预测光谱反射比。
在新的打印环境中,对于原反向光谱特征化模型的修正其主要步骤如下:
5)选取30个基本覆盖打印机色域的修正样本的光谱反射比,作为反向模型的目标光谱反射比R target,反,尽量避免选择颜色相近的样本。通过原反向光谱特征化模型预测相应的墨水配方。
6)根据预测的墨水配方通过打印机输出得到修正样本的硬拷贝,并经分光光度计测得其复现光谱反射比R rep
7)与正向波长修正系数α λ类似,分别用截距为0的线性方程拟合所有修正样本的R target,反R rep在各波长的数据,并计算得出反向波长修正系数β λ,如以波长400nm为例,其对应的反向波长修正系数β 400R rep,400R target,反,400的关系可表示为
Figure DEST_PATH_IMAGE006
式中,R target,反,400为所有修正样本在波长400nm处的反向目标光谱反射比,R pre,400为所有修正样本在波长400nm处的复现光谱反射比,β 400为波长400nm处的反向波长修正系数。然后,采用最小二乘法求得最优的β 400
对于其他波长也做同样处理,便可得到对应的反向波长修正系数。
8)对于任一输入目标光谱反射比,首先以反向波长修正系数β λ来修正目标光谱反射比R target,反,获得修正后的目标光谱反射比Rtarget,反,以波长400nm为例,可表示为
Figure DEST_PATH_IMAGE008
式中,Rtarget,反,400为修正后在波长400nm处的反向目标光谱反射比,R target,反,400为在波长400nm处的原反向目标光谱反射比。
对于其他波长也做同样处理,可得到对应的反向波长修正系数β λ
将修正后的目标光谱反射比Rtarget,反作为原反向光谱特征化模型的输入目标光谱反射比,经原反向光谱特征化模型即可预测在新的打印环境下的墨水配方,最终实现原反向光谱特征化模型在新的打印环境下的准确复现。

Claims (1)

1.一种打印机光谱特征化模型的修正方法,其特征在于包括正向打印机光谱特征化模型的修正和反向打印机光谱特征化模型的修正,其中:
正向打印机光谱特征化模型的修正的步骤如下:
1)选取20-50个基本覆盖打印机色域的不同于训练样本的修正样本,将其墨水配方作为输入墨水配方,经原正向光谱特征化模型预测相应的光谱反射比R pre
2)在新的打印环境下,通过打印机输出得到修正样本的硬拷贝,并测得其正向目标光谱反射比R target,正
3)采用最小二乘法,分别用截距为0的线性方程,拟合所有修正样本的R preR target,正在各波长的数据,并计算得出正向波长修正系数α λ; 
4)对于任一输入墨水配方,先利用原正向光谱特征化模型进行预测,然后用正向波长修正系数α λ修正原正向光谱特征化模型所预测的光谱反射比,得到基于新打印环境的预测光谱反射比;
反向打印机光谱特征化模型的修正的步骤如下:
5)选取20-50个基本覆盖打印机色域且不同于训练样本的修正样本,将其光谱反射比作为目标光谱反射比R target,反,通过原反向光谱特征化模型预测相应的墨水配方;
6)在新的打印环境下,通过打印机输出预测的墨水配方而得到修正样本的硬拷贝,并测得复现光谱反射比R rep
7)采用最小二乘法,分别用截距为0的线性方程,拟合所有修正样本的R repR target,反在各波长的的数据,并计算得出反向波长修正系数β λ; 
8)对于任一输入目标光谱反射比,以反向波长修正系数β λ来修正目标光谱反射比R target,反,求得修正后的目标光谱反射比Rtarget,反,并作为原反向光谱特征化模型的输入目标光谱反射比,再经原反向光谱特征化模型预测在新打印环境下的墨水配方,实现原反向光谱特征化模型在新的打印环境下的准确复现。
CN 201110042430 2011-02-22 2011-02-22 一种打印机光谱特征化模型的修正方法 Expired - Fee Related CN102180003B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110042430 CN102180003B (zh) 2011-02-22 2011-02-22 一种打印机光谱特征化模型的修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110042430 CN102180003B (zh) 2011-02-22 2011-02-22 一种打印机光谱特征化模型的修正方法

Publications (2)

Publication Number Publication Date
CN102180003A CN102180003A (zh) 2011-09-14
CN102180003B true CN102180003B (zh) 2013-04-24

Family

ID=44566375

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110042430 Expired - Fee Related CN102180003B (zh) 2011-02-22 2011-02-22 一种打印机光谱特征化模型的修正方法

Country Status (1)

Country Link
CN (1) CN102180003B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103538385B (zh) * 2012-11-14 2015-12-23 慈溪市宝日电器有限公司 基于胞元式Neugebauer光谱模型预测超墨量限制多墨样本光谱的方法
CN103530494A (zh) * 2012-12-10 2014-01-22 岑夏凤 一种反向cynsn模型中的子空间搜索算法
CN103544334B (zh) * 2012-12-10 2016-06-22 慈溪市桥头京姬电子元件厂 基于cii和色差的多墨打印机配方选择算法
CN103057292B (zh) * 2013-01-18 2014-11-05 慈溪市附海镇思维工业产品工作设计室 基于人眼特性的打印机定标算法
CN103540187B (zh) * 2013-03-29 2015-03-18 慈溪市校杰电器有限公司 打印机墨水配方的改进方法
CN106464775B (zh) * 2014-01-30 2019-06-14 惠普发展公司,有限责任合伙企业 颜色模型
CN103862858B (zh) * 2014-03-11 2015-12-30 武汉大学 一种多色打印印刷系统光谱分色方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2222341A (en) * 1988-06-16 1990-02-28 Hallmark Cards Producing print of colour monitor display
CN1417035A (zh) * 2001-11-06 2003-05-14 佳能株式会社 喷墨记录方式中的图像修正方法
CN1539646A (zh) * 2003-04-25 2004-10-27 柯尼卡美能达影像株式会社 图像形成装置的输出修正方法及图像形成装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3364249D1 (en) * 1982-03-30 1986-07-31 Crosfield Electronics Ltd Video retouching systems
JPS63254887A (ja) * 1987-04-10 1988-10-21 Konica Corp カラ−パツチの作成方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2222341A (en) * 1988-06-16 1990-02-28 Hallmark Cards Producing print of colour monitor display
CN1417035A (zh) * 2001-11-06 2003-05-14 佳能株式会社 喷墨记录方式中的图像修正方法
CN1539646A (zh) * 2003-04-25 2004-10-27 柯尼卡美能达影像株式会社 图像形成装置的输出修正方法及图像形成装置

Also Published As

Publication number Publication date
CN102180003A (zh) 2011-09-14

Similar Documents

Publication Publication Date Title
CN102180003B (zh) 一种打印机光谱特征化模型的修正方法
KR100278567B1 (ko) 3 차원 색소를 보다 높은 차원의 색소로 변환시키기 위한 시스템, 방법 및 프로그램
US8280161B2 (en) Spot color table compensation
USRE42673E1 (en) State-space based modeling of pixel elements of a dynamically varying color marking device
JP6888507B2 (ja) プロファイル調整方法、プロファイル調整プログラム、プロファイル調整装置、及び、プロファイル調整システム
CN103873736A (zh) 一种基于多维表的颜色闭环修正方法和系统
US20120200867A1 (en) Image processing device, image processing method and program for same, and printing system
KR100278430B1 (ko) 외부적으로 정의된 4 차원 색소를 등가의 4차원 색소로 변환시키기 위한 시스템, 방법 및 프로그램
US7782490B2 (en) Sensitivity matrix determination for adaptive color control
CN111866303B (zh) 颜色转换信息生成方法、以及颜色转换信息生成装置
US20090161183A1 (en) Optimal node placement for multi-dimensional profile luts for arbitrary media and halftones using parameterized minimization
CN102238297B (zh) 一种生成色彩对照文件的方法及系统
CN1406757A (zh) 自动彩色映像变换
JP5736382B2 (ja) 最適化されたプリンタ・キャリブレーション方法
JP2003125210A (ja) 色処理方法、記憶媒体、色処理装置、色変換装置、および画像形成装置
CN103978783A (zh) 一种用于静电数字印刷网点阶调的控制方法
Sun et al. Modified spectral Neugebauer model for printer characterization
US20070070367A1 (en) Linearization and characterization process of an image producing device
JP2011116129A (ja) 色処理装置およびその方法
US20130148140A1 (en) Systems and methods for inkjet printing
CN103057292B (zh) 基于人眼特性的打印机定标算法
US20180034980A1 (en) Calibration system that ensure improved accuracy while reducing calculation cost, calibration method, and recording medium therefor
CN111800557B (zh) 颜色转换信息生成方法、以及颜色转换信息生成装置
CN103540187B (zh) 打印机墨水配方的改进方法
JP2007259203A (ja) 測色方法および画像処理装置

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130424

Termination date: 20160222