一种基于多尺度梯度差分熵的图像质量评价方法

一种基于多尺度梯度差分熵的图像质量评价方法
【技术领域】
[0001] 本发明属于图像质量评价技术领域，设及一种全参考图像质量评价方法，具体设 及一种基于多尺度梯度差分滴的图像质量评价方法，适用于图像编码质量评价中。
【背景技术】
[0002] 据统计，人类接收到的70%信息来自于视觉。观察图像是人类认识世界的有效途 径。尤其是随着互联网技术的快速发展，图像给本发明日常生活带来了巨大的变化，并在医 疗、新闻、娱乐等多个领域占据着重要地位。而在图像获取、处理及传输的过程中，图像质量 不可避免的会受到影响，出现失真（模糊失真、噪声失真、编码失真等），从而无法达到应用 的要求。因此，如何评价图像质量的好坏是非常重要的。
[0003] 常用的图像质量评价方法包括主观图像质量评价方法和客观图像质量评价方法。 主观图像质量评价方法主要是W人的主观感受为主，具体方法是采用专口的观察人员，根 据事先制定的准则对一组图像进行判断，最后对评判结果进行加权求和得出最终的评价结 果。主观评价方法耗时长，成本高，容易受到人的情绪、自身条件等因素的影响，因此，目前 应用较为广泛的是客观图像质量评价方法。
[0004] 客观图像质量评价方法包括全参考图像质量评价方法、部分参考图像质量评价方 法和无参考图像质量评价方法。全参考图像质量评价方法对原始图像和失真图像之间的 差异性进行判断，可W取得较好的评价效果。常用的全参考图像质量评价方法包括峰值信 噪比（PeakSi即alNoiseRatio,PSNR)、结构相似度（Struc1:uralSIMilarity，SSIM)， 视觉信噪比（VisualSi即al-to-NoiseRatio,VSNR)、信息保真度（VisualInformation Fidelity,VI巧等。不同的质量评价方法适用于不同失真类型图像的质量评价中。
[0005] 客观图像质量评价方法是否有效的依据是看其客观评价值与人类主观观测意 见的一致性。该一致性可W通过相关系数的大小来表示，分别是；Spearman相关系数 (SpearmanRank-orderCorrelationCoefficient，SROCC)，Kendall相关系数化endall Rank-orderCorrelationCoefficient，KROCC)，Pearson相关系数(PearsonLinear CorrelationCoefficient，CC)。
[0006] 客观图像质量评价方法的优点是速度快、费用低、应用广泛，不受主观因素的影 响。不足之处在于其只能在某些方面有限度的模仿人眼的主观视觉特性，常常会出现与主 观评价结果不一致的情况，不同的评价方法适用于不同的应用领域，难W建立适用于任何 领域的评价方法。

【发明内容】

[0007] 本发明针对现有技术的不足，提供一种基于多尺度梯度差分滴的图像质量评价方 法。
[0008] 本发明所采用的技术方案是：一种基于多尺度梯度差分滴的图像质量评价方法， 其特征在于，包括W下步骤：
[0009] 步骤1 ;将原始图像和失真图像均转换成灰度图像；
[0010] 步骤2 ;基于步骤1的结果，计算不同尺度空间内的梯度差分滴，并将不同尺度空 间内的梯度差分滴进行加权平均；
[0011] 步骤3 ;计算原始图像和失真图像对应的灰度图像均值的比值；
[0012] 步骤4 ;计算原始图像和失真图像对应的灰度图像方差的比值；
[0013] 步骤5 ;基于步骤2、步骤3和步骤4的结果，进行综合评价。
[0014] 作为优选，步骤2的具体实现包括W下子步骤：
[0015] 步骤2. 1 ;基于步骤1的结果，构建不同尺度空间；高斯卷积核是尺度变换的唯一 线性核，原始图像和失真图像的尺度空间分别通过对应的灰度图像与高斯卷积核进行卷积 获得，不同的尺度由5参数来确定；
[0016] 步骤2. 2 ;对不同尺度下的原始图像和失真图像通过下采样的方式生成不同采样 间隔的图像；
[0017] 步骤2.3 ;计算不同尺度及采样间隔下的原始图像和失真图像的梯度大小，并进 行差分运算；
[0018] 步骤2. 4 ;基于步骤2. 3的结果，进行量化及绝对值运算；
[001引步骤2. 5 ;基于步骤2. 4的结果，计算滴；
[0020] 步骤2. 6 ;基于步骤2. 5的结果，将不同尺度不同采样间隔下的滴进行加权平均。
[0021] 作为优选，步骤3的具体实现包括W下子步骤：
[0022] 步骤3. 1 ;计算原始图像对应的灰度图像均值；
[002引步骤3. 2 ;计算失真图像对应的灰度图像均值；
[0024] 步骤3. 3 ;基于步骤3. 2和步骤3. 1的结果，计算两者的比值。
[00巧]作为优选，步骤4的具体实现包括W下子步骤：
[002引步骤4. 1 ;计算原始图像对应的方差；
[0027] 步骤4. 2 ;计算失真图像对应的方差；
[002引步骤4. 3 ;基于步骤4. 1和步骤4. 2的结果，计算两者的比值。
[0029] 本专利中主要针对图像编码质量评价的应用需求，从人眼观测特性出发，引入了 SIFT(ScaleInvariantFea1:ureTransform,SIFT)算法中的尺度空间的概念，通过提取尺 度空间下的轮廓特征差异，建立一种与人的主观感受较为吻合的图像质量评价方法。
[0030] 本发明引入尺度空间的基本思想，尺度空间中各尺度图像的模糊程度逐渐变大， 能够模拟人在不同距离下目标在视网膜上的形成过程，因此通过提取尺度空间下的轮廓特 征，并通过滴来反映原始图像和失真图像在尺度空间下轮廓特征的差异，与主观评价结果 具有较好的一致性。
【附图说明】
[0031] 图1 ;是本发明实施例的流程图。
[0032] 图2 ;是本发明实施例1的原始图像。
[0033] 图3 ;是本发明实施例1的失真图像。
[0034] 图4 ;是本发明实施例1的采样方法示意图。
[00巧]图5 ;是本发明实施例1的高斯金字塔示意图。
【具体实施方式】
[0036] 为了便于本领域普通技术人员理解和实施本发明，下面结合附图及实施例对本发 明作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本发明，并不 用于限定本发明。
[0037] 请见图1、图2和图3,图2是原始图像，图3是失真图像，图像空间分辨率为 512*512,本实施例1提供的一种基于多尺度梯度差分滴的图像质量评价方法，包括W下步 骤：
[003引步骤1 ;将原始图像和失真图像均转换成灰度图像；
[00測转换方法如下；
[0040] I(X,y) = 0. 2989*R(X,y) +0. 5870*G(X,y) +0. 1140地（X,y);
[0041] I(x，y)代表了转换后的灰度图像在空间坐标（x，y)处对应的像素大小，R(x，y)、 G(x，y)和B(x，y)分别代表原始图像（失真图像）在空间坐标（X，y)处对应的R、G和B分 量大小。
[0042] 步骤2 ;基于步骤1的结果，计算不同尺度空间内的梯度差分滴，并将不同尺度空 间内的梯度差分滴进行加权平均；
[0043] 本步骤进一步包括子步骤：
[0044] 步骤2. 1 ;基于步骤1的结果，构建不同尺度空间；高斯卷积核是尺度变换的唯一 线性核，原始图像和失真图像的尺度空间可W分别通过对应的灰度图像与高斯卷积核进行 卷积获得，不同的尺度由5参数来确定。
[0045] G(x,y, 5 )高斯卷积核的定义如下；
[0046]
[0047] 卷积方法定义如下：
[0048] L(X，y，5 ) =G(X，y，5 ) *1 (X，y);
[0049] 在卷积运算时，需要确定卷积运算的窗口大小，在本实施例中，窗口大小为3*3。尺 度参数5分别为1,1. 5,2,2. 5。调整5,本发明可W获得4个不同的高斯卷积核，分别为：
[0052] 通过上述4个高斯卷积核，可W获得4幅不同尺度空间的原始图像和4幅不同尺 度空间的失真图像。
[0053] 步骤2. 2 ;下采样采用隔点采样的方式，在本实施例中，采样间隔分为设置为3、4、 5、6、7、8、9、10。采样方法如下：假设采样间隔为T，也就是说在原图的水平方向和垂直方 向每隔T取一个像素，将所有取得的像素按照行和列的顺序依次排列后获得下采样后的图 像，采样间隔为4时的采样方法示意图如图4所示。通过步骤2. 1和步骤2. 2,可W构建原 始图像和失真图像的高斯金字塔，其示意图如图5所示，每个高斯金字塔共包含4*8 = 32 幅图像。
[0054] 步骤2. 3;计算不同尺度及采样间隔下的原始图像和失真图像的梯度大小，并进 行差分运算。
[00巧]梯度大小Ga(x,y)计算公式如下；
[0059] 其中，Img(x，y)代表不同尺度及采样间隔下的图像在空间坐标为（X，y)处的像素 灰度大小。
[0060] 差分运算的计算方法如下：
[0061]Diresult(X，y) =GaOriginal(x，y)-GaDistortion(x，y);
[0062]Ga化iginal(X，y)代表不同尺度及采样间隔下的原始图像对应的梯度大小， GaDistcxrtion(x,y)代表不同尺度及采样间隔下的失真图像对应的梯度大小。
[0063] 步骤2. 4;基于步骤2. 3的结果，进行量化及绝对值运算。
[0064] 考虑到人眼分辨率的有限性，本发明仅仅考虑差异达到一定数值的梯度带来的影 响，因此，本发明首先对梯度大小的差异进行量化取整，在此实施例中，量化因子设为8。进 一步忽略梯度差异的符号，对上述结果取其绝对值。其计算方法如下：
[0065]Diresult(X,y)=油S(round(Diresult(X,y)/8));
[0066] 油s( ?)代表绝对值运算，rouncK?)代表取整运算。
[0067] 步骤2. 5;基于步骤2. 4的结果，计算滴。
[0068] 滴的计算方法如下：
[0069]Entropy=-PilogaPi;
[0070]Pi代表不同数值在Diresult中出现的概率。
[0071] 步骤2. 6 ;基于步骤2. 5的结果，将不同尺度及采样间隔下的滴进行加权平均。
[0072] 在本实施例中，不同尺度及采样间隔下对应的滴分别如表1所示：
[0073] 表1不同尺度及采样间隔下对应的滴
[0074]
[00巧]对上述结果进行加权平均后的结果为0. 2393。
[0076] 步骤3 ;计算原始图像和失真图像对应的灰度图像均值的比值；
[0077] 本步骤进一步包括子步骤：
[0078] 步骤3. 1 ;计算原始图像对应的灰度图像均值；
[0079] 在本实施例中，原始图像对应的灰度图像均值为131. 03。
[0080] 步骤3. 2 ;计算失真图像对应的灰度图像均值；
[0081] 在本实施例中，失真图像对应的灰度图像均值为131. 08。
[008引步骤3. 3 ;基于步骤3. 2和步骤3. 1的结果，计算两者的比值。
[0083] 在本实施例中，两者的比值为0.9996.
[0084] 步骤4 ;计算原始图像和失真图像对应的灰度图像方差的比值；
[0085] 本步骤进一步包括子步骤：
[008引步骤4. 1 ;计算原始图像对应的方差；
[0087] 在本实施例中，原始图像对应的灰度图像方差为52. 2336。
[008引步骤4. 2 ;计算失真图像对应的方差；
[0089] 在本实施例中，失真图像对应的灰度图像方差为51. 6626。
[0090] 步骤4. 3 ;基于步骤4. 1和步骤4. 2的结果，计算两者的比值。
[0091] 在本实施例中，两者的比值为1. 0111。
[0092] 步骤5 ;基于步骤2、步骤3和步骤4的结果，进行综合评价。
[0093] 综合评价方法为=者相乘取负号，其结果为-0. 2393*0. 9996*1. 0111 = -0. 2419。
[0094] 本实施例2采用国际通用CSIQ图像质量评价数据（Mosta卵arent distortion:full-referenceimagequalityassessmentandtheroleof strategy,JournalofElectronicImaging),它由30张无失真参考图像和866张失真图 像组成，图像质量采用差分平均主观分DM0S组成。该里本发明主要测量了JPEG压缩失 真、JPEG2000压缩失真、加性高斯白噪声和模糊失真四种类型。采用S种评价指标来评价 不同图像质量评价算法的性能；斯皮尔曼相关系数（Spearmanrankordercorrelation coefficient,SROCC)，坎德尔相关系数化endallrankordercorrelationcoefficient, KROCC)和皮尔森相关系数（Pearsoncorrelationcoefficient)。相关系数的数值越大， 代表相关性越强。为了进行对比，本发明也对目前应用较为广泛的SSIM评价指标进行了 实验。相关实验结果如表2和表3所示。对比表2和表3可W发现，本方法与主观评价方 法的相关性在JPEG压缩失真、JPEG2000压缩失真及模糊失真时均强于SSIM评价指标。因 此，本发明可W得出如下结论；本方法与主观评价方法具有较强的一致性，可W被广泛应用 于各种图像质量评价系统中，尤其是视频/图像编码质量评价等应用场合。
[0095] 表2本文评价方法与主观评价方法的相关性
[0096]
[0099] 本发明引入尺度空间的基本思想，尺度空间中各尺度图像的模糊程度逐渐变大， 能够模拟人在不同距离下目标在视网膜上的形成过程，因此与主观评测过程具有相似性。 本发明通过提取不同尺度和不同分辨率下原始图像和失真图像梯度差分滴的信息，通过其 对图像的失真程度进行评测，达到了较好的评价效果。
[0100] 应当理解的是，本说明书未详细阐述的部分均属于现有技术。
[0101] 应当理解的是，上述针对较佳实施例的描述较为详细，并不能因此而认为是对本 发明专利保护范围的限制，本领域的普通技术人员在本发明的启示下，在不脱离本发明权 利要求所保护的范围情况下，还可W做出替换或变形，均落入本发明的保护范围之内，本发 明的请求保护范围应W所附权利要求为准。
【主权项】
1. 一种基于多尺度梯度差分熵的图像质量评价方法，其特征在于，包括以下步骤： 步骤1:将原始图像和失真图像均转换成灰度图像； 步骤2 :基于步骤1的结果，计算不同尺度空间内的梯度差分熵，并将不同尺度空间内 的梯度差分熵进行加权平均； 步骤3 :计算原始图像和失真图像对应的灰度图像均值的比值； 步骤4 :计算原始图像和失真图像对应的灰度图像方差的比值； 步骤5 :基于步骤2、步骤3和步骤4的结果，进行综合评价。2. 根据权利要求1所述的基于多尺度梯度差分熵的图像质量评价方法，其特征在于： 步骤2的具体实现包括以下子步骤： 步骤2. 1 :基于步骤1的结果，构建不同尺度空间；高斯卷积核是尺度变换的唯一线 性核，原始图像和失真图像的尺度空间分别通过对应的灰度图像与高斯卷积核进行卷积获 得，不同的尺度由S参数来确定； 步骤2. 2 :对不同尺度下的原始图像和失真图像通过下采样的方式生成不同采样间隔 的图像； 步骤2. 3 :计算不同尺度及采样间隔下的原始图像和失真图像的梯度大小，并进行差 分运算； 步骤2. 4 :基于步骤2. 3的结果，进行量化及绝对值运算； 步骤2. 5 :基于步骤2. 4的结果，计算熵； 步骤2. 6 :基于步骤2. 5的结果，将不同尺度不同采样间隔下的熵进行加权平均。3. 根据权利要求1所述的基于多尺度梯度差分熵的图像质量评价方法，其特征在于： 步骤3的具体实现包括以下子步骤： 步骤3. 1 :计算原始图像对应的灰度图像均值； 步骤3. 2 :计算失真图像对应的灰度图像均值； 步骤3. 3 :基于步骤3. 2和步骤3. 1的结果，计算两者的比值。4. 根据权利要求1所述的基于多尺度梯度差分熵的图像质量评价方法，其特征在于： 步骤4的具体实现包括以下子步骤： 步骤4. 1 :计算原始图像对应的方差； 步骤4. 2 :计算失真图像对应的方差； 步骤4. 3 :基于步骤4. 1和步骤4. 2的结果，计算两者的比值。
【专利摘要】本发明公开了一种基于多尺度梯度差分熵的图像质量评价方法，包括如下步骤：步骤1，将原始图像和失真图像均转换成灰度图像；步骤2，基于步骤1的结果，计算不同尺度空间内的梯度差分熵，并将不同尺度空间内的梯度差分熵进行加权平均；步骤3，计算原始图像和失真图像对应的灰度图像均值的比值；步骤4，计算原始图像和失真图像对应的灰度图像方差的比值；步骤5，基于步骤2，步骤3和步骤4的结果，进行综合评价。本发明引入尺度空间的基本思想，基于人在不同距离下目标在视网膜上的形成过程，通过提取尺度空间下的轮廓特征，并通过熵来反映原始图像和失真图像在尺度空间下轮廓特征的差异，与主观评价结果具有较好的一致性。
【IPC分类】G06T7/00
【公开号】CN104881878
【申请号】CN201510325238
【发明人】田昕, 李松, 郑国兴, 周辉, 杨晋陵, 高俊玲
【申请人】武汉大学
【公开日】2015年9月2日
【申请日】2015年6月12日