发明名称 一种综合多特征图像检索方法
摘要 本发明涉及一种综合多特征图像检索方法,包括进行图像特征的提取、索引和特征匹配,所述图像特征包括颜色特征、纹理特征和形状特征,所述图像的颜色特征提取包括:(1)将图像规格化为128×128像素;(2)将一幅图像分成m×n个小块;(3)对于每一小块图像的每个像素计算得出C’值,选定主C’值,由各主C’值形成一个对应的二维矩阵A。本发明通过对传统图像颜色特征的提取方法做出了改进,改进了传统局部颜色直方图法,比一般的基于颜色的图像检索方法的查准率有了较大的提高;采用了融合颜色、纹理、形状多种图像特征的图像检索方法,能够有效地提高方法的查准率。
申请公布号 CN101551823B 申请公布日期 2011.06.08
申请号 CN200910097741.5 申请日期 2009.04.20
申请人 浙江师范大学 发明人 朱信忠;赵建民;徐慧英;袁杰
分类号 G06F17/30(2006.01)I 主分类号 G06F17/30(2006.01)I
代理机构 杭州赛科专利代理事务所 33230 代理人 王桂名
主权项 1.一种综合多特征图像检索方法,包括进行图像特征的提取、索引和特征匹配,所述图像特征包括颜色特征、纹理特征和形状特征,其特征在于:所述图像检索方法包括下述步骤:(1)将例图和当前检索图像大小规格化为128×128像素;(2)将规格化后的图像颜色空间转换到HSV空间,同时进行量化,量化后共67级;(3)将图像分成32×32大小的块共16块,得出每块的颜色直方图,同时将这些块继续分成8×8大小的小块,提取各小块的主C′值,得到整幅图形的二维矩阵;每一小块图像的主C′值按下述方法提取:对于每一小块图像的每个像素计算得出C’值,选定主C’值,由各主C’值形成一个对应的二维矩阵A,大小为m×n,表示为A={a<sub>ij</sub>} i=1,2,…m;j=1,2,…n;主C’值的选定原则是:将像素个数最多的C’值作为该小块图像的主C’值,如果多个C’值的像素个数相同且都最大时则取平均值后取整做为主C’值;图像C’值的计算方法如下:对图像选用HSV颜色空间,HSV各分量的范围为:h∈[0,360],s∈[0,1],v∈[0,1],规定h、s、v为RGB空间颜色转换到HSV空间的连续值,H、S、V为h、s、v三个量量化后的离散值,量化方法如下式:<img file="FSB00000431642500021.GIF" wi="1575" he="746" /><maths num="0001"><![CDATA[<math><mrow><mi>S</mi><mo>,</mo><mi>V</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mi>s</mi><mo>,</mo><mi>v</mi><mo>&lt;</mo><mn>0.2</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>0.2</mn><mo>&le;</mo><mi>s</mi><mo>,</mo><mi>v</mi><mo>&lt;</mo><mn>0.7</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>s</mi><mo>,</mo><mi>v</mi><mo>&GreaterEqual;</mo><mn>0.7</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>量化后合成一维的颜色矢量C,合成方法如下:①对于v≤0.2的颜色认为是黑色,C=0;②对于s≤0.1,且v>0.2的颜色按亮度v划分为三种灰度,分别为:v∈(0.2,0.5]深灰,v∈(0.5,0.8]浅灰和v∈(0.8,1]白色,C的相应值分别为C=1,2,3;③其它情况下由式(3)算出C:C=4+9H+3S+V                                            (3)进而,C’取值如下式:<maths num="0002"><![CDATA[<math><mrow><msup><mi>C</mi><mo>,</mo></msup><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mi>H</mi></mtd><mtd><mi>C</mi><mo>&NotEqual;</mo><mn>0,1,2,3</mn></mtd></mtr><mtr><mtd><mn>7</mn></mtd><mtd><mi>C</mi><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mn>8</mn></mtd><mtd><mi>C</mi><mo>=</mo><mn>1,2</mn></mtd></mtr><mtr><mtd><mn>9</mn></mtd><mtd><mi>C</mi><mo>=</mo><mn>3</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>(4)按分块颜色直方图相交法算出两幅图像的相似度Similarity<sub>11</sub>,公式如下:<maths num="0003"><![CDATA[<math><mrow><msub><mi>Similarity</mi><mn>11</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>16</mn></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>16</mn></munderover><mo>{</mo><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mn>66</mn></munderover><mi>min</mi><mrow><mo>(</mo><msub><mi>a</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>b</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>]</mo><mo>/</mo><mrow><mo>(</mo><mn>32</mn><mo>&times;</mo><mn>32</mn><mo>)</mo></mrow><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中a<sub>i,j</sub>,b<sub>i,j</sub>分别表示两幅图像的第i(1≤i≤16)块的C值为j(0≤j≤66)的像索个数;(5)按顺序扫描步骤(3)得到的二维矩阵,得出图像的颜色相邻矩阵,并算出两幅图像的相似度Similarity<sub>12</sub>,具体步骤如下:(5.1)建立一个k×k,k=10的矩阵P,各元素初始值为0;将上步得到的矩阵A按Z字扫描顺序进行扫描得到扫描序列E,设E<sub>i</sub>与E<sub>i+1</sub>是扫描序列中的一对相继出现的颜色,E<sub>i</sub>在E<sub>i+1</sub>的前面,则P中相应元素P[E<sub>i</sub>,E<sub>i+1</sub>]自增1,反复直到扫描完成;(5.2)建立k×k的矩阵D,D中元素的计算公式如下式:<maths num="0004"><![CDATA[<math><mrow><msub><mi>D</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>则矩阵D就是该图像的颜色相邻矩阵;(5.3)采用下述方法进行相似度表征:设D<sub>A</sub>,D<sub>B</sub>分别表示两幅图像A,B的颜色相邻矩阵,大小皆为k×k,k=10,则两幅图像间相似度定义如下式:<maths num="0005"><![CDATA[<math><mrow><msub><mi>Similarity</mi><mn>12</mn></msub><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>B</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><mrow><mo>(</mo><msub><mi>D</mi><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></msub><mo>-</mo><msub><mi>&mu;</mi><mi>A</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>D</mi><msub><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></msub><mo>-</mo><msub><mi>&mu;</mi><mi>B</mi></msub><mo>)</mo></mrow></mrow><msqrt><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msup><mrow><mo>(</mo><msub><mi>D</mi><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></msub><mo>-</mo><msub><mi>&mu;</mi><mi>A</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msup><mrow><mo>(</mo><msub><mi>D</mi><msub><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></msub><mo>-</mo><msub><mi>&mu;</mi><mi>B</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow></msqrt></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,<maths num="0006"><![CDATA[<math><mrow><msub><mi>&mu;</mi><mi>A</mi></msub><mo>=</mo><mfrac><mn>1</mn><msup><mi>k</mi><mn>2</mn></msup></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>D</mi><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msub><mi>&mu;</mi><mi>B</mi></msub><mo>=</mo><mfrac><mn>1</mn><msup><mi>k</mi><mn>2</mn></msup></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>D</mi><msub><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>(6)按下式(10)对similarity<sub>11</sub>和similarity<sub>12</sub>进行合成得到图像颜色特征相似度Similarity<sub>1</sub>:Similarity<sub>1</sub>=Similarity<sub>11</sub>×0.6+Similarity<sub>12</sub>×0.4        (10)(7)首先将图像转换成灰度图,再由CANNY算子得出图像的边界曲线,然后得到图像链码直方图表示的形状特征,最后再按直方图相交法计算形状特征相似度Similarity<sub>2</sub>;(8)将步骤(7)中得到的图像边界曲线,提取出基于游程长度统计的图像纹理特征参数T<sub>1</sub>,T<sub>2</sub>,T<sub>3</sub> T<sub>4</sub>,按下式(11)计算出两幅图像的纹理特征相似度Similarity<sub>3</sub>:<maths num="0008"><![CDATA[<math><mrow><msub><mi>Similarity</mi><mn>3</mn></msub><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>4</mn></munderover><mo>[</mo><mo>|</mo><msub><mi>T</mi><mrow><mi>a</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>T</mi><mrow><mi>b</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>|</mo><mo>/</mo><mi>max</mi><mrow><mo>(</mo><msub><mi>T</mi><mrow><mi>a</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>,</mo><msub><mi>T</mi><mrow><mi>b</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中T<sub>a,i</sub>,T<sub>b,i</sub>分别表示两幅图像的T<sub>i</sub>值;(9)将上面得到的图像颜色特征相似度Similarity<sub>1</sub>、形状特征相似度Similarity<sub>2</sub>、纹理特征相似度Similarity<sub>3</sub>按下式(12)加权相加合成为两幅图像总的相似度:<maths num="0009"><![CDATA[<math><mrow><mi>similarity</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><msub><mi>w</mi><mi>i</mi></msub><mo>*</mo><msub><mi>similarity</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中w<sub>1</sub>=0.6,w<sub>2</sub>=0.2,w<sub>3</sub>=0.2;(10)对图像库中的每幅图像重复以上步骤,计算出每幅图像与例图的相似度;(11)将图像按相似度的逆序排列,返回给用户。
地址 321004 浙江省金华市婺城区迎宾大道688号