基于图像的数字油画画布的自动生成方法
【专利摘要】基于图像的数字油画画布的自动生成方法，能够将任意图像自动转换为数字油画画布，并合理提示笔刷的走向和大小的方法。该方法包括以下步骤：首先,提出一种统计滤波方法，去除Mean-Shift图像分割产生的边界锯齿，获得光滑的边界线条；其次，采用K均值聚类选出代表性颜色，并提出一种自适应数字显示调整方法，保证相应数字的显示清晰可见、不与边界重叠；然后，计算图像空间的方向场，生成能够提示笔刷大小和走向的流线；最后，合成上述边界线条、颜色数字和笔刷流线，生成与输入图像内容一致的数字油画画布。
【专利说明】基于图像的数字油画画布的自动生成方法
【技术领域】
[0001] 本发明涉及一种利用计算机技术的数字油画画布的自动生成方法，能够将任意的 一幅图像自动转换为有效的数字油画画布，在显示区域轮廓和数字编号的同时，清晰地提 示出笔刷的大小和走向。
【背景技术】
[0002] 传统的油画是在一张空白的画布上完成的，需要轮廓勾勒、颜料调配、笔刷绘制等 许多步骤，专业性极强，需要长期的学习才能有所掌握。数字油画是一种新的艺术形式，它 提供给绘画者带有轮廓、颜色提示信息的画布和预先调配好的颜料，从而使所有喜爱绘画 但又缺乏绘画基础的人都能直接上手，轻松绘制出自己的绘画作品。数字油画继承了儿童 填色画的简单直观，避免了传统油画的专业限制，具有较高的观赏性和体验性，深受广大 业余绘画者的喜爱，自上世纪五十年代发明以来，迅速流行于世界各地。
[0003] 然而，目前的数字油画画布都需要专业人员的精心设计，因此市场上在售的数字 油画大都以特定的风景画为主，难以个性化定制。事实上，数字油画画布的设计极具挑战 性：需要用线条仔细逐一勾勒出整个场景中颜色相近的区域，精心挑选出具有代表性的若 干种颜色，并为每个区域标记上代表这些颜色的数字符号。一个画面丰富的场景往往具有 成百上千个区域，设计过程不仅难度极较大且费时费力。因此，亟需一种数字油画画布的自 动生成方法，打破其专业设计的局限，进一步提高个性化定制的体验性。
[0004] 艺术家在绘制油画过程中，除了颜料的颜色，每一次运笔的笔刷大小和走向也颇 为讲究。一般而言，在背景等面积较大的区域往往使用较大的笔刷，而在需要表现细节的区 域中使用较小的笔刷；笔刷的走向则是主要沿着区域的边界或区域内纹理的走向。目前的 数字油画画布只提示了颜料的使用区域和颜色，忽略了笔刷大小和走向等其他重要因素， 缺乏绘画经验的人容易产生笔触杂乱的问题或如同儿童填色画一样缺乏层次感。因此，亟 需一种在数字油画画布上提示笔刷大小和走向的方法，进一步改善数字油画的艺术效果。
【发明内容】
[0005] 为克服数字油画画布的设计困难，本发明提供了一种利用计算机技术的数字油画 画布的自动生成方法，能够将任意的一副图像自动转换为有效的数字油画画布，在显示区 域轮廓和数字编号的同时，清晰地提示出笔刷的大小和走向。
[0006] 该方法的具体步骤如下：
[0007] 步骤1、用Mean-Shift图像分割算法对输入图像I进行分割，并用每个区域的平均 颜色代替原有颜色，得到分割图像Ir。
[0008] 步骤2、将Ir的区域边界通过统计滤波进行边界光顺，得到边界图像Ie。
[0009] 步骤3、采用K均值算法将Ir的每个区域的平均颜色进行聚类分析。
[0010] 步骤4、根据Ir中每个区域的轮廓，自适应调整与其对应的颜色编号的大小和位 置。
[0011] 步骤5、根据输入图像I和边界图像Ie计算方向场。
[0012] 步骤6、通过方向场的线积分卷积，得到提示笔刷大小和走向的流线。
[0013]步骤7、将边界线条、数字编号和提示流线合成到一起，得到最终的数字油画画布。[0014]进一步，所述步骤1中，对于处在位置（x，y)的像素，记其颜色为（R，G，B)，用位 置-颜色空间的五维向量$ = 构建一个初始样本，其中a是将调整空间 尺度的归一化量，本发明中取255/图像长度。那么，Mean-Shift分割可以表示为一个寻找 中心位置的迭代过程：
【权利要求】
1. 基于图像的数字油画画布的自动生成方法，包括以下步骤： 步骤1、用Mean-Shift图像分割算法对输入图像I进行分割，并用每个区域的平均颜色 代替原有颜色，得到分割图像Ir。 步骤2、将Ir的区域边界通过统计滤波进行边界光顺，得到边界图像Ie。 步骤3、采用K均值算法将Ir的每个区域的平均颜色进行聚类分析。 步骤4、根据Ir中每个区域的轮廓，自适应调整与其对应的颜色编号的大小和位置。 步骤5、根据输入图像I和边界图像Ie计算方向场。 步骤6、通过方向场的线积分卷积，得到提示笔刷大小和走向的流线。 步骤7、将边界线条、数字编号和提示流线合成到一起，得到最终的数字油画画布。
2. 如权利要求1所述的基于图像的数字油画画布的自动生成方法，其特征在于：所述 步骤1中，对于处在位置（X，y)的像素，记其颜色为（R，G，B)，用位置-颜色空间的五维向 量殳=(〇^，〇|/,/?.0./^构建一个初始样本,其中〇是将调整空间尺度的归一化量,本发 明中取255/图像长度。那么，Mean-Shift分割可以表示为一个寻找中心位置的迭代过程：
其中，k表示当前迭代的次数，ffc表示第k次迭代的均值位置，nk是该局部邻域内初始 样本的数量，每次迭代都是对所有满足S A局部邻域条件的初始样本求平均，S用 于控制局部邻域的大小。该迭代是一个根据初始样本的分布不断移动局部邻域中心的过 程，当该中心位置几乎不移动时终止迭代。此外，Mean-Shift算法从每一个初始样本出发 分别进行上述迭代。迭代终止后，将终止中心位置较近的初始样本归为一类，从而完成图像 的分割，最后用每个区域的平均颜色IcJ取代原有的像素颜色。
3. 如权利要求1所述的基于图像的数字油画画布的自动生成方法，其特征在于：所述 步骤2中，对于当前像素 p，记其所在的区域编号为λ p，统计其局部邻域内所有像素的区域 编号I |p-q| I彡s}，并将当前像素 p的区域编号赋值为出现次数最多的区域编号入，。 该滤波保持了图像分割整体结果的同时，有效去除了边界的锯齿。其中，参数s是控制边界 光滑程度的一个整数，s越大，边界线条越光滑。
4. 如权利要求1所述的基于图像的数字油画画布的自动生成方法，其特征在于：所述 步骤3中，将每个区域的平均颜色{cj (I < j < N，N为区域总数）作为输入，采用K均值算 法进行聚类分析。它是一个迭代求解的过程：在第k次迭代中，计算所有&与代表颜色￠:? 的距离Ilci - €|丨|，并将&归类到距离最近的代表颜色所在的类中。接着，根据上述归 类，重新计算每一个代表颜色Cfl的值，即对归类到Cf丨的所有颜色求平均，并进入下一次迭 代，直至颜色误差变化小于某一阈值或迭代次数超过一定的值。需要指出的是，为了提高 算法稳定性和收敛速度，选取面积最大的K个区域的平均颜色作为代表颜色的初始值Cf。
5. 如权利要求1所述的基于图像的数字油画画布的自动生成方法，其特征在于：步骤 4所述的自适应调整数字的大小和位置，保证数字编号在最终数字油画画布中清晰可见，所 述的步骤4由以下步骤组成： (4. 1)我们沿着该中心位置的上、下、左、右四个方向进行搜索，假设在其中一个方向首 先遇到该区域内的像素 P1，继续沿着该方向搜索，直至找到离开该区域的位置或图像边界 P2，那么我们将数字的中心位置设为Ok = (Pl+P2)/2。 (4. 2)为保证数字显示的可见性，数字在区域内部应尽可能的大，即不妨记最大允许的 大小为R。然而，位置在区域边界附近的数字，由于不能超出区域的范围，显示区域却非常受 限。为此，可以采用类似于Mean-Shift分割的迭代方法，逐步移动中心位置0 k，直至不与边 界重叠或迭代次数超过10次为止：
该迭代通过计算局部邻域内区域编号相同的像素的平均位置，使得中心位置Ok逐渐远 离区域边界。其中，表示以Ok为中心，以R为半径的邻域范围，而nk是该邻域 范围内的编号相同的像素的数量。 (4.3)对于面积很小或形状细长的区域，可能放不下大小为R的数字符号。为此，我们 只能适当缩小数字符号的大小，以保证其显示在区域内部。
6. 如权利要求1所述的基于图像的数字油画画布的自动生成方法，其特征在于：所述 步骤5中，首先，计算方向场。假设输入图像为I，统计滤波后的图像为F，分别计算它们的 Sobel梯度:和丨VF，然后构建张量场VJt ·▽/ + VFt · VF，高斯模糊该张量场后， 分析其特征值，对应的特征向量即为梯度和切向方向。
7. 如权利要求1所述的基于图像的数字油画画布的自动生成方法，其特征在于：步骤 6所述的以一种合理的方式，为数字油画画布添加提示笔刷走向的提示线，所述的步骤6由 以下步骤组成： (6. 1)我们在图像空间撒1000个随机点，逐一从这些点出发进行卷积。生成切向流线： Xk+1 一 Xk+t (Xk) 其中，t(xk)是在Xk处的单位切向量。由于Xk往往不在整数像素位置上，我们采用双线 性插值的方法计算t (Xk)。每一次卷积，我们都可以获得表示为多边形的一条流线{xk}。 (6.2) 假设笔刷的长度和宽度分别是t和ω，当卷积的长度超过t，或者遇到区域边 界、数字符号或已经卷积过的位置时，卷积终止。 (6.3) 此外，每次卷积后，都需要在图像空间作相应的记录，即以流线上的点{xk}为圆 心、以ω为半径的圆逐一进行标记，以保证流线之间的距离大于ω。需要强调的是，笔刷的 长度和宽度与所处区域的面积有关，即I = ω = ω_*β，其中1_和（〇_分 别是最大允许的笔刷长度和宽度，j s取2. 5, A为归一化的区域面积。
【文档编号】G06T7/60GK104392433SQ201410616110
【公开日】2015年3月4日 申请日期:2014年11月5日 优先权日:2014年11月5日
【发明者】陈佳舟, 何栋, 秦绪佳 申请人:浙江工业大学