1.一种基于子空间细分插值的打印机配方计算方法，其特征在于包括以下步骤：

1)对C、M、Y每一个通道，从设备驱动最小值到设备驱动最大值以一定的间隔选定n个点，共组成n×n×n个训练样本；

2)通过打印机打印这n×n×n个训练样本，并利用分光光度计测量所有训练样本得到CIELAB值，从而建立了n×n×n的CMY到CIELAB查找表；

3)对于目标CIELAB值(LAB)T，在查找表中找到距离目标点最近的点(LAB)C，其对应的C、M和Y通道驱动值分别记作dC,k，dM,1和dY,m，其中k、1、m分别表示C、M和Y通道采样点的序号，其范围均在1到n之间；(LAB)T与点(LAB)C之间的距离可以通过CIELAB空间中的欧几里德距离或CIE1976L*a*b*的色差公式计算得到；

4)根据点(LAB)C与点(LAB)T之间的距离判断(LAB)C是否满足要求(小于预设的阈值)，如果满足要求，则(LAB)C对应的驱动值(dC,k，dM,1，dY,m)即为所求的驱动值；否则，以(dC,k，dM,1，dY,m)为中心点，以(dC,k+1-dC,k-1)/2为边长建立新的子空间，即(dC,k-1+dC，k)/2、(dC，k+dC，k+1)/2，(dM,1-1+dM，1)/2、(dM,1+dM,1+1)/2，(dY,m-1+dY,m)/2、(dY,m+dY,m+1)/2六个平面组成新子空间的8个顶点。注意要确保(dC,k-1+dC，k)/2、(dC，k+dC，k+1)/2、(dM,1-1+dM,1)/2、(dM,1+dM,1+1)/2、(dY,m-1+dY,m)/2、(dY,m+dY,m+1)/2均在0～255范围内，如果dC,k为0，则应选取dC,k、(dC,k+dC,k+1)/2与其它通道采样点组成新的子空间；如果dC，k为255，则应选取dC,k、(dC,k-1+dC，k)/2与其它通道采样点组成新的子空间。同样地，M和Y通道也做类似处理。同时，确保新的子空间顶点对应的C、M和Y通道驱动值为整数，如果不是整数，则按照四舍五入原则将其强制变换为整数；

5)通过基于立方体插值的正向色度特征化方法预测新的子空间顶点对应的CIELAB值，建立新的查找表；

6)查找该子空间顶点中与(LAB)T最近的点，若该最近的点与(LAB)T之间小于预设的阈值，则其对应的驱动值即为(LAB)T所求的驱动值；若大于预设的阈值，则判断新的子空间的所有边长是否均为1，若均为1，则该最近的点对应的驱动值即为(LAB)T所求的驱动值，若不均为1，则以该最近点为中心点，以原子空间边长的一半建立新的子空间，再重新通过第4)步开始迭代，直到得到所求的C、M和Y通道驱动值为止。

2.根据权利要求1所述的基于子空间细分插值的打印机配方计算方法，其特征在于所

述步骤2)中通过打印机打印这n×n×n个训练样本，并利用分光光度计测量所有训练样本得到CIELAB值，从而建立了n×n×n的CMY到CIELAB查找表。本专利也适用于除CIELAB色彩空间外的其他色彩空间，例如CIEXYZ色彩空间，CIECAM02中的色彩空间等。

3.根据权利要求1所述的基于子空间细分插值的打印机配方计算方法，其特征在于所

述步骤4)中根据点(LAB)C与点(LAB)T之间的距离判断(LAB)C是否满足要求(小于预设的阈值)，如果满足要求，则(LAB)C对应的驱动值(dC,k，dM,1，dY,m)即为所求的驱动值；否则，以(dC,k，dM,1，dY,m)为中心点，以(dC，k+1-dC,k-1)/2为边长建立新的子空间，即(dC,k-1+dC，k)/2、(dC，k+dC，k+1)/2，(dM,1-1+dM,1)/2、(dM,1+dM,1+1)/2，(dY,m-1+dY,m)/2、(dY,m+dY,m+1)/2六个平面组成新子空间的8个顶点；本专利也可适用于以(dC,k，dM,1，dY,m)为中心点，以其他边长建立新的子空间，例如以3×(dC,k+1-dC,k-1)/4为边长建立新的子空间。