本发明属于数码喷墨印花颜色管理领域，针对数码喷墨印花颜色管理中印染颜色与显示器颜色存在较大误差以及图像丢失纹理层次等问题，提出了一种基于空间网络用于纺织品喷墨印染的颜色色域映射方法。
背景技术：
现有的数码喷墨印花颜色管理系统中，实现显示器端rgb颜色到印染机端cmyk颜色的转换过程包含三个主要步骤：首先，利用显示器色彩特性文件iccprofile将rgb颜色值转换到标准颜色空间；然后利用公共色域映射算法实现显示器色域内颜色到印染机色域内颜色的映射；最后，利用印染机色彩特性文件iccprofile，将映射后的颜色值转换为cmyk颜色值。
在颜色管理软件进行颜色转换过程中，一个核心的部分为色域空间的映射，色域空间映射算法的好坏直接影响着颜色转换的一致性结果。目前很多研究者提出了不同的色域映射算法，众多的色域映射算法可以被不同的分类准则划分为不同的类别。
根据颜色再现流程中源色域与目标色域之间的位置关系，颜色色域映射算法可以分为：基于设备-设备(device-to-device)的色域映射算法和基于图像-设备(image-to-device)的色域映射算法。基于设备-设备的色域映射方法不需要考虑图像的特点，并且可以利用查找表等方法计算色域间的映射关系，其原理相对简单，易于理解且计算方便。但该方法与原稿图像的特征无关，源色域与目标色域的描述本身都是一次性确定，在同一流程中对任何图片都采用相同的映射方法映射，对未超出目标设备色域的原稿图像颜色信息常会造成过度裁切或压缩现象，造成图像细节信息过度丢失。基于图像-设备的色域映射算法要根据原稿图像的颜色特性进行，除了考虑颜色本身与目标色域的相对位置和大小关系外，还要考虑周围颜色信息的影响。此类色域映射方法计算复杂度高，但可以保留源图像中的空间细节。
根据对源色域中颜色信息的处理关系，可以将目前已有的色域映射算法分为两类：逐点色域映射算法(pointwisegammutmappingalgorithms,pgma)和空间域色域映射算法(spatialgamutmappingalgorithms,sgma)。逐点色域映射算法主要考虑在标准连接空间中将输入的颜色信息逐点映射到输出点上，而不考虑图像的类型和内容。广为流行和业内所接受的icc色彩管理流程就是基于pgma算法而建立的。空间色域映射算法不仅考虑输入图像的颜色特征，还考虑图像相邻像素之间的影响，根据图像的内容和空间位置关系，对图像的颜色分别进行处理，以获得高质量的复制图像，是目前的研究的难点之一，现阶段并未被工业界实际采用。
逐点色域映射算法中的大多数算法属于基于设备-设备的色域映射算法，对于逐点色域映射算法，英国学者morovic和luo进行了大量的研究和总结，并将其分为两大类：色域裁切算法和色域压缩算法。色域裁切算法的主要思想是将源色域中位于目标色域之外的颜色映射到目标色域的边缘，而位于目标色域之内的颜色保持不变。对于目标色域之外的颜色，此类算法指定某种映射准则，使其沿着特定的方向映射到目标色域的边界。常用的色域裁切算法有sclip算法,hpminδe算法和cusp算法等。色域压缩方法是将源色域内的颜色信息整体压缩到目标色域内，该方法不仅将目标色域之外的颜色压缩，目标色域之内的颜色也进行压缩。色域压缩方法主要包含线形压缩算法、sgck算法、s-clipping算法、llin、lnlin、lclip、slingclsp和cllin算法等。
icc色彩管理流程就是基于上述设备-设备色域映射算法建立的，该颜色管理流程在颜色映射一致性方面主要存在两方面的不足：颜色的转换存在较大的偏差、经转换后的图像存在纹理层次信息的丢失。图像裁切算法因进行了裁切，对于目标色域之外的图像颜色细节细信息丢失较多。图像压缩算法能较好的保持图像细节部分，但颜色信息失真较多，颜色饱和度也降低了。
技术实现要素：
本发明要克服现有技术的上述缺点，我们提出了一种基于空间网络用于纺织品喷墨印染的颜色色域映射方法。为了减少颜色转换过程中转换颜色空间产生的颜色损失，根据显示器显示整数rgb值的特点，本发明在rgb颜色空间以显示器颜色整数点为顶点建立颜色网络n＝(v,e,w)，每个顶点与其26个邻域利用无向边进行连接，边的权重w利用ciede2000距离公式计算。在建立好的颜色网络中，利用新提出的多源最短路径算法求得印染机色域外颜色到印染机色域边缘颜色的映射关系。本发明不仅可以保持色域内颜色的精准映射而且可以尽最大可能地保留图像的层次纹理信息。
基于空间网络用于纺织品喷墨印染的颜色色域映射方法，包含步骤：
1)实验样本设计；
印染机端，任意cmyk颜色值，可以通过cmk、cyk、myk三种中的某一种颜色值替代。对于cmyk四色通道，分别以0％、10％、25％、40％、55％、70％、85％、100％为间隔进行采样，获得每个通道的采样点后，以cmk、cyk、myk三种方式进行组合，通过去除冗余颜色块，获得1352个印染颜色样本。显示器端，对rgb值每个通道以10％为间隔选取11个等级的样本值，通过三个通道组合获得显示样本集；
2)rgb与cmyk之间样本对映射关系的建立；
rgb与cmyk之间样本对映射关系的建立主要包含三个步骤：首先通过测色仪分别获取显示器样本rgb与印染机样本cmyk的颜色值vrgb、vcmyk，将采集到的颜色值vrgb、vcmyk转换到同一观测环境下颜色值v′rgb、v′cmyk。然后，利用rbf神经网络建立v′rgb与rgb颜色之间的转换函数f：
rgb＝f(v),v∈v′rgb,rgb∈rgb
最后，利用转换函数f实现v′cmyk到rgb之间的转换，建立cmyk与rgb样本对之间的映射关系；
3)利用ciede2000距离在rgb颜色空间建立26邻域网络；
在rgb颜色空间建立颜色空间网络n＝(v,e,w)，其中v为顶点的集合，e为顶点间边的集合。网络以rgb整数点为顶点，每个顶点与其26个邻域点利用无向边连接，边的权重w利用ciede2000距离公式计算。
4)利用多源最短路径建立色域外颜色的映射关系；
基于新建立的颜色网络n＝(v,e,w)，新提出的多源最短路径算法可以为印染机色域外的每个点搜索到对应颜色最相近的色域边缘点，建立色域外颜色点到色域边缘颜色点的多对一映射关系。
多源最短路径算法中符号及公式的定义：
sb:印染机色域边缘颜色点集合；
q:存储已搜索到颜色点的队列；
q(ri,gi,bi)：点q(r,g,b)邻近的第i个点，其中i∈[1,26]，i为正整数；
w(q(ri,gi,bi),q(r,g,b))：点q(ri,gi,bi)和点q(r,g,b)之间的权重；
m:显示器色域中的每个点与其26邻域的ciede2000距离查找表；
t:三维查找表，输入某个rgb颜色值作为下标，输出为1表示该点在印染机色域内，输出为0则反之；
d:三维查找表，输入色域外某个rgb颜色值，输出距离该点最近色域边缘点的距离；
p:三维查找表，输入色域外某个rgb颜色值，输出距离该点最近的色域边缘点；
b:三维查找表，输入色域外某个rgb颜色值，输出该点26个邻域中距离和该点相近的后续节点；
l：三维查找表，输入某个色域外某个rgb颜色值，输出为0表示该点未被访问过，输出为1则反之。
u:三维查找表，输入色域外某个rgb颜色值，输出为1表示队列q中对应该点已失效，否则反之。
多源最短路径算法的求解过程如下：
输入：t，sb，m
①初始化：
初始化距离查找表d为0；
将所有sb存储到队列q中；
初始化访问查找表l为0；
初始化查找表u为0；
利用色域边缘点集合sb初始化最近点查找表p中边缘点距离
其最近的点为本身，其他点为空；
②判断队列q是否为空，空则跳转到
③获取队列q中队首的点q(r,g,b),并从队列中删除该点
④如果u(r,g,b)＝1跳转回②
⑤遍历点q(r,g,b)的26个邻域点，遍历结束后跳转到②
⑥如果点q(ri,gi,bi)在显示器色域外或t(ri,gi,bi)＝1，跳转回⑤继续遍历
⑦如果l(ri,gi,bi)＝1，跳转到⑨
⑧更新l(ri,gi,bi)为1，将q(ri,gi,bi)插入到队列q队尾，将q(ri,gi,bi)插入到b(r,g,b)相邻的后续节点，更新p(ri,gi,bi)的值为p(r,g,b)，更新d(ri,gi,bi)的值d(ri,gi,bi)＝d(r,g,b)+w(q(ri,gi,bi),q(r,g,b))，跳转到⑤
⑨如果d(ri,gi,bi)<d(r,g,b)+w(q(ri,gi,bi),q(r,g,b)),则跳转到⑤
⑩将q(ri,gi,bi)插入到队列q队尾，将q(ri,gi,bi)插入到b(r,g,b)相邻的后续节点，更新p(ri,gi,bi)的值为p(r,g,b)，更新d(ri,gi,bi)的值d(ri,gi,bi)＝d(r,g,b)+w(q(ri,gi,bi),q(r,g,b))
输出：p
本发明的优点：在减少颜色映射中存在损失的同时，实现色域内颜色精准映射，保留图像的层次纹理等信息。
具体实施方式
1)实验样本设计及颜色数据采集；
现印染机端，任意cmyk颜色值，都可以通过cmk、cyk、myk三种中的某一种组合进行替代。对于印染机的四个通道cmyk,分别以0％、10％、25％、40％、55％、70％、85％、100％为间隔进行采样，获得每个通道的采样点后，以cmk、cyk、myk三种方式进行组合，通过去除冗余颜色块，获得1352个颜色样本。在无色彩管理的介入的情况下，直接印染到纺织品上，利用测色仪获得cmyk样本对应颜色值vcmyk。显示器端，对rgb值每个通道以10％为间隔选取11个等级的样本值，通过显示器端显示样本值，利用测色仪获得对应颜色值vrgb。测色仪在采集颜色过程中，显示器端采集颜色与印染布匹上采集的颜色光照条件不一致，对于采集得到的颜色值vrgb、vcmyk,首先将其转换为视锥细胞刺激值，调整视锥细胞刺激值到标准观测环境，获得标准视锥细胞刺激值，将其转换为固定光照环境d50下颜色值v′rgb、v′cmyk。通过上述步骤，可将显示器端与印染布匹上采集得到的颜色值标准化到0～1之间。
2)rgb与cmyk之间样本对映射关系的建立；
获得颜色值v′rgb、v′cmyk后，利用rbf神经网络建立v′rgb与rgb颜色之间的转换函数f：
rgb＝f(v),v∈v′rgb,rgb∈rgb
然后，利用转换函数f将v′cmyk转换到rgb颜色值，获得cmyk与rgb样本的映射关系。
3)利用ciede2000距离在rgb颜色空间建立26邻域网络；
在rgb颜色空间建立颜色空间网络n＝(v,e,w)，其中v为顶点的集合，e为顶点间边的集合。网络以rgb整数点为顶点，每个顶点与其26个邻域点利用无向边连接，边的权重w利用ciede2000距离公式计算。
4)色域内与色域外颜色映射关系的建立；
显示器的色域比印染机色域大，且显示器色域包含印染机色域。对于印染机色域内的颜色，可以通过一对一建立映射关系；对于印染机色域外的颜色，现有裁切映射算法会丢失图片的层次纹理信息。上述多源最短路径算法可以建立色域外颜色映射关系，该方法建立的色域映射关系可以保证色域内颜色的精准映射、色域外颜色色差最小映射地同时保留更多的图像层次纹理信息。
通过上述五个步骤，可以在rgb颜色空间建立rgb到cmyk之间的映射查找表，实现显示器颜色到印染机颜色精准映射，最大程度保留图像的层次纹理信息。
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举，本发明的保护范围的不应当被视为仅限于实施例所陈述的具体形式，本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。