《一种基于色彩聚类的乱针绣铺底绣稿生成方法.pdf》由会员分享，可在线阅读，更多相关《一种基于色彩聚类的乱针绣铺底绣稿生成方法.pdf（13页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 104063562 A (43)申请公布日 2014.09.24 C N 1 0 4 0 6 3 5 6 2 A (21)申请号 201410334507.0 (22)申请日 2014.07.14 G06F 17/50(2006.01) (71)申请人南京大学 地址 210023 江苏省南京市栖霞区仙林大道 163号南京大学 (72)发明人孙正兴 胡佳高 李博 宋有成 (74)专利代理机构江苏圣典律师事务所 32237 代理人胡建华 (54) 发明名称 一种基于色彩聚类的乱针绣铺底绣稿生成方 法 (57) 摘要 本发明公开了一种基于色彩聚类的乱针绣铺 底绣稿生成方法，。
2、包括以下步骤：将输入图像转 换到合适的色彩空间，以便更准确的度量颜色之 间的差异程度；采用绣线库中可选绣线颜色对图 像进行颜色量化，统计量化图像的直方图，并对直 方图中占比最多的几种颜色进行聚类，从而得到 输入图像的构成颜色；利用图像构成颜色初始化 聚类算法，聚类得到每个像素的类别及绣制每一 类所需的绣线颜色；对聚类结果做优化处理，消 除聚类结果中的孤立的像素点及零散的小面积区 域，得到平滑的图像色块区域划分，生成铺底绣 稿；利用得到的铺底绣稿进行铺底绣制。 (51)Int.Cl. 权利要求书2页 说明书7页 附图3页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页。
3、 说明书7页 附图3页 (10)申请公布号 CN 104063562 A CN 104063562 A 1/2页 2 1.一种基于色彩聚类的乱针绣铺底绣稿生成方法，其特征在于，包括以下步骤： 步骤一，色彩空间转换：将输入图像Image转换到CIE Lab色彩空间； 步骤二，图像构成颜色分析：基于绣线颜色库，计算输入图像的颜色直方图，通过对直 方图进行分析得到图像构成颜色； 步骤三，图像颜色聚类：根据图像构成颜色对图像像素进行聚类，得到每个像素的类 别，及绣制每一类别所需的绣线颜色； 步骤四，铺底绣稿生成：对聚类结果进行优化，消除聚类结果中的孤立的像素点及零散 的小面积区域，得到平滑的划分区域，。
4、生成乱针绣铺底绣稿； 步骤五，铺底绣制：在步骤四得到的铺底绣稿上进行铺底绣制。 2.根据权利要求1所述的一种基于色彩聚类的乱针绣铺底绣稿生成方法，其特征在 于，步骤二中所述图像构成颜色分析包括以下步骤： 步骤21，图像颜色量化：找出与每个像素最接近的绣线颜色，并用该绣线颜色代替像 素原始颜色，得到量化图像QuaImg； 步骤22，计算颜色直方图：统计量化图像QuaImg中各种颜色所占的比例，并按占比从 大到小排序，得到量化图像QuaImg的颜色直方图Hist，取出占比最高的前N种颜色，N为满 足的最小N值，Hist i 表示直方图中的第i种颜色的占比，i取值1N； 步骤23，得到图像构成颜色：。
5、根据用户输入的颜色数K对步骤22得到的N种颜色使用 层次聚类算法进行聚类，得到K类的中心色即为图像的K个构成颜色。 3.根据权利要求2所述的一种基于色彩聚类的乱针绣铺底绣稿生成方法，其特征在 于，步骤三中所述图像颜色聚类包括以下步骤： 步骤31，设置初始聚类中心，即将步骤二中得到的图像构成颜色设置为聚类中心； 步骤32，对每个像素p，计算其与每个聚类中心色的色差，并将该像素划分到距离最近 的聚类中心，得到该像素p的类别l p 1,2,.,K，即将像素p划分到第l p 类； 步骤33，根据步骤32得到的所有像素的类别，重新计算聚类中心，新的聚类中心为绣 线颜色库中与属于该类别的所有像素的平均色色。
6、差最小的绣线颜色； 步骤34，判断步骤33迭代过程中聚类中心与步骤32比是否发生变化，若未发生变化， 则结束迭代过程，将当前聚类中心颜色定为绣制所需绣线颜色，转步骤四，否则转步骤32。 4.根据权利要求3所述的一种基于色彩聚类的乱针绣铺底绣稿生成方法，其特征在 于，步骤四中所述的铺底绣稿生成部分包括以下步骤： 步骤41，利用图的最小割算法进行平滑处理，重新分配各像素的类别，通过最小化如下 能量函数实现： E pImage D p (l p )+ p,q V(p,q)(l p ,l q )， 其中，D p (l p )为区域数据项，表示像素p与其类别l p 的匹配程度，V(p,q)为边界平滑 项。
7、，表示相邻像素p和像素q之间的不一致性，表示输入图像8-邻域系统中的无序像素 对集合，为权重因子，取值范围010，(l p ,l q )为指示函数，用于指示像素p的类 别l p 和像素q的类别l q 是否相同，定义如下： 权 利 要 求 书CN 104063562 A 2/2页 3 步骤42，计算优化后图像中各区域的像素个数，将区域像素个数小于阈值S的区域 合并到与其邻接边界最长的相邻区域中，最终生成包含区域划分信息和所选绣线颜色的铺 底绣稿，其中S表示输入图像Image的像素总数，取值范围02。 5.根据权利要求4所述的一种基于色彩聚类的乱针绣铺底绣稿生成方法，其特征在 于，步骤41中所述区。
8、域数据项D p (l p )计算公式如下： D p (l p )-log(Prob(p,l p )； 其中，Prob(p,l p )表示像素p属于类别l p 的概率，利用像素p的实际颜色与步骤三得 到的每一类的聚类中心的色差计算得出，计算公式如下： 其中，Image(p)表示输入图像中像素p的颜色，ClusterColor(j)表示步骤三中得到的 第j类的聚类中心颜色，j1,.,K。 6.根据权利要求5所述的一种基于色彩聚类的乱针绣铺底绣稿生成方法，其特征在 于，步骤41中所述边界平滑项V(p,q)计算公式如下： 其中，Image(q)表示输入图像中像素q的颜色，variance为集合中所有像。
9、素对之间 色差的方差，d(p,q)为像素p和像素q在输入图像中的坐标距离。 权 利 要 求 书CN 104063562 A 1/7页 4 一种基于色彩聚类的乱针绣铺底绣稿生成方法 技术领域 0001 本发明涉及一种乱针绣铺底绣稿的计算机自动生成方法，属于计算机辅助设计技 术领域，具体地说是一种基于色彩聚类的乱针绣铺底绣稿生成方法。 背景技术 0002 乱针绣是由常州艺人杨守玉教授于二十世纪三十年代创立的一种刺绣艺术，具 有精美的艺术价值和深厚的文化价值。目前江苏省和常州市已将其列入首批非物质文化 遗产，中国政府也已将其列入世界非物质文化遗产申请计划。乱针绣的绣制过程主要分 为三个阶段：填充大面。
10、积色块的铺底、勾勒对象区域及其轮廓的粗绣和刻画对象区域显 著细节的精绣，且乱针绣的艺术价值主要体现在绣品粗绣和精绣阶段。目前已有多种计 算机辅助乱针绣制作方法，例如：中国专利“一种计算机辅助乱针绣制作方法”，专利号： ZL201010175217.8；中国专利“一种采用遗传算法的计算机辅助乱针绣制作方法”，公开号： CN201310211288.2；中国专利“一种基于反应扩散方程的计算机辅助乱针绣制作方法”，公 开号：CN201310456178.2和中国专利“采用模糊聚类和随机行走的计算机辅助乱针绣制作 方法”，公开号：CN201310539428.9；等，但这些方法均侧重于面向整个绣制过程。
11、中的绣线选 择和针法生成，以期实现乱针绣绣品各阶段的完全数字化制作，并未单独考虑每一阶段绣 制。然而，计算机辅助乱针绣制作的核心并非在于实现绣制过程的全自动化，而是在于保持 其艺术价值的有效人机协作方式。由于乱针绣的铺底主要进行衬托主题的大面积色块填 充，其针法相对简单而又含大量重复工作，因此，数字化乱针绣铺底绣制与人工化粗绣和精 绣的有机结合是实现计算机辅助乱针绣制作的最实用途径。而数字化乱针绣铺底绣制的关 键在于乱针绣铺底绣稿的自动生成，其主要涉及绣线颜色选择和色块区域划分两个问题。 0003 由于在绣制过程中，可选的绣线颜色数量有限，通常只有几百种，例如绣线生产厂 商DMC的绣线颜色表包。
12、含454种颜色，而待绣制的图像通常具有大量的颜色，如：常见的24 位的真彩色图像最多可能有2 24 16777216种颜色。另外，受到实际可操作性的约束，实 际绣制过程中不可能选择所有的几百种绣线颜色，通常只是从其中选择少量几种颜色进行 绣制，比如，常见的电脑绣花机可以同时使用的颜色数只有6、9或12种，如果使用的绣线颜 色数较多，绣制时需要频繁地停机换线，影响机器绣制的速度。因此，需要根据输入图像的 色彩分布特征从绣线颜色表中选择一个子集，使用少量的颜色来重建输入图像。这种采用 较少数量的颜色来重建具有丰富色彩的图像是数字图像处理基本技术色彩量化的主要研 究内容。目前已有多种算法，比较传统的。
13、算法有统一量化法、频度序列法、中位切割法、中值 分裂法、八叉树法和聚类法等。近年来，部分研究者又提出了一些新的方法或改进方法以改 善重建图像的视觉效果，例如文献1：周兵,沈钧毅,彭勤科.CoQuWeiP使用加权积的 色彩量化算法J.计算机辅助设计与图形学学报,2004,16(4):470-475.中提出的加权积 色彩量化算法CoQuWeiP；文献2：Hu X,Wang T,Li D.A new approach of color quantization based on ant colony clustering algorithmC.Information Technology:Codin。
14、g and Computing,2005.ITCC2005.International Conference on.IEEE,2005,1:102-108.中 说 明 书CN 104063562 A 2/7页 5 提出的基于蚁群聚类的色彩量化算法；文献3：Ozdemir D,Akarun L.Fuzzy algorithms for combined quantization and ditheringJ.Image Processing,IEEE Transactions on,2001,10(6):923-931.中提出的结合色彩抖动的色彩量化；文献4：Baolong G,Xiang F.。
15、A modied octree color quantization algorithmC.Communications and Networking in China,2006.ChinaCom06.First International Conference on.IEEE,2006:1-3.中的 改进八叉树等。这些方法应用于乱针绣铺底绣线颜色选择时存在一个共同的问题，即它们 选择出的颜色集并不是从可选的绣线颜色集中选出的，这些算法按照某种优化条件求解得 到全局和局部最优解，如果从绣线颜色表寻找相似颜色来替代这些算法的全局最优或局部 最优解，则得到的解可能与最优解相差很大。 0004 另一。
16、方面，对于区域划分问题，文献5：田绪安,李丙社,徐一兵.刺绣图像 轮廓图的计算机生成系统J.西北纺织工学院学报,1998,2.使用边缘检测方法，得 到区域划分结果。然而该方法更侧重于区域轮廓，要得到闭合的区域需要更复杂的操 作。文献6：陈兆乾,谢俊元.刺绣打版系统中彩色空间聚类的方法与应用J.软件学 报,1996,7(A00):85-90.中使用颜色聚类方法得到区域划分结果。但是该方法仅利用了像 素的颜色信息，而未考邻域像素的信息，导致聚类结果中会有大量零散的面积很小的区域， 给绣制带来不便。文献7：陈圣国,孙正兴,项建华,等.计算机辅助乱针绣制作技术研 究J.计算机学报,2011,34(3)。
17、:526-538.中先对图像进行平滑处理，然后用绣线颜色替 代平滑后的像素颜色。但这样并未完全解决大量零散小区域存在的问题，且平滑过程会改 变图像原有颜色，导致绣制结果会产生一定色差。 发明内容 0005 发明目的：本发明所要解决的技术问题是针对现有技术的不足，提供一种适用于 乱针绣铺底绣制使用的绣稿，用于支持计算机辅助乱针绣铺底绣制。 0006 为了解决上述技术问题，本发明公开了一种基于色彩聚类的乱针绣铺底绣稿生成 方法，包括以下步骤： 0007 步骤一：色彩空间转换：为了能够准确的将颜色在数值上的距离映射到人眼感知 色差上，将输入图像Image转换到CIE Lab色彩空间。 0008 步骤。
18、二，图像构成颜色分析：基于绣线颜色库，计算输入图像的颜色直方图，通过 直方图分析得到图像的构成颜色。 0009 步骤三，图像颜色聚类：在图像构成颜色的指导下对图像像素进行聚类，得到每个 像素的类别，及绣制每一类所需的绣线颜色。 0010 步骤四，铺底绣稿生成：对聚类结果进行优化，消除聚类结果中的孤立的像素点及 零散的小面积区域，得到平滑的划分区域，生成乱针绣铺底绣稿。 0011 步骤五，铺底绣制：在步骤四得到的铺底绣稿上进行铺底绣制。 0012 本发明步骤二中所述图像构成颜色分析包括以下步骤：步骤21，图像颜色量化： 找出与每个像素最接近的绣线颜色，并用该绣线颜色代替像素原始颜色，得到量化图像。
19、 QuaImg。 0013 步骤22，计算颜色直方图：统计量化图像QuaImg中各种颜色所占的比例，并按占 比从大到小排序，得到量化图像QuaImg的颜色直方图Hist，取出占比最高的前N种颜色，N 说 明 书CN 104063562 A 3/7页 6 为满足的最小N值，Hist i 表示直方图中的第i种颜色的占比，i取值1 N。 0014 步骤23，得到图像构成颜色：根据用户输入的颜色数K对步骤22得到的N种颜色 使用层次聚类算法进行聚类，得到K类的中心色即为图像的K个构成颜色，即为图像的构成 颜色。 0015 本发明步骤三中所述图像颜色聚类包括以下步骤：步骤31，设置初始聚类中心， 即将步。
20、骤二中得到的图像构成颜色设置为聚类中心。 0016 步骤32，对每个像素p，计算其与每个聚类中心色的色差，并将该像素划分到距离 最近的聚类中心，得到该像素p的类别l p 1,2,.,K，即将像素p划分到第l p 类。 0017 步骤33，根据步骤32得到的所有像素的类别，重新计算聚类中心，新的聚类中心 为绣线颜色库中与属于该类别的所有像素的平均色色差最小的绣线颜色。 0018 步骤34，判断步骤33迭代过程中聚类中心与步骤32比是否发生变化，若未发生变 化，则结束迭代过程，将当前聚类中心颜色定为绣制所需绣线颜色，转步骤四，否则转步骤 32。 0019 本发明步骤四所述铺底绣稿生成包括以下步骤：。
21、步骤41，利用图的最小割算法进 行平滑处理，通过最小化如下能量函数实现： 0020 E pImage D p (l p )+ p,q V(p,q)(l p ,l q )， 0021 其中，D p (l p )为区域数据项，表示像素p与其类别l p 的匹配程度，V(p,q)为边界平 滑项，表示相邻像素p和像素q之间的不一致性，表示输入图像8-邻域系统中的无序像 素对集合，为权重因子，取值范围010，(l p ,l q )为指示函数，用于指示像素p的 类别l p 和像素q的类别l q 是否相同，定义如下： 0022 0023 步骤42，计算优化后图像中各区域的像素个数，将区域像素个数小于阈值S的 。
22、区域合并到与其邻接边界最长的相邻区域中，最终生成包含区域划分信息和所选绣线颜色 的铺底绣稿，其中S表示输入图像Image的像素总数，取值范围02。 0024 步骤41中所述区域数据项D p (l p )计算公式如下： 0025 D p (l p )-log(Prob(p,l p )； 0026 其中，Prob(p,l p )表示像素p属于类别l p 的概率，利用像素p的实际颜色与步骤 三得到的每一类的聚类中心的色差计算得出，计算公式如下： 0027 0028 其中，Image(p)表示输入图像中像素p的颜色，ClusterColor(j)表示步骤三中得 到的第j类的聚类中心，j1,.,K。 0。
23、029 步骤41中所述边界平滑项V(p,q)计算公式如下： 说 明 书CN 104063562 A 4/7页 7 0030 0031 其中，Image(q)表示输入图像中像素q的颜色，variance为集合中所有像素对 之间色差的方差，d(p,q)为像素p和像素q在输入图像中的坐标距离。 0032 本申请基于绣线色彩库，对图像进行构成颜色分析，在构成颜色的指导下对图像 进行颜色聚类，得到绣线颜色，并在聚类后使用平滑算法消除聚类结果中零散的面积很小 的区域，从而得到平滑的色块区域划分，生成乱针绣铺底绣稿，以便进行计算机辅助铺底绣 制。 0033 有益效果：本发明具有以下优点：首先，本发明可以在仅。
24、需用户输入颜色数目的 前提下自动的对输入图像按颜色进行区域划分。其次，本发明能够自动的选择出绣制各个 区域所需的绣线颜色。 附图说明 0034 下面结合附图和具体实施方式对本发明做更进一步的具体说明，本发明的上述和 /或其他方面的优点将会变得更加清楚。 0035 图1是本发明的处理流程示意图。 0036 图2a是实施例的输入图像示例示意图。 0037 图2b是对图2a输入图像采用本申请所述方法进行处理，生成的铺底绣稿示意图。 0038 图3是按图2b提供的划分及绣线颜色进行铺底绣制的结果。 0039 图4a是对图2a输入图像进行构成颜色分析过程中，利用绣线颜色库进行颜色量 化的结果。 0040。
25、 图4b是对图2a输入图像进行构成颜色分析，得到的图像构成颜色。 0041 图5是对图2a输入图像进行颜色聚类，得到的聚类结果示意图。 具体实施方式 ： 0042 如图1所示，本发明公开的是一种基于色彩聚类的乱针绣铺底绣稿生成方法，具 体包括以下步骤： 0043 步骤一，色彩空间转换：为了能够更准确的度量颜色之间的差异，需要将输入图像 转换到合适的色彩空间。 0044 步骤二，图像构成颜色分析：基于绣线颜色库，计算输入图像的颜色直方图，通过 直方图分析得到图像的构成颜色。 0045 步骤三，图像颜色聚类：在图像构成颜色的指导下对图像像素进行聚类，得到每个 像素的类别，及绣制每一类所需的绣线颜色。
26、。 0046 步骤四，铺底绣稿生成：对聚类结果进行优化，消除聚类结果中的孤立的像素点及 零散的小面积区域，得到平滑的划分区域，生成乱针绣铺底绣稿。 0047 步骤五，铺底绣制：在步骤四得到的铺底绣稿上进行铺底绣制。 0048 本发明的输入图像为自然界的可供绣制的任意彩色图像。 0049 下面具体介绍各个步骤的主要流程： 说 明 书CN 104063562 A 5/7页 8 0050 1.色彩空间转换 0051 由于常用的RGB色彩空间模型中，各种颜色在数值上的差距并不能很好的反 映人眼直观上对色彩差异的感觉。针对这个问题，可以先将图像色彩空间转换到CIE Lab(CIE法语：Commissio。
27、n internationale de lclairage，国际照明委员会。Lab： CIE1976(L*,a*,b*)色彩空间)或CIE Luv(CIE1976(L*,u*,v*)色彩空间)等能够直观的 反映人眼感知色差的颜色空间上。本申请将输入图像转换到CIE Lab色彩空间上进行计算。 CIE Lab色彩空间规范中将颜色C 1 与C 2 之间的色差定义为C 1 、C 2 在L、a、b三个通道上的欧 氏距离： 0052 0053 2.图像构成颜色分析 0054 基于绣线颜色库ThreadColor，分析输入图像Image的构成颜色，需用户事先指定 颜色聚类类别数K。绣线颜色库ThreadC。
28、olor内包括现有所有绣线的颜色。包括以下几个 步骤： 0055 步骤21，图像颜色量化：即用绣线颜色表示输入图像。对每个像素p，计算其与 每种绣线颜色的距离。建立量化图像QuaImg，令量化图像中每个像素的颜色为与输入图 像中对应位置像素颜色最接近的绣线颜色，即：QuaImg(p)ThreadColor(x)，其中x满足 Dis(ThreadColor(x),Image(p)Dis(ThreadColor(y),Image(p)，Dis()表示 在步骤1中的色彩空间上计算两种颜色之间色差的函数。 0056 步骤22，计算颜色直方图：步骤21得到的量化图像QuaImg仅由少量有限的颜色 构成，。
29、先统计QuaImg中各种颜色所占的比例，并按占比从大到小排序，得到QuaImg的颜色 直方图Hist。取出其中前N种颜色，N的为满足的最小N值，Hist i 表示直 方图中的第i种颜色的占比，i取值1N，即保证所有取出的颜色总占比不小于80的同 时，使取出的颜色最少。 0057 步骤23，得到图像构成颜色：对步骤22得到的N种颜色使用文献8： Defays D.An efficient algorithm for a complete link methodJ.The Computer Journal,1977,20(4):364-366.中提出的层次聚类算法进行聚类，得到K类的中心色即为 图。
30、像的K个构成颜色MainColor1,.,K。 0058 3.图像颜色聚类 0059 在步骤二中得到的图像构成颜色的指导下，基于绣线颜色库ThreadColor在输入 图像Image上进行K均值聚类，得到每个像素的类别，及绣制每一类所需的绣线颜色。包括 以下几个步骤： 0060 步骤31，将步骤二中得到的图像构成颜色设置为聚类中心： 0061 ClusterColor(m)MainColor(m),m1,K； 0062 ClusterColor(m)表示第m个聚类中心的颜色。 0063 步骤32，对每个像素p，计算其与每个聚类中心的距离，并将该像素划分到距离最 近的那个聚类中心上，得到该像素的。
31、类别l p 1,2,.,K，即将像素p划分到第l p 类。 0064 步骤33，根据所有像素的划分，重新计算聚类中心。新的聚类中心为绣线颜色库中 说 明 书CN 104063562 A 6/7页 9 与类别平均色最接近的绣线颜色。即对每一类别t，先计算出该类所包含所有像素的平均 色，然后计算该平均色与所有绣线颜色的距离，将距离最小的绣线颜色ThreadColor(n)设 置为此类新的聚类中心，即ClusterColor(t)ThreadColor(n)。 0065 步骤34，判断步骤33迭代过程中聚类中心与步骤32比是否发生变化，若未发生变 化，则结束迭代过程，将当前聚类中心颜色定为绣制所需绣。
32、线颜色，转步骤四，否则转步骤 32。 0066 4.铺底绣稿生成 0067 对步骤三中得到的聚类结果进行优化，消除聚类结果中孤立的像素点及零散的小 面积区域，得到均匀平滑的划分区域，生成铺底绣稿。包括以下几个步骤 0068 步骤41，利用图的最小割算法进行平滑处理，重新分配各像素的类别。本发明通过 最小化如下能量函数实现： 0069 E pImage D p (l p )+ p,q V(p,q)(l p ,l q ) 0070 其中，D p (l p )为区域数据项，表示像素p与其类别l p 的匹配程度，D p (l p )越小，表 示两者越匹配，像素p的属于类别l p 的可能性越大；V(p,。
33、q)为边界平滑项，表示像素p和像 素q之间的不一致性，V(p,q)越小，表示像素p与像素q越不一致，两者划分到不同区域的 可能性越大；表示原图像8-邻域系统中的无序像素对集合，为权重因子，取值范围 010，(l p ,l q )为指示函数，指示像素p的类别l p 和像素q的类别l q 是否相同，定义如 下： 0071 0072 区域数据项D p (l p )设置为像素p的属于类别l p 的概率Prob(p,l p )的负对数，即： 0073 D p (l p )-log(Prob(p,l p )； 0074 其中，像素p属于类别l p 的概率Prob(p,l p )利用像素p的实际颜色与步骤三。
34、得到 的每一类的聚类中心的色差计算得出，计算公式如下： 0075 0076 利用如下公式计算出边界平滑项V(p,q)： 0077 0078 其中，variance为中所有像素对之间色差的方差，d(p,q)为像素p和像素q在 图像中的坐标距离。权重因子可设置为8。 0079 本申请使用文献9：Boykov Y,Veksler O,Zabih R.Fast approximate energy minimization via graph cutsJ.IEEE Transactions on Pattern Analysis and Machine I ntelligence,2001,23(11。
35、):1222-1239.中描述的方法求解该能量最小化问题。 0080 步骤42，计算优化后图像中各区域的像素个数，将区域像素个数小于阈值S的 说 明 书CN 104063562 A 7/7页 10 区域合并到与其邻接边界最长的相邻区域中。通常可设为0.5，S表示输入图像Image 的像素总数。最终生成的铺底绣稿包含区域划分信息和步骤三得到的绣线颜色。 0081 5.铺底绣制 0082 在步骤四得到的铺底绣稿上进行铺底绣制。采用“中国专利ZL201010175217.8一 种计算机辅助乱针绣制作方法”中的绣制方法，使用所选出绣线颜色对划分得到的平滑色 块区域进行铺底绣制。 0083 实施例 00。
36、84 本实施例中，绣线颜色库为绣线生产厂商DMC的454种绣线颜色。如图2a为输入 的待处理图像，通过本发明所述的一种基于色彩聚类的乱针绣铺底绣稿生成方法，可从图 2a中的输入图像得到如图2b的铺底绣稿，图2b中包含区域划分和所选绣线颜色，每个区域 中括号里面的数值表示从绣线颜色库中选择出来的用来铺底绣制该区域的绣线颜色的RGB 值。最终能够得到图3所示的铺底绣制效果。由于本发明的特殊性，本实施例中的附图，需 要使用灰度效果显示。 0085 具体实施过程如下： 0086 输入如图2a的图像后，用户输入图像颜色聚类数K2。 0087 步骤一中，将输入图像转换到CIE Lab色彩空间。 0088 。
37、步骤二中，对输入图像进行主色调分析。先进行颜色量化，计算出每个像素颜色最 接近的绣线颜色，用对应的绣线颜色替换原像素颜色，得到如图4a所示的由有限的离散颜 色组成的图像。然后统计量化后图像的颜色直方图，对占比最大的前N种颜色进行聚类，得 到如图4b所示的图像构成颜色。 0089 步骤三中，利用图像构成颜色作为初始聚类中心，对输入图像进行颜色聚类，聚类 收敛后得到新的中心色和每个像素对应的中心，如图5。这里的聚类中心色是绣线颜色库中 的某一种颜色。 0090 步骤四中，对图5的聚类结果进行优化，使用图的最小割算法优化后，合并像素个 数少于0.5S的区域后，得到乱针绣铺底绣稿，结果如图2b。 00。
38、91 步骤五中，利用步骤四得到的铺底绣稿中给出的绣线颜色和色块区域划分结果， 进行模拟的铺底绣制，绣制结果如图3。 0092 本发明提供了一种基于色彩聚类的乱针绣铺底绣稿生成方法，具体实现该技术方 案的方法和途径很多，以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的 普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进 和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以 实现。 说 明 书CN 104063562 A 10 1/3页 11 图1 图2a 说 明 书 附 图CN 104063562 A 11 2/3页 12 图2b 图3 说 明 书 附 图CN 104063562 A 12 3/3页 13 图4a 图4b 图5 说 明 书 附 图CN 104063562 A 13 。