CN106671594B - 一种确定印刷三原色密度的方法 - Google Patents
一种确定印刷三原色密度的方法 Download PDFInfo
- Publication number
- CN106671594B CN106671594B CN201611248758.2A CN201611248758A CN106671594B CN 106671594 B CN106671594 B CN 106671594B CN 201611248758 A CN201611248758 A CN 201611248758A CN 106671594 B CN106671594 B CN 106671594B
- Authority
- CN
- China
- Prior art keywords
- value
- green
- yellow
- density
- red
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 66
- 238000007639 printing Methods 0.000 title claims abstract description 60
- 239000003086 colorant Substances 0.000 title claims abstract description 31
- 230000004075 alteration Effects 0.000 claims abstract description 40
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000002310 reflectometry Methods 0.000 claims abstract description 10
- 230000001419 dependent effect Effects 0.000 claims description 9
- 238000012886 linear function Methods 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 6
- LEYJJTBJCFGAQN-UHFFFAOYSA-N chembl1985378 Chemical compound OC1=CC=C2C=CC=CC2=C1N=NC(C=C1)=CC=C1N=NC1=CC=C(S(O)(=O)=O)C=C1 LEYJJTBJCFGAQN-UHFFFAOYSA-N 0.000 claims description 3
- 238000005286 illumination Methods 0.000 claims description 3
- 238000012887 quadratic function Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims 1
- 230000003595 spectral effect Effects 0.000 abstract description 15
- 238000005516 engineering process Methods 0.000 abstract description 10
- 238000003908 quality control method Methods 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 description 12
- 238000000205 computational method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 241001062009 Indigofera Species 0.000 description 2
- 206010054949 Metaplasia Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000015689 metaplastic ossification Effects 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B41—PRINTING; LINING MACHINES; TYPEWRITERS; STAMPS
- B41F—PRINTING MACHINES OR PRESSES
- B41F33/00—Indicating, counting, warning, control or safety devices
- B41F33/0036—Devices for scanning or checking the printed matter for quality control
Abstract
本发明属于印刷过程质量控制技术领域,尤其涉及一种确定印刷三原色密度的方法。本发明的方法是设计一条每一墨区包含青、品红、黄、红、绿、蓝等六种色块的测控条,印刷后检测每一色块的光谱反射率,计算出每一色块的三刺激值和其中青、品红、黄色块的主密度值,利用最小二乘法分别拟合青、品红、黄、红、绿、蓝的三刺激值和青、品红、黄主密度值的关系;再分别针对青、品红、黄三色,根据三刺激值和主密度之间的函数关系,循环计算主密度从低到高变化时,青、品红、黄、红、绿、蓝的三刺激值和色度值,再计算与标准色度值之间的色差。本发明所述的方法只需以扫描方式检测一张或以上印张上各墨区的光谱反射率,对样本数量要求低,通过程序计算确定最合理的密度,效率高,对于使用者的技术要求低。
Description
技术领域
本发明属于印刷过程质量控制技术领域,尤其涉及一种确定印刷三原色密度的方法。
技术背景
针对一种给定的纸张、油墨组合,多大的印刷墨量才是最合理的问题,传统的方法是通过调整印刷时输出墨量大小,获得墨色深浅不同的印张,检测不同印张上的相对反差,以相对反差最大的印张的密度为最佳实地密度来控制青、品红、黄、黑四色的印刷墨量。文献《一种胶印油墨最佳印刷参数的确定方法研究》(2010年中国印刷与包装学术会议会刊,蒋文燕,王志宏)通过建立相对反差和实地密度的关系曲线,求出最佳实地密度和相对反差。与直接选择相对反差最大的样张的密度作为最佳密度相比,通过建立相对反差和实地密度的关系曲线,再确定的最佳实地密度更加准确。依据相对反差来确定最佳密度的方法可以确定当前生产条件下的最适合的墨量要求,但需要印刷并检测较多的测试样张,成本高、效率低,并且由于不同印刷企业各自依据本厂生产条件确定最佳密度,同一印前文件在不同印刷厂的印刷效果一致性差。
由于印刷生产专业分工协助程度的提高,为了使印前制作的文件由不同的印刷厂印刷后,能达到相同的、并与原稿尽可能一致的效果,不再根据相对反差来控制墨量,而是要求按照标准化的墨色要求来控制墨量,如《ISO 12647印刷过程控制标准:2004》是平版印刷生产过程的控制标准,该标准把纸张分为五类,针对每一类纸张,给出青、品红、黄、黑、红、绿、蓝色块的标准色度值,要求通过调整印刷墨量,使印张上的青、品红、黄、黑、红、绿、蓝色块的色度值尽可能接近标准色度值。一种方法通过调整印刷输出墨量,印刷大量墨色深浅不同的印张,或者一张印张上不同位置墨色深浅不同,检测深浅不同的青、品红、黄、黑色块的色度值,针对每一色,找出与标准色度值之间色差最小的色块,该色块的密度作为该色的最佳密度用于控制印刷墨量。这种方法需要印刷大量的测试样张,需要大量的人工检测与计算,成本高,效率低。有人提出采用比较不同印张上青、品红、黄、红、绿、蓝色的色差之和或平均色差,取色差和或平均色差最小的样张来确定青、品红、黄的密度,但由于红、绿、蓝色的色度值同时受二种原色的密度影响,需要检测计算大量的密度大小分布不同的样本才能获得较为合理的结果,效率太低。
中国专利公开号为CN105716995A,申请日为2016年3月18日的专利申请文件公开了一种确定最佳青、品红、黄、黑四色印刷密度的方法,针对青、品红、黄、黑四色中每一色,通过建立该色三刺激值和主密度之间的函数关系,再确定与标准色色差最小时的密度为最佳密度值。该方法通过建立的函数关系,自动地计算出最佳密度值,能大幅度节省材料和时间浪费,提高了生产效率。但是该方法在确定青、品红、黄的最佳密度时,是依据青、品红、黄三种原色的色度值各自尽量接近标准值来确定密度,没有顾及红、绿、蓝三种二次色也需要尽量接近标准值,而二次色的色度值会随着形成该二次色的二种基本色的密度而变化。当原色色度值尽量接近标准时,二次色色度值可能偏离标准较大,或者原色和二次色的综合颜色效果并不合理。
发明内容
本发明所要解决的技术问题是高效地确定平版印刷时青、品红、黄三色最合理密度值的方法,所述的最合理青、品红、黄密度能够使青、品红、黄、红、绿、蓝六色的色度值尽可能地接近规定的标准色度值。为此,本发明采用以下技术方案:
本发明的方法是设计一条每一墨区包含青、品红、黄、红、绿、蓝等六种色块的测控条,印刷后检测每一色块的光谱反射率,计算出每一色块的三刺激值和其中青、品红、黄色块的主密度值,利用最小二乘法分别拟合青、品红、黄、红、绿、蓝的三刺激值和青、品红、黄主密度值的关系;再分别针对青、品红、黄三色,根据三刺激值和主密度之间的函数关系,循环计算主密度从低到高变化时,青、品红、黄、红、绿、蓝的三刺激值和色度值,再计算与标准色度值之间的色差,能够使得六色总色差值最小的青、品红、黄的主密度值确定为最合理密度值。
在采用上述技术方案的基础上,本发明还可采用以下进一步的技术方案:
具体包括以下步骤:
a.以印刷机的一个墨区宽度为单位,当印刷机墨区数为奇数时,设计一条长度小于纸张宽度,并且为一个墨区宽度奇数倍的数字测控条;当印刷机墨区数为偶数时,设计一条长度小于纸张宽度,并且为一个墨区宽度偶数倍的数字测控条。测控条的每一墨区宽度内从左到右依次排列青、蓝、品红、绿、红、黄六种色块。印前制作时,在大版文件的拖梢边边缘,沿印刷滚筒轴向添加所述数字测控条,数字测控条的轴向中心与大版文件的轴向中心对齐。
b.印刷附着所述测控条的印张,利用扫描式分光光度计测量印张上的测控条,记录测控条的每一墨区的青、蓝、品红、红、绿、黄色色块的光谱反射率数据,利用青、品红、黄色块的光谱反射率数据计算青、品红、黄三色的XYZ值和主密度值,利用红、绿、蓝色块的光谱反射率数据计算每一色块的XYZ值。其中XYZ值指国际照明委员会创立的CIE 1931 XYZ颜色空间的坐标值,其中青色块的主密度指青色块的ISO标准红光密度,品红色块的主密度指品红色块的ISO标准绿光密度,黄色块的主密度指黄色块的ISO标准蓝光密度。计算XYZ值采用成熟的技术,如GB/T 19437-2004;计算青、品红、黄色的主密度值采用成熟的技术,如GB/T 23649-2009。
c.根据每一青色块的主密度值DC和X值XC,利用最小二乘法进行线性拟合,得到以DC为自变量,XC为因变量的线性函数采用同样的方法,得到以青色块的Y值和Z值随着主密度而变化的函数和
采用同样的方法,分别得到品红色的XYZ值随品红色主密度而变化的线性关系函数黄色的XYZ值随黄色主密度而变化的线性关系函数
以每一墨区内的青色块主密度DC、同一墨区内的品红色块主密度DM为自变量,同一墨区内的蓝色块的X值XB为因变量,根据二元二次函数利用最小二乘法进行拟合:
其中a0、a1、a2、a3、a4为经拟合得到的系数;采用同样的方法,得到蓝色块的Y值和青色块主密度、品红色块主密度之间的关系函数蓝色块的Z值和青色块主密度、品红色块主密度之间的关系函数
采用同样的方法,得到绿色块的XYZ值和青色主密度、黄色主密度之间的关系函数
采用同样的方法,得到红色块的XYZ值和品红色主密度、黄色主密度之间的关系函数
d.预设初始的最小色差ΔEmin为600;
根据青色的XYZ值随着主密度值而变化的函数 循环计算当青色主密度值从0.70开始,以0.01为间距,直至1.80时对应的XYZ值,再根据XYZ值计算L*a*b*值,再计算青色当前L*a*b*值与标准L*a*b*值之间的色差ΔEC,其中L*a*b*值指国际照明委员会推荐的CIE 1976L*a*b*空间的颜色坐标值。计算L*a*b*值和色差采用成熟的技术,如GB/T 19437-2004。
针对青色主密度值每一次取值,根据品红色的XYZ值随着品红色主密度值而变化的函数蓝色的XYZ值随青色主密度和品红色主密度而变化的函数 循环计算当品红色主密度值从0.70开始,以0.01为间距,直至1.80时对应的品红色XYZ值和蓝色XYZ值,再根据品红色的XYZ值计算品红色的L*a*b*值,根据蓝色的XYZ值计算蓝色的L*a*b*值,再计算品红色当前L*a*b*值与品红色标准L*a*b*值之间的色差ΔEM,计算蓝色当前L*a*b*值与蓝色标准L*a*b*值之间的色差ΔEB。
针对品红色主密度值每一次取值,根据黄色的XYZ值随着黄色主密度值而变化的函数绿色的XYZ值随青色主密度和黄色主密度而变化的函数红色的XYZ值随品红色主密度和黄色主密度而变化的函数 循环计算当黄色主密度值从0.60开始,以0.01为间距,直至1.60时对应的黄色XYZ值、绿色XYZ值和红色XYZ值,再根据黄色XYZ值计算黄色的L*a*b*值,根据绿色XYZ值计算绿色的L*a*b*值,根据红色XYZ值计算红色的L*a*b*值,再计算黄色当前L*a*b*值与黄色标准L*a*b*值之间的色差ΔEY,计算绿色当前L*a*b*值与绿色标准L*a*b*值之间的色差ΔEG,计算红色当前L*a*b*值与红色标准L*a*b*值之间的色差ΔER。计算青、品红、黄、红、绿、蓝六色的总色差ΔE6,计算公式表示为:
ΔE6=ΔEC+ΔEM+ΔEB+ΔEY+ΔEG+ΔER;
如果当前六色总色差值ΔE6小于最小色差值ΔEmin,令最小色差值ΔEmin等于当前色差值ΔE6,表示为:
ΔEmin=ΔE6;
取当前的青、品红、黄密度为最合理密度DC_best、DM_best、DY_best,表示为:
DC_best=DC;
DM_best=DM;
DY_best=DY。
进一步地,步骤c中,利用最小二乘法根据二元二次多项式拟合蓝色块的值和同一墨区内的青、品红主密度值之间的关系,绿色块的值和同一墨区内的青、黄主密度值之间的关系,红色块的值和同一墨区内的品红、黄主密度值之间的关系;步骤d中,随青、品红、黄三色主密度的变化,以能够使青、品红、黄、红、绿、蓝六色总色差最小的密度确定为青、品红、黄三色最合理密度。
本发明的有益效果:本发明所述的方法所确定的青、品红、黄三色密度能够使得青、品红、黄、红、绿、蓝六色色度值与标准色度值之间的色差之和最小,精度高,有利于满足印刷生产标准化的要求;本发明所述的方法只需以扫描方式检测一张或以上印张上各墨区的光谱反射率,对样本数量要求低,通过程序计算确定最合理的密度,效率高,对于使用者的技术要求低。
附图说明:
图1为本发明的青、品红、黄色最合理密度确定方法的实现流程图。
图2为本发明设计的数字测控条的其中一个墨区示意图。
具体实施方式
结合附图,对本发明提供的一种确定印刷三原色密度的方法作进一步说明。
印刷生产过程中,为了通过分工协助提高生产效率,同时又能正确地还原原稿的颜色,需要整个生产过程按照标准化要求进行质量控制。标准化的印刷生产要求青、品红、黄等原色和红、绿、蓝等二次色的色度值尽可能接近标准色度值要求。对于给定的纸张和油墨,青、品红、黄等原色的色度值随着青、品红、黄的密度而变化;而红、绿、蓝等二次色由于是二种原色叠印而成,其色度值会随着形成该二次色的二种颜色的密度而变化,由于二次色色度值随二种原色密度而变化的关系相对复杂,目前生产中,一般只根据使青、品红、黄等原色度值尽可能接近标准要求的原则来确定青、品红、黄的最合理密度,用于控制印刷过程中的输出墨量。但在这种青、品红、黄密度条件下,红、绿、蓝可能并不符合标准要求,或者综合比较青、品红、黄、红、绿、蓝六种颜色与相应标准的总色差并不是最小的。本发明所要解决的技术问题是提供一种确定平版印刷时青、品红、黄三色最合理密度值的方法,所述的最合理青、品红、黄密度能够使青、品红、黄、红、绿、蓝六色与相应标准之间色差之和达到最小,更好符合标准化生产的要求。
本发明的方法是设计一条每一墨区包含青、品红、黄、红、绿、蓝色块的测控条,如图2所示。印刷后检测每一色块的光谱反射率,计算出每一色块的三刺激值和每一青、品红、黄色块的主密度值,利用最小二乘法分别拟合青、品红、黄、红、绿、蓝的三刺激值和青、品红、黄主密度值的关系;再分别针对青、品红、黄三色,根据三刺激值和主密度之间的函数关系,循环计算主密度从低到高的变化过程时,青、品红、黄、红、绿、蓝的三刺激值和色度值,再计算与标准色度值之间的色差,能够使得六色总色差值最小时的青、品红、黄的主密度值确定为最合理密度值,如图1所示。具体包括以下步骤:
步骤a以印刷机的一个墨区宽度为单位,当印刷机墨区数为奇数时,设计一条长度略小于纸张宽度,并且为一个墨区宽度奇数倍的数字测控条;当印刷机墨区数为偶数时,设计一条长度略小于纸张宽度,并且为一个墨区宽度偶数倍的数字测控条。测控条的每一墨区宽度内从左到右依次排列青、蓝、品红、绿、红、黄色色块。印前制作时,在大版文件的拖梢边边缘,沿印刷滚筒轴向添加所述数字测控条,数字测控条的轴向中心与大版文件的轴向中心对齐。
步骤b以印刷机的一个墨区宽度为单位,当印刷机墨区数为奇数时,设计一条长度小于印刷纸张宽度2厘米或以上,并且为一个墨区宽度奇数倍的数字测控条;当印刷机墨区数为偶数时,设计一条长度小于印刷纸张宽度2厘米或以上,并且为一个墨区宽度偶数倍的数字测控条。测控条的每一墨区宽度内从左到右依次排列青、蓝、品红、绿、红、黄色色块。
例如,某一型号的印刷机有34个墨区,每一墨区的宽度为30毫米,印刷纸张的宽度为1000毫米,可以设计一条包含32个墨区宽度,即宽度为960毫米的测控条。测控条的每一墨区内依次排列青、蓝、品红、绿、红、黄色色块。印前制作时,在大版文件的拖梢边边缘,沿印刷滚筒轴向添加上述制作的数字测控条,数字测控条的轴向中心与大版文件的轴向中心对齐,使得印刷时,控制条轴向中心与印刷机轴向中心对齐,控制条上的每一墨区和印刷机的墨区位置相对应。
步骤b印刷附着所述测控条的印张,利用扫描式分光光度计测量印张上的测控条,记录测控条的每一墨区的青、蓝、品红、红、绿、黄色色块的光谱反射率数据,如步骤a所述测控条,扫描一张样张上的测控条,可快速获得青、蓝、品红、红、绿、黄色各32个色块的光谱反射率数值;如果扫描检测n张样张,则获得青、蓝、品红、红、绿、黄色各32×n个色块的光谱反射率数值。利用每一青、品红、黄色色块的光谱反射率数据计算该色块的XYZ值和主密度值,利用每一红、绿、蓝色块的光谱反射率数据计算该色块的XYZ值。
利用物体色的光谱反射率计算XYZ值有成熟的技术,计算方法如式(1)-(4)所示。利用物体色的光谱反射率计算主密度值有成熟的技术,计算方法如式(5)所示。
X=K·∑(Sλ·ρλ·x2) (2)
Y=K·∑(Sλ·ρλ·y2) (3)
Z=K·∑(Sλ·ρλ·z2) (4)
其中Sλ表示标准照明体的光谱功率分布,ρλ表示各个色块的波长为λ对应的光谱反射率,x2、y2、z2表示国际照明委员会规定的CIE 1931标准色度观察者光谱三刺激值,K为系数,X、Y、Z表示测试图上色块在标准照明体下的XYZ值。
其中D为物体的密度值,Πλ表示波长为λ对应的光谱乘积,根据所选择的响应状态和模拟滤色片的不同,计算青色主密度时选择红光滤色片,计算品红色主密度时选择绿光滤色片,计算黄色主密度时选择蓝光滤色片,ISO标准有规定的Πλ值,可以从已有资料中查到。
步骤c根据每一色块青色色块的主密度值DC和X值XC,利用最小二乘法进行线性拟合,得到以DC为自变量,XC为因变量的线性函数根据每一青色色块的主密度值DC和Y值YC,利用最小二乘法进行线性拟合,得到以DC为自变量,YC为因变量的线性函数根据每一青色色块的主密度值DC和Z值ZC,利用最小二乘法进行线性拟合,得到以DC为自变量,ZC为因变量的线性函数
采用同样的方法,分别得到品红色的XYZ值随品红色主密度而变化的线性关系函数黄色的XYZ值随黄色主密度而变化的线性关系函数
以每一墨区内的青色色块主密度DC、同一墨区内的品红色块主密度DM为自变量,同一墨区内的蓝色块的X值XB为因变量,根据式(6)所示的二元二次函数利用最小二乘法进行拟合:
其中a0、a1、a2、a3、a4为经拟合得到的系数。采用同样的方法,得到蓝色块的Y值和青色块主密度、品红色块主密度之间的关系函数蓝色块的Z值和青色块主密度、品红色块主密度之间的关系函数
采用同样的方法,得到绿色块的X、Y、Z值和青色主密度、黄色主密度之间的关系函数
采用同样的方法,得到红色块的X、Y、Z值和品红色主密度、黄色主密度之间的关系函数
步骤d根据充分的生产实践经验,青、品红、黄、红、绿、蓝六色油墨在各种纸张上印刷时,每一实地色块的色度值与标准色度值之间的色差不可能大于100,所以预设初始的青、品红、黄、红、绿、蓝六色色差之和的最小值ΔEmin为600;
根据充分的生产实践经验,青、品红色实地色块的主密度不可能小于0.70或大于1.80,黄色实地色块的主密度不可能小于0.60或大于1.60,同时由于人眼所能察觉的密度变化阀值约为0.05,而目前的密度计检测精度一般为0.01。因此,本方法根据青色的XYZ值随着主密度值而变化的函数循环计算当青色主密度值从0.70开始,以0.01为间距,直至1.80时对应的XYZ值,再根据XYZ值计算L*a*b*值,计算L*a*b*值采用成熟的技术,计算方法如式(7)-(15)所示。
当
当
当
当
当
当
其中Xe、Ye、Ze表示白场的XYZ值,L*、a*、b*分别表示某一XYZ值对应的L*a*b*值中的L*值、a*值、b*值。
再计算该L*a*b*值与青色标准L*a*b*值之间的色差ΔEC,计算色差采用成熟的技术,计算方法如式(16)所示。
其中分别表示标准L*a*b*值,ΔE表示当前L*a*b*值与标准L*a*b*值之间的色差。
在循环计算青色色差随着青色主密度值逐渐增加而变化的过程中,针对青色主密度值每一次取值,根据品红色的XYZ值随着品红色主密度值而变化的函数蓝色的XYZ值随青色主密度和品红色主密度而变化的函数循环计算当品红色主密度值从0.70开始,以0.01为间距,直至1.80时对应的品红色XYZ值和蓝色XYZ值,再根据品红色的XYZ值计算品红色的L*a*b*值,根据蓝色的XYZ值计算蓝色的L*a*b*值,再计算品红色当前L*a*b*值与品红色标准L*a*b*值之间的色差ΔEM,计算蓝色当前L*a*b*值与蓝色标准L*a*b*值之间的色差ΔEB。
在循环计算品红色和蓝色的色差随着品红色主密度值逐渐增加而变化的过程中,针对品红色主密度值每一次取值,根据黄色的XYZ值随着黄色主密度值而变化的函数绿色的XYZ值随青色主密度和黄色主密度而变化的函数红色的XYZ值随品红色主密度和黄色主密度而变化的函数循环计算当黄色主密度值从0.60开始,以0.01为间距,直至1.60时对应的黄色XYZ值、绿色XYZ值和红色XYZ值,再根据黄色的XYZ值计算黄色的L*a*b*值,根据绿色的XYZ值计算绿色的L*a*b*值,根据红色的XYZ值计算红色的L*a*b*值,再计算黄色当前L*a*b*值与黄色标准L*a*b*值之间的色差ΔEY,计算绿色当前L*a*b*值与绿色标准L*a*b*值之间的色差ΔEG,计算红色当前L*a*b*值与红色标准L*a*b*值之间的色差ΔER。再计算青、品红、黄、红、绿、蓝六色的总色差ΔE6,计算方法如式(17)所示:
ΔE6=ΔEC+ΔEM+ΔEB+ΔEY+ΔEG+ΔR (17)
如果当前六色总色差值ΔE6小于最小色差值ΔEmin,令最小色差值ΔEmin等于当前色差值ΔE6,如式(18)所示:
ΔEmin=ΔE6 (18)
取当前的青、品红、黄密度为最合理密度DC_best、DM_best、DY_best,如式(19)、(20)、(21)所示:
DC_best=DC (19)
DM_best=DM (20)
DY_best=DY (21)
当全部循环计算完成后,程序中当前贮存的DC_best、DM_best、DY_best就分别是青色、品红色和黄色的最合理密度值,当前贮存的ΔEmin就是在青色、品红色和黄色的最合理密度值时,青、品红、黄、红、绿、蓝六色的总色差。
综上所述,本发明的方法是通过扫描检测色块的光谱反射率,计算获得所需的三刺激值、色度值、色差、主密度值,利用最小二乘法建立青、品红、黄、红、绿、蓝的三刺激值和青、品红、黄主密度值的关系,循环计算青、品红、黄三色主密度从低到高变化时,青、品红、黄、红、绿、蓝六色总色差值,取色差最小时的青、品红、黄的主密度值确定为最合理密度值。
虽然本发明已通过参考优选的实施例进行了图示和描述,但是,本专业普通技术人员应当了解,在权利要求书的范围内,可作形式和细节上的各种各样变化。
Claims (8)
1.一种确定印刷三原色密度的方法,其特征在于包括以下步骤:
a.以印刷机的一个墨区宽度为单位,设计一条长度小于纸张宽度的数字测控条,该测控条的每一墨区宽度内从左到右依次排列青、蓝、品红、绿、红、黄色色块,并在印前制作时在大版文件添加所述数字测控条;
b.印刷附着所述测控条的印张,测量印张上的测控条以记录测控条的每一墨区的光谱反射率数据,并计算青、品红、黄色色块的XYZ值和主密度值,计算红、绿、蓝色色块的XYZ值;
c.根据青、品红、黄色色块的主密度值D和X值X,利用最小二乘法进行线性拟合得到以D为自变量、X为因变量的线性函数f,并得到蓝、绿、红色块的XYZ值分别和青、品红、黄色色块主密度之间的关系函数f;
所述步骤c中利用最小二乘法根据二元二次多项式拟合蓝色块的XYZ值和同一墨区内的青、品红主密度值之间的关系,绿色块的XYZ值和同一墨区内的青、黄主密度值之间的关系,红色块的XYZ值和同一墨区内的品红、黄主密度值之间的关系;
d.预设初始的最小色差ΔEmin为600,根据每一个色块的XYZ随主密度值而变化的函数循环计算从0.70开始、以0.01为间距直至1.80时对应的XYZ值,在根据XYZ值计算L*a*b*值,再计算当前L*a*b*值与标准L*a*b*值之间的色差;
若六色总色差值小于最小色差值ΔEmin,则令最小色差值ΔEmin等于当前色差值;
其中,XYZ值指国际照明委员会创立的CIE 1931 XYZ颜色空间的坐标值,X值指颜色空间的X坐标值,Y值指颜色空间的Y坐标值,Z值指颜色空间的Z坐标值。
2.根据权利要求1所述的一种确定印刷三原色密度的方法,其特征在于步骤a还包括:当印刷机墨区数为奇数时,所述数字测控条为一个墨区宽度的奇数倍;当印刷机墨区数为偶数时,所述数字测控条为一个墨区宽度的偶数倍。
3.根据权利要求1所述的一种确定印刷三原色密度的方法,其特征在于所述数字测控条设于大版文件的拖梢边边缘并沿印刷滚筒轴向添加,所述数字测控条的轴向中心与大版文件的轴向中心对齐。
4.根据权利要求1所述的一种确定印刷三原色密度的方法,其特征在于所述测量印张的测控条的方式为利用扫描式分光光度计。
5.根据权利要求1所述的一种确定印刷三原色密度的方法,其特征在于所述步骤c还包括:
根据每一青色块的主密度值DC和X值XC,利用最小二乘法进行线性拟合,得到以DC为自变量,XC为因变量的线性函数并采用同样的方法,得到以青色块的Y值和Z值随着主密度而变化的函数和
采用同样的方法,分别得到品红色的XYZ值随品红色主密度而变化的线性关系函数黄色的XYZ值随黄色主密度而变化的线性关系函数
以每一墨区内的青色块主密度DC、同一墨区内的品红色块主密度DM为自变量,同一墨区内的蓝色块的X值XB为因变量,根据二元二次函数利用最小二乘法进行拟合:
其中a0、a1、a2、a3、a4为经拟合得到的系数;
采用同样的方法,得到蓝色块的Y值和青色块主密度、品红色块主密度之间的关系函数蓝色块的Z值和青色块主密度、品红色块主密度之间的关系函数
采用同样的方法,得到绿色块的XYZ值和青色主密度、黄色主密度之间的关系函数
采用同样的方法,得到红色块的XYZ值和品红色主密度、黄色主密度之间的关系函数
6.根据权利要求1所述的一种确定印刷三原色密度的方法,其特征在于所述步骤d还包括:
根据青色的XYZ值随着主密度值而变化的函数 循环计算当青色主密度值从0.70开始,以0.01为间距,直至1.80时对应的XYZ值,再根据XYZ值计算L*a*b*值,再计算青色当前L*a*b*值与标准L*a*b*值之间的色差ΔEC;
针对青色主密度值每一次取值,根据品红色的XYZ值随着品红色主密度值而变化的函数蓝色的XYZ值随青色主密度和品红色主密度而变化的函数循环计算当品红色主密度值从0.70开始,以0.01为间距,直至1.80时对应的品红色XYZ值和蓝色XYZ值,再根据品红色的XYZ值计算品红色的L*a*b*值,根据蓝色的XYZ值计算蓝色的L*a*b*值,再计算品红色当前L*a*b*值与品红色标准L*a*b*值之间的色差ΔEM,计算蓝色当前L*a*b*值与蓝色标准L*a*b*值之间的色差ΔEB;
针对品红色主密度值每一次取值,根据黄色的XYZ值随着黄色主密度值而变化的函数绿色的XYZ值随青色主密度和黄色主密度而变化的函数红色的XYZ值随品红色主密度和黄色主密度而变化的函数 循环计算当黄色主密度值从0.60开始,以0.01为间距,直至1.60时对应的黄色XYZ值、绿色XYZ值和红色XYZ值,再根据黄色XYZ值计算黄色的L*a*b*值,根据绿色XYZ值计算绿色的L*a*b*值,根据红色XYZ值计算红色的L*a*b*值,再计算黄色当前L*a*b*值与黄色标准L*a*b*值之间的色差ΔEY,计算绿色当前L*a*b*值与绿色标准L*a*b*值之间的色差ΔEG,计算红色当前L*a*b*值与红色标准L*a*b*值之间的色差ΔER;
计算青、品红、黄、红、绿、蓝六色的总色差ΔE6,计算公式表示为:
ΔE6=ΔEC+ΔEM+ΔEB+ΔEY+ΔEG+ΔER;
如果当前六色总色差值ΔE6小于最小色差值ΔEmin,令最小色差值ΔEmin等于当前色差值ΔE6,表示为:
ΔEmin=ΔE6。
7.根据权利要求1所述的一种确定印刷三原色密度的方法,其特征在于所述步骤d还包括:
取当前的青、品红、黄密度为最合理密度DC_best、DM_best、DY_best,表示为:
DC_best=DC;
DM_best=DM;
DY_best=DY。
8.根据权利要求7所述的一种确定印刷三原色密度的方法,其特征在于随青、品红、黄三色主密度的变化,以能够使青、品红、黄、红、绿、蓝六色总色差最小的密度确定为青、品红、黄三色最合理密度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611248758.2A CN106671594B (zh) | 2016-12-29 | 2016-12-29 | 一种确定印刷三原色密度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611248758.2A CN106671594B (zh) | 2016-12-29 | 2016-12-29 | 一种确定印刷三原色密度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106671594A CN106671594A (zh) | 2017-05-17 |
CN106671594B true CN106671594B (zh) | 2018-10-30 |
Family
ID=58873406
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611248758.2A Active CN106671594B (zh) | 2016-12-29 | 2016-12-29 | 一种确定印刷三原色密度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106671594B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110021048B (zh) * | 2019-03-22 | 2021-04-02 | 杭州电子科技大学上虞科学与工程研究院有限公司 | 一种报纸印刷的墨色评价方法 |
CN114987076A (zh) * | 2021-05-11 | 2022-09-02 | 广东新宏泽包装股份有限公司 | 一种确定印刷四色的最佳实地密度的方法 |
CN115122793B (zh) * | 2022-07-25 | 2024-03-22 | 印工社(青岛)数字科技有限公司 | 一种扫描式印刷测控条及其制作方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101365042A (zh) * | 2008-09-01 | 2009-02-11 | 西安理工大学 | 高保真打印输出L*a*b*图像的方法 |
CN103714245A (zh) * | 2013-12-19 | 2014-04-09 | 杭州电子科技大学 | 一种获取铜版纸印刷最佳实地密度的方法 |
CN105716995A (zh) * | 2016-03-18 | 2016-06-29 | 杭州电子科技大学 | 一种确定最佳印刷密度的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9259953B2 (en) * | 2013-09-27 | 2016-02-16 | Eastman Kodak Company | Tactile images having coefficient of friction differences |
-
2016
- 2016-12-29 CN CN201611248758.2A patent/CN106671594B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101365042A (zh) * | 2008-09-01 | 2009-02-11 | 西安理工大学 | 高保真打印输出L*a*b*图像的方法 |
CN103714245A (zh) * | 2013-12-19 | 2014-04-09 | 杭州电子科技大学 | 一种获取铜版纸印刷最佳实地密度的方法 |
CN105716995A (zh) * | 2016-03-18 | 2016-06-29 | 杭州电子科技大学 | 一种确定最佳印刷密度的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106671594A (zh) | 2017-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105716995B (zh) | 一种确定最佳印刷密度的方法 | |
CN106541702B (zh) | 一种印刷质量校正系统及其方法 | |
CN101518980B (zh) | 综合目测、密度测量与色度测量的印刷色彩品质检测方法 | |
CN105144680B (zh) | 光谱打印映射 | |
EP2651112B1 (en) | Color processing method, color processing device, and color processing system | |
CN106585090B (zh) | 一种印刷网点面积增大值的检测与校正方法 | |
CN106671594B (zh) | 一种确定印刷三原色密度的方法 | |
CN100589521C (zh) | 高保真打印输出L*a*b*图像的方法 | |
CN107395911B (zh) | 图像形成装置以及颜色变换控制方法 | |
NZ525412A (en) | Spectral color control method for printing ink supply | |
Milder et al. | Inkjet printers linearization using 3D gradation curves | |
CN107856409A (zh) | 一种基于吸收光谱的印刷专色配色方法 | |
CN102092206A (zh) | 一种用复频谱颜色特征数值检测印刷颜色的方法 | |
CN102729612A (zh) | 一种胶印专色油墨多目标优化配色方法 | |
CN107796773A (zh) | 一种凹版印刷专色半色调图像质量的控制方法 | |
CN107009734A (zh) | 装饰纸数码印刷色彩控制与还原方法 | |
CN102756549B (zh) | 基于观测光源的打样颜色一致性控制方法 | |
CN102238297B (zh) | 一种生成色彩对照文件的方法及系统 | |
CN107650508B (zh) | 一种胶印过程控制的机器视觉方法 | |
WO2002072355A9 (en) | Method and apparatus for controlling color of a printing press | |
CN201456554U (zh) | 一种调整油墨控制信息的系统 | |
CN105303527A (zh) | 一种光致褪色喷绘艺术品的图像复原方法 | |
FI111140B (fi) | Laatutietojen keruu offsetrotaatiopainossa | |
CN104869282B (zh) | 色谱方面的特性化数据的换算 | |
CN107728966A (zh) | 用于凹版印刷的数码彩喷配色工艺 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: Hangzhou City, Zhejiang province 310000 dragon Qianjiang Economic Development Zone Hollywood Road No. 103 building 1-4 Applicant after: Of safety science and Technology Co Ltd Address before: Hangzhou City, Zhejiang province 310000 dragon Qianjiang Economic Development Zone Hollywood Road No. 103 building 1-4 Applicant before: HANGZHOU HAOBO PRINTING CO., LTD. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |