发明名称 织物图案创意要素自动提取方法
摘要 本发明公开了一种织物图案创意要素自动提取方法,包括几何元素提取流程和色彩元素提取流程,几何元素提取流程包括消除纹理噪声、建立区域作用域图和区域合并;色彩元素提取流程包括对输入图像进行显著性分析,得到反映该图像的显著性图,运用颜色量化中的量化标准分别针对划分的两区域进行颜色值的量化。与现有技术相比,本发明术可以提取出更加符合人的主观视觉理解、更具美学价值的色彩主题。
申请公布号 CN104778703A 申请公布日期 2015.07.15
申请号 CN201510181282.4 申请日期 2015.04.15
申请人 天津大学 发明人 刘世光;姜亚茜;罗华荣
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 李素兰
主权项 一种织物图案创意要素自动提取方法,包括几何元素提取流程和色彩元素提取流程,其特征在于,所述几何元素提取流程和所述色彩元素提取流程,具体包括以下步骤:步骤1、实现所述几何元素提取流程,包括:步骤1‑1、通过多次均值漂移算法迭代得到织布图案中一个给定球体范围区域内点的均值,作为密度梯度的估计,直到得到局部稳定状态中心点的聚类中心点,完成纹理噪声的消除;S是p维空间下一个有限数目点的集合,是一个半径为r的球面区域,中心点为x,漂移特征向量为y,满足|y‑x|<r,期望的新的点的坐标μ=y‑x,通过如下的公式(1)计算得到:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&mu;</mi><mo>=</mo><msub><mo>&Integral;</mo><msub><mi>S</mi><mi>x</mi></msub></msub><mrow><mo>(</mo><mi>y</mi><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mi>p</mi><mrow><mo>(</mo><mi>y</mi><mo>|</mo><msub><mi>S</mi><mi>x</mi></msub><mo>)</mo></mrow><mi>dy</mi><mo>=</mo><msub><mo>&Integral;</mo><msub><mi>S</mi><mi>x</mi></msub></msub><mrow><mo>(</mo><mi>y</mi><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mfrac><mrow><mi>p</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><mi>p</mi><mrow><mo>(</mo><mi>y</mi><mo>&Element;</mo><msub><mi>S</mi><mi>x</mi></msub><mo>)</mo></mrow></mrow></mfrac><mi>dy</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000699711360000011.GIF" wi="1413" he="143" /></maths>步骤1‑2、基于图的图像分割方法先对一幅织物图案图像建立一个无向图G=(V,E),每个像素点对应织物图像中每个节点,连接节点的边对应图像像素点之间三维色彩空间下的差值,块与块之间为节点的形式,预先定义每个块内部的差异Int(C):<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>Int</mi><mrow><mo>(</mo><mi>C</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>max</mi><mrow><mi>e</mi><mo>&Element;</mo><mi>MST</mi><mrow><mo>(</mo><mi>C</mi><mo>,</mo><mi>E</mi><mo>)</mo></mrow></mrow></munder><mi>&omega;</mi><mrow><mo>(</mo><mi>e</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000699711360000012.GIF" wi="732" he="106" /></maths>其中,<img file="FDA0000699711360000014.GIF" wi="143" he="72" />且MST(C,E)表示最小生成树算法形成的图的子集,ω(e)表示无向图中边e的权重;预先定义块与块之间的差异Dif(C<sub>1</sub>,C<sub>2</sub>):<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>Dif</mi><mrow><mo>(</mo><msub><mi>C</mi><mn>1</mn></msub><mo>,</mo><msub><mi>C</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><munder><mi>min</mi><mrow><msub><mi>v</mi><mi>i</mi></msub><mo>&Element;</mo><msub><mi>C</mi><mn>1</mn></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>&Element;</mo><msub><mi>C</mi><mn>2</mn></msub><mo>,</mo><mrow><mo>(</mo><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>&Element;</mo><mi>E</mi></mrow></munder><mi>&omega;</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000699711360000013.GIF" wi="973" he="118" /></maths>其中v<sub>j</sub>,v<sub>i</sub>为无向图G=(V,E)中的顶点,ω(v<sub>i</sub>,v<sub>j</sub>)为两顶点路径的权值;如果C<sub>1</sub>和C<sub>2</sub>没有交集,则将Dif(C<sub>1</sub>,C<sub>2</sub>)置为无穷大;将边以递增的顺序进行排序,对区域C<sub>1</sub>和C<sub>2</sub>进行合并,如果两个区域相连的边值Dif(C<sub>1</sub>,C<sub>2</sub>)小于公式(5)则合并,否则不合并;MInt(C<sub>1</sub>,C<sub>2</sub>)=min(Int(C<sub>1</sub>)+τ(C<sub>1</sub>),Int(C<sub>2</sub>)+τ(C<sub>2</sub>))   (5)这里τ(C)=k/|C|,其中k是合并的尺度系数来调节合并的块的尺寸;步骤1‑3、对织物图案“连续性”的处理,具体为块与块的各种几何连续性合并;计算块与块的几何关系即块与块的包含性关系、相邻性关系和不相连性关系;其中:令C<sub>1</sub>、C<sub>2</sub>是过分割图像中的两个块,则二者存在包含性关系的概率ρ<sub>in</sub>(C<sub>1</sub>,C<sub>2</sub>)描述为公式(6):<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&rho;</mi><mi>in</mi></msub><mrow><mo>(</mo><msub><mi>C</mi><mn>1</mn></msub><mo>,</mo><msub><mi>C</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>cross</mi><mrow><mo>(</mo><mi>box</mi><mrow><mo>(</mo><msub><mi>C</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>,</mo><mi>box</mi><mrow><mo>(</mo><msub><mi>C</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mrow><mi>min</mi><mrow><mo>(</mo><mi>arc</mi><mrow><mo>(</mo><msub><mi>C</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>,</mo><mi>arc</mi><mrow><mo>(</mo><msub><mi>C</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000699711360000021.GIF" wi="1113" he="195" /></maths>其中,box()代表块的最小包围盒,cross()代表两个最小包围盒的交叉面积,arc()表示块的最小包围盒面积,min()两个数的最小值;令E<sub>1</sub>、E<sub>2</sub>为过分割图像中两个块的轮廓,则二者存在相邻性关系的概率描述为公式(7):<img file="FDA0000699711360000022.GIF" wi="1399" he="186" />其中p表示像素,count()计算几何中元素个数,len()代表周长,θ未设定的阈值,通过求两个相交块的边缘相交距离除以较小周长的边界来衡量相邻性的概率;令C<sub>1</sub>、C<sub>2</sub>为两个过分割块,则二者存在不相邻性关系的概率表述为公式(8):<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&rho;</mi><mi>disc</mi></msub><mrow><mo>(</mo><msub><mi>C</mi><mn>1</mn></msub><mo>,</mo><msub><mi>C</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><munder><mi>&Sigma;</mi><mrow><mi>p</mi><mo>&Element;</mo><msub><mi>C</mi><mn>1</mn></msub></mrow></munder><mrow><mo>(</mo><msub><mi>I</mi><mi>p</mi></msub><mo>/</mo><msub><mi>N</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><munder><mi>&Sigma;</mi><mrow><mi>p</mi><mo>&Element;</mo><msub><mi>C</mi><mn>2</mn></msub></mrow></munder><mrow><mo>(</mo><msub><mi>I</mi><mi>p</mi></msub><mo>/</mo><msub><mi>N</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>/</mo><msub><mi>I</mi><mi>rand</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000699711360000023.GIF" wi="1481" he="190" /></maths>其中,N<sub>1</sub>和N<sub>2</sub>分别代表C<sub>1</sub>、C<sub>2</sub>中像素点的个数,I<sub>p</sub>代表像素点p的强度值,I<sub>rand</sub>表示色彩空间中强度的范围;上述与合并相关的几何连续性的计算的前提是两个块的中心点距离小于设定的值;步骤2、实现所述织物色彩主题提取流程,包括:步骤2‑1、针对一幅图像I计算坐标(x',y')处显著性值S(x',y'),如公式(9)所示:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>S</mi><mrow><mo>(</mo><msup><mi>x</mi><mo>&prime;</mo></msup><mo>,</mo><msup><mi>y</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>=</mo><mo>|</mo><msub><mi>I</mi><mi>&mu;</mi></msub><mo>-</mo><msub><mi>I</mi><msub><mi>&omega;</mi><mi>hc</mi></msub></msub><mrow><mo>(</mo><msup><mi>x</mi><mo>&prime;</mo></msup><mo>,</mo><msup><mi>y</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000699711360000031.GIF" wi="1311" he="123" /></maths>其中,I<sub>μ</sub>表示图像全局的平均特征向量,<img file="FDA0000699711360000032.GIF" wi="217" he="93" />表示坐标(x',y')处像素值经过核大小为5×5的高斯差分滤波变换的特征向量值;实验中ω<sub>hc</sub>的值设置为π/2.75;步骤2‑2、将颜色转换到基于视觉的HSV色彩空间,对HSV三个分量进行不等间隔的划分,实现颜色主题的量化提取,颜色提取算法如公式(10)所示:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mi>h</mi><mo>&Element;</mo><mo>(</mo><mo>-</mo><mn>45,20</mn><mo>]</mo></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mi>h</mi><mo>&Element;</mo><mo>(</mo><mn>20,50</mn><mo>]</mo></mtd></mtr><mtr><mtd><mn>2</mn></mtd><mtd><mi>h</mi><mo>&Element;</mo><mo>(</mo><mn>50,75</mn><mo>]</mo></mtd></mtr><mtr><mtd><mn>3</mn></mtd><mtd><mi>h</mi><mo>&Element;</mo><mo>(</mo><mn>75,155</mn><mo>]</mo></mtd></mtr><mtr><mtd><mn>4</mn></mtd><mtd><mi>h</mi><mo>&Element;</mo><mo>(</mo><mn>155,195</mn><mo>]</mo></mtd></mtr><mtr><mtd><mn>5</mn></mtd><mtd><mi>h</mi><mo>&Element;</mo><mo>(</mo><mn>195,275</mn><mo>]</mo></mtd></mtr><mtr><mtd><mn>6</mn></mtd><mtd><mi>h</mi><mo>&Element;</mo><mo>(</mo><mn>275</mn><mo>,</mo><mo>-</mo><mn>45</mn><mo>]</mo></mtd></mtr></mtable></mfenced><mo>,</mo><mi>S</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mi>s</mi><mo>&Element;</mo><mo>[</mo><mn>0,0.2</mn><mo>]</mo></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mi>s</mi><mo>&Element;</mo><mo>(</mo><mn>0.2,0.7</mn><mo>]</mo></mtd></mtr><mtr><mtd><mn>2</mn></mtd><mtd><mi>s</mi><mo>&Element;</mo><mo>(</mo><mn>0.7,1</mn><mo>]</mo></mtd></mtr></mtable></mfenced><mo>,</mo><mi>V</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mi>v</mi><mo>&Element;</mo><mo>[</mo><mn>0,0.2</mn><mo>]</mo></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mi>v</mi><mo>&Element;</mo><mo>(</mo><mn>0.2,0.7</mn><mo>]</mo></mtd></mtr><mtr><mtd><mn>2</mn></mtd><mtd><mi>v</mi><mo>&Element;</mo><mo>(</mo><mn>0.7,1</mn><mo>]</mo></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000699711360000033.GIF" wi="1802" he="673" /></maths>
地址 300072 天津市南开区卫津路92号