本发明属于图像相似性计算技术领域，特别是一种基于频域视觉显著性分析的图像相似性计算方法。
背景技术
图像相似性计算的目的是衡量两幅图像内容的相似性程度，该技术广泛应用于基于内容的图像检索和图像分类中。相比于传统的基于标签的图像检索和分类，基于内容的方法直接考察图像包含的具体内容，更符合人们进行图像检索和分类的习惯，同时避免了繁琐的人工标注图像标签的过程。基于内容的图像检索或图像分类不需要对相关图像进行精确匹配，而是通过计算查询图像和候选图像之间视觉特征的相似度来完成，通过提取图像的纹理、颜色或者形状等特征形成特征向量，来代表对应的图像。在图像检索中，判断图像之间是否相似主要通过比较不同图像特征向量之间的距离。距离越大说明图像的相似度越小，距离越小说明图像相似度越大，即内容相似的图像间应该具有较小的距离，而内容不相似的图像间应具有较大的距离。按照人们的视觉习惯，图像内容主要指图像中的显著性前景部分，而图像中的背景信息在比较图像相似性时起到的作用比较小。
目前大部分图像相似性计算方法，都是通过计算图像之间的距离，从而对整幅图像的相似性进行计算，如果两幅图像背景颜色很相似，那么即使各自包含的前景对象完全不一样，也会得到较高的相似性，因此对整幅图像的相似性进行计算的方法难以满足基于内容的图像检索和图像分类的需求。如果先对图像进行精细的前景和背景分割，那么由于现有的前景分割算法计算量大，效率低，难以满足海量图像在线检索和分类的应用需求。
技术实现要素：
本发明的目的在于提供一种计算效率高、抗背景颜色一致能力强的基于频域视觉显著性分析的图像相似性计算方法。
实现本发明目的的技术解决方案为：一种基于频域视觉显著性分析的图像相似性计算方法，包括以下步骤：
步骤1、频域视觉显著性检测：将两张待比较的彩色数字图像转换为灰度图像，转换后分别进行二维离散傅里叶变换，对变换得到的幅度谱进行滤波，然后再分别进行傅里叶反变换，得到两张待比较图像的显著性图；
步骤2、显著性前景和非显著性背景权值比率计算：根据显著性图检测显著性前景区域和非显著性背景区域，分别计算显著性前景权重比率和非显著性背景权重比率；
步骤3、图像相似性计算：按照距离计算公式，计算待比较图像之间的距离，距离越小则图像相似度越高。
进一步地，步骤1所述的频域视觉显著性检测，具体如下：
步骤1.1、图像灰度变换：两张待比较的彩色数字图像i1和i2均为rgb彩色图像，根据式(1)将i1和i2分别转换为灰度图像，得到i1g和i2g：
i1g＝i1-r×0.299+i1-g×0.587+i1-b×0.114(1)
式(1)中i1-r、i1-g和i1-b分别表示i1的r、g、b三个颜色分量；同样得到i2g；
步骤1.2、二维离散傅里叶变换：将i1g和i2g分别进行二维离散傅里叶变换得到i1g-f和i2g-f，如式(2)
式(2)中h、w分别为图像i1g的高和宽，(x1g,y1g)、(u1g-f,v1g-f)分别为二维离散傅里叶变换前、后的像素坐标，j为虚数单位，r(u1g-f,v1g-f)为傅里叶变换后的实数部分，i(u1g-f,v1g-f)为傅里叶变换后的虚数部分；
式(2)中：
步骤1.3、幅度谱滤波：将式(2)中的幅度谱置为1，即|f(u1g-f,v1g-f)|＝1，则i1g-f成为式(5)所示i'1g-f的形式：
同样得i'2g-f；
步骤1.4、傅里叶逆变换：使用式(6)对i'1g-f和i'2g-f分别进行二维离散傅里叶逆变换得到i'1g和i'2g，即
i'1g和i'2g即根据频域计算得到的i1和i2的显著性图像，i'1g和i'2g的尺寸与各自原图像一致，均为灰度图像，灰度值越高的像素点代表原图像i1和i2中该位置对应的视觉显著性值高。
进一步地，步骤2所述的显著性前景和非显著性背景权值比率计算，具体如下：
步骤2.1、显著性前景和非显著性背景检测：显著性图像i'1g和i'2g中，显著性值大于等于θ的像素点分别定义为显著性前景集合
步骤2.2、权值计算：使用式(7)分别计算i1和i2的显著性前景权值比率rf
式(7)中size(·)表示显著性前景集合
使用式(8)分别计算i1和i2的非显著性背景权值比率：
式(8)中size(·)表示背景集合
进一步地，步骤3所述的图像相似性计算，具体如下：
步骤3.1、显著性前景和非显著性背景距离计算：分别统计
步骤3.2、图像距离计算：根据式(11)计算i1和i2之间的距离d(i1,i2)：
计算得到的d(i1,i2)值越大，则i1和i2的相似度越低；d(i1,i2)值越小，则i1和i2的相似度越高。
本发明与现有技术相比，其显著优点为：(1)抗背景颜色一致能力强：首先通过频域分析得到原图像的显著性前景，然后计算前景和背景的权重比率，计算图像间距离时前景被赋予较大权重，而背景被赋予较小权重，减小背景对图像相似性计算的影响，即使两幅图像背景颜色一致，只要显著性前景差异较大，依然能计算得到较大的图像距离值；(2)计算效率高：计算图像显著性时，使用的图像快速傅里叶变换和反变换，避免了在空间域进行大量的图像特征提取计算，具有很高的计算效率。
下面结合附图和具体实施方式对本发明作进一步的详细描述。
附图说明
图1为本发明基于频域视觉显著性分析的图像相似性计算方法的主流程图。
图2为本发明中频域视觉显著性检测步骤的流程图。
图3为本发明中显著性前景和非显著性背景权值比率计算步骤的流程图
图4为本发明中图像相似性计算步骤的流程图。
图5为本发明实施例中图像相似性计算示例图，其中(a)、(b)、(c)为3张待比较图像转换得到的灰度图，(d)、(e)、(f)为分别对(a)、(b)、(c)进行显著性检测得到的显著性图。
具体实施方式
如图1所示，本发明基于频域视觉显著性分析的图像相似性计算方法，包括以下步骤：
步骤1、频域视觉显著性检测：将两张待比较的彩色数字图像转换为灰度图像，转换后分别进行二维离散傅里叶变换，对变换得到的幅度谱滤波，然后再分别进行傅里叶反变换，得到两张待比较图像的显著性图；
如图2所示，所述的频域视觉显著性检测，具体如下：
步骤1.1、图像灰度变换：两张待比较的彩色数字图像i1和i2均为rgb彩色图像，根据式(1)将其分别转换为灰度图像，得到i1g和i2g：
i1g＝i1-r×0.299+i1-g×0.587+i1-b×0.114(1)
式(1)中i1-r、i1-g和i1-b分别表示i1的r、g、b三个颜色分量；
同样可得i2g；下面涉及到对i1和i2分别做同样的运算时均以i1为例说明，i2用同样方法计算即可。
步骤1.2、二维离散傅里叶变换：将i1g和i2g分别进行二维离散傅里叶变换得到i1g-f和i2g-f，如式(2)
式(2)中h，w为图像i1g的高和宽，u1g-f、v1g-f、x1g、y1g为像素坐标，j为虚数单位，r(u1g-f,v1g-f)为傅里叶变换后的实数部分，i(u1g-f,v1g-f)为傅里叶变换后的虚数部分；
式(2)中：
步骤1.3、幅度谱滤波：将式(2)中的幅度谱置为1，即|f(u1g-f,v1g-f)|＝1，则i1g-f成为式(5)所示i'1g-f的形式，即
同样可得i'2g-f
步骤1.4、傅里叶逆变换：使用式(6)对i'1g-f和i'2g-f分别进行二维离散傅里叶逆变换得到i'1g和i'2g，即
i'1g和i'2g即根据频域计算得到的i1和i2的显著性图像，i'1g和i'2g与各自原图像的分辨率一致，形式均为灰度图像，灰度值越高的像素点代表原图像i1和i2中该位置对应的显著性值高。
步骤2、显著性前景和非显著性背景权值比率计算：根据显著性图检测显著性前景区域和非显著性背景区域，分别计算显著性前景权重比率非显著性和背景权重比率；
如图3所示，所述的显著性前景和非显著性背景权值比率计算，具体如下：
步骤2.1、显著性前景和非显著性背景检测：显著性图像i'1g和i'2g中显著性值大于等于θ的像素点分别定义为显著性前景集合
步骤2.2、权值计算：使用式(7)分别计算i1和i2的显著性前景权值比率rf
式(7)中size(·)表示显著性前景集合
使用式(8)分别计算i1和i2的非显著性背景权值比率：
式(8)中size(·)表示背景集合
步骤3、图像相似性计算：按照定义的距离计算公式，计算待比较图像之间的距离，距离越小则图像相似度越高。
如图5所示，所述的图像相似性计算，具体如下：
步骤3.1、前景和背景距离计算：分别统计
步骤3.2、图像距离计算：根据式(11)计算i1和i2之间的距离d(i1,i2)：
计算得到的d(i1,i2)值越大，则i1和i2的相似度越低；d(i1,i2)值越小，则i1和i2的相似度越高。
实施例1
图5为图像相似性计算示例。其中，图5(a)，图5(b)和图5(c)为3张待比较rgb彩色图像转换的灰度图像，其中图5(a)和图5(b)前景相似，原图背景颜色差异较大，图5(b)和图5(c)背景颜色相似，前景差异较大。图5(d)，图5(e)和图5(f)为分别对图5(a)，图5(b)和图5(c)进行显著性检测得到的显著图。用本方法计算得到的图5(a)和图5(b)的相似性为0.1645，图5(c)和图5(b)的相似性为0.1925，表明图5(a)和图5(b)更加相似。直接统计整张图像的直方图，使用常见的全图比较的方法，依然使式(9)计算图像间的距离，计算得到图5(a)和图5(b)的相似性为0.2081，图5(c)和图5(b)的相似性为0.1856，表明图5(c)和图5(b)更加相似。该例说明本发明抗背景颜色一致能力较强，而全图比较的方法易受背景颜色一致影响。