CN104751819A - 一种显示设备反向特征化方法 - Google Patents
一种显示设备反向特征化方法 Download PDFInfo
- Publication number
- CN104751819A CN104751819A CN201510180007.0A CN201510180007A CN104751819A CN 104751819 A CN104751819 A CN 104751819A CN 201510180007 A CN201510180007 A CN 201510180007A CN 104751819 A CN104751819 A CN 104751819A
- Authority
- CN
- China
- Prior art keywords
- value
- tristimulus values
- rgb
- motivation
- primary colours
- 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
Abstract
本发明涉及一种显示设备的反向特征化方法,能有效解决现有技术中显示设备反向特征化时需要储存大量数据,特征化精度不够高的问题。其技术方案是:对于需要显示的颜色,利用其色度值XYZ初步计算其对应的驱动值RGB值,随后,将该RGB值代入显示设备正向特征化模型,计算其对应的三刺激值X′Y′Z′,并计算XYZ和X′Y′Z′之间的色差;若该色差大于指定的阈值,则重新调整驱动值RGB并重复上述过程,直到色差满足需求为止。本发明可以明显提高显示设备反向特征化精度,加快计算速度,提高显示设备再现颜色的精确度,减小显示设备颜色复现误差,可以显著提高显示设备显示图像的精度,适用于商用LCD显示器以及各类移动终端显示屏幕的反向特征化。
Description
技术领域
本发明涉及颜色显示、复制以及色彩管理领域,属于色彩复制再现技术领域中的显示设备反向特征化方法。具体涉及显示设备的反向特征化方法,尤其是一种基于三刺激值反向计算驱动值以及基于正向特征化模型计算的反向特征化方法。
背景技术
近年来,显示技术的发展非常迅速,彩色图像显示设备不断更新换代,传统的CRT显示设备几乎已经完全被各类LCD显示设备所代替。随着科学技术的不断进步,智能手机、平板电脑等移动终端的功能越来越强大,使用者经常会在智能手机、平板电脑等移动终端上观看视频、图片等信息,电脑不再是唯一的选择。
随着使用者对颜色复现与显示的精度要求越来越高,为保证颜色的正确显示,做好各类显示设备的特征化成为必要。再现一个已知色度值的颜色时,首先由该颜色的色度值经由反向特征化模型计算出其对应的驱动值,随后显示设备才能依据该驱动值显示颜色,一般情况下使用基于色度的特征化方法。LCD屏的显示设备与传统的CRT显示设备相比,其色品恒定性及通道可加性都较差,因此研究者提出了多种针对LCD显示屏的特征化方法,但这些模型大部分是针对专业液晶显示器提出的正向特征化模型,即已知驱动值,可以通过特征化模型计算出其对应的色度值。这些模型无法由已知的色度值计算出其对应的驱动值,而且基本上没有考虑环境光对显示屏显示颜色的影响;徐艳芳等人的研究结果表明,环境光照对显示设备的颜色再现的影响不容忽视。目前常用的显示设备反向特征化方法是查找表法,该方法需要存储大量数据,计算量大,计算速度慢,没有考虑环境光的影响。
发明内容
针对以上情况,本发明的目的在于提供一种适用于在常规使用环境下的,具有RGB三基色的普通商用液晶显示器或移动终端类显示屏的反向特征化方法。本发明考虑了环境光对显示设备再现颜色的影响,反向特征化过程通过构建数学模型,不需要存储大量的样本数据;本发明中反向特征化方法是在正向特征化方法的基础上提出的,与正向特征化使用的是同一个计算模型,正向特征化是由已知的驱动值求解其对应的三刺激值,反向特征化则是由已知的三刺激值求其对应的驱动值。正向模型分为两步,第一步假设通道可加,初步预测给定驱动值对应的三刺激值,第二步,分空间补偿第一步的预测误差;对于反向特征化模型,首先由目标颜色的三刺激值XYZ初步计算输入驱动值RGB,然后由正向模型预测该输入驱动值对应的三刺激值X1Y1Z1,计算预测值X1Y1Z1与已知三刺激值XYZ之间的色差,若色差小于设定的阈值,则计算值RGB就是已知三刺激值XYZ的颜色所对应的输入驱动值;否则,在修正计算得到的输入驱动值RGB后,再次用正向模型预测其对应三刺激值X′Y′Z′,直到预测值X′Y′Z′与已知三刺激值XYZ之间的色差满足设定的色差要求。
该方法是采取以下技术方案实现的。
一种显示设备的反向特征化方法,其特征在于其对显示设备实现反向特征化的步骤如下。
第一步:选取样本点,分别拟合显示设备各基色驱动值与主刺激值及色度值比值的关系,利用以上两类关系,可计算出各基色任意驱动值对应的三刺激值,利用通道可加性可以预测任意颜色的三刺激值,计算对该预测值进行补偿的补偿系数矩阵,建立正向特征化模型,其具体过程按以下具体步骤执行。
步骤1-1:拟合显示设备各基色驱动值与主刺激值的关系,其具体步骤如下。
步骤1-1-1:分别对显示设备各基色驱动值dR、dG、dB按选定间隔取19个样本点,测量各基色样本点的驱动值对应的三刺激值XmYmZm。
步骤1-1-2:分段拟合各基色的驱动值与主刺激值之间的关系,得到关系式XR=f 1(R)、YG=f 2(G)和ZB=f 3(B),其中,OP(O=X,Y,Z;P=R,G,B)表示各基色的三刺激值XYZ,XR、YG和ZB 分别是R基色、G基色和B基色的主刺激值;R、G、B分别表示R基色,G基色和B基色的单通道驱动值。
步骤1-2:拟合各基色驱动值与色度值比值的关系,其具体步骤如下。
步骤1-2-1:利用步骤1-1中测量样本点所得三刺激值,计算各样本点的色品值比值。对于R基色,其色品值比值分别是 和,G基色的色品值比值是和,B基色的色品值比值是和。
步骤1-2-2:分别拟合各基色驱动值与色度值比值的关系,R基色为(y R/x R)=F1(R)和(z R/x R)=F2(R),G基色为(x G/y G)=F3(G)和(z G/y G)=F4(G),B基色为(x B/z B)=F5(B)和(y B/z B)=F6(B)。
步骤1-3:计算各基色相等时,即灰色样本点对应的三刺激值,并将RGB颜色空间划分为多个子空间,其具体步骤如下。
步骤1-3-1:假设通道可加,将等驱动值下各基色的三刺激值相加,得到理论灰色在不同驱动值下的三刺激值XgrayYgrayZgray。
步骤1-3-2:对步骤1-3中的灰色样本点,测量其三刺激值Xgray1Ygray1Zgray1。
步骤1-3-3:计算灰色样本点的测量值 Xgray1Ygray1Zgray1与理论计算值XgrayYgrayZgray的差值,根据差值随驱动值的变化情况确定分界点,将对应各基色的驱动值分为多段,即将RGB空间划分为多个子空间。
步骤1-4:在上述步骤基础上,对任意已知驱动值RGB的颜色,能初步预测其驱动值对应的三刺激值,计算该预测值的补偿系数,其具体步骤如下。
步骤1-4-1:在每一个子空间中,选择样本点,测量这些样本点的驱动值RGB对应的三刺激值XRGBYRGBZRGB。
步骤1-4-2:根据通道可加性原理,利用步骤1-1和1-2中得到的关系式,初步预测这些样本点的驱动值RGB对应的三刺激值XRGB1YRGB1ZRGB1。
步骤1-4-3:计算样本点测量值XRGBYRGBZRGB与模型初步预测值XRGB1YRGB1ZRGB1的差值ΔXΔYΔZ。
步骤1-4-4:令A_RGB表示驱动值矢量的齐次坐标,B_Δ表示其对应的误差矢量,即A_RGB=[R G B 1],B_Δ=[ΔX ΔY ΔZ],T表示误差补偿系数矩阵T,则在每个子空间中的都有B_Δ=A_RGB*T。由于B_Δ和A_RGB的具体数据均已知,可采用伪逆法计算出误差系数矩阵T。
第二步:获取需要显示的颜色的色度值,如果该色度值是L*a*b*值,则将其转换为三刺激值XYZ。
第三步:由色度值三刺激值XYZ初步预测其对应的输入驱动值RGB,其计算过程按以下具体步骤执行。
步骤3-1:令ZB=Z,将ZB代入正向特征化方法中拟合的B值和其对应的主刺激值ZB的关系式ZB=f 3(B)中,计算驱动值B;利用正向特征化模型依次计算B值对应的三刺激值XB和YB。
步骤3-2:令YG=Y-YB,将YG代入正向特征化方法中拟合的G值和其对应的主刺激值YG值的关系式YG=f 2(G)中,计算驱动值G;利用正向特征化模型依次计算G值对应的三刺激值XG和ZG。
步骤3-3:令XR=X-XB-XG,将XR代入正向特征化方法中拟合的R值和其对应的主刺激值XR值的关系式XR=f 1(R)中,计算驱动值R;利用正向特征化模型依次计算R值对应的三刺激值YR和ZR。
第四步:将驱动值RGB代入正向特征化模型,计算其对应的三刺激值X′Y′Z′,计算过程按以下步骤进行。
步骤4-1:将R值代入关系式XR=f 1(R),G和B值分别代入关系式YG=f 2(G)和ZB=f 3(B)中,计算各基色的主刺激值XR,YG和ZB。
步骤4-2:将R值代入到关系式(y R/x R)=F1(R)、(z R/x R)=F2(R)中,G值代入到关系式(x G/y G)=F3(G)、(z G/y G)=F4(G)中,B值代入到关系式(x B/z B)=F5(B)、(y B/z B)=F6(B)中,计算各基色驱动值对应的色度值比值。
步骤4-3:依据三刺激值和色度值之间的关系,利用步骤4-1中的主刺激值和步骤4-2中的色品值比值,计算各基色的驱动值对应的另外两个刺激值。
步骤4-4:假设通道可加,计算模型初步预测的三刺激值X1Y1Z1。
步骤4-5:根据RGB值的大小,判断该颜色点所属的子空间,选择相应的误差补偿系数矩阵T。
步骤4-6:利用RGB值以及误差补偿系数矩阵T计算误差值,即三刺激值的补偿误差值ΔXΔYΔZ。
步骤4-7:将步骤4-4的计算结果和步骤4-6的计算结果相加,得到正向模型预测的三刺激值X′Y′Z′。
第五步:计算已知色度值XYZ和正向模型预测的X′Y′Z′之间的CIE1976色差,并用ΔE* 2表示。
第六步:判断ΔE* 2是否达到指定阈值要求,若ΔE* 2不大于指定阈值,则第三步求出的RGB值就是所求颜色色度值XYZ对应的驱动值,转第八步中止;若ΔE* 2大于指定阈值,则转第七步。
第七步:按照如下步骤调整RGB值,并重复第四步到第六步。
步骤7-1:比较计算值X′Y′Z′与已知值XYZ的差值,令ΔXX=(X′-X)/XRMAX,ΔYY=(Y′-Y)/YGMAX,ΔZZ=(Z′-Z)/ZBMAX,其中,X′Y′Z′是第四步所得计算值,XYZ是所求已知颜色对应驱动值的色度三刺激值,XRMAX是基色单通道R的驱动值为最大值时对应的主刺激值XR,YGMAX是基色单通道G的驱动值为最大值时对应的主刺激值YG,ZBMAX是基色单通道B的驱动值为最大值时对应的主刺激值ZB。
步骤7-2:比较|ΔXX|,|ΔYY|和|ΔZZ|,若|ΔXX|最大,转到步骤7-3,若|ΔYY|最大,转到步骤7-4,若|ΔZZ|最大,转到步骤7-5;若|ΔXX|=|ΔYY|,且都大于|ΔZZ|,或者|ΔXX|=|ΔZZ|,且都大于|ΔYY|,转到步骤7-3;若|ΔYY|=|ΔZZ|,且都大于|ΔXX|,转到步骤7-4。
步骤7-3:若ΔXX>0,用R-ΔR代替原来的R值;若ΔXX<0,用R+ΔR代替原来的R值。
步骤7-4:若ΔYY>0,用G-ΔG代替原来的G值;若ΔYY<0,用G+ΔG代替原来的G值。
步骤7-5:若ΔZZ>0,用B-ΔB代替原来的B值;若ΔZZ<0,用B+ΔB代替原来的B值。
第八步:结束。
本发明设计科学、推广性强、预测精度高,反向特征化过程中采用三刺激值初步计算驱动值,简单易用,结合正向特征化模型,用简单的迭代递推法逐步修正驱动值,拟合效果好,提高了模型的预测精度,正反向特征化采用相同的数学模型结构,便于工业生产的实际应用,利于推广应用。
附图说明
图1是一种显示设备反向特征化方法的流程示意图。
图2是显示设备正向特征化模型的流程示意图。
图3是由三刺激值XYZ初步计算驱动值RGB的流程示意图。
图4是由驱动值RGB值计算三刺激值X′Y′Z′的正向特征化流程示意图。
图5是调整驱动值RGB值的流程示意图。
具体实施例
下面结合附图,通过具体实施例对本发明作进一步详述,本实施例以IPAD2为研究对象,测量环境为在标准光源灯箱下,打开标准日光光源模拟常规使用环境的亮度,照度计测量垂直IPAD2表面的照度是320LX。以下实施例只是描述性的,不是限定性的,因此不能以此限定本发明的保护范围。
图1是一种显示设备反向特征化方法的流程示意图。
如图1所示,一种显示设备反向特征化方法的具体步骤应主要包括以下几步。
第一步:建立正向特征化模型,其具体过程按以下具体步骤执行。
图2是显示设备正向特征化模型的流程示意图。如图2所示,显示设备正向特征化应包括以下几步。
步骤1-1:拟合显示设备各基色驱动值与其主刺激值的关系,其具体步骤如下。
步骤1-1-1:分别对显示设备各基色驱动值dR、dG、dB按选定间隔取19个样本点,测量各基色样本点的驱动值对应的三刺激值XmYmZm。
步骤1-1-2:分段对各基色的驱动值与主刺激值之间的关系进行拟合,为方便计算,本实施例选择128作为驱动值的分界点,将RGB基色的驱动值分为两段,分段拟合驱动值与主刺激值之间的关系,得到关系式XR=f 1(R),YG=f 2(G)和ZB=f 3(B)。其中,XR是R基色的主刺激值,YG是G基色的主刺激值,ZB是B基色的主刺激值;R、G、B分别表示R基色,G基色和B基色的单通道驱动值。
步骤1-2:拟合各基色驱动值与色度值比值的关系,其具体步骤如下。
步骤1-2-1:利用步骤1-1中测量样本点所得三刺激值,计算不同驱动值下的色品值比值。对于R基色,其色品值比值分别是 和 ,G基色的色品值比值是和,B基色的色品值比值是和,色品值的计算公式如式1:
式1
其中,xyz是色品值,XmYmZm是基色样本点三刺激值的测量值。
步骤1-2-2:分别拟合各基色驱动值与色度值比值的关系,R基色为(y R/x R)=F1(R)和(z R/x R)=F2(R),G基色为(x G/y G)=F3(G)和(z G/y G)=F4(G),B基色为(x B/z B)=F5(B)和(y B/z B)=F6(B)。
步骤1-3:计算各基色相等时,即灰色样本点对应的三刺激值,并将RGB颜色空间划分为多个子空间,其具体步骤如下。
步骤1-3-1:假设通道可加,将等驱动值下各基色的三刺激值相加,得到理论灰色在不同驱动值下的三刺激值XgrayYgrayZgray,其计算公式如式2:
式2
其中,XgrayYgrayZgray是各基色驱动值相等时,即理论灰色的三刺激值。
步骤1-3-2:对步骤1-3中的灰色样本点,测量其三刺激值Xgray1Ygray1Zgray1。
步骤1-3-3:计算灰色样本点的测量值 Xgray1Ygray1Zgray1与理论计算值XgrayYgrayZgray的差值,根据差值随驱动值的变化情况确定分界点,将对应各基色的驱动值分为多段,即将RGB空间划分为多个子空间。
步骤1-4:对任意颜色,初步预测其驱动值RGB对应的三刺激值,计算该预测值的补偿系数,其具体步骤如下。
步骤1-4-1:在每一个子空间中,选择多个样本点,测量这些样本点的驱动值RGB对应的三刺激值XRGBYRGBZRGB。
步骤1-4-2:根据通道可加性原理,利用步骤1-1和1-2中得到的关系式,初步预测这些样本点的驱动值RGB对应的三刺激值XRGB1YRGB1ZRGB1。
步骤1-4-3:计算样本点测量值XRGBYRGBZRGB与初步预测值XRGB1YRGB1ZRGB1的差值ΔXΔYΔZ,计算公式如式3:
式3
其中,XRGBYRGBZRGB和XRGB1YRGB1ZRGB1分别由步骤1-4-1和1-4-2所得。
步骤1-4-4:令A_RGB表示驱动值矢量的齐次坐标,B_Δ表示其对应的误差矢量,即A_RGB=[R G B 1],B_Δ=[ΔX ΔY ΔZ],T表示误差补偿系数矩阵T,则在每个子空间中都有:
B_Δ=A_RGB*T 式4
由于B_Δ和A_RGB的具体数据均已知,可采用伪逆法计算出误差系数矩阵T。
第二步:获取需要显示的颜色的色度值,如果该色度值是L*a*b*值,则用色度学中L*a*b*色空间到XYZ色空间的转换公式,将其转换为三刺激值XYZ。
第三步:由色度值三刺激值XYZ初步预测其对应的输入驱动值RGB。
图3是由三刺激值XYZ初步计算驱动值RGB的流程示意图。
如图3所示,由三刺激值XYZ初步计算驱动值RGB的计算过程按以下具体步骤执行。
步骤3-1:令ZB=Z,将ZB代入正向特征化方法中拟合的B值和其对应的主刺激值ZB的关系式ZB=f 3(B)中,计算驱动值B;利用正向特征化模型依次计算B值对应的三刺激值XB和YB。
步骤3-2:令YG=Y-YB,将YG代入正向特征化方法中拟合的G值和其对应的主刺激值YG值的关系式YG=f 2(G)中,计算驱动值G;利用正向特征化模型依次计算G值对应的三刺激值XG和ZG。
步骤3-3:令XR=X-XB-XG,将XR代入正向特征化方法中拟合的R值和其对应的主刺激值XR值的关系式XR=f 1(R)中,计算驱动值R;利用正向特征化模型依次计算R值对应的三刺激值YR和ZR。
第四步:将计算所得驱动值RGB代入正向特征化模型,计算其对应的三刺激值X′Y′Z′。
图4是由驱动值RGB值计算三刺激值X′Y′Z′的正向特征化流程示意图。
如图4所示,由驱动值RGB值计算三刺激值X′Y′Z′的计算步骤如下。
步骤4-1:将R值代入关系式XR=f 1(R),G和B值分别代入关系式YG=f 2(G)和ZB=f 3(B)中,计算各基色的主刺激值XR,YG和ZB。
步骤4-2:将R值代入关系式(y R/x R)=F1(R)和(z R/x R)=F2(R)中,G值代入关系式(x G/y G)=F3(G)和(z G/y G)=F4(G)中,B值代入关系式(x B/z B)=F5(B)和(y B/z B)=F6(B)中,计算各基色驱动值对应的色度值比值。
步骤4-3:依据三刺激值和色度值之间的关系,利用步骤4-1中的主刺激值和步骤4-2中的色品值比值,计算各基色的驱动值对应的另外两个刺激值,计算公式如式5:
式5
其中,XYZ和xyz分别是下标RGB表示的各基色对应的三刺激值和色品值。
步骤4-4:假设通道可加,计算模型初步预测的三刺激值,计算公式如式6:
式6
其中,X1Y1Z1即是依据正向模型,对由已知色度值的颜色反向求得的驱动值RGB所对应的三刺激值进行预测,所得的预测值。
步骤4-5:根据RGB值的大小,判断该颜色点所属的子空间,选择相应的误差补偿系数矩阵T。
步骤4-6:利用RGB值以及误差补偿系数矩阵T计算误差值,即三刺激值的补偿误差值ΔXΔYΔZ。
步骤4-7:将步骤4-4的计算结果和步骤4-6的计算结果相加,得到正向模型预测的三刺激值X′Y′Z′,计算公式如式7:
式7
其中,X1Y1Z1是步骤4-4的计算值,ΔXΔYΔZ是步骤4-6计算的误差补偿值。
第五步:计算已知色度值XYZ和正向模型预测值X′Y′Z′之间的CIE1976色差,并用ΔE* 2表示。
第六步:判断ΔE* 2是否达到指定阈值要求,若ΔE* 2不大于指定阈值,则第三步求出的RGB值就是所求颜色色度值XYZ对应的驱动值,转第八步中止;若ΔE* 2大于指定阈值,转第七步。本实施例中指定阈值为3。
第七步:调整RGB值,并重复第四步到第六步。
图5是调整驱动值RGB值的流程示意图。
如图5所示,按如下步骤调整RGB值。
步骤7-1:比较计算值X′Y′Z′与已知值XYZ的差值,计算公式如式8:
式8
其中,X′Y′Z′是第四步所得计算值,XYZ是所求已知颜色对应驱动值的色度三刺激值,XRMAX是基色单通道R的驱动值为最大值时对应的主刺激值XR,YGMAX是基色单通道G的驱动值为最大值时对应的主刺激值YG,ZBMAX是基色单通道B的驱动值为最大值时对应的主刺激值ZB;本实施例中,各单通道基色的最大值均取为255。
步骤7-2:比较|ΔXX|,|ΔYY|和|ΔZZ|,若|ΔXX|最大,转到步骤7-3,若|ΔYY|最大,转到步骤7-4,若|ΔZZ|最大,转到步骤7-5;若|ΔXX|=|ΔYY|,且都大于|ΔZZ|,或者|ΔXX|=|ΔZZ|,且都大于|ΔYY|,转到步骤7-3;若|ΔYY|=|ΔZZ|,且都大于|ΔXX|,转到步骤7-4。
步骤7-3:若ΔXX>0,用R-ΔR代替原来的R值;若ΔXX<0,用R+ΔR代替原来的R值,本实施例中取ΔR=1。
步骤7-4:若ΔYY>0,用G-ΔG代替原来的G值;若ΔYY<0,用G+ΔG代替原来的G值,本实施例中取ΔG=1。
步骤7-5:若ΔZZ>0,用B-ΔB代替原来的B值;若ΔZZ<0,用B+ΔB代替原来的B值,本实施例中取ΔB=1。
第八步:结束。
Claims (4)
1.一种显示设备的反向特征化方法,其特征在于,其对于显示设备实现反向特征化的步骤如下:
第一步:选取样本点,分别拟合显示设备各基色驱动值与主刺激值及各基色驱动值与色度值比值的关系,利用以上两类关系式,可计算出各基色任意驱动值对应的三刺激值,利用通道可加性可以预测任意颜色的三刺激值,计算对该预测值进行补偿的补偿系数矩阵,建立正向特征化模型;
第二步:获取需要显示的颜色的色度值,如果该色度值是L*a*b*值,则将其转换为三刺激值XYZ;
第三步:由色度值三刺激值XYZ初步预测其对应的输入驱动值RGB;
第四步:将驱动值RGB代入正向特征化模型,计算其对应的三刺激值X′Y′Z′;
第五步:计算已知色度值XYZ和正向模型预测的X′Y′Z′之间的CIE1976色差ΔE* 2;
第六步:判断ΔE* 2是否达到指定阈值要求,若ΔE* 2不大于指定阈值,则第三步求出的RGB值就是所求颜色色度值XYZ对应的驱动值,转第八步中止;若ΔE* 2大于指定阈值,转第七步;
第七步:调整RGB值,并重复第四步到第六步;
第八步:结束。
2.根据权利要求1所述的一种移动终端显示屏幕的特征化方法,其特征在于第三步,由色度值三刺激值XYZ初步预测其对应的输入驱动值RGB,其具体步骤如下:
步骤3-1:令ZB=Z,将ZB代入正向特征化方法中拟合的B值和其对应的主刺激值ZB的关系式ZB=f 3(B)中,计算驱动值B;利用正向特征化模型依次计算B值对应的三刺激值XB和YB;
步骤3-2:令YG=Y-YB,将YG代入正向特征化方法中拟合的G值和其对应的主刺激值YG值的关系式YG=f 2(G)中,计算驱动值G;利用正向特征化模型依次计算G值对应的三刺激值XG和ZG;
步骤3-3:令XR=X-XB-XG,将XR代入正向特征化方法中拟合的R值和其对应的主刺激值XR值的关系式XR=f 1(R)中,计算驱动值R;利用正向特征化模型依次计算R值对应的三刺激值YR和ZR。
3.根据权利要求1所述的一种移动终端显示屏幕的特征化方法,其特征在于第四步,将驱动值RGB代入正向特征化模型,计算其对应的三刺激值X′Y′Z′,其具体步骤如下:
步骤4-1:将R值代入关系式XR=f 1(R),G和B值分别代入关系式YG=f 2(G)和ZB=f 3(B)中,计算各基色的主刺激值XR,YG和ZB;
步骤4-2:将R值代入关系式(y R/x R)=F1(R)和(z R/x R)=F2(R)中,G值代入关系式(x G/y G)=F3(G)和(z G/y G)=F4(G)中,B值代入关系式(x B/z B)=F5(B)和(y B/z B)=F6(B)中,计算各基色驱动值RGB对应的色度值比值;
步骤4-3:依据三刺激值和色度值之间的关系,利用步骤4-1中的主刺激值和步骤4-2中的色品值比值,计算各基色的驱动值对应的另外两个刺激值,计算公式如式1:
式1
其中,XYZ和xyz分别是下标RGB表示的各基色对应的三刺激值和色品值;
步骤4-4:假设通道可加,计算模型初步预测的三刺激值,计算公式如式2:
式2
其中,X1Y1Z1即是依据正向模型,对由已知色度值的颜色反向求得的驱动值RGB所对应的三刺激值进行预测,所得的预测值;
步骤4-5:根据RGB值的大小,判断该颜色点所属的子空间,并选择相应的误差补偿系数矩阵T;
步骤4-6:利用RGB值以及误差补偿系数矩阵T计算误差值,即三刺激值的补偿误差值ΔXΔYΔZ;
步骤4-7:将步骤4-4的计算结果和步骤4-6的计算结果相加,得到正向模型预测的三刺激值X′Y′Z′,计算公式如式3:
式3
其中,X1Y1Z1是步骤4-4的计算值,ΔXΔYΔZ是步骤4-6计算的误差补偿值。
4.根据权利要求1所述的一种移动终端显示屏幕的特征化方法,其特征在于第七步,调整RGB值,其具体步骤如下:
步骤7-1:比较计算值X′Y′Z′与已知值XYZ的差值,计算公式如式4:
式4
其中,X′Y′Z′是第四步所得计算值,XYZ是所求已知颜色对应驱动值的色度三刺激值,XRMAX是基色单通道R的驱动值为最大值时对应的主刺激值XR,YGMAX是基色单通道G的驱动值为最大值时对应的主刺激值YG,ZBMAX是基色单通道B的驱动值为最大值时对应的主刺激值ZB;
步骤7-2:比较|ΔXX|,|ΔYY|和|ΔZZ|,若|ΔXX|最大,转到步骤7-3,若|ΔYY|最大,转到步骤7-4,若|ΔZZ|最大,转到步骤7-5;若|ΔXX|=|ΔYY|,且都大于|ΔZZ|,或者|ΔXX|=|ΔZZ|,且都大于|ΔYY|,转到步骤7-3;若|ΔYY|=|ΔZZ|,且都大于|ΔXX|,转到步骤7-4;
步骤7-3:若ΔXX>0,用R-ΔR代替原来的R值,;若ΔXX<0,用R+ΔR代替原来的R值;
步骤7-4:若ΔYY>0,用G-ΔG代替原来的G值;若ΔYY<0,用G+ΔG代替原来的G值;
步骤7-5:若ΔZZ>0,用B-ΔB代替原来的B值;若ΔZZ<0,用B+ΔB代替原来的B值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510180007.0A CN104751819B (zh) | 2015-04-16 | 2015-04-16 | 一种显示设备反向特征化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510180007.0A CN104751819B (zh) | 2015-04-16 | 2015-04-16 | 一种显示设备反向特征化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104751819A true CN104751819A (zh) | 2015-07-01 |
CN104751819B CN104751819B (zh) | 2017-04-12 |
Family
ID=53591386
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510180007.0A Expired - Fee Related CN104751819B (zh) | 2015-04-16 | 2015-04-16 | 一种显示设备反向特征化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104751819B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105606543A (zh) * | 2016-03-29 | 2016-05-25 | 温州大学 | 一种便携式香蕉成熟度检测仪 |
CN106997751A (zh) * | 2016-01-25 | 2017-08-01 | 曲阜师范大学 | 一种显示设备的反向光谱特征化方法 |
CN108462863A (zh) * | 2018-02-11 | 2018-08-28 | 上海健康医学院 | 一种基于复合模型的显示设备色空间转换方法 |
CN109147683A (zh) * | 2017-06-28 | 2019-01-04 | 上海健康医学院 | 一种复合的显示设备的反向特征化方法 |
CN110880306A (zh) * | 2019-11-01 | 2020-03-13 | 南京图格医疗科技有限公司 | 一种医用显示器色彩还原校正方法 |
CN112071257A (zh) * | 2020-07-31 | 2020-12-11 | 佛山市青松科技股份有限公司 | 一种led屏校正方法、装置、存储介质及led屏 |
CN112286604A (zh) * | 2020-09-28 | 2021-01-29 | 广州视源电子科技股份有限公司 | 显示画面的颜色调整方法、装置和计算机设备 |
CN112735352A (zh) * | 2020-12-31 | 2021-04-30 | Oppo广东移动通信有限公司 | 显示设备色温自适应的调整方法、装置、设备及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006038980A (ja) * | 2004-07-23 | 2006-02-09 | Sony Corp | 色温度補正用データ形成装置 |
CN101188096A (zh) * | 2007-12-18 | 2008-05-28 | 上海广电集成电路有限公司 | 一种液晶显示器的色域映射算法和颜色管理系统 |
CN102511050A (zh) * | 2009-09-30 | 2012-06-20 | 伊斯曼柯达公司 | 用于优化显示器配置文件的方法 |
CN103354087A (zh) * | 2013-07-01 | 2013-10-16 | 广东威创视讯科技股份有限公司 | 拼接显示装置颜色校正的方法和系统 |
CN104219512A (zh) * | 2014-09-30 | 2014-12-17 | 曲阜师范大学 | 一种显示设备色域边界描述方法 |
-
2015
- 2015-04-16 CN CN201510180007.0A patent/CN104751819B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006038980A (ja) * | 2004-07-23 | 2006-02-09 | Sony Corp | 色温度補正用データ形成装置 |
CN101188096A (zh) * | 2007-12-18 | 2008-05-28 | 上海广电集成电路有限公司 | 一种液晶显示器的色域映射算法和颜色管理系统 |
CN102511050A (zh) * | 2009-09-30 | 2012-06-20 | 伊斯曼柯达公司 | 用于优化显示器配置文件的方法 |
CN103354087A (zh) * | 2013-07-01 | 2013-10-16 | 广东威创视讯科技股份有限公司 | 拼接显示装置颜色校正的方法和系统 |
CN104219512A (zh) * | 2014-09-30 | 2014-12-17 | 曲阜师范大学 | 一种显示设备色域边界描述方法 |
Non-Patent Citations (1)
Title |
---|
徐艳芳等: "基于色度比特性的显示器色度特性化", 《液晶与显示》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106997751A (zh) * | 2016-01-25 | 2017-08-01 | 曲阜师范大学 | 一种显示设备的反向光谱特征化方法 |
CN105606543A (zh) * | 2016-03-29 | 2016-05-25 | 温州大学 | 一种便携式香蕉成熟度检测仪 |
CN109147683A (zh) * | 2017-06-28 | 2019-01-04 | 上海健康医学院 | 一种复合的显示设备的反向特征化方法 |
CN108462863A (zh) * | 2018-02-11 | 2018-08-28 | 上海健康医学院 | 一种基于复合模型的显示设备色空间转换方法 |
CN108462863B (zh) * | 2018-02-11 | 2020-06-05 | 上海健康医学院 | 一种基于复合模型的显示设备色空间转换方法 |
CN110880306A (zh) * | 2019-11-01 | 2020-03-13 | 南京图格医疗科技有限公司 | 一种医用显示器色彩还原校正方法 |
CN112071257A (zh) * | 2020-07-31 | 2020-12-11 | 佛山市青松科技股份有限公司 | 一种led屏校正方法、装置、存储介质及led屏 |
CN112071257B (zh) * | 2020-07-31 | 2021-12-10 | 佛山市青松科技股份有限公司 | 一种led屏校正方法、装置、存储介质及led屏 |
CN112286604A (zh) * | 2020-09-28 | 2021-01-29 | 广州视源电子科技股份有限公司 | 显示画面的颜色调整方法、装置和计算机设备 |
CN112286604B (zh) * | 2020-09-28 | 2024-04-16 | 广州视源电子科技股份有限公司 | 显示画面的颜色调整方法、装置和计算机设备 |
CN112735352A (zh) * | 2020-12-31 | 2021-04-30 | Oppo广东移动通信有限公司 | 显示设备色温自适应的调整方法、装置、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN104751819B (zh) | 2017-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104751819A (zh) | 一种显示设备反向特征化方法 | |
CN107945729B (zh) | 转换方法及电路、显示装置及驱动方法和电路、存储介质 | |
CN107665684B (zh) | 一种色彩Mura补偿方法 | |
CN104981861B (zh) | 信号转换装置和方法 | |
CN111243550B (zh) | 调节装置、显示屏的伽马电压和白平衡的调节方法与介质 | |
CN107784975B (zh) | Amoled显示装置的亮度和色度自动调整方法及系统 | |
CN103218988B (zh) | 一种rgb信号到rgbw信号的图像转换方法及装置 | |
CN103080999B (zh) | 用于具有背光调制的显示器的自适应色彩校正 | |
CN102394041B (zh) | 白平衡调整方法 | |
US10291892B2 (en) | White balance method of four-color pixel system | |
CN107863083B (zh) | 一种显示装置的驱动方法及驱动装置 | |
WO2015196608A1 (zh) | 一种rgb信号到rgbw信号的图像转换方法及装置 | |
CN102509541B (zh) | 色彩调整装置、色彩调整方法以及显示器 | |
CN104969285B (zh) | 电子设备、显示控制装置及其方法 | |
KR102587865B1 (ko) | 표시장치 및 그의 영상 처리방법 | |
WO2019071777A1 (zh) | 显示驱动方法、装置及设备 | |
CN101866621A (zh) | 一种led显示屏校正系统和方法 | |
WO2019061677A1 (zh) | 一种OLED显示模组的3Gamma校正方法及系统 | |
KR102601853B1 (ko) | 표시장치 및 그의 영상 처리방법 | |
US20180336811A1 (en) | Method and system for regulating brightness and chromaticity of display panel | |
CN110223650A (zh) | 四色像素的白平衡参数生成方法和显示装置 | |
CN106023929A (zh) | 显示装置的白平衡调整方法及其系统 | |
CN112292851B (zh) | 色域校正方法及装置 | |
WO2015172336A1 (zh) | 图像的显示方法以及显示系统 | |
KR20100073648A (ko) | 표시 장치의 색 보정 방법 및 장치 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170412 |
|
CF01 | Termination of patent right due to non-payment of annual fee |