CN102521650A - 一种基于粒子群优化的专色油墨配色方法 - Google Patents
一种基于粒子群优化的专色油墨配色方法 Download PDFInfo
- Publication number
- CN102521650A CN102521650A CN2011104053907A CN201110405390A CN102521650A CN 102521650 A CN102521650 A CN 102521650A CN 2011104053907 A CN2011104053907 A CN 2011104053907A CN 201110405390 A CN201110405390 A CN 201110405390A CN 102521650 A CN102521650 A CN 102521650A
- Authority
- CN
- China
- Prior art keywords
- value
- concentration
- formula
- expression
- particle
- 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.)
- Granted
Links
Images
Landscapes
- Spectrometry And Color Measurement (AREA)
Abstract
Description
技术领域
本发明属于印刷技术领域,具体涉及一种基于粒子群优化的专色油墨配色算法。
背景技术
随着人们生活品质的提高,越来越多的人开始追求印刷色彩对原稿的精确再现。由于专色不通过青、品红、黄、黑四色叠印的方法来合成颜色,而是采用特制的专色油墨来印刷该颜色,且每一种专色都有其本身固定的色相,从而能够保证印刷中颜色的准确性。因此,如何精确地实现专色油墨配色是印刷色彩精确再现的关键所在。
目前,主要的专色配色算法有三刺激值匹配和光谱匹配两种方法,三刺激值配色通过指定照明条件,以达到目标色和配方样品之间的三刺激值误差最小,但是由于存在同色异谱现象,因此在不同照明环境下两种原先颜色一致的印刷品所呈现出来的颜色也会出现差异,同时三刺激值配色算法的预测结果与实际配方相差很大,而配方修正算法的修正结果与初始化条件密切相关,而且经常发散,能收敛到确切值的结果也与实际配方有比较大的差异。与三刺激值匹配方法相比,直接光谱匹配方法寻求的是目标色和配方样品的光谱反射率曲线完全一致,原则上可以得到光谱异构性很低的配方,针对此配色算法,求解最终配方值实际上转换成了解约束优化问题,传统的优化方法容易陷入局部极值。因此,如何避免局部极值,收敛慢,使得目标色和配方样品的光谱反射率曲线完全一致是有待解决的技术问题。
发明内容
本发明的目的是针对现有技术存在的不足,提出了一种基于粒子群优化的专色油墨配色算法,利用该方法既能避免优化算法陷入局部极值,又能避免同色异谱对预测结果的干扰,解决了专色油墨实现精确配色的问题。本发明解决上述技术问题所采取的方法是:
但是在实验中发现这实际上是一条凹向上方的曲线,这主要是由于印刷时,纸张表面比较光滑,当油墨浓度较低时,纸张表面会产生明显的反射,而伴随着油墨浓度的增大,纸张表面反射所带来的影响会不断降低,因此值会随之上升。为了使值保持不变,通常采用多项式拟合,一般采用如下三阶表达式:
其中,, 表示底材的值矩阵,表示目标的值矩阵, , 表示油墨单位浓度的值,, 为修正系数。配色算法的精度取决于值的计算准确性,通过式(2)可以将值与浓度的关系曲线修正为直线。但是在实际计算过程中发现,由于已知,因此式(2)所表示的多项式曲线必须要保证过指定点,这给实际的拟合运算带来不便,为了解决这一问题,采用三次样条曲线拟合方法,经过改进,式( 1) 可写成:
步骤(3)粒子个体的位置修正。针对步骤(2)中生成的配方,根据粒子群优化算法公式进行调整,并生成新的浓度配方组合。具体是:
本发明与现有技术相比具有以下有益效果:通过基于粒子群优化的专色油墨配色方法可以使配色精度更加精确,配色效率大大提高,并且此方法所得的配色样品与标准样品之间的光谱分布非常接近,说明同色异谱异构性很低。
附图说明
图1是基色青(C)在不同梯度浓度下的光谱反射率曲线;
图2是基色品(M)在不同梯度浓度下的光谱反射率曲线;
图3是基色黄(Y)在不同梯度浓度下的光谱反射率曲线;
图4是No.1的标准(STD)与对应配方样品(Sample)的光谱反射率曲线;
图5是No.2的标准(STD)与对应配方样品(Sample)的光谱反射率曲线;
图6是No.3的标准(STD)与对应配方样品(Sample)的光谱反射率曲线;
图7是No.4的标准(STD)与对应配方样品(Sample)的光谱反射率曲线;
图8是No.5的标准(STD)与对应配方样品(Sample)的光谱反射率曲线。
具体实施方式
下面结合附图对本发明作进一步描述。
步骤(1)定标着色,通过多项式拟合方法求解修正系数,建立基础数据库。
先对每一种油墨分别以一定的浓度等级进行梯度着色,即单独印刷,一般采用6-7个浓度梯度,太少会影响式(2)的拟合精度,太多则会造成计算复杂度增加,之后确定底材和定标油墨的值。具体操作是选择Hangzhou Toka Ink Chemical Co., Ltd.生产的黄(Y)、品(M)、青(C)三种油墨作为基色,利用IGT C1胶印油墨打样机,按照不同的梯度浓度与冲淡剂混合,按每一浓度调拌均匀并在550N压力下进行打样,得到基色青(C)在不同梯度浓度下的光谱反射率曲线(参见图1)、基色品(M)在不同梯度浓度下的光谱反射率曲线(参见图2)、基色黄(Y)在不同梯度浓度下的光谱反射率曲线(参见图3),并以此作为基础数据库进行配色。
但是在实验中发现这实际上是一条凹向上方的曲线,这主要是由于印刷时,纸张表面比较光滑,当油墨浓度较低时,纸张表面会产生明显的反射,而伴随着油墨浓度的增大,纸张表面反射所带来的影响会不断降低,因此值会随之上升。为了使值保持不变,通常采用多项式拟合,一般采用如下三阶表达式:
其中,, 表示底材的值矩阵,表示目标的值矩阵, , 表示油墨单位浓度的值,, 为修正系数。配色算法的精度取决于值的计算准确性,通过式(2)可以将值与浓度的关系曲线修正为直线。但是在实际计算过程中发现,由于已知,因此式(2)所表示的多项式曲线必须要保证过指定点,这给实际的拟合运算带来不便,为了解决这一问题,采用三次样条曲线拟合方法,经过改进,式( 1) 可写成:
步骤(3)粒子个体的位置修正。针对步骤(2)中生成的配方,根据粒子群优化算法公式进行调整,并生成新的浓度配方组合。具体是:
具体操作是采用5个已知配方的印刷标样,将光谱反射率值输入到本文提出的配色算法中,执行步骤(2)到步骤(4)的程序,将计算得到的配方与实际配方做比较如下。
将上表中算法得到的配方在相同条件下进行打样,通过X-Rite 530分光光度计测量其色度值得到结果如下表所示。
Claims (1)
1.一种基于粒子群优化的专色油墨配色方法,其特征在于该方法包括以下步骤:
步骤(3)粒子个体的位置修正;针对步骤(2)中生成的配方,根据粒子群优化算法公式进行调整,并生成新的浓度配方组合;具体是:
步骤(4)检查结束条件;设定结束条件为或达到最大迭代次数,检查是否满足结束条件,若不满足则转回步骤(3)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110405390.7A CN102521650B (zh) | 2011-12-08 | 2011-12-08 | 一种基于粒子群优化的专色油墨配色方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110405390.7A CN102521650B (zh) | 2011-12-08 | 2011-12-08 | 一种基于粒子群优化的专色油墨配色方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102521650A true CN102521650A (zh) | 2012-06-27 |
CN102521650B CN102521650B (zh) | 2014-04-02 |
Family
ID=46292561
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110405390.7A Expired - Fee Related CN102521650B (zh) | 2011-12-08 | 2011-12-08 | 一种基于粒子群优化的专色油墨配色方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102521650B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102975490A (zh) * | 2012-11-27 | 2013-03-20 | 台州森林彩印包装有限公司 | 一种考虑成本因素的胶印专色油墨配色方法 |
CN103612483A (zh) * | 2013-12-06 | 2014-03-05 | 中国印刷科学技术研究所 | 一种基于光谱匹配的油墨配色方法 |
CN104730025A (zh) * | 2015-04-01 | 2015-06-24 | 河南工业大学 | 一种基于太赫兹光谱的混合物定量分析方法 |
CN108501556A (zh) * | 2018-02-01 | 2018-09-07 | 西安理工大学 | 一种薄膜凹版印刷专色配色方法 |
CN108875792A (zh) * | 2018-05-25 | 2018-11-23 | 魔金真彩网络科技(长沙)有限公司 | 一种基于机器学习的涂料修色方法及系统 |
CN112285095A (zh) * | 2020-09-15 | 2021-01-29 | 中国科学院上海技术物理研究所 | 基于弹性粒子群优化的火星物质分析仪在轨定标方法 |
CN113910796A (zh) * | 2021-11-15 | 2022-01-11 | 南京林业大学 | 一种基于k-m理论的印刷油墨配色方法 |
CN114316667A (zh) * | 2022-01-05 | 2022-04-12 | 上海紫江彩印包装有限公司 | 一种专色油墨及其调色方法和应用 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1946814A (zh) * | 2004-04-26 | 2007-04-11 | 昭和电工株式会社 | 涂料及其用途 |
CN1993698A (zh) * | 2004-08-03 | 2007-07-04 | 纳幕尔杜邦公司 | 预测化学混合物的性质的方法和仪器 |
CN101604151A (zh) * | 2009-06-29 | 2009-12-16 | 杭州电子科技大学 | 一种用于印刷色彩的检测控制方法 |
-
2011
- 2011-12-08 CN CN201110405390.7A patent/CN102521650B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1946814A (zh) * | 2004-04-26 | 2007-04-11 | 昭和电工株式会社 | 涂料及其用途 |
CN1993698A (zh) * | 2004-08-03 | 2007-07-04 | 纳幕尔杜邦公司 | 预测化学混合物的性质的方法和仪器 |
CN101604151A (zh) * | 2009-06-29 | 2009-12-16 | 杭州电子科技大学 | 一种用于印刷色彩的检测控制方法 |
Non-Patent Citations (1)
Title |
---|
林剑等: "基于改进遗传算法的胶印质量控制方法研究", 《包装工程》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102975490A (zh) * | 2012-11-27 | 2013-03-20 | 台州森林彩印包装有限公司 | 一种考虑成本因素的胶印专色油墨配色方法 |
CN103612483A (zh) * | 2013-12-06 | 2014-03-05 | 中国印刷科学技术研究所 | 一种基于光谱匹配的油墨配色方法 |
CN103612483B (zh) * | 2013-12-06 | 2016-03-09 | 中国印刷科学技术研究所 | 一种基于光谱匹配的油墨配色方法 |
CN104730025A (zh) * | 2015-04-01 | 2015-06-24 | 河南工业大学 | 一种基于太赫兹光谱的混合物定量分析方法 |
CN104730025B (zh) * | 2015-04-01 | 2017-04-26 | 河南工业大学 | 一种基于太赫兹光谱的混合物定量分析方法 |
CN108501556A (zh) * | 2018-02-01 | 2018-09-07 | 西安理工大学 | 一种薄膜凹版印刷专色配色方法 |
CN108875792A (zh) * | 2018-05-25 | 2018-11-23 | 魔金真彩网络科技(长沙)有限公司 | 一种基于机器学习的涂料修色方法及系统 |
CN112285095A (zh) * | 2020-09-15 | 2021-01-29 | 中国科学院上海技术物理研究所 | 基于弹性粒子群优化的火星物质分析仪在轨定标方法 |
CN113910796A (zh) * | 2021-11-15 | 2022-01-11 | 南京林业大学 | 一种基于k-m理论的印刷油墨配色方法 |
CN113910796B (zh) * | 2021-11-15 | 2022-11-11 | 南京林业大学 | 一种基于k-m理论的印刷油墨配色方法 |
CN114316667A (zh) * | 2022-01-05 | 2022-04-12 | 上海紫江彩印包装有限公司 | 一种专色油墨及其调色方法和应用 |
Also Published As
Publication number | Publication date |
---|---|
CN102521650B (zh) | 2014-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102521650A (zh) | 一种基于粒子群优化的专色油墨配色方法 | |
CN102729612A (zh) | 一种胶印专色油墨多目标优化配色方法 | |
CN102799895B (zh) | 基于最小二乘支持向量机的胶印油墨配色方法 | |
CN101100807A (zh) | 纺织品色牢度色差评级方法 | |
US20100189348A1 (en) | Spot color table compensation | |
US8953219B2 (en) | Methods for determining printer adjustments to accommodate proof/print substrate color differences | |
CN104902153A (zh) | 一种多光谱相机色彩校正方法 | |
CN102511050B (zh) | 用于优化显示器配置文件的方法 | |
CN107766681A (zh) | 一种印刷品黑色油墨预测光谱的校正技术 | |
CN104309300A (zh) | 一种彩印油墨计算机配色方法 | |
CN102238297B (zh) | 一种生成色彩对照文件的方法及系统 | |
CN105049841B (zh) | 利用单通道预滤光片提高彩色相机显色能力的方法 | |
CN105667069A (zh) | 一种基于比尔朗伯定律的光谱配色方法 | |
CN113487681B (zh) | 基于光谱灵敏度曲线和光源光谱优化的相机颜色标定方法 | |
CN105681624A (zh) | 一种基于光谱的专色叠印颜色预测方法 | |
CN105818518A (zh) | 一种用于热转印技术的凹印四色油墨配色方法 | |
LU501026B1 (en) | Dosing ink for digital printing on reflective substrates | |
CN106671594A (zh) | 一种确定印刷三原色密度的方法 | |
CN113791746B (zh) | 一种胶印专色余墨再配色方法、装置、系统及介质 | |
CN103264578B (zh) | 基于密度测量的印刷专色配色方法 | |
CN108389236A (zh) | 一种荧光喷墨油墨色光灰平衡控制方法 | |
TWI267448B (en) | Simulation of profile gives prediction of specific colors | |
KR20110003527A (ko) | 이펙트 안료를 포함하는 컬러 포뮬러의 제조 방법 | |
CN104869282A (zh) | 色谱方面的特性化数据的换算 | |
CN110880306A (zh) | 一种医用显示器色彩还原校正方法 |
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 |
Granted publication date: 20140402 Termination date: 20161208 |
|
CF01 | Termination of patent right due to non-payment of annual fee |