基于吸光度的PET薄膜专色配方预测方法
张研1,2, 周世生2,*, 曹从军2, 任鹏刚2
1. 西安理工大学机械与精密仪器工程学院, 陕西 西安 710048
2. 西安理工大学印刷包装工程重点实验室, 陕西 西安 710048
*通讯联系人 E-mail: zhoushisheng@xaut.edu.cn

作者简介: 张 研, 1987年生, 西安理工大学机械与精密仪器工程学院博士研究生 E-mail: 573300946@qq.com

摘要

基于吸光度与物质浓度的相关关系建立了一种PET薄膜印刷的混合油墨组分预测模型。 该方法首先依据光的传播理论和光在墨层与承印材料中的多重内反射原理, 利用基色和混合专色薄膜印品在黑白两色基底上的反射光谱, 建立具有高透低反特性薄膜印品透射光谱的获取方法, 并利用透射光谱求解吸收光谱; 然后对专色吸光度与基色浓度进行回归分析, 确定吸光度与基色浓度线性相关度高的特征波段, 再根据郎伯比尔定律, 利用特征波长处的专色样本吸光度与基色样本吸光度建立组分预测模型, 求解专色样本的配方比例; 最后, 分析预测配方与预设浓度的偏差, 并按照预测配方浓度重新打样, 采用色差与光谱均方根误差对计算配方色与目标专色的匹配程度进行分析, 验证模型精度。 以凹印PET薄膜双基色混合样品为实验对象, 对提出的方法进行测试。 光谱分析表明, 基色样本和专色样本光谱反射率曲线在黑白两种基底上有明显差异, 但都随基色油墨浓度变化具有相同的变化规律, 专色样本的光谱透射曲线随着基色浓度的变化接近浓度比例大的基色曲线, 在400~580 nm区间, 专色的吸收光谱随着基色油墨A浓度降低而升高, 在580~700 nm区间随着基色油墨A浓度降低而降低。 除570~590 nm吸收曲线相交区间以外, 在可见光范围内专色油墨吸光度与基色油墨浓度之间的判定系数 R2均大于0.95, 平均值为0.990 0, 两者具有强线性相关关系。 分别选取基色A和B波长为520 nm( R2为0.994 2)和700 nm( R2为0.998 5)处的吸光度代入配色模型, 并利用最小二乘法求解目标专色的配方比例。 实验结果表明, 六组预测浓度与预设浓度相比较, 配方平均偏差为2.5%, 无显著偏差。 配方色样本与专色样本色差最大值为1.98, 最小值为0.30, 平均值为0.85, 均满足GMI对专色复制的要求, 其中5组专色色差小于1.5, 符合专色忠实复制的要求; 6组专色的光谱均方根误差最大值为2.93%, 最小值为0.49%, 平均值为1.40%, 基本实现了同色同谱的高精度颜色复制。 验证了该方法对于提高专色油墨配色精度和改善PET薄膜印刷质量具有显著效果, 可为PET薄膜专色配方预测提供科学方法。

关键词: 吸光度; 光谱透射率; PET薄膜印刷; 计算机配色
中图分类号:O433.5+1 文献标志码:A
An Absorbance-Based Prediction Method of Spot Color Formula for PET Film
ZHANG Yan1,2, ZHOU Shi-sheng2,*, CAO Cong-jun2, REN Peng-gang2
1. School of Mechanical and Precision Instrument Engineering, Xi’an University of Technology, Xi’an 710048, China;
2. Shaanxi Provincial Key Laboratory of Printing and Packaging Engineering, Xi’an University of Technology, Xi’an 710048, China;
Abstract

A model for predicting spot color ink formula for PET film printing is proposed based on the correlation between absorbance and the concentrations of primary inks. Firstly, according to the methodology with the concept of light propagation and the multiple internal reflectance of ink layer and substrate, a method for obtaining the spectral transmittance of the thin film prints with high transmission characteristic is established by making use of the reflectance of the film prints with back and white substrates. The absorption spectrum of the film prints in the visible spectrum is obtainedby the spectral transmittance. Then, the regression equation between the absorbance of spot color inks and the concentrations of a primary color ink is established to determine characteristic wavelengths with strong linear correlation. After that, according to Lambert-Beer law, the prediction model is established by using the absorbance of the spot color samples and the absorbance of the primary color samples at the characteristic wavelengths, in order to obtain the proportions of primary color inks. Finally, the deviation between the predicted formulas and the actual concentrations is analyzed. In order to evaluate the accuracy of the prediction method, the remade samples are compared with target samples in terms of chromatic difference and the spectralroot mean square error. Some dual-componentspot color samples produced on PET filmby gravureare used as experimental subjects to verify the proposed method. The analysis of spectrum shows that the reflectance on the black substrate is significantly different from the reflectance on the black substrate, but both have the same trend with the variation of primary ink concentration. The transmittance curves of spot colorsamples are situated between the curves of primary color samples and movecloser to the primary color with higher concentration. The absorption spectrum of samples increase in the regionof 400~580 nm with the decrease of concentration of primary ink A, and decrease in the region of 580~700 nm as the concentration of A decreases. Except for the region of 570~590 nm, the linear correlation coefficients R2 between absorbance of spot color inks and the concentration of primary color A are higher than 0.95 and the average value is 0.990 0, which means a strong linear correlation in the visible spectrum range. The absorbancevalues at the wavelength 520 nm ( R2 of 0.994 2) and 700 nm ( R2 of 0.998 5) of the primaryinks A and B are selected to predict the formulas of spot colors by using the least squares method. The results show that the 6 groups of predicted concentrations are with 2.5% deviation from the actual concentrations of target samples, which means no significant difference. The maximum chromatic difference between the targetsamples and remade samples is estimated to be 1.98, the minimum to be 0.30, and the average value to be 0.85, which satisfies the requirement of spot color reproduction. 5 of the 6 groups are smaller than 1.5, which satisfies the requirement of faithful reproduction. The maximum RMSE is 2.95%, and the minimum is 0.49%, and the average value is 1.40%, which means a high precision color reproduction in the visible spectrum. It is confirmed that the proposed method could effectively improve the printing quality and the spot color matching precision, which may provide a scientific method for the predictive study of spot color ink formula of PET film printing.

Keyword: Absorbance; Spectral transmittance; PET film print; Computer color matching
引 言

在包装印刷领域广泛使用的专色, 不仅可以凸显产品的品质, 建立品牌专属性, 并且由于其颜色的特殊性, 是常用的防伪手段。 专色一般是通过基色油墨混合调配制得。 随着计算机技术的发展, 专色配色技术更加自动化、 智能化, 但如何保证配色的准确性和感官的一致性一直是色彩学领域的研究难点。 在纸张印刷配色研究中有基于密度配色方法、 基于三刺激值配色方法与基于神经网络等配色方法[1], 在纤维混色中有Stearns-Noechel配色模型和Friele配色模型[2, 3], 相关的研究均以反射光谱为基础, 应用最为广泛的是基于Kubelka-Munk理论的配色模型[4, 5], 它将反射光谱转换为物质成分的吸收系数与散射系数, 并根据线性加和性原理建模。 但已有的研究均以不透明材料为前提, 此类配色模型显然不适用于具有高透光学特性的薄膜印品。 PET聚酯薄膜[6, 7]作为常用的包装印刷薄膜基材, 是以聚对苯二甲酸乙二醇酯为原料的薄膜材料, 其配色仍以人工经验为主, 既浪费生产时间与材料, 又降低生产效率, 针对透明薄膜印品配色相关的研究鲜有报道。 另外, 常规测量方法获得的光谱反射率中包含部分基底表面反射的光能, 影响模型预测配方的精度。 据此, 基于薄膜印品的吸收光谱, 建立排除基底影响的计算机专色配方预测模型。 吸收光谱是物质分子、 原子或原子核对某一波长范围的光选择性吸收所呈现的光谱曲线, 并随着浓度的变化吸收强度也不同, 可以对物质进行定性、 定量及结构分析[8]

1 光谱配色模型
1.1 配色模型

物体对入射光线的选择性吸收和反射是物体表面产生颜色的本质原因, 对于半透明薄膜印品, 根据朗伯比尔定律[9]

A=lg1T=εcl(1)

其中, A为吸光度; T为光谱透射率; ε 为摩尔吸收系数; c吸光物质的摩尔浓度, mol· L-1; l为吸收层厚度, cm。

在光谱范围内各波长处的吸收程度是不相同的。 朗伯比尔定律描述了当一束平行单色光垂直通过某一均匀非散射的吸光物质时, 其吸光度与吸光物质的浓度c及吸收层厚度l成正比。 当吸收物质为多种物质混合而成时, 混合物质的吸光度为各组成物质吸光度之和, 可表示为

A=i=1NAi=i=1Nεi0lci(z)dz(2)

当待测样品满足朗伯比尔定律时, 吸光度与物质浓度呈现极强的线性关系。 根据这种线性关系, 利用不同波长下混合物质的吸光度, 可以求解其各个组分的浓度。

基于吸光度的油墨专色配色方法如式(3)所示。 该模型表示当某波长处混合油墨的吸光度与基色油墨浓度符合强线性相关时, 混合油墨的吸光度可用各基色油墨的浓度与其吸光度乘积之和表示。 求解过程为典型的最小值优化问题, 可用最小二乘法进行计算, 使式(3)最小的浓度c即为所求的专色配方。 因浓度具有明确的物理意义, 要求浓度必须满足非负性和有界性。

minf(c1, c2, , cn)=i=1nciAX, i(λ)-AM(λ)(3)s.t. c1+c2++cn=1c1, c2, , cn0

其中, ci为第i个基色油墨在混合油墨M中的浓度比例; n为基色油墨数量; AX, i(λ )为第i个基色油墨在波长λ 处的吸光度; AM(λ )为混合油墨在波长λ 处的吸光度。

为求解专色配方, 首先需获得基色油墨和目标专色油墨在可见光范围的吸收光谱, 由式(1)可知, 吸收光谱可通过油墨的光谱透射率换算获得, 再通过对基色油墨与专色样本进行浓度与吸光度的相关性分析以确定线性度高的特征波段。

1.2 光谱透射率计算

在上层介质为油墨下层介质为承印材料的印刷品中, 当光线入射到油墨表面, 部分光线在油墨表面发生反射, 部分光线穿透油墨层到达承印基底, 这部分光线继续在承印材料表面及内部传播, 产生反射和透射现象, 如图1所示。 此外, 光线还在油墨层内、 承印基材层内以及油墨与承印材料的界面中发生多重内反射[10]。 通常利用分光光度计、 光谱仪等设备可直接测得光谱反射率, 油墨光谱透射率则可通过计算获得。

图1 墨层和承印材料的反射和透射示意图Fig.1 Reflectance and transmittance of ink and substrate

定义Rx为油墨的光谱反射率, Tx为油墨的光谱透射率, Ru为承印材料的光谱反射率, Tu为承印材料的光谱透射率, 则油墨叠加承印材料的光谱反射率Rxu可以表示为

Rxu=Rx+RuTx2(1+RuRx+)(4)

数列求和得

Rxu=Rx+RuTx21-RxRu(5)

油墨叠加承印材料的光谱透射率Txu表示为

Txu=TxTu(1+RuRx+)(6)

求和

Txu=TxTu1-RxRu(7)

结合式(5)和式(7), 分离出变量RxTx

Rx=Tu2Rxu-RuTxu2Tu2-Ru2Txu2(8)Tx=(Rxu-Rx)1Ru-Rx(9)

为求解RxTx, 将印制好油墨的PET薄膜(理论透射率100%)样条分别放置在黑白两种不透明漫反射基材上利用分光光度计进行光谱反射率测量, 得到两组油墨叠加基材的光谱反射率数据Rxw(白底)和Rxb(黑底)。 将RxwRxb以及两种基材的光谱反射率数据Ruw(白底)和Rub(黑底)分别代入式(5)可以推导出油墨的光谱反射率

Rx=RuwRxb-RubRxwRuw(1+RubRxb)-Rub(1+RuwRxw)(10)

RxwRuw代入式(9)即求出油墨的光谱透射率

Tx=(Rxw-Rx)1Ruw-Rx(11)

1.3 模型性能评价方法及指标

评价配色模型精度最直观的方法是按照模型预测配方重新进行打样, 将按照预测配方重新打样的样本定义为配方色样本, 按照预设配方印制的样本定义为目标专色样本, 对配方色与目标专色采用色差和光谱均方根误差两个指标进行差异分析。 色差评价是利用分光光度计测量配方色与目标专色的色度值, 采用CIEDE2000色差公式[11, 12, 13]进行计算。 计算公式如式(12)

ΔE00={ΔL'KLSL2+ΔC'KCSC2+ΔH'KHSH2+RTΔC'KCSC2ΔH'KHSH2}12(12)

在包装印刷领域, 根据GMI(国际图形测量公司)认证标准, 专色色差小于1.5则可视为无色差的忠实复制, 色差小于2.8表示有小量色差, 但满足专色复制要求。

2 实验与结果分析
2.1 样本制备

利用RK printing proofer凹印打样机按照预先设置的配方浓度制备若干样本。 具体过程是利用METTLER TOLEDO ME204电子天平称取300g凹印蓝墨和300g凹印红墨分别与稀释剂按照3: 1的比例混合作为调配专色的基色油墨A和B。 实验中所用稀释剂为乙酸乙酯和乙醇按1: 1比例混合得到的溶液。 基色油墨配制完成后, 按照表1所示的预设浓度准确称量基色油墨A和B并充分调配, 然后利用凹印打样机在0.15 mm PET薄膜上印刷基色与专色, 获得建模所需的基色油墨样本(S01, S02)和专色油墨样本(S03, S04, S05, S06, S07, S08)。 打样完成后静置样品至油墨充分干燥。

表1 实验制备样本的预设配方 Table 1 Recipes of printing samples
2.2 光谱数据预处理

采用爱色丽X-Rite CI64分光光度计对所有样本进行光谱反射率测量。 测量时将样本平铺在基材上方, 每个样本分别以白色纸张和黑色纸张两种基材为背景, 黑白纸张为均匀漫反射不透明材质, 同时黑白两种基材的光谱反射率也将被记录下来, 测量过程要求在平整的台面上进行。 每个样本测量3次, 取其平均值作为各样本最终的测量值。 分光光度计的采样波长为400~700 nm, 采样间隔为10 nm。

实验共有基色样本2组, 专色样本6组, 基色样本与专色样本的光谱反射率测量结果如图2所示。 图2(a)是白色基底uw与白色基底上各样本的光谱反射率, 图2(b)是黑色基底ub与黑色基底上各样本的光谱反射率。 由于背景色不同, 样本的光谱反射率有明显差异, 但都随基色油墨浓度变化具有相同的变化规律。

图2 黑白基底与样本的光谱反射率
(a): 白色基底; (b): 黑色基底
Fig.2 Reflectance of samples on different backgrounds
(a): White background; (b): Black background

测量完成后, 将测得的光谱反射率代入式(10)与式(11)计算所有样本的光谱透射率并绘制曲线, 结果如图3(a)所示。 计算所得的光谱透射率排除了背景底色的影响, 仅表征PET薄膜上墨层的透射特性。 图3(a)中的透射率曲线显示出专色样本的透射性能介于两组基色样本之间, 且随着基色A浓度降低, B浓度升高, 专色样本的光谱透射曲线接近浓度比例大的基色。

图3 样本光谱透射率曲线(a)和吸收曲线(b)Fig.3 Transmittance curves (a) and Absorbance curves (b) of samples

将光谱透射率代入式(1)计算样本的吸收光谱并绘制曲线得到如图3(b)所示结果。 该图表明: 专色样本对光谱的吸收强度随组分浓度比例改变而变化, 在400~580 nm区间, S03— S08的吸收光谱随着基色A浓度降低而升高; 在580~700 nm区间, S03— S08的吸收光谱随着基色A浓度降低而降低, 所有样本吸收曲线相交于570~590 nm波长区域, 另外所有样本的光谱反射率曲线(白底)、 光谱透射率曲线也均相交于此处。

2.3 配方预测

在求解配方比例之前, 首先将样本的吸光度与基色A浓度进行回归分析, 验证吸光度与基色浓度的线性关系, 实验结果如图4所示。 在可见光400~700 nm范围内除波长570~590 nm区域外, 其他各波长处的吸光度与基色A浓度的判定系数R2均大于0.95, 平均值为0.990 0, 线性程度高。 570~590 nm区间为各浓度样本光谱曲线相交的位置, 吸光度与浓度不具有典型的线性关系。

图4 可见波长范围内吸光度与基色A浓度的判定系数Fig.4 Coefficient of determination between absorbance and concentration of primary ink A in visible spectral regin

由于专色样本制作时仅使用了两种基色, 从图4中选用两组线性程度高的特征波段用于配方预测, 特征波段的选择应避开570~590 nm波长区间, 这样既可保证预测精度, 又可使运算过程简洁。 特征波段选取的计算结果如图5所示。 520 nm波长处相关系数R2为0.994 2, 700 nm波长处的相关系数R2为0.998 5, 均呈强线性相关关系, 符合配色模型计算要求。

图5 特征波长吸光度与基色A浓度的回归分析Fig.5 Regression analysis between absorbance and concentration of primary ink A at characteristic wavelengths

将基色样本S01, S02在波长为520与700 nm处的吸光度代入配色模型, 利用最小二乘法分别求解目标样本S03, S04, S05, S06, S07和S08的配方浓度比例, 计算结果如表2所示。 从预测浓度与预设浓度的偏差角度分析, 两者非常接近, 平均偏差为2.5%。

表2 目标样本配方预测 Table 2 Predicted recipes of target samples
2.4 配方精度验证

表2中结果可以观察到预测配方与预设浓度十分接近, 为进一步评价配方预测精度, 将基色A与B按照表2所示的预测浓度比例混合并打样, 配方样本命名为C03— C08, 分别对应专色样本S03— S08, C03— C08制备过程与S03— S08一致。 待配方样本完全干燥后, 将专色样本与配方样本平铺放置在标准光源箱中D50光源下进行主观评价实验, 选取15名视觉正常人员参与评价颜色色相、 明度和饱和度差异, 得出的结论是C07与S07亮度与色相有微小差异, 其余样本均无显著差别。

除主观评价外, 将图6中所有样本平放在白色基底上, 利用分光光度计测量光谱反射率与色度值, 并用式(12)计算相应样本的色差。 测量结果如表3所示。 从表中可以看出, 6组样本中, 目标样本与配方样本之间色差最大值为1.98, 最小值为0.30, 平均值为0.85, 全部满足专色复制要求, 其中S03, S04, S05, S06和S08色差均小于1.5, 满足无色差忠实复制要求。

图6 目标样本与配方样本主观评价Fig.6 Visual evaluation of target samples and remade samples

表3 目标样本与配方样本的的精度评价(测量条件: D50光源/2° 视场) Table 3 Chromatic difference and RMSE of samples (D50/2° )

图7描述的是C03— C08配方样本与相应目标样本S03— S08的光谱反射率曲线的匹配程度, 从图7可以看出, 配方样本C03, C04, C05, C08与目标样本的光谱反射率曲线几乎重合, C06与C07两个样本在红色波段与目标样本存在少许差异。 结合表3可知, 目标样本与配方样本之间均方根误差RMSE最大值为2.93%, 最小值为0.49%, 平均值为1.40%, 基本实现了同色同谱的高精度颜色复制。

图7 目标样本与配方样本光谱反射率曲线Fig.7 Spectral reflectance curves of target samples and remade samples

3 结 论

针对PET薄膜专色配色, 利用薄膜高透低反性质和黑白两种基底建立了获得薄膜印品光谱反射率和光谱透射率的方法。 根据郎伯比尔定律以及吸光度与浓度在特征波长处的强线性相关关系, 提出了基于吸光度的薄膜印品配色方法, 并以凹版印刷专色配方预测为例, 从颜色色差与光谱均方根误差两方面验证了配色的有效性和准确性。 该方法无需制作基色油墨梯尺, 建模样本少, 计算过程简洁, 一次预测就可以获得与真实比例极为接近的预测结果, 解决了现有薄膜印品配色精度无法保证的问题, 也为其他印刷方式或承印材料印刷提供了新的配色思路。 本方法仅用二元混合的专色进行了验证, 后续研究将对三元及三元以上混和的专色进行验证。

The authors have declared that no competing interests exist.

参考文献
[1] Furferi R, Governi L. Color Research & Application, 2011, 36(3): 179. [本文引用:1]
[2] Shen J, Zhou X, Ma H, et al. Textile Research Journal, 2016, 87(3): 296. [本文引用:1]
[3] Li Rong, Gu Feng. Coloration Technology, 2010, 122(2): 74. [本文引用:1]
[4] Schabbach L M, Bondioli F, Fredel M C. Journal of the European Ceramic Society. 2011, 33(5): 659. [本文引用:1]
[5] Neuman M, Coppel L, Edstroem P. Nordic Pulp & Paper Research Journal, 2012, 27(2): 426. [本文引用:1]
[6] Stinson R M, Obendorf S K. Journal of Applied Polymer Science, 2015, 62(12): 2121. [本文引用:1]
[7] GUO Mei, XIONG Ying, XU Peng, et al(郭梅, 熊英, 徐鹏, ). Polymer Materials Science and Engineering(高分子材料科学与工程), 2010, 26(2): 54. [本文引用:1]
[8] LIU Hui-jun, TAO Shao-hua, YANG Bing-chu, et al(刘慧军, 陶少华, 杨兵初, ). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2015, 35(12): 3490. [本文引用:1]
[9] Mayerhoefer Thomas G, Muschke Harald, Popp Juergen. ChemPhysChem, 2016, 17: 1948. [本文引用:1]
[10] Klein G A. Industrial Color Physics. Springer-Verlag, Berlin, 2004. [本文引用:1]
[11] Luo M R, Cui G, Rigg B. Color Research & Application, 2001, 26(5): 340. [本文引用:1]
[12] Sharma G, Wu W, Dalal E N. Color Research & Application, 2005, 30: 21. [本文引用:1]
[13] Pecho O E, Ghinea R, Alessand retti R, et al. Dental Materials, 2016, 32(1): 82. [本文引用:1]