技术领域
[0001] 本
发明涉及一种计算机辅助乱针绣制作方法，属于计算机
图像处理和
计算机图形学技术领域。
背景技术
[0002] 乱针绣是把画理和绣理结合在一起创造出来的一种新的刺绣艺术，她是当代刺绣艺术家杨守玉女士创始于二十世纪三十年代，她的最大特点就是打破了几千年来传统刺绣“密接其针，排比其线”的平面刺绣，而是让作者自由地表达自己的思想感情，用纵横交叉、长短不一、疏密重叠、灵活多变的针法刺绣出来的立体画面，实现“画理”与“绣理”的有机结合。乱针绣的制作者多以现实生活为依据，选题广泛，视觉空间意识强，视觉与心理像契合，用色明快，语言丰富，利用照相术，强调绘画性，将具象写实与装饰相结合。目前江苏省已将其列入首批非物质文化遗产，中国政府也已将其列入世界非物质文化遗产
申请计划。然而，乱针绣因其技术难度相对较大，制作周期相对较长，故从事乱针绣艺术工作的人数相对较少，成功的乱针绣艺术精品更是凤毛麟
角，难得一见。精美的乱针绣作品成本昂贵，即使是普通的乱针绣作品在市场上也难觅其踪。因此，采用计算机技术直接从油画、摄影等数字图像生成乱针绣作品不仅可以缩短乱针绣绣品的制作周期，降低制作成本，减小绣者的劳动强度，同时对乱针绣艺术的传承和保护也具有重要的意义。
[0003] 已有诸如文献1李俊,张华,王崇骏：智能化刺绣CAD系统中的工作流技术研究.计算机科学,2005,32(3):97-100.所述的日本的Tajima、德国的ZSK、美国的Wilcom等刺绣打版系统，往往需要人工交互编辑每个针迹的坐标序列，且生成的针迹排布规整，只能处理简单的图案，不满足乱针绣工艺的要求。由于乱针绣的特殊性，国内外很少有计算机辅助乱针绣制作的研究，主要的研究包括中国
专利《一种计算机辅助乱针绣制作方法》，公开号：CN101859335A，该专利以彩色图像作为输入，采用乱针绣绣线
颜色选择和针迹参数计算生成乱针绣绣品的针迹序列，并将其输入刺绣机得到乱针绣绣品。然而，实际的乱针绣制作过程中不仅要考虑绣线颜色的排布，还要考虑绣线方向的排布、长度的排布等多种因素；在针迹参数计算过程中，其采用试探法仅根据图像局部区域的颜色差异、梯度方向确定绣线方向而没有考虑到图像区域的纹理性、重要性等特征对针迹排布的影响。另外一个中国专利《一种采用遗传
算法的计算机辅助乱针绣制作方法》，申请号201310211288.2：该专利同样以彩色图像作为输入，其定义了相应的针法模型参数，并通过
遗传算法优化适应度函数得到针法模型中的高层参数，并利用高层参数得到最终的针迹序列。然而，其定义的针法模型并没有考虑相邻针迹之间的关系，而本专利通过建立笔画邻域图建立了针迹与针迹之间的拓扑关系，并通过5个扩散因子有效的控制针迹的方向、长度及颜色的排布策略；另外，该专利采用遗传算法进行的
迭代过程，运行时间比较长。
[0004] 乱针绣的制作过程，主要是将输入的彩色图像转化为乱针绣制作所需要的针迹序列的过程，该过程属于基于笔画的绘制过程。基于笔画的绘制过程主要考虑笔画的排布问题，如文献2Hertzmann A.Painterly rendering with curved brush strokes of thmultiple sizes,In Proceedings of the25 Annual Conference on Computer Graphics and Interactive Techniques,SIGGRAPH’98,Orlando,Florida,USA,1998;453-460. 在中间结果图像与参考图像中颜色差值大于一定
阈值的区域防止笔画的起始点。文献
3Litwinowicz P.Proceessing images and video for an impressionist effect,In th
Proceedings of the24 Annual Conference on Computer Graphics and Interactive Techniques,SIGGRAPH’97,Los Angeles,California,USA,1997;407-414.在交替
像素位置防止笔画以保证输出结果中没有空洞。文献4Mao X,Nagasaka Y,Imamiya A.Automatic generation of pencil drawing from2D image using line integral convolution,In CAD/Graphics,volume9,Kunming,China,2001;240-248.通过生成一定
密度的噪音确定笔画的密度，噪音密度由图像
色调影响，利用低通
滤波器生成铅笔笔画。然而，乱针绣主要利用针线排布方向的多变，针线的长短参差以及用色的灵活多变表现物象，而这些方法均不能直接应用于乱针绣制作过程中。
发明内容
[0005] 发明目的：本发明所要解决的技术问题是针对
现有技术的不足，提供一种基于反应扩散过程的计算机辅助乱针绣制作方法。
[0006] 技术方案：本发明公开了一种基于反应扩散过程的计算机辅助乱针绣制作方法，其特征在于，包括以下步骤：
[0007] 步骤一，图像特征提取：以彩色图像作为输入，提取包括对象区域，区域矢量场和图像显著性在内的图像特征；
[0008] 步骤二，建立乱针绣针法模型并计算高层参数：建立包含单根绣线、单个交叉针和交叉针排布的低、中、高三层参数表示的乱针绣针法模型，并根据图像特征计算参数表示的乱针绣针法模型中控制交叉针排布的高层参数；
[0009] 步骤三，乱针绣针迹序列反应扩散过程：对每个对象区域，依据步骤二中获取的高层参数生成初始的针迹序列集合，建立笔画邻域图，并结合步骤一中提取的图像特征，对初始的针迹序列集合中每个交叉针的方向、长度、颜色色调、颜色
饱和度及颜色
亮度分别进行五个反应扩散过程，然后为每个针迹从绣线库中选择对应的绣线颜色，得到更新后的乱针绣针迹序列；
[0010] 步骤四，将步骤三中得到的乱针绣针迹序列坐标转换为电脑刺绣机使用的以毫米为单位的坐标格式，将其输入电脑刺绣机，从而得到实际的乱针绣绣品。
[0011] 本发明步骤一中包括以下步骤：
[0012] 步骤11：对输入的彩色图像，采用文献5Zeng K，Zhao M T，Xiong C M，et al.From image parsing to painterly rendering[J].ACM Transactions on Graphics，2009，29(1)：2:1-2:11中所述image-parsing
图像分割方法将图像分割成不同的对象区域；
[0013] 步骤12：对于分割后的对象区域，采用文献6Kang H，Lee S Y，Chui C thK.Coherent line drawing[C]//Proceedings of the5 International Symposium on Non-photorealistic Animation and Rendering.New York：ACM Press，2007：43-50中所述Edge Tangent Flow 方法提取区域矢量场；该特征为对象区域的每个像素点赋予一个值，该值的范围为[0,π)；
[0014] 步骤13：对输入的彩色图像，采用文献7Cheng M M,Zhang G X,Mitra N J,et al.Global contrast based salient region detection[C]//Computer Vision and Pattern Recognition(CVPR),2011IEEE Conference on.IEEE,2011:409-416中所述基于区域
对比度的显著性提取方法提取图像显著性特征；该特征就是一幅灰度图像，亮度越大，显著性越高，其值范围为[0,255]。
[0015] 本发明步骤二中包括以下步骤：
[0016] 建立乱针绣针法模型：所述乱针绣针法模型包括低、中、高三个层次的参数表示，低层为单根绣线，中层为单个交叉针、高层为交叉针排布；其中，低层参数集合为η={ηk:k=1～2N}，其中， 表示单根绣线k对应的参数集，为绣线起点坐标、 为绣线终点坐标，rk、gk、bk分别为绣线在
RGB颜色空间的R、G、B通道值，N表示交叉针总数，k为绣线序号；中层参数集合为χ={χi:i=1～N}，其中，χi=(μi,ξi,θi,βi,hi,si,vi)为单个交叉针i对应的参数集，μi为交叉针中心点，ξi为交叉针长度，θi为交叉针方向，βi为交叉针夹角，hi、si、vi分别为在HSV颜色空间中的交叉针色调、交叉针饱和度、交叉针亮度，i为交叉针序号；高层参数集合为Ω={ρ,λθ,λξ,λh,λs,λv}，其中，ρ为交叉针排布的密度，λθ、λξ、λh、λs、λv分别为交叉针排布的方向扩散因子、长度扩散因子、颜色色调扩散因子、颜色饱和度扩散因子、颜色亮度扩散因子；
[0017] 计算高层参数：高层参数集合Ω中ρ的取值范围为(0.0,1.0]，ρ取0时为最稀疏，没有针迹，ρ取1.0时为最稠密，实际上机过程中针迹的排布在保持较大
覆盖率的前提下同时保证不会排布过密，本发明经过计算取值0.80～0.98，优选0.89，λθ由步骤12中提取的区域矢量场计算区域矢量场方差值进一步计算得到，λξ由步骤13中提取的图像显著性特征方差值进一步计算得到，λh、λs、λv分别由对象区域内所有像素颜色色调方差、饱和度方差、亮度方差进一步计算得到。
[0018] 本发明步骤二中计算高层参数，λθ、λξ、λh、λs、λv具体的计算公式为：
[0019] λθ、λξ、λh、λs、λv均可取正值和负值，取正值会分别使得对象区域内的交叉针的方向、长度、颜色色调、颜色饱和度、颜色亮度这五个参数值相似，而且，正值越大，相似性越大，取负值则会使其值差异比较大，且负值越小，差异越大。
[0020] λθ=0.5-σθ，其中， 为区域矢量场方差，通过对大量对象区域的实验统计发现，该方差的分布范围大致为(0.1,0.92)，方差小于0.5时区域矢量场分布比较均匀，λθ取值为正，大于0.5时则比较乱，λθ取值为负，因此公式中以
0.5为
阀值。obj为对象区域的像素点集合，Num(obj)为对象区域所有像素点的个数，pixel为对象区域中的任意一个像素点，θpixel为对象区域内pixel像素点对应的区域矢量场值，为对象区域内所有像素点对应的区域矢量场值的平均值，计算公式为
[0021]
[0022] λξ=50-σξ，其中， 为区域显著性值方差，通过实验统计分析发现，该方差的分布范围大致为(0,153)，该方差小于50时，对象区域的显著性比较接近，λξ取值为正，反之则为负，Salpixel为对象区域内pixel像素点对应的显著性值， 为对象区域内所有像素点对应的显著性值的平均值，计算公式为
[0023]
[0024] λh=0.8-σh，其中， 为区域颜色色调方差，通过实验统计分析发现，该方差的分布范围大致为(0,3.05)，该方差小于0.8时，对象区域的色调分布比较集中，λh取值为正，反之则为负，hpixel为对象区域内pixel像素点在原图中的颜色在HSV颜色空间的H通道值， 为对象区域内所有像素点颜色在HSV颜色空间的H通道值的平均值，计算公式为
[0025] λs=0.5-σs，其中， 为区域颜色饱和度方差，通过实验统计分析发现，该方差的分布范围大致为(0,1.05)，该方差小于0.5时，对象区域的饱和度分布比较集中，λh取值为正，反之则为负，spixel为对象区域内pixel像素点在原图中的颜色在HSV颜色空间的S通道值，为对象区域内所有像素点颜色在HSV颜色空间的S通道值的平均值，计算公式为
[0026] λv=0.5-σv，其中， 为区域颜色亮度方差，通过实验统计分析发现，该方差的分布范围大致为(0,1.14)，该方差小于0.5时，对象区域的亮度分布比较集中，λh取值为正，反之则为负，vpixel为对象区域内pixel像素点在原图中的颜色在HSV颜色空间的V通道值， 为对象区域内所有像素点颜色在HSV颜色空间的V通道值的平均值，计算公式为
[0027] 本发明步骤三中包括以下步骤：
[0028] 步骤31：对每个对象区域，根据步骤二中获取的交叉针排布的密度ρ确定
采样点个数N，采用均匀泊松盘采样方法选取N个像素点位置作为初始的针迹序列集合Λ={Λi:i=1～N}，其中，N为交叉针总数，与采样点个数相同，Λi={χi,τi}为每个交叉针χi与其邻域交叉针队列τi的组合，τi初始为空，按照序号i分别设置交叉针χi的交叉针中心点μi为第i个采样点的坐标值，并分别为交叉针χi的交叉针长度ξi、交叉针方向θi、交叉针夹角βi、交叉针色调hi、交叉针饱和度si、交叉针亮度vi取随机值作为初始值，ξi的取值范围为[4,25]，θi的取值范围为[0,π)，βi的取值范围为 hi的取值范围为[0,360)，si和vi的取值范围为[0.0,1.0]；
[0029] 步骤32：对针迹序列集合Λ中每个组合Λi执行步骤321～步骤322；
[0030] 步骤321：依据步骤12中提取的区域矢量场获取组合Λi的交叉针中心点μi处的矢量场方向，以交叉针中心点μi为原点，获取的矢量场方向及其垂直方向为X轴和Y轴，建立局部笛卡尔
坐标系，并划分为四个象限；
[0031] 步骤322：依据交叉针中心点μi及四个象限，分别在每个象限内选择一个与交叉针中心点μi距离最小的交叉针，并将该交叉针序号加入到邻域交叉针队列τi中，邻域交叉针队列τi中包含四个象限中的交叉针的序号；
[0032] 步骤33：由步骤31得到针迹序列集合Λ中每个组合Λi的交叉针χi，由步骤32得到每个组合Λi的邻域交叉针队列τi，将二者组合成完整的针迹序列集合Λ；
[0033] 步骤34：对针迹序列集合Λ执行步骤341～步骤345；
[0034] 步骤341：执行方向反应扩散过程；
[0035] 步骤342：执行长度反应扩散过程；
[0036] 步骤343：执行颜色色调反应扩散过程；
[0037] 步骤344：执行颜色饱和度反应扩散过程；
[0038] 步骤345：执行颜色亮度反应扩散过程；
[0039] 步骤35：对针迹序列集合Λ中的每个针迹从绣线库中选择对应的绣线颜色，得到更新后的乱针绣针迹序列。
[0040] 本发明步骤341中所述方向反应扩散过程，具体步骤为：
[0041] 步骤3411：设定初始迭代次数K1=1；
[0042] 步骤3412：对针迹序列集合Λ中的每个组合Λi，执行步骤34121～步骤34124：
[0043] 步骤34121：获取组合Λi的交叉针方向θi，获取步骤12提取的区域矢量场中位*于组合Λi的交叉针中心点μi的矢量场方向θi ；
[0044] 步骤34122：计算组合Λi的交叉针中心点μi与邻域交叉针队列τi中每个交叉针的交叉针中心点的距离，并将该距离的反比设置为邻域交叉针队列τi中每个交叉针与序号为i的交叉针的权重；
[0045] 步骤34123：计算方向偏量值 其中，τi为组合Λi的邻域交叉针队列，j为τi中每个交叉针的序号，ωij为序号为j的交叉针与序号为i的交叉针的权重，θj为序号为j的交叉针的交叉针方向，εθ=rand(0,0.05)为方向随机因子；
[0046] 步骤34124：计算组合Λi新的交叉针方向值θi=θi+Δθ；
[0047] 步骤3413：根据针迹序列集合Λ中每个组合Λi更新后的交叉针方向值，计算方向
马尔科夫
能量方程值
[0048] 步骤3414：如果当前第K1次迭代的能量值与第K1-1次迭代的能量值差值小于阈值0.5（通过实验统计分析得到，对方向扩散过程来说，该能量值差值阈值取0.5时结果趋于稳定），则迭代结束，否则，转步骤3412继续执行，且迭代次数K1加1。
[0049] 本发明步骤342中所述长度反应扩散过程，具体步骤为：
[0050] 步骤3421：设定初始迭代次数K2=1；
[0051] 步骤3422：对针迹序列集合Λ中的每个组合Λi，执行步骤34221～步骤34224：
[0052] 步骤34221：获取组合Λi的交叉针长度ξi，依据步骤13提取的图像显著性特征计算长度ξi*，计算公式为 其中，Isal(μi)为组合Λi的交叉针中心点μi处的图像显著性值；
[0053] 步骤34222：计算组合Λi的交叉针中心点μi与邻域交叉针队列τi中每个交叉针的交叉针中心点的距离，并将该距离的反比设置为邻域交叉针队列τi中每个交叉针与序号为i的交叉针的权重；
[0054] 步骤34223：计算长度偏量值 其中，ξj为序号为j的交叉针的交叉针长度，εξ=rand(0,0.2)为长度随机因子；
[0055] 步骤34224：计算组合Λi新的交叉针长度值ξi=ξi+Δξ；
[0056] 步骤3423：根据针迹序列集合Λ中每个组合Λi更新后的交叉针长度值，计算长度马尔科夫能量方程值
[0057] 步骤3424：如果当前第K2次迭代的能量值与第K2-1次迭代的能量值差值小于阈值1.0（通过实验统计分析得到，对长度扩散过程来说，该能量值差值阈值取1.0时结果趋于稳定），则迭代结束，否则，转步骤3422继续执行，且迭代次数K2加1。
[0058] 本发明步骤343中所述颜色色调反应扩散过程，具体步骤为：
[0059] 步骤3431：设定初始迭代次数K3=1；
[0060] 步骤3432：对针迹序列集合Λ中的每个组合Λi，执行步骤34321～步骤34324：
[0061] 步骤34321：获取组合Λi的交叉针色调hi，获取输入彩色图像中位于组合Λi的*交叉针中心点μi的像素在HSV颜色空间的H通道值hi ；
[0062] 步骤34322：计算组合Λi的交叉针中心点μi与邻域交叉针队列τi中每个交叉针的交叉针中心点的距离，并将该距离的反比设置为邻域交叉针队列τi中每个交叉针与序号为i的交叉针的权重；
[0063] 步骤34323：计算色调偏量值 其中，hj为序号为j的交叉针的交叉针色调，εh=rand(0,1.0)为色调随机因子；
[0064] 步骤34324：计算组合Λi新的交叉针色调值hi=hi+Δh；
[0065] 步骤3433：根据针迹序列集合Λ中每个组合Λi更新后的交叉针色调值，计算色调马尔科夫能量方程值
[0066] 步骤3434：如果当前第K3次迭代的能量值与第K3-1次迭代的能量值差值小于阈值0.2（通过实验统计分析得到，对长度扩散过程来说，该能量值差值阈值取0.2时结果趋于稳定），则迭代结束，否则，转步骤3432继续执行，且迭代次数K3加1。
[0067] 本发明步骤344中所述颜色饱和度反应扩散过程，具体步骤为：
[0068] 步骤3441：设定初始迭代次数K4=1；
[0069] 步骤3442：对针迹序列集合Λ中的每个组合Λi，执行步骤34421～步骤34424：
[0070] 步骤34421：获取组合Λi的交叉针饱和度si，获取输入彩色图像中位于组合Λi*的交叉针中心点μi的像素在HSV颜色空间的S通道值si ；
[0071] 步骤34422：计算组合Λi的交叉针中心点μi与邻域交叉针队列τi中每个交叉针的交叉针中心点的距离，并将该距离的反比设置为邻域交叉针队列τi中每个交叉针与序号为i的交叉针的权重；
[0072] 步骤34423：计算饱和度偏量值 其中，sj为序号为j的交叉针的交叉针饱和度，εs=rand(0,0.001)为饱和度随机因子；
[0073] 步骤34424：计算组合Λi新的交叉针饱和度值si=si+Δs；
[0074] 步骤3443：根据针迹序列集合Λ中每个组合Λi更新后的交叉针饱和度值，计算饱和度马尔科夫能量方程值
[0075] 步骤3444：如果当前第K4次迭代的能量值与第K4-1次迭代的能量值差值小于阈值0.05（通过实验统计分析得到，对长度扩散过程来说，该能量值差值阈值取0.05时结果趋于稳定），则迭代结束，否则，转步骤3442继续执行，且迭代次数K4加1。
[0076] 本发明步骤345中所述颜色亮度反应扩散过程，具体步骤为：
[0077] 步骤3451：设定初始迭代次数K5=1；
[0078] 步骤3452：对针迹序列集合Λ中的每个组合Λi，执行步骤34521～步骤34524：
[0079] 步骤34521：获取组合Λi的交叉针亮度vi，获取输入彩色图像中位于组合Λi的*交叉针中心点μi的像素在HSV颜色空间的V通道值vi ；
[0080] 步骤34522：计算组合Λi的交叉针中心点μi与邻域交叉针队列τi中每个交叉针的交叉针中心点的距离，并将该距离的反比设置为邻域交叉针队列τi中每个交叉针与序号为i的交叉针的权重；
[0081] 步骤34523：计算亮度偏量值 其中，vj为序号为j的交叉针的交叉针亮度，εv=rand(0,0.001)为亮度随机因子；
[0082] 步骤34524：计算组合Λi新的交叉针亮度值vi=vi+Δv；
[0083] 步骤3453：根据针迹序列集合Λ中每个组合Λi更新后的交叉针亮度值，计算亮度马尔科夫能量方程值
[0084] 步骤3454：如果当前第K5次迭代的能量值与第K5-1次迭代的能量值差值小于阈值0.05（通过实验统计分析得到，对长度扩散过程来说，该能量值差值阈值取0.05时结果趋于稳定），则迭代结束，否则，转步骤3452继续执行，且迭代次数K5加1。
[0085] 本发明步骤四将以像素为单位的针迹坐标乘以绣线宽度即可得到以毫米为单位的电脑刺绣机坐标，最后将其输入电脑刺绣机得到乱针绣绣品。
[0086] 本发明的特色在于根据乱针绣的特点对乱针绣绣制过程中线条方向、长度、颜色的排布问题进行了处理，定义了针法模型，尤其是针法模型高层参数中的五个扩散因子，并通过对针迹序列建立笔画邻域图构造针迹与针迹之间的拓扑关系，并在此
基础上利用五个反映扩散过程有效的控制针迹的排布方向、长短及颜色搭配。
[0087] 有益效果：本发明具有以下优点：1、允许用户任意选取彩色图像作为输入，得到其相应的乱针绣绣品，减少了用户的工作量，提高了乱针绣制作效率；2、定义了一种乱针绣针法模型，结合了乱针绣技法的特点并将其绣法模式化；3、本发明通过建立笔画邻域图，以五个扩散因子为参数并通过反应扩散过程自动控制绣制过程中针迹的方向、长度、颜色等因素，使得绣制过程更加灵活多变。综上所述，本发明以彩色图像为输入并采用计算机辅助手段自动完成乱针绣绣品的制作过程，有助于乱针绣艺术的传承和发展。
附图说明
[0088] 下面结合附图和具体实施方式对本发明做更进一步的具体说明，本发明的上述和/或其他方面的优点将会变得更加清楚。
[0090] 图2为本发明的对象区域提取结果。
[0091] 图3为本发明的区域矢量场提取结果。
[0092] 图4为本发明的图像显著性提取结果。
[0093] 图5为本发明的笔画邻域图。
[0094] 图6为本发明最终乱针绣绣制效果图。具体实施方式：
[0095] 本发明公开了一种基于反应扩散过程的计算机辅助乱针绣制作方法，如图1所示，下面分别介绍各
实施例部分的主要流程：
[0096] 1、图像特征提取
[0097] 本发明步骤一中提取的图像特征包括对象区域、区域矢量场、图像显著性三种，具体描述如下。
[0098] 1.1对象区域提取
[0099] 在乱针绣中，不同的对象其绣法也不同，所以需要将输入的彩色图像分割成不同的对象区域，并对不同的对象区域采用不同的针法参数以得到最终的绣制结果。如本发明步骤11所述，本发明采用文献5Zeng K，Zhao M T，Xiong C M，et al.From image parsing to painterly rendering[J].ACM Transactions on Graphics，2009，29(1)：2:1-2:11中所述image-parsing图像分割方法将图像分割成M不同的对象区域{R1,R2,...,RM}，每个对象区域Rm为一组像素点集合，m为对象区域序号，集合中任意像素点至少与该集合中另一像素点相邻，分割结果如图2所示，对附图2而言，实施例中的对象区域为{R1,R2}。
[0100] 1.2区域矢量场提取
[0101] 在实际的乱针绣绣制过程中，走线的方向要考虑不同对象区域的纹理走向，区域的矢量场可以用来指导走线的方向。本发明采用文献6Kang H，Lee S Y，Chui C thK.Coherent line drawing[C]//Proceedings of the5 International Symposium on Non-photorealistic Animation and Rendering.New York：ACM Press，2007：
43-50中所述Edge Tangent Flow方法为每个对象区域Rm提取相应的区域矢量场
其中，Θm为Rm的区域矢量场，g为Rm中每个像素点的序号，
Num(Rm)为Rm中像素点的个数， 为Rm中第g个像素点的矢量场值，该值范围为[0,π)，矢量场提取结果如图3所示，对附图3而言，实施例中包含附图2中2个对象区域的矢量场Θ1和Θ2。
[0102] 1.3图像显著性提取
[0103] 在乱针绣中，绣制线条的长短与对象的重要性有关系，例如简单的背景区域重要性比较低，一般采用较长的线进行绣制，而比较重要的前景对象为了刻画更加细腻，绣制时线条长度较短为好，因此，可以依据提取图像显著性特征指导绣制线条长度信息。本发明采用文献7Cheng M M,Zhang G X,Mitra N J,et al.Global contrast based salient region detection[C]//Computer Vision and Pattern Recognition(CVPR),2011IEEE Conference on.IEEE,2011:409-416中所述基于区域对比度的显著性提取方法提取图像显著性特征，该特征为包含所有对象区域每个像素点的显著性值的集合，某个像素点的显著性值越高，则该像素点的重要性越高，其显著性值的范围为[0,255]，显著性提取结果如图4所示，对附图4而言，越亮的像素点位置显著性值越高，反之，越暗责显著性值越低。
[0104] 2、乱针绣针法模型
[0105] 本发明步骤二中定义了一种乱针绣针法模型的三层参数集合，并且包括乱针绣针法模型中高层参数的获取方法。
[0106] 2.1针法模型建立
[0107] 乱针绣以交叉针作为基本的单元，针法模型包括低、中、高三个层次的参数集合，低层为单根绣线，中层为单个交叉针、高层为交叉针排布。
[0108] 低层参数集合为η={ηk:k=1～2N}，其中， 表示单根绣线k对应的参数集， 为绣线起点坐标、 为绣线终点坐标，rk、gk、bk分别为绣线在RGB颜色空间的R、G、B通道值，N表示交叉针总数，k为绣线序号；
[0109] 中层参数集合为χ={χi:i=1～N}，其中，χi=(μi,ξi,θi,βi,hi,si,vi)为单个交叉针i对应的参数集，μi为交叉针中心点，ξi为交叉针长度，θi为交叉针方向，βi为交叉针夹角，hi、si、vi分别为在HSV颜色空间中的交叉针色调、交叉针饱和度、交叉针亮度，i为交叉针序号；
[0110] 高层参数集合为Ω={ρ,λθ,λξ,λh,λs,λv}，其中，ρ为交叉针排布的密度，λθ、λξ、λh、λs、λv分别为交叉针排布的方向扩散因子、长度扩散因子、颜色色调扩散因子、颜色饱和度扩散因子、颜色亮度扩散因子；
[0111] 2.1高层参数计算
[0112] 高层参数集合Ω中每个高层参数的获取如下：
[0113] ρ根据实验统计值取0.89为优；
[0114] λθ=0.5-σθ，其中， 为区域矢量场方差，obj为对象区域的像素点集合，Num(obj)为对象区域所有像素点的个数，pixel为对象区域中的任意一个像素点，θpixel为对象区域内pixel像素点对应的区域矢量场值， 为对象区域内所有像素点对应的区域矢量场值的平均值，计算公式为
[0115] λξ=50-σξ，其中， 为区域显著性值方差，Salpixel为对象区域内pixel像素点对应的显著性值， 为对象区域内所有像素点对应的显著性值的平均值，计算公式为
[0116] λh=0.8-σh，其中， 为区域颜色色调方差，hpixel为对象区域内pixel像素点在原图中的颜色在HSV颜色空间的H通道值， 为对象区域内所有像素点颜色在HSV颜色空间的H通道值的平均值，计算公式为
[0117]
[0118] λs=0.5-σs，其中， 为区域颜色饱和度方差，spixel为对象区域内pixel像素点在原图中的颜色在HSV颜色空间的S通道值，为对象区域内所有像素点颜色在HSV颜色空间的S通道值的平均值，计算公式为
[0119]
[0120] λv=0.5-σv，其中， 为区域颜色亮度方差，vpixel为对象区域内pixel像素点在原图中的颜色在HSV颜色空间的V通道值，为对象区域内所有像素点颜色在HSV颜色空间的V通道值的平均值，计算公式为
[0121]
[0122] 3、反应扩散过程
[0123] 结合步骤一中提取的对象区域、区域矢量场、图像显著性特征及步骤二中建立的针法模型及计算的高层参数Ω，本发明步骤三中对每个对象区域Rm执行以下步骤：生成初始的针迹序列集合，建立笔画邻域图，分别对初始的针迹序列集合中每个交叉针的方向、长度、颜色色调、颜色饱和度及颜色亮度五个参数进行五个反应扩散过程，然后为针迹序列集合中的每根绣线从绣线库中选择最相近的颜色，得到更新后的乱针绣针迹序列。
[0124] 3.1初始针迹序列生成
[0125] 步骤311：对对象区域Rm，根据ρ确定采样像素点个数 采用文献8Gamito M N,Maddock S C.Accurate multidimensional Poisson-disk sampling[J].ACM Transactions on Graphics(TOG),2009,29(1):8中的方法随机选取N个像素点，可知N
[0126] 步骤312：在这个N个像素点处放置交叉针组成初始的针迹序列集合Λ={Λi:i=1～N}，其中，Λi={χi,τi}为每个交叉针χi与其邻域交叉针队列τi的组合，τi初始为空，i为采样像素点的序号，按照序号i分别设置χi中的μi值为第i个采样像素点的坐标，并分别为χi中的ξi、θi、βi、hi、si、vi取随机值作为初始值，ξi的取值范围为[4,25]，θi的取值范围为[0,π)，βi的取值范围为 hi的取值范围为[0,360)，si和vi的取值范围为[0.0,1.0]；
[0127] 3.2笔画邻域图建立
[0128] 步骤321：将当前对象区域Rm的针迹序列集合Λ中每个Λi的邻域交叉针队列τi设置为空，对每个Λi执行步骤322～步骤324；
[0129] 步骤322：从Θm中得到像素点位置μi=(xi,yi)的方向θi，以μi为原点，以θi的方向为X轴，θi+90°的方向为Y轴建立局部笛卡尔坐标系，分别计算得到4个象限的角度范围：[θi,θi+90°)，[θi+90°，θi+180°)，[θi+180°，θi+270°)，[θi+270°，θi+360°)，分别为每个象限建立4个队列list1，list2，list3，list4，初始每个队列为空；
[0130] 步骤323：对Λ中除Λi以外的所有交叉针的位置点μj=(xj,yj)，其中，j=1～N,j≠i，依次计算3个变量值：dx=xj-xi，dy=yj-yi， 对每个象限的角度范围[αmin,αmax)，其中，αmin和αmax分别为象限角度范围的下限值和上限值，首先计算v1x=cos(αmin)，v1y=sin(αmin)，v2x=cos(αmax)，v2y=sin(αmax)，然后计算以下3个 值
δ3=arccos((v1x*v2x+v1y*v2y))，如果δ1≤δ3且δ2≤δ3，则当前点μj=(xj,yj)属于[αmin,αmax)表示的象限，将μj=(xj,yj)所在的交叉针Λj的序号j加入到对应的队列中。
[0131] 步骤324：对每个队列listt，t=1,2,3,4，计算listt中每个序号对应的交叉针的位置点μj=(xj,yj)与μi=(xi,yi)的欧式距离选择距离最小的交叉针并将其序号加到τi中。
[0132] 建立的笔画邻域图如图5所示，对附图5而言，实施例中包含附图2中2个对象区域的2个笔画邻域图，为了更清晰地表现该笔画邻域图的拓扑关系，实施例中的2个笔画邻域图在步骤311中采用的参数ρ均为0.72，且将每个交叉针中心点与邻域交叉针队列中的交叉针中心点相连。
[0133] 3.3交叉针方向扩散
[0134] 步骤331：设定初始迭代次数K1=1；
[0135] 步骤332：对针迹序列集合Λ中的每个组合Λi，执行步骤3321～步骤3324：
[0136] 步骤3321：获取组合Λi的交叉针方向θi，获取步骤12提取的区域矢量场中位于*组合Λi的交叉针中心点μi的矢量场方向θi ；
[0137] 步骤3322：计算组合Λi的交叉针中心点μi与邻域交叉针队列τi中每个交叉针的交叉针中心点的距离，并将该距离的反比设置为邻域交叉针队列τi中每个交叉针与序号为i的交叉针的权重；
[0138] 步骤3323：计算方向偏量值 其中，τi为组合Λi的邻域交叉针队列，j为τi中每个交叉针的序号，ωij为序号为j的交叉针与序号为i的交叉针的权重，θj为序号为j的交叉针的交叉针方向，εθ=rand(0,0.05)为方向随机因子；
[0139] 步骤3324：计算组合Λi新的交叉针方向值θi=θi+Δθ；
[0140] 步骤333：根据针迹序列集合Λ中每个组合Λi更新后的交叉针方向值，计算方向马尔科夫能量方程值
[0141] 步骤334：如果当前第K1次迭代的能量值与第K1-1次迭代的能量值差值小于0.5，则迭代结束，否则，转步骤332继续执行，且迭代次数K1加1。
[0142] 3.4交叉针长度扩散
[0143] 步骤341：设定初始迭代次数K2=1；
[0144] 步骤342：对针迹序列集合Λ中的每个组合Λi，执行步骤3421～步骤3424：
[0145] 步骤3421：获取组合Λi的交叉针长度ξi，依据步骤13提取的图像显著性特征计*算长度ξi，计算公式为 其中，Isal(μi)为组合Λi的交叉针中心点μi
处的图像显著性值；
[0146] 步骤3422：计算组合Λi的交叉针中心点μi与邻域交叉针队列τi中每个交叉针的交叉针中心点的距离，并将该距离的反比设置为邻域交叉针队列τi中每个交叉针与序号为i的交叉针的权重；
[0147] 步骤3423：计算长度偏量值 其中，ξj为序号为j的交叉针的交叉针长度，εξ=rand(0,0.2)为长度随机因子；
[0148] 步骤3424：计算组合Λi新的交叉针长度值ξi=ξi+Δξ；
[0149] 步骤343：根据针迹序列集合Λ中每个组合Λi更新后的交叉针长度值，计算长度马尔科夫能量方程值
[0150] 步骤344：如果当前第K2次迭代的能量值与第K2-1次迭代的能量值差值小于1.0，则迭代结束，否则，转步骤342继续执行，且迭代次数K2加1。
[0151] 3.5颜色色调扩散
[0152] 步骤351：设定初始迭代次数K3=1；
[0153] 步骤352：对针迹序列集合Λ中的每个组合Λi，执行步骤3521～步骤3524：
[0154] 步骤3521：获取组合Λi的交叉针色调hi，获取输入彩色图像中位于组合Λi的交*叉针中心点μi的像素在HSV颜色空间的H通道值hi ；
[0155] 步骤3522：计算组合Λi的交叉针中心点μi与邻域交叉针队列τi中每个交叉针的交叉针中心点的距离，并将该距离的反比设置为邻域交叉针队列τi中每个交叉针与序号为i的交叉针的权重；
[0156] 步骤3523：计算色调偏量值 其中，hj为序号为j的交叉针的交叉针色调，εh=rand(0,1.0)为色调随机因子；
[0157] 步骤3524：计算组合Λi新的交叉针色调值hi=hi+Δh；
[0158] 步骤353：根据针迹序列集合Λ中每个组合Λi更新后的交叉针色调值，计算色调马尔科夫能量方程值
[0159] 步骤354：如果当前第K3次迭代的能量值与第K3-1次迭代的能量值差值小于0.2，则迭代结束，否则，转步骤352继续执行，且迭代次数K3加1。
[0160] 3.6颜色饱和度扩散
[0161] 步骤361：设定初始迭代次数K4=1；
[0162] 步骤362：对针迹序列集合Λ中的每个组合Λi，执行步骤3621～步骤3624：
[0163] 步骤3621：获取组合Λi的交叉针饱和度si，获取输入彩色图像中位于组合Λi的交叉针中心点μi的像素在HSV颜色空间的S通道值si*；
[0164] 步骤3622：计算组合Λi的交叉针中心点μi与邻域交叉针队列τi中每个交叉针的交叉针中心点的距离，并将该距离的反比设置为邻域交叉针队列τi中每个交叉针与序号为i的交叉针的权重；
[0165] 步骤3623：计算饱和度偏量值 其中，sj为序号为j的交叉针的交叉针饱和度，εs=rand(0,0.001)为饱和度随机因子；
[0166] 步骤3624：计算组合Λi新的交叉针饱和度值si=si+Δs；
[0167] 步骤363：根据针迹序列集合Λ中每个组合Λi更新后的交叉针饱和度值，计算饱和度马尔科夫能量方程值
[0168] 步骤364：如果当前第K4次迭代的能量值与第K4-1次迭代的能量值差值小于0.05，则迭代结束，否则，转步骤362继续执行，且迭代次数K4加1。
[0169] 3.7颜色亮度扩散
[0170] 步骤371：设定初始迭代次数K5=1；
[0171] 步骤372：对针迹序列集合Λ中的每个组合Λi，执行步骤3721～步骤3724：
[0172] 步骤3721：获取组合Λi的交叉针亮度vi，获取输入彩色图像中位于组合Λi的交*叉针中心点μi的像素在HSV颜色空间的V通道值vi ；
[0173] 步骤3722：计算组合Λi的交叉针中心点μi与邻域交叉针队列τi中每个交叉针的交叉针中心点的距离，并将该距离的反比设置为邻域交叉针队列τi中每个交叉针与序号为i的交叉针的权重；
[0174] 步骤3723：计算亮度偏量值 其中，vj为序号为j的交叉针的交叉针亮度，εv=rand(0,0.001)为亮度随机因子；
[0175] 步骤3724：计算组合Λi新的交叉针亮度值vi=vi+Δv；
[0176] 步骤373：根据针迹序列集合Λ中每个组合Λi更新后的交叉针亮度值，计算亮度马尔科夫能量方程值
[0177] 步骤374：如果当前第K5次迭代的能量值与第K5-1次迭代的能量值差值小于0.05，则迭代结束，否则，转步骤372继续执行，且迭代次数K5加1。
[0178] 3.8绣线颜色选择
[0179] 输入：绣线颜色库 γ为绣线颜色的序号，W为绣线颜色库中绣线颜色的总数， 分别为绣线颜色在RGB颜色空间的R、G、B通
道值，对象区域Rm的针迹序列集合Λ={Λi:i=1～N}；
[0180] 输出：对象区域Rm的绣线集合η={ηk:k=1～2N},每根绣线为
[0181]
[0182] 步骤381：设置对象区域Rm的绣线集合η为空；
[0183] 步骤382：对针迹序列集合Λ中的每个Λi执行步骤3821～步骤3824；
[0184] 步骤3821：提取其χi中的参数ξi、θi、βi、hi、si、vi，设置两个新的绣线η2i、η2i+1加入到η中，2i和2i+1分别为绣线的序号；
[0185] 步骤3822：分别计算η2i和η2i+1的起点和终点坐标，
[0186] 步骤3823：利用颜色空间转换方法将χi中HSV颜色空间的hi、si、vi值转化为RGB颜色空间的三通道值rgbr，rgbg，rgbb；
[0187] 步骤3824：绣线颜色库中的颜色值为RGB值，从绣线颜色库中选择与RGB三通道值rgbr，rgbg，rgbb距离最小的绣线颜色的RGB值 作为η2i和η2i+1的颜色值，
[0188] 最后的绣制效果如图6所示，对附图6而言，实施例中一个是前景的花朵区域，另一个是背景区域，其中，花朵区域的高层绘制参数Ω={ρ,λθ,λξ,λh,λs,λv}的取值分别为{ρ=0.89，λθ=0.6，λξ=2.4，λh=0.09，λs=0.13，λv=0.05}，背景区域的高层绘制参数的取值分别为{ρ=0.89，λθ=-0.15，λξ=5.3，λh=-2.19，λs=-0.5，λv=0.4}。
[0189] 本发明实现了一种基于反应扩散过程的计算机辅助乱针绣制作方法，具体实现该技术方案的方法和途径很多，以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
