CN115366236B - 一种实现数智化颜色实时调整的方法 - Google Patents
一种实现数智化颜色实时调整的方法 Download PDFInfo
- Publication number
- CN115366236B CN115366236B CN202211079861.4A CN202211079861A CN115366236B CN 115366236 B CN115366236 B CN 115366236B CN 202211079861 A CN202211079861 A CN 202211079861A CN 115366236 B CN115366236 B CN 115366236B
- Authority
- CN
- China
- Prior art keywords
- color
- data
- ink
- production
- printing
- 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
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B28—WORKING CEMENT, CLAY, OR STONE
- B28B—SHAPING CLAY OR OTHER CERAMIC COMPOSITIONS; SHAPING SLAG; SHAPING MIXTURES CONTAINING CEMENTITIOUS MATERIAL, e.g. PLASTER
- B28B11/00—Apparatus or processes for treating or working the shaped or preshaped articles
- B28B11/001—Applying decorations on shaped articles, e.g. by painting
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B28—WORKING CEMENT, CLAY, OR STONE
- B28B—SHAPING CLAY OR OTHER CERAMIC COMPOSITIONS; SHAPING SLAG; SHAPING MIXTURES CONTAINING CEMENTITIOUS MATERIAL, e.g. PLASTER
- B28B17/00—Details of, or accessories for, apparatus for shaping the material; Auxiliary measures taken in connection with such shaping
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B28—WORKING CEMENT, CLAY, OR STONE
- B28B—SHAPING CLAY OR OTHER CERAMIC COMPOSITIONS; SHAPING SLAG; SHAPING MIXTURES CONTAINING CEMENTITIOUS MATERIAL, e.g. PLASTER
- B28B17/00—Details of, or accessories for, apparatus for shaping the material; Auxiliary measures taken in connection with such shaping
- B28B17/0063—Control arrangements
- B28B17/0081—Process control
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/001—Texturing; Colouring; Generation of texture or colour
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Ceramic Engineering (AREA)
- Mechanical Engineering (AREA)
- Structural Engineering (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Facsimile Image Signal Circuits (AREA)
Abstract
本发明公开了一种实现数智化颜色实时调整的方法,包括以下步骤:S1、提前将其做好环境色彩管理;S2、将瓷砖样板生产文件A0录入数据库,生成其样板制造的发色情况及范围,打印并且烧制;S3、对烧制的样板色卡在生产环境M0的发色情况下进行色彩数据采集与分析;S4、再次生产时,选取生产文件A0在现下的生产环境M1生产相应发色的生产色卡,打印烧制,读取信息;S5、将数据建立色彩信息样本库和新环境数据信息库使用计算机通过计算公式样本库的色彩信息比对,将数据进行校准;S6、通过再次生产文件数据生成生产环境M1的样板;该实现数智化颜色实时调整的方法可防止瓷砖的不同批次的颜色出现色差,提高陶瓷的色彩质量的同时,能够提高陶瓷的生产效率。
Description
技术领域
本发明涉及陶瓷砖色彩处理的技术领域,具体的说,尤其是一种实现数智化颜色实时调整的方法。
背景技术
陶瓷的喷墨印刷使得瓷砖装饰化、个性化、小批量、多花色的发展趋向。由于瓷砖在生产时通常都是有多条线进行生产,虽然第一批生产的瓷砖颜色相同,但同一款产品在同一生产线,在不同时间再次生产时,生产条件已经发生变化,例如喷墨打印设备所使用的墨水品牌、窑炉温度等变化,打印出来的颜色也会发生变化,这种情况同样发生在同一款产品更换生产线再次生产时,颜色也会与上次生产完全不同,从而导致瓷砖的不同批次的颜色出现色差,而且跟第一批次相比,会越来越偏离。
陶瓷中的色彩数据模型与常规的印刷色是不一样的。在陶瓷印刷系统中,基本印刷色和平面艺术中使用的青涩、品红、黄色和黑色相差甚远。陶瓷基本颜色墨水色有:蓝色-棕色-黄色-黑色-橘色-粉色-绿色,特殊颜色墨水:包裹色红-包裹色黄,一般常用于生产线的组合有:蓝棕黄黑橘、蓝棕黄黑橘包裹红、蓝棕黄黑橘包裹红黄、蓝棕橘黑黄、蓝棕橘黑黄包裹红、蓝棕橘黑黄包裹红黄、蓝棕包裹黄黑包裹红等组合,不同于传统印刷行业,陶瓷印刷中墨水的顺序不同都会行程完全不同的发色情况,传统的墨水色阶一般是5-10为色阶,灰度信息控制全凭人为经验。
因为陶瓷的色彩是受多方面影响,如坯体的干湿程度及用料配方,釉面的干湿程度,平整度,窑炉的温度高低波动,喷墨机喷墨的压力力度,甚至是天气的温度影响等等,在生产时通常都是有多条线进行生产,虽然第一批生产的瓷砖颜色相同,但同一款产品在同一生产线,在不同时间再次生产时,生产条件已经发生变化,这种情况同样发生在同一款产品更换生产线再次生产时,颜色也会与上次生产完全不同,从而导致瓷砖的不同批次的颜色出现色差,而且跟第一批次相比,会越来越偏离。
发明内容
为了克服现有技术存在的缺陷,本发明提供一种实现数智化颜色实时调整的方法,旨在解决上述背景技术中的问题。
本发明解决其技术问题所采用的技术方案是:一种实现数智化颜色实时调整的方法,包括以下步骤:
S1、针对拍摄的相机,提前将其做好环境色彩管理;
S2、将瓷砖样板生产文件A0录入数据库,生成其样板制造的发色情况及范围,打印并且进行烧制;
S3、使用相机对烧制的样板色卡在生产环境M0的发色情况下进行色彩数据采集与分析,其中,采集的色彩数据为Lab,通过颜色转换转换关系Tp(·),分析转换至相对应生产线的打印设备的陶瓷喷墨印刷墨水通道数据:
A(L0,a0,b0)=Tp[A(C0,M0,Y0,K0,......Sn0)]作为样板色彩数据,并录入形成数据建立色彩信息样本库,作为下次在生产时的数据比对依据;
S4、再次生产时,选取生产文件A0在现下的生产环境M1生产相应发色的生产色卡,打印烧制,读取信息,其中,采集色卡数据CIELab,CIEXYZ后,通过颜色转换转换关系TN(·),分析转换至相对应生产线的打印设备的陶瓷喷墨印刷墨水通道数据:
Ax(Lx,ax,bx)=TN[Ax(Cx,Mx,Yx,Kx,......Snx)]作为再次生产当下发色校准数据,保存生成新环境数据信息库,作为再次生产文件Ax;
S5、将数据建立色彩信息样本库和新环境数据信息库使用计算机通过计算公式样本库的色彩信息比对,将数据进行校准,生成新环境生产的校准曲线,其校准公式为:
ΔP=Snx-Sn0;
S6、通过校准曲线校准处理生成再次生产文件数据,通过再次生产文件数据生成生产环境M1的样板。
进一步的,所述步骤S1中的环境色彩管理还包括方法:
灯光色温光源校准;
调节LCC白平衡;
使用Krome camera软件的Colorcheck模块测量数码相机色彩复原性针对此环境制作相应的色彩配置文件,以用于后面拍摄瓷砖色卡时统一色彩配置文件。
进一步的,所述打印设备是直接接收软件导入的CMYK的数值,并按该数值进行控制墨量相应的发色,输出与现下生产波动因素导致的墨量发色相匹配的CMYK的墨量值,通过计算机的计算公式样本库的色彩信息比对计算输出相应的墨量,控制墨量输出。
进一步的,所述生产文件A0和再次生产文件A1进行比对,根据印刷色度数据匹配转换关系公式转换关系:f-1(x)=f(x+p),并且对其进行校准至瓷砖样板发色情况A0。
进一步的,所述步骤S5中的校准的方法为印刷模式与印刷模式之间的转换,根据校准原理模型的创建的规则,将不同的生产线墨水数量及顺序,将其转换至相匹配的生产线墨水,设置是已改变的CMYK+Sn至样板的文件A0打印数据CMYK+Sn,匹配关系为:
TNC(·),获得样板制造的发色多通道印刷色与再次生产的发色多通道印刷色,进行对比及分色转换:
Ax(C1,M1,Y1,K1)=TNc[Ax(C1,M1,Y1,K1)]
针对生产样板使用的S种油墨分量[S1,S2,…,SM]则可写成[CMYK+S1,S2…,Sn],已生产样板色彩数据根据再次生产印刷色度数据匹配转换关系公式:f-1(x)=f(x+p),获得“变量CMYK+Sn”匹配至“样板的CMYK+S”。
具体地,所述校准原理模型的创建的规则是通过相机读取每种油墨CMYK+Sn不同的灰度值,在不同的环境下所表现的色彩值Lab&XYZ,并得到相关的曲线特性和数据分析以及存储。
更具体地,所述油墨为Sn组时,就会有n组的Lab的曲线特征图,每种油墨每一等级的灰度都对应了一定的Lab值,N组油墨的组合数对应了/>组的Lab数值,而且这些数据均可被表达在CIELab和XYZ色域空间里,记做f(y)=f(L,a,b);
所述油墨为多组组成的N维空部分色域间时,颜色矫正公式为:
f(y)=f-1(x),则有fx(Lx,ax,bx)=TN[Ax(Cx,Mx,Yx,Kx,,,Snx)],油墨打印输入灰度值均为正整数。
更具体地,所述校准原理模型的建立方法包括:
设一组色彩C0可以表示为M0环境下数值为[C0,M0,Y0,K0...Sn0]=Tnp[L0,a0,b0],另一M1环境下两组数值变为色彩C1:[C0,M0,Y0,K0...Sn0]=Tnp[L1,a1,b1],[C1,M1,Y1,K1...Sn1]=Tnp[L0,a0,b0],要使得M1环境下Tnp[L1,a1,b1]≈Tnp[L0,a0,b0],那么在CIELab空间C1[L1,a1,b1]点要往C0[L0,a0,b0]点位移,同理对应的N维空间点C1[C1,M1,Y1,K1...Sn1]要往C0[C0,M0,Y0,K0...Sn0]位移,假设在N维[C,M,Y,K...Snp]空间里存在一增量p=[△c,△m,△y,△k,...△Sn]值,使得[C1,M1,Y1,K1...Sn1]=[C0+△c,M0+△m,Y0+△y,K0+△k,...Sn0+△Sn]=Tnp[L0,a0,b0],就表示色彩C0在M1环境下被矫正或者说色彩C0被复制到了M1环境里,用函数式表达即为:
f(Lx,ax,bx)=npf-1(Cx,Mx,Yx,Kx...Snx)=npf-1(C0+△c,M0+△m,Y0+△y,K0+△k,...Sn0+△Sn)
其中,X代表环境状况,不是序数词,0代表原始条件,目的是两条件下颜色得到复制f(Lx,ax,bx)=f(L0,a0,b0)。
更具体地,所述计算机的数据计算包括单组油墨色彩的灰度阶调层次的处理以及空间数据查找表的建立与使用,函数关系为:
其中,x为原灰度值,y为层次转换曲线后的灰度值,YL,Ya,Yb灰度值与x值为非线性曲线关系,YL,Ya,Yb灰度值与x值满足对数函数关系:
y=uevx+k,进行变换后Y=vlnx+u;
其中Y=lny,x=lnx,Y可以形成一组lnv为斜率和u为截距的曲线,Y=vlnx+u可简化转换成线性方程y=vx+u;
同时,对油墨可以进行数据测量获得数值(x,y),同时绘制出散点图,对散点图进行最小二乘法原理求线性方程和方差分析,其中,y=f(lnx,w),其中w=[w1,w2,.....wn]T为待定系数,为寻找y=f(lnx,w)的参数w的最优估计值,对于给定m组为m>n,观测数据(xi,yi)(i=1,2,3...m),求目标函数:
由最小二乘法拟合可以得出待定系数:
w=[w1,w2,.....wn]T,
其中,y=vx+u满足关系式
y=w0+w1X1+w2X2,.....+wnXn,
并形成方程组:
YL=wL0+wL1X1+wL2X2,.....+wLnXn
Ya=wa0+wa1X1+wa2X2,.....+wanXn
Yb=wb0+wb1X1+wb2X2,.....+wbnXn
记
从而形成回归方程E(yi)=w0+w1x1+w2x2,.....+wnxn,
得到样本多元线性回归方程
回归模型方差值
拟合优度
R2代表了数据的拟合线性的好坏程度,值越小代表拟合度好。
更具体地,所述灰度阶调层次的处理以及空间数据查找表的建立的建立方法包括:
采用大数据数字查找表的方法进行数字式样本库集数据比对转换,对每组油墨可以根据空间转换函数关系式,建立一个一维的数据查找表,原图像的灰度值x为表项索引数据值,y为对应矫正目标值存储于表中,利用每种油墨灰度值查找表,给与索引x值后可以找到对应每种油墨对应的目标灰度值y数组,对原图像所有目标灰度值进行转换,即可获得转换后的灰度值构成的新图像。
本发明的有益效果是:通过组合形成的校准曲线进行对比,以实时通过曲线的波动变化进行反作用调整,针对已经生产过的产品,进行再次复产的颜色及信息校准至样板的过程,从而令到当前生产的产品的颜色保持与第一次生产的产品颜色相同,防止瓷砖的不同批次的颜色出现色差,提高陶瓷的色彩质量的同时,能够提高陶瓷的生产效率,而且整个色彩调节过程是将数据生成校准曲线来进行,陶瓷的色彩不受坯体的干湿程度及用料配方、釉面的干湿程度、平整度、窑炉的温度高低波动等因素的影响。
附图说明
图1为本发明的颜色调整方法的流程图。
图2为变量CMYK+Sn匹配至样板的CMYK+S的关系图。
图3为本发明实施例中的Lab曲线图。
图4为本发明实施例中的色彩样板数据模块CMYK+X的曲线图。
图5为本发明实施例中的蓝色阶的组合校准曲线图。
图6为本发明实施例中的棕色阶的组合校准曲线图。
图7为本发明实施例中的黄色阶的组合校准曲线图。
图8为本发明实施例中的黑色阶的组合校准曲线图。
图9为本发明实施例中的油墨X的Lab值分析图。
图10为本发明实施例中的原图原始灰度级和处理后图像灰度级的曲线关系图。
图11为本实施例中的灰度值查表示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明。在此需要说明的是,对于这些实施方式的说明用于帮助理解本发明,但并不构成对本发明的限定。此外,下面所描述的本发明各个实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示,一种实现数智化颜色实时调整的方法,包括以下步骤:
S1、针对拍摄的相机,提前将其做好环境色彩管理,其中环境色彩管理方法包括:
灯光色温光源校准;
调节LCC白平衡;
使用Krome camera软件的Colorcheck模块测量数码相机色彩复原性针对此环境制作相应的色彩配置文件,以用于后面拍摄瓷砖色卡时统一色彩配置文件。
S2、将瓷砖样板生产文件A0录入数据库,生成其样板制造的发色情况及范围,打印并且进行烧制;
S3、使用相机对烧制的样板色卡在生产环境M0的发色情况下进行色彩数据采集与分析,其中,采集的色彩数据为Lab,通过颜色转换转换关系Tp(·),分析转换至相对应生产线的打印设备的陶瓷喷墨印刷墨水通道数据:
A(L0,a0,b0)=Tp[A(C0,M0,Y0,K0,......Sn0)]作为样板色彩数据,并录入形成数据建立色彩信息样本库,作为下次在生产时的数据比对依据;
S4、再次生产时,选取生产文件A0在现下的生产环境M1生产相应发色的生产色卡,打印烧制,读取信息,其中,采集色卡数据CIELab,CIEXYZ后,通过颜色转换转换关系TN(·),分析转换至相对应生产线的打印设备的陶瓷喷墨印刷墨水通道数据:
Ax(Lx,ax,bx)=TN[Ax(Cx,Mx,Yx,Kx,......Snx)]作为再次生产当下发色校准数据,保存生成新环境数据信息库,作为再次生产文件Ax;
S5、将数据建立色彩信息样本库和新环境数据信息库使用计算机通过计算公式样本库的色彩信息比对,将数据进行校准,生成新环境生产的校准曲线,其校准公式为:
ΔP=Snx-Sn0;
其中,校准曲线所使用的校准的方法为印刷模式与印刷模式之间的转换,根据校准原理模型的创建的规则,将不同的生产线墨水数量及顺序,将其转换至相匹配的生产线墨水,设置是已改变的CMYK+Sn至样板的文件A0打印数据CMYK+Sn,匹配关系为:
TNC(·),获得样板制造的发色多通道印刷色与再次生产的发色多通道印刷色,进行对比及分色转换:
Ax(C1,M1,Y1,K1)=TNc[Ax(C1,M1,Y1,K1)],
参见图2所示,针对生产样板使用的S种油墨分量[S1,S2,…,SM]则可写成[CMYK+S1,S2…,Sn],已生产样板色彩数据根据再次生产印刷色度数据匹配转换关系公式:f-1(x)=f(x+p),获得“变量CMYK+Sn”匹配至“样板的CMYK+S”。
S6、通过校准曲线校准处理生成再次生产文件数据,通过再次生产文件数据生成生产环境M1的样板。
本实施例通过组合形成的校准曲线进行对比,以实时通过曲线的波动变化进行反作用调整,针对已经生产过的产品,进行再次复产的颜色及信息校准至样板的过程,从而令到当前生产的产品的颜色保持与第一次生产的产品颜色相同,防止瓷砖的不同批次的颜色出现色差,提高陶瓷的色彩质量的同时,能够提高陶瓷的生产效率,而且整个色彩调节过程是将数据生成校准曲线来进行,陶瓷的色彩不受坯体的干湿程度及用料配方、釉面的干湿程度、平整度、窑炉的温度高低波动等因素的影响。
需要说明的是,打印设备是直接接收软件导入的CMYK的数值,并按该数值进行控制墨量相应的发色,输出与现下生产波动因素导致的墨量发色相匹配的CMYK的墨量值,通过计算机的计算公式样本库的色彩信息比对计算输出相应的墨量,控制墨量输出,其关系如下:
颜色控制值(CMYK+Xn)→墨水设备打印机→输出色(ClELab,ClEXYZ,CMYK+Sn墨量);
本实施例中的生产文件A0和再次生产文件A1进行比对,根据印刷色度数据匹配转换关系公式转换关系:f-1(x)=f(x+p),并且对其进行校准至瓷砖样板发色情况A0。
如图3所示,本实施例中的校准原理模型的创建的规则是通过相机读取每种油墨CMYK+Sn不同的灰度值,在不同的环境下所表现的色彩值Lab&XYZ,并得到相关的曲线特性和数据分析以及存储;本实施例主要采用了油墨的Lab值的数据分析,如图3所示的一组不同油墨在某一特定环境条件下的Lab的数据特征图(纵向表示L,a,b的数据值,横向表示墨水输入量)。
当油墨为Sn组时,就会有n组的Lab的曲线特征图,每种油墨每一等级的灰度都对应了一定的Lab值,N组油墨的组合数对应了/>组的Lab数值,而且这些数据均可被表达在CIELab和XYZ色域空间里,记做f(y)=f(L,a,b);
当油墨为多组组成的N维空部分色域间时,颜色矫正公式为:
f(y)=f-1(x),则有fx(Lx,ax,bx)=TN[Ax(Cx,Mx,Yx,Kx,,,Snx)],油墨打印输入灰度值均为正整数。
如图4所示,根据每个需要用的墨水进行1:1的解析相应的发色情况,录入信息处理制成色彩样板数据模块,形成色彩样板数据模块CMYK+X。
如图5至8所示,蓝棕黄黑色阶的每个墨水在实际生产时的波动情况,蓝色参数是原始的,橙色参数是因生产情况(生产线窑炉烧成曲线的变化,打印设备的更换或参数的变化,油墨的变化,介质物化性质的变更,气候氛围等等)变更,发色的数据也发生实时变化,而我们的校准则是检测到数据变动后,与我们的彩样板数据模块样板的标准值数据进行一一比对后,计算出数据的变化值,根据数据的波动变化反馈作用到原图像计算出新的灰度值,此生成新的灰度值图像在变化后的环境里生产可以保证和样板颜色一致。
本实施例的校准原理模型的建立方法包括:设一组色彩C0可以表示为M0环境下数值为[C0,M0,Y0,K0...Sn0]=Tnp[L0,a0,b0],另一M1环境下两组数值变为色彩C1:[C0,M0,Y0,K0...Sn0]=Tnp[L1,a1,b1],[C1,M1,Y1,K1...Sn1]=Tnp[L0,a0,b0],要使得M1环境下Tnp[L1,a1,b1]≈Tnp[L0,a0,b0],那么在CIELab空间C1[L1,a1,b1]点要往C0[L0,a0,b0]点位移,同理对应的N维空间点C1[C1,M1,Y1,K1...Sn1]要往C0[C0,M0,Y0,K0...Sn0]位移,假设在N维[C,M,Y,K...Snp]空间里存在一增量p=[△c,△m,△y,△k,...△Sn]值,使得[C1,M1,Y1,K1...Sn1]=[C0+△c,M0+△m,Y0+△y,K0+△k,...Sn0+△Sn]=Tnp[L0,a0,b0],就表示色彩C0在M1环境下被矫正或者说色彩C0被复制到了M1环境里,用函数式表达即为:
f(Lx,ax,bx)=npf-1(Cx,Mx,Yx,Kx...Snx)=npf-1(C0+△c,M0+△m,Y0+
△y,K0+△k,...Sn0+△Sn)
其中,X代表环境状况,不是序数词,0代表原始条件,目的是两条件下颜色得到复制f(Lx,ax,bx)=f(L0,a0,b0)。
如图9所示的油墨X的Lab值分析图,当把色彩记录空间表达为RGB或CIEXYZ时一样可以的转换成CMYK+Sn的公式,Tsp[X,Y,Z]=[C,M,Y,K...Sn],此处不再赘述。
由于N组油墨所组成的色域空间,在做信息空间转换,查找差值信息比对时,当N值很大时,数据库数据量会非常庞大,运算量较大,时间较长,为了简化计算,把CMYK+Sn色域空间N+4组每种油墨单组单独矫正,即有f(Lxi,axi,bxi)=f(Xni),因为当每组单元油墨,也有其对应的发色空间CIELab,它只是混合空间中的一特殊部分空间,当每组油墨都被加以矫正后,其合成的色彩Lab值也就被加以纠正。
本实施中的计算机的数据计算包括单组油墨色彩的灰度阶调层次的处理以及空间数据查找表的建立与使用,函数关系为:
其中,x为原灰度值,y为层次转换曲线后的灰度值,同一x值一定的情况下,y点可以看做散落在非线性曲线/>围一定离散的点,Y值取值范围标准化为,YL值范围(0-255),Ya(-127~127),Yb(-127~127),油墨组灰度值取x(0-100),假设取值都无限连续,由Lab曲线图我们的看到YL,Ya,Yb灰度值与x值为非线性曲线关系,对图像曲线做拟合曲线,YL,Ya,Yb灰度值与x值满足对数函数关系:
y=uevx+k,进行变换后Y=vlnx+u;
其中Y=lny,x=lnx,Y可以形成一组lnv为斜率和u为截距的曲线,Y=vlnx+u可简化转换成线性方程y=vx+u;
同时,对油墨可以进行数据测量获得数值(x,y),同时绘制出散点图,对散点图进行最小二乘法原理求线性方程和方差分析,其中,y=f(lnx,w),其中w=[w1,w2,.....wn]T为待定系数,为寻找y=f(lnx,w)的参数w的最优估计值,对于给定m组为m>n,观测数据(xi,yi)(i=1,2,3...m),求目标函数:
由最小二乘法拟合可以得出待定系数:
w=[w1,w2,.....wn]T,
其中,y=vx+u满足关系式
y=w0+w1X1+w2X2,.....+wnXn,
并形成方程组:
YL=wL0+wL1X1+wL2X2,.....+wLnXn
Ya=wa0+wa1X1+wa2X2,.....+wanXn
Yb=wb0+wb1X1+wb2X2,.....+wbnXn
记
从而形成回归方程E(yi)=w0+w1x1+w2x2,.....+wnxn,
得到样本多元线性回归方程
回归模型方差值
拟合优度
R2代表了数据的拟合线性的好坏程度,值越小代表拟合度好。
设定一定值可以做数据的参考:
由YL,Ya,Yb灰度值与x值满足对数函数关系:
y=uevx+k,进行变换后Y=vlnx+u,(其中Y=lny),Y可以形成一组lnv为斜率和u为截距的曲线,可以通过给定的(lnx,y)与(v,u)关系建立函数矩阵式对数据进行标记色彩灰度值:
因此当测到每种油墨的曲线后,可得到每个油墨对应x灰度下的YL,Ya,Yb,vLi,vai,vbi,uLi,uai,ubi数值,任何环境下的数值均可进行参数的采集和计算,那么在环境发生变化时要得到同一YL,Ya,Yb值时,只需对数据参数进行比对,就可以得到另一灰度X值。
如图10和图11所示,本实施例中的灰度阶调层次的处理以及空间数据查找表的建立的建立方法包括:
采用大数据数字查找表的方法进行数字式样本库集数据比对转换,对每组油墨可以根据空间转换函数关系式,建立一个一维的数据查找表,原图像的灰度值x为表项索引数据值,y为对应矫正目标值存储于表中,利用每种油墨灰度值查找表,给与索引x值后可以找到对应每种油墨对应的目标灰度值y数组,对原图像所有目标灰度值进行转换,即可获得转换后的灰度值构成的新图像。
本实施例还包括有数据准确性验证及补偿步骤,由于每组油墨在打印设备里参数多为正整数值(y取值范围0-100),而且每通道墨水在叠加过程中也会有溢墨,出色域超值等问题,造成计算很多的误差值等问题,甚至还会有查不到数据等问题,因为超色域时得到的矫正值y>100,或y<0,所以需要对数值做取整,以及数据准确性验证及补偿。
使用的是逐组油墨,逐个目标像素点按其像素空间矩阵目标函数进行的逐一计算查找,忽略多油墨组合物造成的油墨量堆积以及油墨堆积叠序对色彩呈现的影响,建立一多维数据集来规避这一问题,由于当目标油墨组数较多时会造成色彩数据集数据量非常庞大,造成运算非常缓慢,图像处理效率会非常的低,运行缓慢,计算机内存等硬件要求会很高,但是多维空间数据集的建立对颜色的矫正还是最全面的。
CIELab和CIEXYZ可看做一个三维的数据集,N组油墨打印机的组合色彩可以看做为一N维空间(可看作是N维部分区域点空间组合,因为油墨灰度值范围0-100的整数值)。若有函数fx(Lx,ax,bx)=TN[Ax(Cx,Mx,Yx,Kx,...Snx)],按颜色的转换模型,计算并建立多维空间数据集,以及CIELab和CIEXYZ可看做一个三维的数据集。在数据集里原始颜色空间空间结构排列一定数量的“节点”,在“节点”对应的存储单元内,节点数量的计算如,当有N维空间,每一维度坐标分为m等分,转换处理以后的目标空间维度是S,每一维度坐标分为p等分,以及CIELab和CIEXYZ,每一维度坐标分为255等分。则对于N维空间节点数Q=(m+1)N,S维空间节点数Q=(p+1)s,CIELab空间节点数Q=2563.存储了转换后的目标颜色数据。
即不同情况下的输出的结果依然是Lab色彩信息的表现是一致的,把不同环境下的色彩颜色做全部的采集建立一个Lab三维色彩数据集,此时两种环境的不同维度的数据节点都同CIELab&XYZ空间做一个颜色数据集,两个维度N,S空间的数据集的数据量分别为DN=3(m+1)N,DS=3Q=(p+1)s(这里取一个数据占一个字节)。节点的数据计算,符合颜色模型函数,同时其节点数均在不同维度空间数量有所不同,但其控制构成的空间结构是可以对应包含的查找对比的。
若目标数据集里的节点不包含原始颜色数据所对应的目标数据,则需要利用数据集里已有的其他节点数据,经颜色插值计算,获得需要的目标颜色。即当给定一图像信息[C0M0Y0K0...Sn0=L0a0b0],经信息数据集建立后录用查找,即可得到环境变化后[L0a0b0=C1M1Y1K1...Sn1],此值就是我们查找到环境变化后油墨输入组数据C1M1Y1K1...Sn1。
由于在多组油墨混合时,在一多维空间数据集里会存在有一组Lab会有多组CMYK...Sn的多种配比组合,也有一定的误差性,但单一组油墨Lab差值计算的单组油墨灰度的误差极小,基本可以对应查找,需要对单组油墨查找数据组[YC,YM,YY,YK...YSn]和多维空间数据集查找油墨灰度值[YC1,YM1,YY1,YK1...YSn1]进行校对与补偿,当ΔY=|YSn-YSn1|min时,两组目标值的y值的平均值才是需要的数值,即y=YSn+Y补偿,而
两组目标值的y值信息需要做一个准确性验证,可以通过误差大小做一个检验,及数据要满足最小乘法原理和曲线的高度拟合,设y值与参与色彩油墨混合数Ki和查找信息误差绝对值ΔY存在函数条件关系式:
y=f(Ki,ΔY),
可设定条件取值,当Ki=1,无关ΔY值或ΔY=0,y=YSn,使用单组油墨计算值,
当Ki<4,若ΔY>2时,y=YSn1,使用多空间数据集查找值;
当Ki<4,若0<ΔY<2时,y=YSn+Y补偿,两者的数值均值;
当Ki>4,若无关ΔY值时,y=YSn1,使用多空间数据集查找值;
当对目标色误差值可以接受的更大的话,还可以通过调整ΔY的大小进行修改。
以上结合附图对本发明的实施方式作了详细说明,但本发明不限于所描述的实施方式。对于本领域的技术人员而言,在不脱离本发明原理和精神的情况下,对这些实施方式进行多种变化、修改、替换和变型,仍落入本发明的保护范围内。
Claims (10)
1.一种实现数智化颜色实时调整的方法,其特征在于,包括以下步骤:
S1、针对拍摄的相机,提前将其做好环境色彩管理;
S2、将瓷砖样板生产文件A0录入数据库,生成其样板制造的发色情况及范围,打印并且进行烧制;
S3、使用相机对烧制的样板色卡在生产环境M0的发色情况下进行色彩数据采集与分析,其中,采集的色彩数据为Lab,通过颜色转换转换关系Tp(·),分析转换至相对应生产线的打印设备的陶瓷喷墨印刷墨水通道数据:
A(L0,a0,b0)=Tp[A(C0,M0,Y0,K0,......Sn0)]作为样板色彩数据,并录入形成数据建立色彩信息样本库,作为下次在生产时的数据比对依据;
S4、再次生产时,选取生产文件A0在现下的生产环境M1生产相应发色的生产色卡,打印烧制,读取信息,其中,采集色卡数据CIELab,CIEXYZ后,通过颜色转换转换关系TN(·),分析转换至相对应生产线的打印设备的陶瓷喷墨印刷墨水通道数据:
Ax(Lx,ax,bx)=TN[Ax(Cx,Mx,Yx,Kx,......Snx)]作为再次生产当下发色校准数据,保存生成新环境数据信息库,作为再次生产文件Ax;
S5、将数据建立色彩信息样本库和新环境数据信息库使用计算机通过计算公式样本库的色彩信息比对,将数据进行校准,生成新环境生产的校准曲线,其校准公式为:
ΔP=Snx-Sn0;
S6、通过校准曲线校准处理生成再次生产文件数据,通过再次生产文件数据生成生产环境M1的样板。
2.根据权利要求1所述的一种实现数智化颜色实时调整的方法,其特征在于,所述步骤S1中的环境色彩管理还包括方法:
灯光色温光源校准;
调节LCC白平衡;
使用Krome camera软件的Colorcheck模块测量数码相机色彩复原性针对此环境制作相应的色彩配置文件,以用于后面拍摄瓷砖色卡时统一色彩配置文件。
3.根据权利要求1所述的一种实现数智化颜色实时调整的方法,其特征在于,所述打印设备是直接接收软件导入的CMYK的数值,并按该数值进行控制墨量相应的发色,输出与现下生产波动因素导致的墨量发色相匹配的CMYK的墨量值,通过计算机的计算公式样本库的色彩信息比对计算输出相应的墨量,控制墨量输出。
4.根据权利要求1所述的一种实现数智化颜色实时调整的方法,其特征在于,所述生产文件A0和再次生产文件A1进行比对,根据印刷色度数据匹配转换关系公式转换关系:f-1(x)=f(x+p),并且对其进行校准至瓷砖样板发色情况A0。
5.根据权利要求1所述的一种实现数智化颜色实时调整的方法,其特征在于,所述步骤S5中的校准的方法为印刷模式与印刷模式之间的转换,根据校准原理模型的创建的规则,将不同的生产线墨水数量及顺序,将其转换至相匹配的生产线墨水,设置是已改变的CMYK+Sn至样板的文件A0打印数据CMYK+Sn,匹配关系为:
TNC(·),获得样板制造的发色多通道印刷色与再次生产的发色多通道印刷色,进行对比及分色转换:
Ax(C1,M1,Y1,K1)=TNc[Ax(C1,M1,Y1,K1)]
针对生产样板使用的S种油墨分量[S1,S2,…,SM]则可写成[CMYK+S1,S2…,Sn],已生产样板色彩数据根据再次生产印刷色度数据匹配转换关系公式:f-1(x)=f(x+p),获得“变量CMYK+Sn”匹配至“样板的CMYK+S”。
6.根据权利要求5所述的一种实现数智化颜色实时调整的方法,其特征在于,所述校准原理模型的创建的规则是通过相机读取每种油墨CMYK+Sn不同的灰度值,在不同的环境下所表现的色彩值Lab&XYZ,并得到相关的曲线特性和数据分析以及存储。
7.根据权利要求6所述的一种实现数智化颜色实时调整的方法,其特征在于,所述油墨为Sn组时,就会有n组的Lab的曲线特征图,每种油墨每一等级的灰度都对应了一定的Lab值,N组油墨的组合数对应了/>组的Lab数值,而且这些数据均可被表达在CIELab和XYZ色域空间里,记做f(y)=f(L,a,b);
所述油墨为多组组成的N维空部分色域间时,颜色矫正公式为:
f(y)=f-1(x),则有fx(Lx,ax,bx)=TN[Ax(Cx,Mx,Yx,Kx,,,Snx)],
油墨打印输入灰度值均为正整数。
8.根据权利要求5所述的一种实现数智化颜色实时调整的方法,其特征在于,所述校准原理模型的建立方法包括:
设一组色彩C0可以表示为M0环境下数值为[C0,M0,Y0,K0...Sn0]=Tnp[L0,a0,b0],另一M1环境下两组数值变为色彩C1:[C0,M0,Y0,K0...Sn0]=Tnp[L1,a1,b1],[C1,M1,Y1,K1...Sn1]=Tnp[L0,a0,b0],要使得M1环境下Tnp[L1,a1,b1]≈Tnp[L0,a0,b0],那么在CIELab空间C1[L1,a1,b1]点要往C0[L0,a0,b0]点位移,同理对应的N维空间点C1[C1,M1,Y1,K1...Sn1]要往C0[C0,M0,Y0,K0...Sn0]位移,假设在N维[C,M,Y,K...Snp]空间里存在一增量p=[△c,△m,△y,△k,...△Sn]值,使得[C1,M1,Y1,K1...Sn1]=[C0+△c,M0+△m,Y0+△y,K0+△k,...Sn0+△Sn]=Tnp[L0,a0,b0],就表示色彩C0在M1环境下被矫正或者说色彩C0被复制到了M1环境里,用函数式表达即为:
f(Lx,ax,bx)=npf-1(Cx,Mx,Yx,Kx...Snx)=npf-1(C0+△c,M0+△m,Y0+△y,K0+△k,...Sn0+△Sn)
其中,X代表环境状况,不是序数词,0代表原始条件,目的是两条件下颜色得到复制f(Lx,ax,bx)=f(L0,a0,b0)。
9.根据权利要求3所述的一种实现数智化颜色实时调整的方法,其特征在于,所述计算机的数据计算包括单组油墨色彩的灰度阶调层次的处理以及空间数据查找表的建立与使用,函数关系为:
其中,x为原灰度值,y为层次转换曲线后的灰度值,YL,Ya,Yb灰度值与x值为非线性曲线关系,YL,Ya,Yb灰度值与x值满足对数函数关系:
y=uevx+k,进行变换后Y=vlnx+u;
其中Y=lny,x=lnx,Y可以形成一组lnv为斜率和u为截距的曲线,Y=vlnx+u可简化转换成线性方程y=vx+u;
同时,对油墨可以进行数据测量获得数值(x,y),同时绘制出散点图,对散点图进行最小二乘法原理求线性方程和方差分析,其中,y=f(lnx,w),其中w=[w1,w2,.....wn]T为待定系数,为寻找y=f(lnx,w)的参数w的最优估计值,对于给定m组为m>n,观测数据(xi,yi)(i=1,2,3...m),求目标函数:
由最小二乘法拟合可以得出待定系数:
w=[w1,w2,.....wn]T,
其中,y=vx+u满足关系式
y=w0+w1X1+w2X2,.....+wnXn,
并形成方程组:
YL=wL0+wL1X1+wL2X2,.....+wLnXn
Ya=wa0+wa1X1+wa2X2,.....+wanXn
Yb=wb0+wb1X1+wb2X2,.....+wbnXn
从而形成回归方程E(yi)=w0+w1x1+w2x2,.....+wnxn,
得到样本多元线性回归方程
回归模型方差值
拟合优度
R2代表了数据的拟合线性的好坏程度,值越小代表拟合度好。
10.根据权利要求9所述的一种实现数智化颜色实时调整的方法,其特征在于,所述灰度阶调层次的处理以及空间数据查找表的建立的建立方法包括:
采用大数据数字查找表的方法进行数字式样本库集数据比对转换,对每组油墨可以根据空间转换函数关系式,建立一个一维的数据查找表,原图像的灰度值x为表项索引数据值,y为对应矫正目标值存储于表中,利用每种油墨灰度值查找表,给与索引x值后可以找到对应每种油墨对应的目标灰度值y数组,对原图像所有目标灰度值进行转换,即可获得转换后的灰度值构成的新图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211079861.4A CN115366236B (zh) | 2022-09-05 | 2022-09-05 | 一种实现数智化颜色实时调整的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211079861.4A CN115366236B (zh) | 2022-09-05 | 2022-09-05 | 一种实现数智化颜色实时调整的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115366236A CN115366236A (zh) | 2022-11-22 |
CN115366236B true CN115366236B (zh) | 2023-09-19 |
Family
ID=84068691
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211079861.4A Active CN115366236B (zh) | 2022-09-05 | 2022-09-05 | 一种实现数智化颜色实时调整的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115366236B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101110892A (zh) * | 2006-07-18 | 2008-01-23 | 佳能株式会社 | 图像处理设备和方法以及信息处理设备和方法 |
WO2014070358A1 (en) * | 2012-11-02 | 2014-05-08 | Electronics For Imaging, Inc. | Generation of a white ink separation |
CN105898105A (zh) * | 2015-02-17 | 2016-08-24 | 柯尼卡美能达株式会社 | 颜色精度验证预测方法以及控制装置 |
CN107009734A (zh) * | 2017-04-13 | 2017-08-04 | 广东省南方彩色制版有限公司 | 装饰纸数码印刷色彩控制与还原方法 |
CN111510580A (zh) * | 2020-05-14 | 2020-08-07 | 佛山市阿瑞斯陶瓷科技有限公司 | 陶瓷色彩校准方法、装置、计算机设备及存储介质 |
CN114338948A (zh) * | 2022-01-17 | 2022-04-12 | 佛山市阿瑞斯数字设备有限公司 | 一种陶瓷生产数字化的颜色实时调整的方法 |
-
2022
- 2022-09-05 CN CN202211079861.4A patent/CN115366236B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101110892A (zh) * | 2006-07-18 | 2008-01-23 | 佳能株式会社 | 图像处理设备和方法以及信息处理设备和方法 |
WO2014070358A1 (en) * | 2012-11-02 | 2014-05-08 | Electronics For Imaging, Inc. | Generation of a white ink separation |
CN104885068A (zh) * | 2012-11-02 | 2015-09-02 | 图像电子公司 | 白色油墨分离的生成 |
CN105898105A (zh) * | 2015-02-17 | 2016-08-24 | 柯尼卡美能达株式会社 | 颜色精度验证预测方法以及控制装置 |
CN107009734A (zh) * | 2017-04-13 | 2017-08-04 | 广东省南方彩色制版有限公司 | 装饰纸数码印刷色彩控制与还原方法 |
CN111510580A (zh) * | 2020-05-14 | 2020-08-07 | 佛山市阿瑞斯陶瓷科技有限公司 | 陶瓷色彩校准方法、装置、计算机设备及存储介质 |
CN114338948A (zh) * | 2022-01-17 | 2022-04-12 | 佛山市阿瑞斯数字设备有限公司 | 一种陶瓷生产数字化的颜色实时调整的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115366236A (zh) | 2022-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR100278567B1 (ko) | 3 차원 색소를 보다 높은 차원의 색소로 변환시키기 위한 시스템, 방법 및 프로그램 | |
EP1267217B1 (en) | Calibration of a multi color imaging system using a predicted color shift | |
US6698860B2 (en) | Spectral color reproduction with six color output | |
US7602532B2 (en) | Highly accurate and rapid scanning by a simple scanning device and color correction technology for a printing device | |
CN107395919B (zh) | 一种通过修改特性文件实现数字印刷颜色校正的方法 | |
CN100589521C (zh) | 高保真打印输出L*a*b*图像的方法 | |
CN109920022B (zh) | 一种制作数码印花小色差校色色卡的方法 | |
CN109120815B (zh) | 颜色选择用比色图表及其生成方法、以及生成装置 | |
CN102082893A (zh) | 数码打样校色方法及装置 | |
CN108621567B (zh) | 用于在喷墨印刷机中进行印刷过程的方法 | |
KR100278430B1 (ko) | 외부적으로 정의된 4 차원 색소를 등가의 4차원 색소로 변환시키기 위한 시스템, 방법 및 프로그램 | |
US7782490B2 (en) | Sensitivity matrix determination for adaptive color control | |
CN102180003B (zh) | 一种打印机光谱特征化模型的修正方法 | |
CN115366236B (zh) | 一种实现数智化颜色实时调整的方法 | |
CN1406757A (zh) | 自动彩色映像变换 | |
US7134737B2 (en) | Closed-loop color correction using factory-measured color cutoffs anchored to field-measured white point | |
JP5736382B2 (ja) | 最適化されたプリンタ・キャリブレーション方法 | |
CN101692690A (zh) | 一种远程校色方法 | |
CN103978783B (zh) | 一种用于静电数字印刷网点阶调的控制方法 | |
CN113553012B (zh) | 多通道打印设备输出效果一致性调整系统与方法 | |
JP2019527151A (ja) | 複数のモジュールを備えたプリンタの較正 | |
JP2002237962A (ja) | カラープロファイルを測定ないし換算する方法および装置 | |
JP2019140502A (ja) | キャリブレーション用色変換テーブルの生成方法、生成プログラム、生成装置、及び、生成システム | |
US11155099B2 (en) | Printer density control mechanism | |
CN107738512B (zh) | 基于icc编辑的改进瓷砖喷墨印刷色彩对样方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |