CN102075667A - 一种基于查表法的反向色彩空间转换方法 - Google Patents
一种基于查表法的反向色彩空间转换方法 Download PDFInfo
- Publication number
- CN102075667A CN102075667A CN 201110028562 CN201110028562A CN102075667A CN 102075667 A CN102075667 A CN 102075667A CN 201110028562 CN201110028562 CN 201110028562 CN 201110028562 A CN201110028562 A CN 201110028562A CN 102075667 A CN102075667 A CN 102075667A
- Authority
- CN
- China
- Prior art keywords
- ciel
- subspace
- look
- color space
- value
- 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
- Color Image Communication Systems (AREA)
- Facsimile Image Signal Circuits (AREA)
- Image Processing (AREA)
Abstract
本发明公开的一种基于查表法的反向色彩空间转换方法,具体按照以下步骤实施:建立设备色彩空间到CIEL*a*b*色彩空间的查找表;分割CIEL*a*b*空间,建立自定义子空间查找表;计算CIEL*a*b*子空间的几何中心;获取CIEL*a*b*子空间的建模数据;建立CIEL*a*b*子空间的多项式转换模型;确定给定待转换点所在的CIEL*a*b*子空间;计算给定待转换点对应的设备空间的颜色值。本发明基于查表法的反向色彩空间转换方法,解决了现有采用三维查表插值法反向变换时的查表不易、转换速度慢的问题,在进行转换时,既能保持转换精度又能使得转换速度满足实用要求。
Description
技术领域
本发明属于色彩管理技术领域,涉及一种从CIEL*a*b*设备无关色彩空间到设备色彩空间的转换方法,具体涉及一种基于查表法的反向色彩空间转换方法。
背景技术
设备色彩空间主要有两种,一种是扫描仪及显示器对应的RGB色彩空间,另一种是彩色打印机或印刷机对应的CMY(或CMYK)色彩空间。CIEL*a*b*色彩空间是国际照明委员会定义的一种与设备无关的,且感知均匀的色彩空间。一般称设备色彩空间向CIEL*a*b*色彩空间的转换为正向转换,而CIEL*a*b*色彩空间向设备色彩空间的转换为反向转换。
基于ICC的色彩管理通过将颜色从一台设备的色彩空间转换到设备无关色彩空间,再从设备无关色彩空间转换到另一台设备的色彩空间,以达到颜色在不同设备之间传输时失真最小的目的。常用的设备无关色彩空间为CIEL*a*b*色彩空间,因此CIEL*a*b*色彩空间到设备色彩空间转换是色彩管理的核心问题之一,转换的精度和速度决定了转换方法能否满足实际需求。
常用的色彩空间转换方法有多项式回归法、三维查表插值法、神经网络法、模糊逻辑法等几类,其中,三维查表插值法是广泛使用的一种方法。使用三维查表插值法做设备色彩空间向CIEL*a*b*设备无关色彩空间转换即正向变换时,一般是将RGB(或CMY)设备色彩空间分割为若干规则排列的小立方体子空间,每个小立方体8个顶点的RGB值和测量得到的对应的CIEL*a*b*值构成查找表的数据项,因此,当任意给定一个RGB颜色点时,能够容易的在查找表中找到该颜色点所在的小立方体子空间,及其8个顶点的RGB值和对应的CIEL*a*b*值,继而通过三维插值得到该RGB颜色点对应的CIEL*a*b*值,该方法转换速度快,且精度可以随着查找表中数据项的增加而提高,从而满足实用需求。但使用三维查表插值法做CIEL*a*b*设备无关色彩空间向设备色彩空间转换即反向变换时,当任意给定一个CIEL*a*b*颜色点时,由于这两类色彩空间的非线性关系,查找表中的CIEL*a*b*测量数据是无规则排列的,因此很难在查找表中找出包含该CIEL*a*b*颜色点的不规则六面体,及其8个顶点的CIEL*a*b*值和对应的RGB值,并进行插值计算得出设备色彩空间的RGB颜色值。针对三维查表插值法的不足,相继提出了一些改进算法,如预测查找方法、立方体迭代方法和反向查找方法等。预测查找方法是将CIEL*a*b*色彩空间用L*C*h*色彩空间来代替,通过对明度、彩度和色调的调整,迭代得出满足色差要求的RGB值,但通过迭代计算出指定精度的RGB值需要大量时间;立方体细分迭代方法是通过迭代细分待转换CIEL*a*b*值对应RGB值所在的立方体集的方法得到在给定色差允值范围之内的RGB值,通过迭代得到指定精度的RGB值,转换速度慢;反向查找方法是通过寻找待转换CIEL*a*b*值对应的8个相邻点,利用一种插值和迭代相结合的方法找到较准确的RGB值,但该方法在寻找CIEL*a*b*对应的8个相邻点和通过迭代计算指定精度的RGB值时需要大量的时间。目前针对查找表插值做反向变换的改进方法转换精度可以满足需求,但由于在进行变换时大都采用了迭代算法而使得色彩空间转换的速度变慢,难以满足实际应用的需要。
发明内容
本发明的目的是提供一种基于查表法的反向色彩空间转换方法,解决了现有采用三维查表插值法反向变换时的查表不易、转换速度慢的问题,在进行转换时,既能保持转换精度又能使得转换速度满足实用要求。
本发明所采用的技术方案是,一种基于查表法的反向色彩空间转换方法,具体按照以下步骤实施:
步骤1:建立设备色彩空间到CIEL*a*b*色彩空间的查找表;
步骤2:根据步骤1得到的查找表中的CIEL*a*b*数据,分割CIEL*a*b*空间,建立自定义子空间查找表;
步骤3:根据步骤2得到的自定义子空间查找表,计算CIEL*a*b*子空间的几何中心;
步骤4:根据步骤1得到的查找表及步骤3得到的CIEL*a*b*子空间的几何中心,获取CIEL*a*b*子空间的建模数据;
步骤5:根据步骤4得到的CIEL*a*b*子空间的建模数据,建立CIEL*a*b*子空间的多项式转换模型;
步骤6:根据步骤2得到的自定义子空间查找表,确定给定待转换点所在的CIEL*a*b*子空间;
步骤7:根据步骤5得到的CIEL*a*b*子空间的多项式转换模型及步骤6得到的给定待转换点所在的CIEL*a*b*子空间,计算给定待转换点对应的设备空间的颜色值,实现CIEL*a*b*色彩空间到设备色彩空间的转换。
本发明的特点还在于,
其中步骤1建立设备色彩空间到CIEL*a*b*色彩空间的查找表,具体按照以下步骤实施:
若设备色彩空间为RGB空间,将RGB空间的每个坐标轴M级分割,使用分割点的R、G、B值在Photoshop中制作RGB模式的色靶文件,色靶中包含M*M*M个色块,显示色靶文件后测量每个色块的CIEL*a*b*值,获得RGB设备色彩空间到CIEL*a*b*色彩空间的查找表数据,建立查找表;
若设备色彩空间为CMYK空间,则在Photoshop中将RGB模式色靶文件转换成CMYK模式,再读取每个色块的C、M、Y、K值,打印色靶文件后测量每个色块的CIEL*a*b*值,获得CMYK设备色彩空间到CIEL*a*b*色彩空间的查找表数据,建立查找表。
其中步骤2分割CIEL*a*b*空间,建立自定义子空间查找表,具体按照以下步骤实施:
根据步骤1得到的查找表中的CIEL*a*b*数据,分别找出最小的L值和最大的L值,并将其整量化,得到Lmin和Lmax,将CIEL*a*b*空间中L轴从Lmin到Lmax的部分以nL个分割点进行分割,即将CIEL*a*b*空间以L轴的Lmin到Lmax的部分为对象顺序分割为nL-1层;
根据步骤1得到的查找表中的CIEL*a*b*数据,找到a和b共同的最小和最大值,并将其整量化,得到ABmin和ABmax,分别将CIEL*a*b*空间中的a轴和b轴从ABmin到ABmax的部分以nAB个分割点进行分割,即将CIEL*a*b*空间以a轴的ABmin到ABmax的部分为对象顺序分割为nAB-1层和以b轴的ABmin到ABmax的部分为对象顺序分割为nAB-1层;
从而将CIEL*a*b*空间规则分割成(nL-1)(nAB-1)(nAB-1)个矩形体子空间,L轴、a轴、b轴上的分割点构成了自定义子空间查找表。
其中步骤3计算CIEL*a*b*子空间的几何中心,具体按照以下步骤实施:
根据步骤2得到的自定义子空间查找表,得到每个CIEL*a*b*矩形体子空间的位置,即子空间的min Li、max Li、min ai、max ai、min bi和max bi,计算每个子空间几何中心点,公式如下:
CLi=(min Li+max Li)/2
Cai=(min ai+max ai)/2,
Cbi=(min bi+max bi)/2
其中CLi、Cai、Cbi为第i个CIEL*a*b*子空间几何中心点的CIEL*a*b*值,0<i≤(nL-1)(nAB-1)(nAB-1)。
其中步骤4获取CIEL*a*b*子空间的建模数据,具体按照以下步骤实施:
计算步骤3得到的每个CIEL*a*b*子空间的几何中心点的CLi、Cai、Cbi值与步骤1得到的查找表中的Lj、aj、bj值之间的色差,色差的计算公式如下:
其中0≤j≤M*M*M,0<i≤(nL-1)(nAB-1)(nAB-1);
取出最小的N个色差值在步骤1查找表中所对应的CIEL*a*b*值和设备空间颜色值作为各CIEL*a*b*子空间的建模数据。
其中步骤5建立CIEL*a*b*子空间的多项式转换模型,具体按照以下步骤实施:
从CIEL*a*b*到设备色彩空间转换的多项式方程,如下式所示,
其中X、Y、Z、…为设备色彩空间颜色值的各个分量;r、s、t=1、2、3;v1=L,v2=a,v3=b;
根据步骤4得到的每个CIEL*a*b*子空间所对应的N组CIEL*a*b*值和设备空间颜色值的建模数据,利用多项式回归的方法计算各CIEL*a*b*子空间的多项式转换方程的系数,即p0~pj、q0~qj、m0~mj、…的值。
其中步骤6确定给定待转换点所在的CIEL*a*b*子空间,具体按照以下步骤实施:对于任意给定的待转换点的CIEL*a*b*值,根据步骤2得到的自定义子空间查找表,查找包含待转换点的矩形体子空间的位置。
其中步骤7计算给定待转换点对应的设备空间的颜色值,具体按照以下步骤实施:根据步骤6确定的给定待转换点所在的CIEL*a*b*子空间以及步骤5建立的该子空间的多项式转换模型,计算待转换点对应的设备空间颜色值的各个分量,实现CIEL*a*b*色彩空间到设备色彩空间的转换。
本发明通过设备色彩空间到CIEL*a*b*色彩空间的查找表,对CIEL*a*b*色彩空间进行分割,构建自定义子空间查找表,方便查找待转换的CIEL*a*b*颜色值所在的CIEL*a*b*子空间的位置;在建模阶段根据精度要求为每个CIEL*a*b*子空间建立一一对应的多项式转换模型,可保证色彩空间反向转换的精度;在进行反向转换时,通过自定义子空间查找表,可快速定位待转换点所在的CIEL*a*b*子空间,以及与其对应的转换模型,避免了在转换过程中使用迭代方法,保证了转换速度,使得该方法在进行CIEL*a*b*色彩空间到设备色彩空间转换时的速度和精度可以同时满足实际应用的需要。
附图说明
图1是本发明方法采用的L轴、a轴和b轴的分割示意图;
图2是本发明方法的自定义子空间查找表;
图3是本发明方法的第i个CIEL*a*b*子空间示意图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明基于查表法的反向色彩空间转换方法,具体按照以下步骤实施:
步骤1:建立设备色彩空间到CIEL*a*b*色彩空间的查找表,
调节彩色打印或显示设备,保证设备正常使用,将RGB空间的每个坐标轴M级分割,使用分割点的R、G、B值在Photoshop中制作RGB模式的色靶文件,色靶中包含了M*M*M个色块;如果设备空间为CMYK,则在Photoshop中将RGB模式色靶文件转换成CMYK模式,再读取每个色块的CMYK值,打印或显示色靶文件后测量每个色块的CIEL*a*b*值,获得设备色彩空间到CIEL*a*b*色彩空间的查找表数据,建立查找表。
步骤2:分割CIEL*a*b*空间,建立自定义子空间查找表,
根据步骤1查找表中的CIEL*a*b*数据,分别找出最小和最大的L值,并将其整量化,得到Lmin和Lmax,将CIEL*a*b*空间中L轴从Lmin到Lmax的部分以nL个分割点进行分割(可以均匀分割,也可以非均匀分割),即将CIEL*a*b*空间以L轴的Lmin到Lmax的部分为对象顺序分割为nL-1层,如图1中的L轴所示。
根据步骤1查找表中的CIEL*a*b*数据,找到a和b共同的最小和最大值,并将其整量化,得到ABmin和ABmax,分别将CIEL*a*b*空间中的a轴和b轴从ABmin到ABmax的部分以nAB个分割点进行分割(可以均匀分割,也可以非均匀分割),即将CIEL*a*b*空间以a轴的ABmin到ABmax的部分为对象顺序分割为nAB-1层和以b轴的ABmin到ABmax的部分为对象顺序分割为nAB-1层,如图1中的a轴和b轴所示。
由此将CIEL*a*b*空间规则分割成(nL-1)(nAB-1)(nAB-1)个矩形体子空间,L轴、a轴、b轴上的分割点构成了自定义子空间查找表,如图2所示。
步骤3:计算CIEL*a*b*子空间的几何中心
根据步骤2得到的自定义子空间查找表,可得每个CIEL*a*b*矩形体子空间的位置,即子空间的min Li、max Li、min ai、max ai、min bi和max bi,如图3所示,计算每个子空间几何中心点,公式如下:
CLi=(min Li+max Li)/2
Cai=(min ai+max ai)/2 (1)
Cbi=(min bi+max bi)/2
其中CLi、Cai、Cbi为第i个CIEL*a*b*子空间几何中心点的CIEL*a*b*值,0<i≤(nL-1)(nAB-1)(nAB-1)。
步骤4:获取CIEL*a*b*子空间的建模数据
将每个CIEL*a*b*子空间的几何中心点CLi、Cai、Cbi值与步骤1查找表中的Lj、aj、bj值计算色差,色差的计算公式如下所示:
其中0≤j≤M*M*M,0<i≤(nL-1)(nAB-1)(nAB-1)。
取出最小的N个色差值在步骤1查找表中所对应的CIEL*a*b*值和设备空间颜色值作为各CIEL*a*b*子空间的建模数据。
步骤5:建立CIEL*a*b*子空间的多项式转换模型
从CIEL*a*b*到设备色彩空间转换的多项式方程,如公式(3)所示,
其中X、Y、Z、…为设备色彩空间颜色值的各个分量;r、s、t=1、2、3;v1=L,v2=a,v3=b。
根据步骤4得出的每个CIEL*a*b*子空间所对应的N组(N大于多项式的系数个数)CIEL*a*b*值和设备空间颜色值的建模数据,利用多项式回归的方法计算各CIEL*a*b*子空间的多项式转换方程的系数,即p0~pj、q0~qj、m0~mj、…的值。
步骤6:确定给定待转换点所在的CIEL*a*b*子空间
对于任意给定的待转换点的CIEL*a*b*值:Lk、ak、bk,根据步骤2得到的自定义子空间查找表,查找包含待转换点Lk、ak、bk的矩形体子空间的位置,即确定min Lk、max Lk、min ak、max ak、min bk和max bk的值。
步骤7:计算给定待转换点对应的设备空间的颜色值
根据步骤6确定的给定待转换点所在的CIEL*a*b*子空间以及步骤5建立的该子空间的多项式转换模型,计算待转换点Lk、ak、bk对应的设备空间颜色值的各个分量,实现CIEL*a*b*色彩空间到设备色彩空间的转换。
实施例
下面以CIEL*a*b*色彩空间到CMYK设备色彩空间转换为例,具体说明本发明基于查表法的反向色彩空间转换方法。
1.建立CIEL*a*b*色彩空间到设备色彩空间的转换模型
步骤1:建立CMYK色彩空间到CIEL*a*b*色彩空间的查找表
调节彩色打印设备,保证设备正常使用,将RGB空间的每个坐标轴以M=11级均匀分割,分割点为0、25、50、75、100、125、150、175、200、225和255;使用分割点的RGB值在Photoshop中制作RGB模式的色靶文件,色靶中包含了11×11×11=1331个色块;在Photoshop自定义CMYK的分色选项中选定分色类型为GCR,黑版产生为中,黑色油墨限制为100%,油墨总量限制为300%,底层颜色添加量为0%后,将RGB模式的色靶文件转换成CMYK模式的色靶文件,并读取每个色块的CMYK值,打印该色靶文件后测量每个色块的CIEL*a*b*值,获得CMYK色彩空间到CIEL*a*b*色彩空间的查找表数据,建立查找表。
步骤2:分割CIEL*a*b*空间,建立自定义子空间查找表
根据步骤1查找表中的CIEL*a*b*数据,分别找出最小和最大的L值,并将其整量化,本实施例中Lmin=6、Lmax=94,将CIEL*a*b*空间中L轴从6到94的部分按照nL=6个分割点6、26、42、58、74、94进行分割,即将CIEL*a*b*空间以L轴从6到94的部分为对象顺序分割为nL-1=5层。
根据步骤1查找表中的CIEL*a*b*数据,找到a和b共同的最小和最大值,并将其整量化,本实施例中ABmin=-50、ABmax=80,分别将CIEL*a*b*空间中的a轴和b轴从-50到80的部分按照nAB=7个分割点-50、-30、-10、10、30、50、80进行分割,即将CIEL*a*b*空间以a轴从-50到80的部分为对象顺序分割为nAB-1=6层和以b轴从-50到80的部分为对象顺序分割为nAB-1=6层。
由此将CIEL*a*b*空间规则分割为(nL-1)(nAB-1)(nAB-1)=5×6×6=180个子空间,L轴、a轴、b轴上的分割点构成自定义子空间查找表。
步骤3:计算CIEL*a*b*子空间的几何中心
根据步骤2得到的自定义子空间查找表,可知每个CIEL*a*b*子空间的位置。如本实施例中第i个子空间的min Li=26、max Li=42、min ai=10、max ai=30、min bi=-30和max bi=-10,根据公式(1),可计算得到该子空间几何中心点的CIEL*a*b*值为CLi=34,Cai=20,Cbi=-20。
步骤4:获取CIEL*a*b*子空间的建模数据
将CIEL*a*b*子空间的几何中心点CLi、Cai、Cbi值(本实施例中0<i<180)与步骤1查找表中的Lj、aj、bj值(本实施例中0<j<1331)按照公式(2)计算色差,取出每个CIEL*a*b*子空间最小的N=120个色差值在步骤1查找表中所对应的CIEL*a*b*值和CMYK值。
步骤5:建立CIEL*a*b*子空间的多项式转换模型
本实施例计算从CIEL*a*b*到CMYK色彩空间转换的多项式方程,如公式(4)所示,根据步骤4得出的每个CIEL*a*b*子空间所对应的N=120组CIEL*a*b*和CMYK值的建模数据,利用多项式回归的方法计算出各CIEL*a*b*子空间的多项式系数,即每个CIEL*a*b*子空间的p0~p10、q0~q10、m0~m10、n0~n10的值。
C=p0+p1L+p2a+p3b+p4La+p5Lb+p6ab+p7L2+p8a2+p9b2+p10Lab
M=q0+q1L+q2a+q3b+q4La+q5Lb+q6ab+q7L2+q8a2+q9b2+q10Lab (4)
Y=m0+m1L+m2a+m3b+m4La+m5Lb+m6ab+m7L2+m8a2+m9b2+m10Lab
K=n0+n1L+n2a+n3b+n4La+n5Lb+n6ab+n7L2+n8a2+n9b2+n10Lab
2.实现CIEL*a*b*色彩空间到CMYK色彩空间的转换
步骤6:确定给定待转换点所在的CIEL*a*b*子空间
给定待转换点的CIEL*a*b*值,本实施例中Lk=56.5、ak=-23.9、bk=41.6;根据自定义子空间查找表,找到包含待转换点Lk、ak、bk的矩形体子空间的位置,即minL=58、maxL=74、mina=-30、maxa=-10、minb=30和maxb=50。
步骤7:计算给定待转换点对应的设备空间CMYK的值
根据步骤6确定的给定待转换点所在的CIEL*a*b*子空间以及步骤5建立的该子空间的多项式转换模型,计算待转换点Lk、ak、bk对应的设备空间CMYK的各个分量值为C=50、M=13、Y=86、K=8,完成CIEL*a*b*色彩空间到CMYK色彩空间的转换。
Claims (8)
1.一种基于查表法的反向色彩空间转换方法,其特征在于,具体按照以下步骤实施:
步骤1:建立设备色彩空间到CIEL*a*b*色彩空间的查找表;
步骤2:根据步骤1得到的查找表中的CIEL*a*b*数据,分割CIEL*a*b*空间,建立自定义子空间查找表;
步骤3:根据步骤2得到的自定义子空间查找表,计算CIEL*a*b*子空间的几何中心;
步骤4:根据步骤1得到的查找表及步骤3得到的CIEL*a*b*子空间的几何中心,获取CIEL*a*b*子空间的建模数据;
步骤5:根据步骤4得到的CIEL*a*b*子空间的建模数据,建立CIEL*a*b*子空间的多项式转换模型;
步骤6:根据步骤2得到的自定义子空间查找表,确定给定待转换点所在的CIEL*a*b*子空间;
步骤7:根据步骤5得到的CIEL*a*b*子空间的多项式转换模型及步骤6得到的给定待转换点所在的CIEL*a*b*子空间,计算给定待转换点对应的设备空间的颜色值,实现CIEL*a*b*色彩空间到设备色彩空间的转换。
2.根据权利要求1所述的基于查表法的反向色彩空间转换方法,其特征在于,所述的步骤1建立设备色彩空间到CIEL*a*b*色彩空间的查找表,具体按照以下步骤实施:
若设备色彩空间为RGB空间,将RGB空间的每个坐标轴M级分割,使用分割点的R、G、B值在Photoshop中制作RGB模式的色靶文件,色靶中包含M*M*M个色块,显示色靶文件后测量每个色块的CIEL*a*b*值,获得RGB设备色彩空间到CIEL*a*b*色彩空间的查找表数据,建立查找表;
若设备色彩空间为CMYK空间,则在Photoshop中将RGB模式色靶文件转换成CMYK模式,再读取每个色块的C、M、Y、K值,打印色靶文件后测量每个色块的CIEL*a*b*值,获得CMYK设备色彩空间到CIEL*a*b*色彩空间的查找表数据,建立查找表。
3.根据权利要求1所述的基于查表法的反向色彩空间转换方法,其特征在于,所述的步骤2分割CIEL*a*b*空间,建立自定义子空间查找表,具体按照以下步骤实施:
根据步骤1得到的查找表中的CIEL*a*b*数据,分别找出最小的L值和最大的L值,并将其整量化,得到Lmin和Lmax,将CIEL*a*b*空间中L轴从Lmin到Lmax的部分以nL个分割点进行分割,即将CIEL*a*b*空间以L轴的Lmin到Lmax的部分为对象顺序分割为nL-1层;
根据步骤1得到的查找表中的CIEL*a*b*数据,找到a和b共同的最小和最大值,并将其整量化,得到ABmin和ABmax,分别将CIEL*a*b*空间中的a轴和b轴从ABmin到ABmax的部分以nAB个分割点进行分割,即将CIEL*a*b*空间以a轴的ABmin到ABmax的部分为对象顺序分割为nAB-1层和以b轴的ABmin到ABmax的部分为对象顺序分割为nAB-1层;
从而将CIEL*a*b*空间规则分割成(nL-1)(nAB-1)(nAB-1)个矩形体子空间,L轴、a轴、b轴上的分割点构成了自定义子空间查找表。
4.根据权利要求1所述的基于查表法的反向色彩空间转换方法,其特征在于,所述的步骤3计算CIEL*a*b*子空间的几何中心,具体按照以下步骤实施:
根据步骤2得到的自定义子空间查找表,得到每个CIEL*a*b*矩形体子空间的位置,即子空间的min Li、max Li、min ai、max ai、min bi和max bi,计算每个子空间几何中心点,公式如下:
CLi=(min Li+max Li)/2
Cai=(min ai+max ai)/2,
Cbi=(min bi+max bi)/2
其中CLi、Cai、Cbi为第i个CIEL*a*b*子空间几何中心点的CIEL*a*b*值,0<i≤(nL-1)(nAB-1)(nAB-1)。
5.根据权利要求1所述的基于查表法的反向色彩空间转换方法,其特征在于,所述的步骤4获取CIEL*a*b*子空间的建模数据,具体按照以下步骤实施:
计算步骤3得到的每个CIEL*a*b*子空间的几何中心点的CLi、Cai、Cbi值与步骤1得到的查找表中的Lj、aj、bj值之间的色差,色差的计算公式如下:
其中0≤j≤M*M*M,0<i≤(nL-1)(nAB-1)(nAB-1);
取出最小的N个色差值在步骤1查找表中所对应的CIEL*a*b*值和设备空间颜色值作为各CIEL*a*b*子空间的建模数据。
6.根据权利要求1所述的基于查表法的反向色彩空间转换方法,其特征在于,所述的步骤5建立CIEL*a*b*子空间的多项式转换模型,具体按照以下步骤实施:
从CIEL*a*b*到设备色彩空间转换的多项式方程,如下式所示,
其中X、Y、Z、…为设备色彩空间颜色值的各个分量;r、s、t=1、2、3;v1=L,v2=a,v3=b;
根据步骤4得到的每个CIEL*a*b*子空间所对应的N组CIEL*a*b*值和设备空间颜色值的建模数据,利用多项式回归的方法计算各CIEL*a*b*子空间的多项式转换方程的系数,即p0~pj、q0~qj、m0~mj、…的值。
7.根据权利要求1所述的基于查表法的反向色彩空间转换方法,其特征在于,所述的步骤6确定给定待转换点所在的CIEL*a*b*子空间,具体按照以下步骤实施:对于任意给定的待转换点的CIEL*a*b*值,根据步骤2得到的自定义子空间查找表,查找包含待转换点的矩形体子空间的位置。
8.根据权利要求1所述的基于查表法的反向色彩空间转换方法,其特征在于,所述的步骤7计算给定待转换点对应的设备空间的颜色值,具体按照以下步骤实施:根据步骤6确定的给定待转换点所在的CIEL*a*b*子空间以及步骤5建立的该子空间的多项式转换模型,计算待转换点对应的设备空间颜色值的各个分量,实现CIEL*a*b*色彩空间到设备色彩空间的转换。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100285623A CN102075667B (zh) | 2011-01-26 | 2011-01-26 | 一种基于查表法的反向色彩空间转换方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100285623A CN102075667B (zh) | 2011-01-26 | 2011-01-26 | 一种基于查表法的反向色彩空间转换方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102075667A true CN102075667A (zh) | 2011-05-25 |
CN102075667B CN102075667B (zh) | 2012-07-04 |
Family
ID=44033998
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011100285623A Expired - Fee Related CN102075667B (zh) | 2011-01-26 | 2011-01-26 | 一种基于查表法的反向色彩空间转换方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102075667B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103106670A (zh) * | 2013-01-16 | 2013-05-15 | 西安理工大学 | 一种计算机输入设备的色彩模型自动建立及转换方法 |
CN103714110A (zh) * | 2012-10-04 | 2014-04-09 | 京瓷办公信息系统株式会社 | 用于创建运行时颜色转换查找表的方法和系统 |
CN106448613A (zh) * | 2016-10-08 | 2017-02-22 | 南京巨鲨医疗科技有限公司 | 一种医用显示器色彩空间均匀化处理方法 |
CN107197119A (zh) * | 2017-05-27 | 2017-09-22 | 华南理工大学 | CIEL*a*b*到CMYK色彩空间的转换方法 |
CN109859702A (zh) * | 2018-12-28 | 2019-06-07 | 南京奥视威电子科技股份有限公司 | 一种3d查找表生成方法、显示器校色方法、显示器校色系统 |
CN110085192A (zh) * | 2019-04-01 | 2019-08-02 | 南京巨鲨显示科技有限公司 | 一种医用显示器颜色曲线的校准方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0581590A2 (en) * | 1992-07-31 | 1994-02-02 | Canon Kabushiki Kaisha | Color processing method |
EP1533994A2 (en) * | 2003-11-20 | 2005-05-25 | Samsung Electronics Co., Ltd. | Apparatus and method for controlling the colors of a color image |
CN1867963A (zh) * | 2003-10-21 | 2006-11-22 | 克雷沃耶提公司 | 色域转换系统和方法 |
CN101621703A (zh) * | 2009-08-06 | 2010-01-06 | 北京中企开源信息技术有限公司 | 一种图像序列色彩空间转换的方法 |
-
2011
- 2011-01-26 CN CN2011100285623A patent/CN102075667B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0581590A2 (en) * | 1992-07-31 | 1994-02-02 | Canon Kabushiki Kaisha | Color processing method |
CN1867963A (zh) * | 2003-10-21 | 2006-11-22 | 克雷沃耶提公司 | 色域转换系统和方法 |
EP1533994A2 (en) * | 2003-11-20 | 2005-05-25 | Samsung Electronics Co., Ltd. | Apparatus and method for controlling the colors of a color image |
CN101621703A (zh) * | 2009-08-06 | 2010-01-06 | 北京中企开源信息技术有限公司 | 一种图像序列色彩空间转换的方法 |
Non-Patent Citations (1)
Title |
---|
《计算机时代》 20051231 胡涛,等 基于局部色彩空间转换的灰平衡校准方法 , 第12期 2 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103714110A (zh) * | 2012-10-04 | 2014-04-09 | 京瓷办公信息系统株式会社 | 用于创建运行时颜色转换查找表的方法和系统 |
CN103714110B (zh) * | 2012-10-04 | 2017-05-17 | 京瓷办公信息系统株式会社 | 用于创建运行时颜色转换查找表的方法和系统 |
CN103106670A (zh) * | 2013-01-16 | 2013-05-15 | 西安理工大学 | 一种计算机输入设备的色彩模型自动建立及转换方法 |
CN103106670B (zh) * | 2013-01-16 | 2015-06-24 | 西安理工大学 | 一种计算机输入设备的色彩模型自动建立及转换方法 |
CN106448613A (zh) * | 2016-10-08 | 2017-02-22 | 南京巨鲨医疗科技有限公司 | 一种医用显示器色彩空间均匀化处理方法 |
CN107197119A (zh) * | 2017-05-27 | 2017-09-22 | 华南理工大学 | CIEL*a*b*到CMYK色彩空间的转换方法 |
CN107197119B (zh) * | 2017-05-27 | 2019-03-05 | 华南理工大学 | CIEL*a*b*到CMYK色彩空间的转换方法 |
CN109859702A (zh) * | 2018-12-28 | 2019-06-07 | 南京奥视威电子科技股份有限公司 | 一种3d查找表生成方法、显示器校色方法、显示器校色系统 |
CN110085192A (zh) * | 2019-04-01 | 2019-08-02 | 南京巨鲨显示科技有限公司 | 一种医用显示器颜色曲线的校准方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102075667B (zh) | 2012-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102075667B (zh) | 一种基于查表法的反向色彩空间转换方法 | |
CN100589521C (zh) | 高保真打印输出L*a*b*图像的方法 | |
CN103144448B (zh) | 一种实现喷墨打印色域最大化的墨量限制方法 | |
EP2717556A2 (en) | Calculating the spectral characteristics of the color resulting from overlaying colorants | |
CN102945556B (zh) | 基于胞元划分纽介堡方程的七色分色算法 | |
US8325394B2 (en) | Hierarchical scanner characterization | |
CN106464775A (zh) | 颜色模型 | |
JP2012049823A (ja) | ガマット形状予測方法、ガマット形状予測プログラム、当該ガマット形状予測方法により予測されたガマットを用いて作成された色変換テーブルを搭載した印刷装置 | |
CN102180003B (zh) | 一种打印机光谱特征化模型的修正方法 | |
JP2021017034A (ja) | 情報処理装置、方法、およびプログラム | |
US20050062992A1 (en) | Image processing apparatus, image processing method, image processing program, and storage medium | |
US8482797B2 (en) | Method and apparatus for calculating metameric colors | |
CN102238297B (zh) | 一种生成色彩对照文件的方法及系统 | |
CN103823943A (zh) | 一种用于色域最大化的多色打印印刷系统拆分式建模方法 | |
CN105667069A (zh) | 一种基于比尔朗伯定律的光谱配色方法 | |
JP2005094160A (ja) | 測色データ補正方法、変換テーブル作成方法及びプログラム | |
JP2003125210A (ja) | 色処理方法、記憶媒体、色処理装置、色変換装置、および画像形成装置 | |
JP4779728B2 (ja) | 色推定モデルの作成方法 | |
CN101902542B (zh) | 颜色变换表调整装置 | |
US7660013B2 (en) | Method of converting a color recipe | |
CN108437440B (zh) | 基于颜色管理的3d打印色彩呈现方法、装置及系统 | |
CN103419511A (zh) | 颜色处理装置和方法 | |
CN104553302A (zh) | 一种灰平衡校准方法及装置 | |
CN106671594A (zh) | 一种确定印刷三原色密度的方法 | |
JP2005051752A (ja) | 視覚的に均一なプリンタ特性サンプル適応生成方法 |
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: 20120704 Termination date: 20150126 |
|
EXPY | Termination of patent right or utility model |