CN103823943A - 一种用于色域最大化的多色打印印刷系统拆分式建模方法 - Google Patents

Abstract

一种用于色域最大化的多色打印印刷系统拆分式建模方法，包括制作并打印各基色的单色梯尺，确定每种单色墨量的超限阈值，测量等比压缩后单色梯尺的光谱值，建立对应关系；对原始多色打印印刷系统拆分打印印刷系统子模型；进行墨色空间全局采样，获得所有打印印刷系统子模型下的采样样本，仿真获取预测光谱反射率信息值；针对所有子模型，使用凸包算法计算得到墨色空间全局采样色域体积；进行K个子模型的组合并使用凸包算法获得每种组合的色域体积，当最大色域覆盖率大于预设阈值时，锁定对应最大色域覆盖率的组合，实现原始打印印刷系统的拆分式建模。该方法具有显著降低系统拆分建模复杂性、避免墨量超限问题、提高输出色彩稳定性等技术优势。

Description

一种用于色域最大化的多色打印印刷系统拆分式建模方法

技术领域

本发明属于多色打印印刷色彩复制技术领域，具体涉及一种用于色域最大化的多色打印印刷系统拆分建模方法。

背景技术

多色打印印刷色彩复制技术是国内外影像复制领域的主流技术之一。此项技术以用户输入的影像色彩控制值为原始数据，依据预先建立的输入输出量之间的关联性模型，实现图像的色彩复制过程。相比传统的打印印刷色彩复制方式，多色打印印刷彩色复制技术具有色彩再现色域范围宽广、复制品色彩层次丰富以及细节清晰等技术优势。

在此项技术实施过程中，为获得具有上述特点的彩色复制品，在正式色彩复制之前需要构建多色打印印刷系统输入输出值之间的关联性模型，以保证影像复制过程的色彩准确性。然而，受多色打印印刷系统本身基色数较多的影响，目前对其系统进行整体建模仍存在采样压力巨大，建模效率较低等问题。此外，考虑到油墨、墨水等着色剂以及承印材料所固有的理化属性，过多种类及过大墨量的多色叠印也往往会造成叠印色彩稳定性较差以及叠印墨量过量等问题。

目前业界解决此类问题的主要方法是对多色打印印刷系统进行拆分建模，该方法的核心思想是将一个多色打印印刷系统拆分为若干个三色或四色子系统模型，通过子系统建模的方式实现多色打印印刷系统的整体建模。此种方法既可以提高系统建模效率，同时也一定程度上避免了墨量超限、叠印色彩稳定性差等问题。目前，在多色打印印刷系统拆分建模领域，现有方法大都是通过枚举、色域分布二维平面分析等方法来实现系统拆分。由于缺少足够的理论方法支撑，前述方法往往存在拆分所得子模型数量较多、建模压力较大以及再现色域损失等问题。对于以上问题，目前学术界及工业界中皆尚未提出从色域最大化角度实现打印印刷系统高效合理的拆分式建模方法。

发明内容

本发明的目的是为了解决背景技术中所述问题，提出一种用于色域最大化的打印印刷系统拆分式建模方法。

本发明的技术方案为一种用于色域最大化的打印印刷系统拆分式建模方法，具体包括以下步骤：

步骤1，根据预设采样间隔值X，制作并打印各基色的单色梯尺；

步骤2，确定每种单色墨量的超限阈值；

步骤3，通过等比压缩方法，依据步骤2所得各单色墨量的相应超限阈值对相应原单色梯尺进行重采样并打印印刷，测量等比压缩后打印印刷所得单色梯尺的光谱值，建立步骤1所得单色梯尺的原始墨量与等比压缩后打印印刷所得单色梯尺的光谱值之间的对应关系；

步骤4，设原始多色打印印刷系统墨色数为P，根据子模型基色墨色数量M，对原始多色打印印刷系统进行枚举式拆分，共得到N个基色墨色数为M的打印印刷系统子模型，其中，P＞M， $C_P^M=N;$

步骤5，对步骤4所得N个打印印刷系统子模型，根据预设采样间隔值Y进行墨色空间全局采样，获得所有打印印刷系统子模型下的采样样本；

步骤6，对所有打印印刷系统子模型下的采样样本仿真获取预测光谱反射率信息值，仿真方式如下，

首先，根据步骤3所得原始墨量与等比压缩后单色梯尺光谱值之间的对应关系，设任一可见光波段标记为λ，按式一进行计算，

${\left(\frac{K}{S}\right)}_{\mathrm{\λ}}=\frac{{(1-P_{\mathrm{\λ}})}^2}{2P_{\mathrm{\λ}}}$    式一

其中，P_λ为墨水颜色对某波段λ的光谱反射率分布，(K/S)_λ为该墨水对某波段λ的光的吸收特性K与散射特性S之间的比值；对于第i种单色墨水，其(K/S)_λ取值记为吸收散射特性比(K/S)_λ,i，按式二进行计算，

${\left(\frac{K}{S}\right)}_{\mathrm{\λ},i}={\left(\frac{K}{S}\right)}_{\mathrm{\λ},i\max }-{\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{paper}}$    式二

其中，(K/S)_λ，imax为第i种单色墨水对某波段λ的实地吸收散射特性比、(K/S)_λ，paper为对某波段λ的纸张吸收散射特性比，i的取值为1,2,…P；

其次，叠印颜色的吸收散射特性比(K/S)_λ，mix由纸张吸收散射特性比(K/S)_λ，paper、各基色吸收散射特性比(K/S)_λ,i计算求得，如式三，

${\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{mix}}={\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{paper}}+\underset{i}{\mathrm{\Σ}}{c}_i{\left(\frac{K}{S}\right)}_{\mathrm{\λ},i}$    式三

其中，c_i为第i种单色墨水墨量值；

再次，求各墨色墨量值在细胞分区形式下的百分制表示c_i,cell如式四，

$c_{i,\mathrm{cell}}=\frac{c_t-c_0}{c_{t,\mathrm{upper}}-c_0}$    式四

其中，c_t表示相应墨色实际墨量值，c₀表示墨量值零点，c_t,upper为墨色实际墨量值所在细胞区间上限值；

随后，求解细胞分区形式下叠印颜色吸收散射特性比(K/S)_λ,mix,cell如式五，

${\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{mix},\mathrm{cell}}={\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{paper}}+\underset{i}{\mathrm{\Σ}}{c}_{i,\mathrm{cell}}{\left(\frac{K}{S}\right)}_{\mathrm{\λ},i,\mathrm{cell}\left(t\right)}$           式五

其中，(K/S)_λ,i,cell(_t)表示各墨色实际墨量值所在细胞区间上限所对应的叠印颜色吸收散射特性比；

最后，预测求解最终叠印光谱值如式六，

$P\left(\mathrm{\λ}\right)=1+{\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{mix},\mathrm{cell}}-\sqrt{{\left({\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{mix},\mathrm{cell}}\right)}^2+2{\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{mix},\mathrm{cell}}}$      式六

其中，P(λ)为某波段λ的预测光谱反射率信息值；

步骤7，针对所有N个打印印刷系统子模型，使用凸包算法计算得到墨色空间全局采样色域体积；针对N个打印印刷系统子模型，进行K个子模型的组合，得到

种组合，并使用凸包算法获得每种组合的色域体积，K的初始值为2；

步骤8，计算步骤7所得每种组合的色域体积与墨色空间全局采样色域体积的比值，得到每种组合的色域覆盖率，当

种组合中最大色域覆盖率大于预设阈值时，锁定对应最大色域覆盖率的组合，实现原始打印印刷系统的拆分式建模；否则，令K=K+1，返回步骤7基于新的K值进行子模型的组合。

而且，步骤8中预设阈值为99%。

而且，M的值取4，Y=2×X。

本发明从多色打印印刷系统拆分建模的思想出发，通过打印单色基色梯尺，确定单色墨量超限阈值并据此对原始墨量信息进行等比压缩；根据单色墨量压缩前后的墨量对应关系，建立原始墨量与墨量压缩后打印印刷测量光谱的对应关系；通过Kubelka-Munk理论对墨色空间全局采样叠印色光谱反射率进行预测，结合色域覆盖率最大化分析方法，在打印印刷系统控制值空间用较少数量的子模型最大程度的逼近打印印刷系统的最大色域，在快速准确的完成多色打印印刷系统拆分式建模的同时，一定程度上避免了墨量超限问题，同时提高了多色复制叠印色彩的稳定性，促进系统更高效的发挥其色彩复制的能力，从而满足高保真彩色复制的需求。本发明在保证多色打印印刷系统模型建模精度的同时，大大降低了建模样本数量，且实施方便，在高保真色彩复制领域具有较高适用性。由于本发明技术方案具有重要应用意义，受到多个项目支持：1.国家自然基金项目61275172，2.国家文物局文物保护领域科学和技术研究一般课题2013-YB-HT-034，3.国家973基础研究子项目2012CB725302。对本发明技术方案进行保护，将对我国相关行业竞争国际领先地位具有重要意义。

附图说明

图1为本发明实施例的流程图。

具体实施方式

本发明技术方案具体实施时可由本领域技术人员采用计算机软件技术运行。结合附图，提供本发明实施例具体描述如下。

如图1所示,实施例提供了一种用于色域最大化的多色打印印刷系统拆分式建模方法，可以在提高建模效率、避免墨量超限以及提高叠印色彩稳定性的同时，用于色域最大化的打印印刷输出效果。实施例采用一台CMYKRGB七色喷墨打印机，其中CMYKRGB分别是青（Cyan）、洋红或品红（Magenta）、黄（Yellow）、黑（Black）、红（Red）、绿（Green）和蓝（Blue）七种颜色的简写。需要说明的是，本发明并不仅仅局限于七色打印机的拆分式建模应用，对于任意多色打印系统和印刷系统的3基色、4基色以及5基色等子模型的拆分式建模同样适用。

实施例主要包括以下步骤：

1）以X为采样间隔值，制作并打印各基色的单色梯尺。

具体实施时，单色梯尺的墨量间隔X可由本领域技术人员自行设定，建议值取5。实施例针对单色油墨，以步长5为采样间隔，生成21级（即取0,5,10,15,……100）单色梯尺并进行打印。

2）确定单色墨量超限的阈值。

如果墨量超限，会出现墨水溢出、墨水起皱等现象。具体实施时，可待打印的单色梯尺干燥之后，利用视觉判断的方法分别找出每种单色墨量超限的临界值，即确定为单色墨量限制的墨量阈值。例如，在对21级均匀采样的C色梯尺进行打印时，发现其墨量超限值为80（墨量值范围为0—100）,80即为C色墨量限制的阈值。

3）通过等比压缩方法，依据2）获得的超限阈值对原单色梯尺进行重采样，打印测量墨量压缩后打印印刷所得单色梯尺的光谱值，建立原始墨量与墨量压缩后打印印刷色彩光谱信息之间的对应关系。

实施例根据确定的各单色梯尺墨量超限阈值，通过等比压缩方法进行重采样，对二次重采样打印测量，获取其光谱反射率信息，由于单色墨量等比压缩前后的墨量信息为一一对应关系，为了后续建模的方便性，将原始墨量值与等比压缩后21级梯尺的光谱值建立一一对应关系。例如，在对C色梯尺进行打印时，发现其墨量超限阈值为80，故将原始0,5,10,15,……,95,100采样的21级梯尺压缩为0,4,8,12,……,76,80，级数不变，并将压缩后梯尺数据导入打印并测量获取光谱，将测量获得的光谱与0,5,10,15,……,95,100采样的21级梯尺墨量值建立一一对应关系。

4）使用数理统计中的排列组合方法，以M为子模型基色墨色数量，对墨色数为P（P＞M）的原始多色打印印刷系统进行枚举式拆分，共得到N个基色墨色数为M的打印印刷系统子模型，其中

M值可事先预设，建议M的值取4。

N个基色墨色数为M的打印印刷系统子模型可记为S1、S2、…、SN。实施例应用多色打印印刷系统拆分思想，将墨色数为7的喷墨打印机进行基色墨色数为4色的子打印机理论拆分，采用数学排列组合方法完成拆分N=

一共得到35个4基色子打印机模型。

5）对所有基色墨色数为M的N个打印印刷系统子模型，根据设定采样间隔值Y进行墨色空间全局采样，获得所有打印印刷系统子模型下的采样样本。

实施例进行墨色空间全局采样，计算获得所有打印印刷系统子模型下的采样样本，组成全局采样样本，建议Y=2×X，实施例中X的取值为5，Y的取值为10。对拆分得到的35个4基色子打印机模型，重新设定采样间隔步长为10（即取0,10,20,……100），进行墨色空间全局采样，计算获得所有字打印机模型下的采样样本，组成全局采样样本，其中样本个数为35*11⁴=512435个，每个4基色子打印机模型下有11⁴个采样样本。

6）以Kubelka-Munk浑浊介质呈色理论为基础，结合其细胞分区优化方法，利用如下系列公式完成各个子模型全局采样样本光谱的仿真计算工作。所有打印机子模型下采样样本即全局采样样本。具体如下：

首先，根据3）建立的原始墨量值与等比压缩后单色梯尺光谱值之间的对应关系，依据公式（1）计算每一单色墨水对应的各单色梯尺每个色样在不同波段的光的吸收特性K与散射特性S之间的比值。一般有多种可见光波段（例如31维或61维，实施例为31维），设任一可见光波段标记为λ，对任一单色墨水：

${\left(\frac{K}{S}\right)}_{\mathrm{\λ}}=\frac{{(1-P_{\mathrm{\λ}})}^2}{2P_{\mathrm{\λ}}}---\left(1\right)$

其中，P_λ为墨水颜色对某波段λ的光谱反射率分布，(K/S)_λ为该墨水对某波段λ的光的吸收特性与散射特性之间的比值。对于第i种单色墨水，其(K/S)_λ取值可记为吸收散射特性比(K/S)_λ,i，可以通过（2）式计算：

${\left(\frac{K}{S}\right)}_{\mathrm{\λ},i}={\left(\frac{K}{S}\right)}_{\mathrm{\λ},i\max }-{\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{paper}}---\left(2\right)$

其中，(K/S)_λ，imax为第i种单色墨水对某波段λ的实地吸收散射特性比、(K/S)_λ，paper为对某波段λ的纸张吸收散射特性比。i的取值为1,2,…P。

其次，对于叠印颜色，其吸收散射特性比(K/S)_λ，mix可由纸张吸收散射特性比(K/S)_λ，paper、各基色吸收散射特性比(K/S)_λ,i计算求得，如式（3）：

${\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{mix}}={\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{paper}}+\underset{i}{\mathrm{\Σ}}{c}_i{\left(\frac{K}{S}\right)}_{\mathrm{\λ},i}---\left(3\right)$

其中，c_i为第i种单色墨水墨量值。

再次，为提高Kubelka-Munk理论预测精度，用细胞分区方法求公式（3）中各墨色墨量值的百分制表示，此处以c_i,cell表示，如（4）式。其中，c_t表示相应墨色实际墨量值，c₀表示墨量值零点，c_t,upper为墨色实际墨量值所在细胞区间上限值。

$c_{i,\mathrm{cell}}=\frac{c_t-c_0}{c_{t,\mathrm{upper}}-c_0}---\left(4\right)$

随后，以公式（5）求解细胞分区形式下叠印颜色吸收散射特性比(K/S)_λ,mix,cell,其中c_i,cell为各墨色墨量值在细胞分区形式下的百分制表示，(K/S)_λ,i,cell(t)表示各墨色实际墨量值所在细胞区间上限所对应的叠印颜色吸收散射特性比

${\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{mix},\mathrm{cell}}={\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{paper}}+\underset{i}{\mathrm{\Σ}}{c}_{i,\mathrm{cell}}{\left(\frac{K}{S}\right)}_{\mathrm{\λ},i,\mathrm{cell}\left(t\right)}---\left(5\right)$

最后，使用（6）式预测求解最终叠印光谱值：

$P\left(\mathrm{\λ}\right)=1+{\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{mix},\mathrm{cell}}-\sqrt{{\left({\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{mix},\mathrm{cell}}\right)}^2+2{\left(\frac{K}{S}\right)}_{\mathrm{\λ},\mathrm{mix},\mathrm{cell}}}---\left(6\right)$

其中，P(λ)为某波段λ的最终叠印光谱值，即预测光谱反射率信息值。按此可得各波段相应的最终叠印光谱值。

最终获得打印机全局空间采样每一样本的预测光谱反射率信息值P(λ)，实施例获得512435个样本的预测光谱反射率信息值，以免于实际打印512435个样本。

7）对35个子打印机预测所得光谱、三维主成分降维后，使用凸包算法计算得到墨色空间全局采样色域体积，然后针对4）中拆分获得的N个子模型，进行K（K<N）个子模型的组合，并使用凸包算法获得其对应子模型组合的色域体积。

实施例获取墨色空间全局采样色域体积的方式为，将所得512435个光谱信息P(λ)，数据形式是512435×31的矩阵，通过主成分分析实现3维降维，使其变为512435×3的矩阵，然后用凸包算法得到总体积。

K的初值取2，每次执行完8）后，如果迭代返回，则K=K+1，再次执行7）时采用新的K值；

实施例中，迭代执行了7）共4次，K的取值从2到5，除第一次执行7）计算了所有35个子打印机模型的色域体积（后续迭代无需重复计算），4次迭代分别计算了使用2个、3个、4个和5个子打印机进行墨色组合的所有子打印机组合情况的色域体积。使用2个、3个、4个和5个子打印机分别有

种组合，每次迭代时对相应的每个组合预测所得光谱利用三维主成分降维后，使用凸包算法计算色域体积。

凸包算法为现有技术，本发明不予赘述。

8）计算色域覆盖率，如果最大色域覆盖率大于阈值，则停止迭代，如果否则K=K+1，返回7）。组合中若干个子打印机的体积除以总体积，即可得到色域覆盖率。通过与墨色空间全局采样色域比较，以最大色域覆盖率为选择依据，当最大色域覆盖率大于阈值时，锁定对应最大色域覆盖率的子模型组合，实现原始打印印刷系统的拆分式建模。

本步骤中，将获得最大色域覆盖率的一组K个子打印印刷系统模型的墨色组合，确立为最终打印印刷系统子模型，实现打印印刷系统拆分式建模。实施例的阈值设为99%。实施例根据之前7）计算的任意一个子打印机墨色组合的色域体积，通过与打印机墨色空间全局采样色域比较，例如K=2时，计算种组合分别的色域体积除以全局采样色域体积的比值，从

个结果里取最大值判断是否大于99%，结果为否，返回7）对K=3的组合进行处理…直到在K=5个子打印机墨色组合案例中，最大值判断是大于99%，RKYM CYGB RYCG RMCBCMYK这5个子打印机墨色组合具有最高色域覆盖率，其对全局采样色域的覆盖率为99.79%，被确立为最终子打印机模型，从而实现了7色喷墨打印机拆分式建模。

本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

Claims (3)

1.一种用于色域最大化的多色打印印刷系统拆分式建模方法，其特征在于，包括以下步骤：

步骤1，根据预设采样间隔值X，制作并打印各基色的单色梯尺；

步骤2，确定每种单色墨量的超限阈值；

步骤3，通过等比压缩方法，依据步骤2所得各单色墨量的相应超限阈值对相应原单色梯尺进行重采样并打印印刷，测量等比压缩后打印印刷所得单色梯尺的光谱值，建立步骤1所得单色梯尺的原始墨量与等比压缩后打印印刷所得单色梯尺的光谱值之间的对应关系；

步骤4，设原始多色打印印刷系统墨色数为P，根据子模型基色墨色数量M，对原始多色打印印刷系统进行枚举式拆分，共得到N个基色墨色数为M的打印印刷系统子模型，其中，P＞M， $C_P^M=N;$

步骤5，对步骤4所得N个打印印刷系统子模型，根据预设采样间隔值Y进行墨色空间全局采样，获得所有打印印刷系统子模型下的采样样本；

步骤6，对所有打印印刷系统子模型下的采样样本仿真获取预测光谱反射率信息值，仿真方式如下，

首先，根据步骤3所得原始墨量与等比压缩后单色梯尺光谱值之间的对应关系，设任一可见光波段标记为λ，按式一进行计算，

${\left(\frac{K}{S}\right)}_{\mathrm{\&lambda;}}=\frac{{(1-P_{\mathrm{\&lambda;}})}^2}{2P_{\mathrm{\&lambda;}}}$    式一

其中，P_λ为墨水颜色对某波段λ的光谱反射率分布，(K/S)_λ为该墨水对某波段λ的光的吸收特性K与散射特性S之间的比值；对于第i种单色墨水，其(K/S)_λ取值记为吸收散射特性比(K/S)_λ,i，按式二进行计算，

${\left(\frac{K}{S}\right)}_{\mathrm{\&lambda;},i}={\left(\frac{K}{S}\right)}_{\mathrm{\&lambda;},i\max }-{\left(\frac{K}{S}\right)}_{\mathrm{\&lambda;},\mathrm{paper}}$    式二

其中，(K/S)_λ，imax为第i种单色墨水对某波段λ的实地吸收散射特性比、(K/S)_λ，paper为对某波段λ的纸张吸收散射特性比，i的取值为1,2,…P；

其次，叠印颜色的吸收散射特性比(K/S)_λ，mix由纸张吸收散射特性比(K/S)_λ，paper、各基色吸收散射特性比(K/S)_λ,i计算求得，如式三，

${\left(\frac{K}{S}\right)}_{\mathrm{\&lambda;},\mathrm{mix}}={\left(\frac{K}{S}\right)}_{\mathrm{\&lambda;},\mathrm{paper}}+\underset{i}{\mathrm{\&Sigma;}}{c}_i{\left(\frac{K}{S}\right)}_{\mathrm{\&lambda;},i}$    式三

其中，c_i为第i种单色墨水墨量值；

再次，求各墨色墨量值在细胞分区形式下的百分制表示c_i,cell如式四，

$c_{i,\mathrm{cell}}=\frac{c_t-c_0}{c_{t,\mathrm{upper}}-c_0}$    式四

其中，c_t表示相应墨色实际墨量值，c₀表示墨量值零点，c_t,upper为墨色实际墨量值所在细胞区间上限值；

随后，求解细胞分区形式下叠印颜色吸收散射特性比(K/S)_λ,mix,cell如式五，

${\left(\frac{K}{S}\right)}_{\mathrm{\&lambda;},\mathrm{mix},\mathrm{cell}}={\left(\frac{K}{S}\right)}_{\mathrm{\&lambda;},\mathrm{paper}}+\underset{i}{\mathrm{\&Sigma;}}{c}_{i,\mathrm{cell}}{\left(\frac{K}{S}\right)}_{\mathrm{\&lambda;},i,\mathrm{cell}\left(t\right)}$       式五

其中，(K/S)_λ,i,cell(t)表示各墨色实际墨量值所在细胞区间上限所对应的叠印颜色吸收散射特性比；

最后，预测求解最终叠印光谱值如式六，

$P\left(\mathrm{\&lambda;}\right)=1+{\left(\frac{K}{S}\right)}_{\mathrm{\&lambda;},\mathrm{mix},\mathrm{cell}}-\sqrt{{\left({\left(\frac{K}{S}\right)}_{\mathrm{\&lambda;},\mathrm{mix},\mathrm{cell}}\right)}^2+2{\left(\frac{K}{S}\right)}_{\mathrm{\&lambda;},\mathrm{mix},\mathrm{cell}}}$     式六

其中，P(λ)为某波段λ的预测光谱反射率信息值；

步骤7，针对所有N个打印印刷系统子模型，使用凸包算法计算得到墨色空间全局采样色域体积；针对N个打印印刷系统子模型，进行K个子模型的组合，得到

种组合，并使用凸包算法获得每种组合的色域体积，K的初始值为2；

步骤8，计算步骤7所得每种组合的色域体积与墨色空间全局采样色域体积的比值，得到每种组合的色域覆盖率，当种组合中最大色域覆盖率大于预设阈值时，锁定对应最大色域覆盖率的组合，实现原始打印印刷系统的拆分式建模；否则，令K=K+1，返回步骤7基于新的K值进行子模型的组合。

2.根据权利要求1所述用于色域最大化的多色打印印刷系统拆分式建模方法，其特征在于：步骤8中预设阈值为99%。

3.根据权利要求1或2所述用于色域最大化的多色打印印刷设备拆分式建模方法，其特征在于：M的值取4，Y=2×X。

Images