CN102075667B - 一种基于查表法的反向色彩空间转换方法 - Google Patents

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值，并将其整量化，得到L_min和L_max，将CIEL^*a^*b空间中L轴从L_min到L_max的部分以nL个分割点进行分割，即将CIEL^*a^*b空间以L轴的L_min到L_max的部分为对象顺序分割为nL-1层；

根据步骤1得到的查找表中的CIEL^*a^*b数据，找到a和b共同的最小和最大值，并将其整量化，得到AB_min和AB_max，分别将CIEL^*a^*b空间中的a轴和b轴从AB_min到AB_max的部分以nAB个分割点进行分割，即将CIEL^*a^*b空间以a轴的AB_min到AB_max的部分为对象顺序分割为nAB-1层和以b轴的AB_min到AB_max的部分为对象顺序分割为nAB-1层；

从而将CIEL^*a^*b空间规则分割成(nL-1)(nAB-1)(nAB-1)个矩形体子空间，L轴、a轴、b轴上的分割点构成了自定义子空间查找表。

其中步骤3计算CIEL^*a^*b子空间的几何中心，具体按照以下步骤实施：

根据步骤2得到的自定义子空间查找表，得到每个CIEL^*a^*b矩形体子空间的位置，即子空间的min L_i、max L_i、min a_i、max a_i、min b_i和max b_i，计算每个子空间几何中心点，公式如下：

CL_i＝(min L_i+max L_i)/2

Ca_i＝(min a_i+max q_i)/2，

Cb_i＝(min b_i+max b_i)/2

其中CL_i、Ca_i、Cb_i为第i个CIEL^*a^*b子空间几何中心点的CIEL^*a^*b值，0＜i≤(nL-1)(nAB-1)(nAB-1)。

其中步骤4获取CIEL^*a^*b子空间的建模数据，具体按照以下步骤实施：

计算步骤3得到的每个CIEL^*a^*b子空间的几何中心点的CL_i、Ca_i、Cb_i值与步骤1得到的查找表中的L_j、a_j、b_j值之间的色差，色差的计算公式如下：

$\mathrm{\Δ}{E}_{{\mathrm{ab}}_{\mathrm{ij}}}=\sqrt{{({\mathrm{CL}}_i-L_j)}^2+{({\mathrm{Ca}}_i-a_j)}^2+{({\mathrm{Cb}}_i-b_j)}^2},$

其中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=p_0+\underset{j=1}{\overset{3}{\mathrm{\Σ}}}{p}_j{v}_r+\underset{j=4}{\overset{9}{\mathrm{\Σ}}}{p}_j{v}_r{v}_s+\underset{j=10}{\overset{36}{\mathrm{\Σ}}}{p}_j{v}_r{v}_s{v}_t$

$Y=q_0+\underset{j=1}{\overset{3}{\mathrm{\Σ}}}{q}_j{v}_r+\underset{j=4}{\overset{9}{\mathrm{\Σ}}}{q}_j{v}_r{v}_s+\underset{j=10}{\overset{36}{\mathrm{\Σ}}}{q}_j{v}_r{v}_s{v}_t,$

其中X、Y、Z、…为设备色彩空间颜色值的各个分量；r、s、t＝1、2、3；v₁＝L，v₂＝a，v₃＝b；

根据步骤4得到的每个CIEL^*a^*b子空间所对应的N组CIEL^*a^*b值和设备空间颜色值的建模数据，利用多项式回归的方法计算各CIEL^*a^*b子空间的多项式转换方程的系数，即p₀～p_j、q₀～q_j、m₀～m_j、…的值。

其中步骤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值，并将其整量化，得到L_min和L_max，将CIEL^*a^*b空间中L轴从L_min到L_max的部分以nL个分割点进行分割(可以均匀分割，也可以非均匀分割)，即将CIEL^*a^*b空间以L轴的L_min到L_max的部分为对象顺序分割为nL-1层，如图1中的L轴所示。

根据步骤1查找表中的CIEL^*a^*b数据，找到a和b共同的最小和最大值，并将其整量化，得到AB_min和AB_max，分别将CIEL^*a^*b空间中的a轴和b轴从AB_min到AB_max的部分以nAB个分割点进行分割(可以均匀分割，也可以非均匀分割)，即将CIEL^*a^*b空间以a轴的AB_min到AB_max的部分为对象顺序分割为nAB-1层和以b轴的AB_min到AB_max的部分为对象顺序分割为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 L_i、max L_i、min a_i、max a_i、min b_i和max b_i，如图3所示，计算每个子空间几何中心点，公式如下：

CL_i＝(min L_i+max L_i)/2

Ca_i＝(min a_i+max a_i)/2                   (1)

Cb_i＝(min b_i+max b_i)/2

其中CL_i、Ca_i、Cb_i为第i个CIEL^*a^*b子空间几何中心点的CIEL^*a^*b值，0＜i≤(nL-1)(nAB-1)(nAB-1)。

步骤4：获取CIEL^*a^*b子空间的建模数据

将每个CIEL^*a^*b子空间的几何中心点CL_i、Ca_i、Cb_i值与步骤1查找表中的L_j、a_j、b_j值计算色差，色差的计算公式如下所示：

$\mathrm{\Δ}{E}_{{\mathrm{ab}}_{\mathrm{ij}}}=\sqrt{{({\mathrm{CL}}_i-L_j)}^2+{({\mathrm{Ca}}_i-a_j)}^2+{({\mathrm{Cb}}_i-b_j)}^2}---\left(2\right)$

其中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=p_0+\underset{j=1}{\overset{3}{\mathrm{\Σ}}}{p}_j{v}_r+\underset{j=4}{\overset{9}{\mathrm{\Σ}}}{p}_j{v}_r{v}_s+\underset{j=10}{\overset{36}{\mathrm{\Σ}}}{p}_j{v}_r{v}_s{v}_t$

$Y=q_0+\underset{j=1}{\overset{3}{\mathrm{\Σ}}}{q}_j{v}_r+\underset{j=4}{\overset{9}{\mathrm{\Σ}}}{q}_j{v}_r{v}_s+\underset{j=10}{\overset{36}{\mathrm{\Σ}}}{q}_j{v}_r{v}_s{v}_t---\left(3\right)$

其中X、Y、Z、…为设备色彩空间颜色值的各个分量；r、s、t＝1、2、3；v₁＝L，v₂＝a，v₃＝b。

根据步骤4得出的每个CIEL^*a^*b子空间所对应的N组(N大于多项式的系数个数)CIEL^*a^*b值和设备空间颜色值的建模数据，利用多项式回归的方法计算各CIEL^*a^*b子空间的多项式转换方程的系数，即p₀～p_j、q₀～q_j、m₀～m_j、…的值。

步骤6：确定给定待转换点所在的CIEL^*a^*b子空间

对于任意给定的待转换点的CIEL^*a^*b值：L_k、a_k、b_k，根据步骤2得到的自定义子空间查找表，查找包含待转换点L_k、a_k、b_k的矩形体子空间的位置，即确定min L_k、max L_k、min a_k、max a_k、min b_k和max b_k的值。

步骤7：计算给定待转换点对应的设备空间的颜色值

根据步骤6确定的给定待转换点所在的CIEL^*a^*b子空间以及步骤5建立的该子空间的多项式转换模型，计算待转换点L_k、a_k、b_k对应的设备空间颜色值的各个分量，实现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值，并将其整量化，本实施例中L_min＝6、L_max＝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共同的最小和最大值，并将其整量化，本实施例中AB_min＝-50、AB_max＝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 L_i＝26、max L_i＝42、min a_i＝10、max a_i＝30、min b_i＝-30和max b_i＝-10，根据公式(1)，可计算得到该子空间几何中心点的CIEL^*a^*b值为CL_i＝34，Ca_i＝20，Cb_i＝-20。

步骤4：获取CIEL^*a^*b子空间的建模数据

将CIEL^*a^*b子空间的几何中心点CL_i、Ca_i、Cb_i值(本实施例中0＜i＜180)与步骤1查找表中的L_j、a_j、b_j值(本实施例中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子空间的p₀～P₁₀、q₀～q₁₀、m₀～m₁₀、n₀～n₁₀的值。

C＝p₀+p₁L+p₂a+p₃b+p₄La+p₅Lb+p₆ab+p₇L²+p₈a²+p₉b²+p₁₀Lab

M＝q₀+q₁L+q₂a+q₃b+q₄La+q₅Lb+q₆ab+q₇L²+q₈a²+q₉b²+q₁₀Lab    (4)

Y＝m₀+m₁L+m₂a+m₃b+m₄La+m₅Lb+m₆ab+m₇L²+m₈a²+m₉b²+m₁₀Lab

K＝n₀+n₁L+n₂a+n₃b+n₄La+n₅Lb+n₆ab+n₇L²+n₈a²+n₉b²+n₁₀Lab

2.实现CIEL^*a^*b色彩空间到CMYK色彩空间的转换

步骤6：确定给定待转换点所在的CIEL^*a^*b子空间

给定待转换点的CIEL^*a^*b值，本实施例中L_k＝56.5、a_k＝-23.9、b_k＝41.6；根据自定义子空间查找表，找到包含待转换点L_k、a_k、b_k的矩形体子空间的位置，即min L＝58、max L＝74、min a＝-30、max a＝-10、min b＝30和max b＝50。

步骤7：计算给定待转换点对应的设备空间CMYK的值

根据步骤6确定的给定待转换点所在的CIEL^*a^*b子空间以及步骤5建立的该子空间的多项式转换模型，计算待转换点L_k、a_k、b_k对应的设备空间CMYK的各个分量值为C＝50、M＝13、Y＝86、K＝8，完成CIEL^*a^*b色彩空间到CMYK色彩空间的转换。

Claims (4)

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色彩空间的查找表数据，建立查找表；

步骤2：根据步骤1得到的查找表中的CIEL^*a^*b数据，分割CIEL^*a^*b空间，建立自定义子空间查找表，具体按照以下步骤实施：

根据步骤1得到的查找表中的CIEL^*a^*b数据，分别找出最小的L值和最大的L值，并将其整量化，得到L_min和L_max，将CIEL^*a^*b空间中L轴从L_min到L_max的部分以nL个分割点进行分割，即将CIEL^*a^*b空间以L轴的L_min到L_max的部分为对象顺序分割为nL-1层；

根据步骤1得到的查找表中的CIEL^*a^*b数据，找到a和b共同的最小和最大值，并将其整量化，得到AB_min和AB_max，分别将CIEL^*a^*b空间中的a轴和b轴从AB_min到AB_max的部分以nAB个分割点进行分割，即将CIEL^*a^*b空间以a轴的AB_min到AB_max的部分为对象顺序分割为nAB-1层和以b轴的AB_min到AB_max的部分为对象顺序分割为nAB-1层；

从而将CIEL^*a^*b空间规则分割成(nL-1)(nAB-1)(nAB-1)个矩形体子空间，L轴、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所述的基于查表法的反向色彩空间转换方法，其特征在于，所述的步骤3计算CIEL^*a^*b子空间的几何中心，具体按照以下步骤实施：

根据步骤2得到的自定义子空间查找表，得到每个CIEL^*a^*b矩形体子空间的位置，即子空间的min L_i、max L_i、min a_i、max a_i、min b_i和max b_i，计算每个子空间几何中心点，公式如下：

CL_i＝(min L_i+max L_i)/2

Ca_i＝(min a_i+max a_i)/2，

Cb_i＝(min b_i+max b_i)/2

其中CL_i、Ca_i、Cb_i为第i个CIEL^*a^*b子空间几何中心点的CIEL^*a^*b值，0＜i≤(nL-1)(nAB-1)(nAB-1)。

3.根据权利要求1所述的基于查表法的反向色彩空间转换方法，其特征在于，所述的步骤6确定给定待转换点所在的CIEL^*a^*b子空间，具体按照以下步骤实施：对于任意给定的待转换点的CIEL^*a^*b值，根据步骤2得到的自定义子空间查找表，查找包含待转换点的矩形体子空间的位置。

4.根据权利要求1所述的基于查表法的反向色彩空间转换方法，其特征在于，所述的步骤7计算给定待转换点对应的设备空间的颜色值，具体按照以下步骤实施：根据步骤6确定的给定待转换点所在的CIEL^*a^*b子空间以及步骤5建立的该子空间的多项式转换模型，计算待转换点对应的设备空间颜色值的各个分量，实现CIEL^*a^*b色彩空间到设备色彩空间的转换。

Images