本发明涉及图像质量处理技术领域,尤其涉及一种图像质量评价方法及系统。
背景技术:
数字图像作为信息表达与交流的元素,获得了广泛的运用。然而,数字图像在采集、压缩、传输、处理、重建过程中伴随着降低质量和失真[1]。因此,图像质量评价(imagequalityassessment,iqa)的研究成为数字图像处理领域的一个重要内容。
根据参考图像是否存在,图像质量评价可以分为全参考(fr)评价、半参考(rr)评价和无参考(nr)评价[2]。mse和psnr是迄今为止运用最为广泛的全参考图像质量评估方法。这两类方法虽然有着计算简单、易于实现等优点,但只计算了像素点之间的绝对误差,没有考虑像素点之间的相关性和人眼视觉系统的感知特性,容易造成在评价图像质量时出现严重不符实际的情况[3][4]。
近年来学者们提出了更加适应人眼视觉特性的评价算法。zhouwang等人[2]在通用图像质量指数(universalimagequalityindex,uqi)[5]的基础上,发展和产生了structuralsimilarityindex(ssim指数)。之后大批学者在ssim的基础上提出了许多改进的方法。后来,linzhang等人[6]基于自然场景分析(nss)的方法,提出引入相位一致性(pc)和梯度相似度(gm)作为评价特征的fsim模型。实验证明这些方法优于mse和psnr,但是ssim评价模型对图像的结构特征还只停留在标量的层面上,导致ssim在图片高度模糊时失去有效性。
自然图像作为二维信号,具有高度结构化的特征,必然具有矢量特性。图像各像素点之间表现出强烈的依赖关系,这种相关性构成了二维图像的结构,而hvs(人类视觉系统,humanvisualsystem)的主要功能就是从视野中获取结构信息。鉴于riesz变换在多维信号处理上的良好表现,linzhang等人[7]构建了rfsim模型,xue-gangluo等人[8]提出了rmfsim模型。但是这些方法只是简单的应用riesz变换构造局部特征,没有充分挖掘单演信号理论的物理含义,因此还有很大的改进空间。
[1]guk,zhaig,yangx,etal.subjectiveandobjectivequalityassessmentforimageswithcontrastchange[c]//imageprocessing(icip),201320thieeeinternationalconferenceon.ieee,2013:383-387.
[2]wangz,bovikac,sheikhhr,etal.imagequalityassessment:fromerrorvisibilitytostructuralsimilarity[j].ieeetransactionsonimageprocessing,2004,13(4):600-612.
[3]wangz,bovikac.meansquarederror:loveitorleaveit?anewlookatsignalfidelitymeasures[j].ieeesignalprocessingmagazine,2009,26(1):98-117.
[4]wangz.applicationsofobjectiveimagequalityassessmentmethods[applicationscorner][j].ieeesignalprocessingmagazine,2011,28(6):137-142.
[5]z.wang,a.c.bovik,h.r.sheik.“auniversalimagequalityindex”,ieeesignalprocessingletters,vol.9,no.3,pp.81-84,2002.
[6]zhangl,zhangl,moux,etal.fsim:afeaturesimilarityindexforimagequalityassessment[j].ieeetransactionsonimageprocessing,2011,20(8):2378-2386.
[7]zhangl,zhangl,moux.rfsim:afeaturebasedimagequalityassessmentmetricusingriesztransforms[c]//imageprocessing(icip),201017thieeeinternationalconferenceon.ieee,2010:321-324.
[8]luoxg,wanghj,wangs.monogenicsignaltheorybasedfeaturesimilarityindexforimagequalityassessment[j].aeu-internationaljournalofelectronicsandcommunications,2015,69(1):75-81.
技术实现要素:
本发明的目的在于解决当前图像质量评价方法对严重失真图像的评价不准确问题,提出了一种结合risez变换和人眼视觉系统特性的图像质量评价方法。该方法充分利用单演信号理论和log-gabor滤波器,同时结合csf(视觉对比度敏感函数,contrastsensitivityfunction)特性曲线来分配不同频段的权重,得到相似度矩阵,然后利用mpc(单演信号相位一致性,monogenicphasecongruency)图来构造池化函数,得到最终的图像质量评价指数,本发明通过在多个数据库上进行测试,结果表明,本发明整体性能优异。
本发明的上述技术问题主要是通过下述技术方案得以解决的:
一种图像质量评价方法,具体步骤如下:
步骤s1:利用log-gabor滤波器分别对参考图像r和失真图像d进行5尺度分解,
步骤s2:基于单演信号理论对分解后的图像进行riesz变换,得到5对单演信号;
步骤s3:对5对单演信号进行局部特征的相似度计算,得到5组单演相似度函数,进而得到5幅单演相似度图;
步骤s4:利用csf曲线设置5幅单演相似度图的权重,然后进行加权求和,得到1幅单演相似度图sm;
步骤s5:计算参考图像r和失真图像d的梯度幅度相似度图sg;
步骤s6:利用步骤s2得到的参考图像r的单演信号进行单演相位一致性计算,得到池化函数;
步骤s7:最后利用池化函数对单演相似度图sm和梯度幅度相似度图sg的乘积进行卷积,得到最终的图像质量评价指数。
其中,所述步骤s1中二维形式的log-gabor滤波器的频域表达式为:
其中,ω0为滤波器的中心频率,σ0决定了滤波器的带宽,θ0为滤波器的方向角,σθ决定了滤波器的角频宽。
其中,由于单个log-gabor滤波器的带宽限制,需要由多个log-gabor滤波器在半径方向和角度方向拼接成一个完整的滤波器组,采用的log-gabor滤波器组的拼接参数为:半径方向拼接数量nr=5,角度方向拼接数量nθ=1。
其中,利用log-gabor滤波器组分别对参考图像r和失真图像d进行5尺度分解,得到分解后的图像rbi(i=1,2,3,4,5)和dbi(i=1,2,3,4,5)log-gabor滤波器组的中心频率ω0分别为
其中,所述步骤s2中riesz变换核的空域表达形式为:
其频域响应为:
(h1(u,v),h2(u,v))=(-ju/||(u,v)||,-jv/||(u,v)||)(3)
以参考图像r为例,对其进行一阶riesz变换,得到:
参考图像r的单演信号rm被定义为:
rm=[r,r1,r2]=[r,h1(x,y)*r,h2(x,y)*r]
=r-(i,j)fr(x,y)(5)
其中,i和j为虚数单位;
riesz变换空间是一个球形坐标系,r、r1、r2分别是空间中的点在球形坐标系三个轴上的投影;在此空间域上,局部幅值a、局部方向θ和局部相位
其中,
其中,所述步骤s2中对分解后的图像进行riesz变换,得到5对单演信号,实现步骤如下:
第一步,对分解后的图像rbi(i=1,2,3,4,5)进行riesz变换,得到5
个单演信号
其中,
第二步,对分解后的图像dbi(i=1,2,3,4,5)进行riesz变换,得到5个单演信号
其中,
其中,步骤s3中对5对单演信号进行局部特征的相似度计算,得到5组单演相似度函数,进而得到5幅单演相似度图,具体为:
第一步,计算riesz变换空间下各点的特征量,包括幅值、相位角、方向角,构造参考图像和失真图像在像素点(x,y)的单演信号相似度:
其中,
第二步,构造参数smi,作为单演信号相似度图:
其中,i=1,2,3,4,5。
其中,步骤s4中利用csf曲线设置5幅单演相似度图的权重,然后进行加权求和,得到1幅单演相似度图sm,具体为:
利用csf曲线设置5幅单演相似度图的权重ki(i=1,2,3,4,5),然后进行加权求和,得到1幅单演相似度图sm:
sm=k1sm1+k2sm2+k3sm3+k4sm4+k5sm5(11)
其中k1=0.3370,k2=0.8962,k3=0.9809,k4=0.9753,k5=0.7411。
其中,所采用的视觉对比敏感度函数(contrastsensitivityfunction,csf)模型:
a(fr)≈2.6[0.0192+0.114fr]exp[-(0.114fr)1.1]
其中,
其中,所述步骤s5计算参考图像r和失真图像d的梯度幅度相似度图sg,实现步骤如下:
第一步,计算参考图像r和失真图像d的各像素点的梯度幅度g;
利用scharr算子来获得像素点的水平梯度gh和垂直梯度gv,然后计算参考图像和失真图像各像素点的梯度幅度g;以参考图像r为例,梯度幅度g计算公式如(13)、(14)所示:
第二步,计算参考图像r梯度幅度图gr和失真图像梯度幅度图gd的梯度幅度相似度图sg:
其中,gr(x,y)和gd(x,y)分别为参考图像和失真图像在像素点(x,y)的梯度幅值,c2和c3为小的正常数,c2=87497.64,c3=65025。
其中,所述步骤s6中利用参考图像r得到的单演信号进行单演相位一致性计算,得到池化函数mpc,实现步骤如下:
根据公式(7),可算出局部能量信息的总和为:
其中,
局部幅度的总和为:
则单演相位一致性模型(也就是池化函数)表示为:
其中,符号
其中,c是过滤响应滤波器的截止值,在此之下相位一致性将会受到抑制;g是一个控制截止明锐度的增益因子;
扩展函数s(x)为:
其中,amax(x,y)是滤波器在(x,y)处的最大响应的振幅,ε是一个小的正常数,其中,g=1.8182,
其中,所述步骤s7中利用池化函数mpc对单演相似度图sm和梯度幅度相似度图sg的乘积进行卷积,得到最终的图像质量评价指数,实现步骤如下:
联合公式(11)(15),同时利用单演相位一致性评价模型作为池化函数mpc,得到最终的质量评价分数:
rvsim=sm·sg*mpc(21)
公式(21)即为提出的图像质量评价指数。
一种采用上述所述的图像质量评价方法的系统,包括:
分解单元,用于利用log-gabor滤波器分别对参考图像r和失真图像d进行5尺度分解,
riesz变换单元,用于基于单演信号理论对分解后的图像进行riesz变换,得到5对单演信号;
第一计算单元,用于对5对单演信号进行局部特征的相似度计算,得到5组单演相似度函数,进而得到5幅单演相似度图;
第二计算单元,用于利用csf曲线设置5幅单演相似度图的权重,然后进行加权求和,得到1幅单演相似度图sm;
第三计算单元,用于计算参考图像r和失真图像d的梯度幅度相似度图sg;
第四计算单元,用于利用riesz变换单元得到的参考图像r的单演信号进行单演相位一致性计算,得到池化函数;
图像质量评价指数获取单元,用于最后利用池化函数对单演相似度图sm和梯度幅度相似度图sg的乘积进行卷积,得到最终的图像质量评价指数。
有益效果:
本发明提出了一种结合risez变换和人眼视觉系统特性的图像质量评价方法。该方法充分利用单演信号理论和log-gabor滤波器,同时结合csf(视觉对比度敏感函数,contrastsensitivityfunction)特性曲线来分配不同频段的权重,得到相似度矩阵,然后利用mpc(单演信号相位一致性,monogenicphasecongruency)图来构造池化函数,得到最终的图像质量评价指数,本发明通过在多个数据库上进行测试,结果表明,本发明整体性能优异。
附图说明
图1是riesz变换空间图。
图2是csf特性曲线图。
图3是图像质量评价方法流程框图。
图4是本发明设计的log-gabor滤波器组和滤波示例。
图5是monarch图、w(x,y)三维曲面图和w(x,y)三维曲旋转图。
图6是c1/c2/c3最优值的确定图。
图7是live数据库中不同iqa算法的评估散点图。
图8是17个iqa模型的gmad竞争图。
图9是图像质量评价系统结构框图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明:
实施例1
如图3所示,本发明提出的结合risez变换和人眼视觉系统特性的图像质量评价方法,(简称rvsim)包含以下步骤:
第一步,将参考图像r和失真图像d分别经过5尺度log-gabor带通滤波,得到图像在5个不同频段的分量rb1/rb2/rb3/rb4/rb5和db1/db2/db3/db4/db5。
第二步,基于单演信号理论对分解后的图像进行riesz变换得到5对单演信号
第三步,利用5对单演信号进行局部特征(包括局部幅度a、局部相位
,进而得到5幅单演相似度图sm1/sm2/sm3/sm4/sm5。
第四步,利用csf曲线设置5幅单演相似度图的权重k1/k2/k3/k4/k5,然后进行加权求和,得到1幅单演相似度图sm。
第五步,计算参考图像r和失真图像d的梯度幅度相似度图sg。
第六步,利用参考图像r得到的单演信号进行单演相位一致性(mpc)计算,得到池化(pooling)函数mpc。
第七步,最后利用池化函数mpc对单演相似度图sm和梯度幅度相似度图sg的乘积进行卷积,得到提出的图像质量评价指数。
本发明提出的结合risez变换和人眼视觉系统特性的图像质量评价方法,第一步中采用的二维形式的loggabor滤波器的频域表达式为:
其中,ω0为滤波器的中心频率,σ0决定了滤波器的带宽,θ0为滤波器的方向角,σθ决定了滤波器的角频宽。
在实际应用中,由于单个log-gabor滤波器的带宽限制,需要由多个log-gabor滤波器在半径方向和角度方向拼接成一个完整的滤波器组。
本发明最终确定的log-gabor滤波器组的拼接参数为:半径方向拼接数量nr=5,角度方向拼接数量nθ=1。为了说明这两个参数选取的合理性,表1列出了取不同拼接参数时,将rvsim算法应用于live数据库得到的srocc/krocc/plcc/rmse性能指标,最高性能用黑体字突出。从表1可以看出,径向拼接数量nr=5、横向拼接数量nθ=1时,rvsim算法取得最佳性能。
表一srocc值与live数据库每个失真类型的比较
利用log-gabor滤波器组分别对参考图像r和失真图像d进行5尺度分解,得到分解后的图像rbi(i=1,2,3,4,5)和dbi(i=1,2,3,4,5)log-gabor滤波器组的中心频率ω0分别为
所述步骤s2中riesz变换核的空域表达形式为:
其频域响应为:
(h1(u,v),h2(u,v))=(-ju/||(u,v)||,-jv/||(u,v)||)(3)
以参考图像r为例,对其进行一阶riesz变换,得到:
参考图像r的单演信号rm被定义为:
rm=[r,r1,r2]=[r,h1(x,y)*r,h2(x,y)*r]
=r-(i,j)fr(x,y)(5)
其中,i和j为虚数单位;
如图1所示,riesz变换空间是一个球形坐标系,r、r1、r2分别是空间中的点在球形坐标系三个轴上的投影;在riesz变换空间域上,局部幅值a、局部方向θ和局部相位
其中,
所述步骤s2中对分解后的图像进行riesz变换,得到5对单演信号,实现步骤如下:
第一步,对分解后的图像rbi(i=1,2,3,4,5)进行riesz变换,得到5个单演信号
其中,
第二步,对分解后的图像dbi(i=1,2,3,4,5)进行riesz变换,得到5个单演信号
其中,
本发明设计的log-gabor滤波器组如图4所示,从左到右的中心频率ω0依次为
从图4中可以看出,
步骤s3中对5对单演信号进行局部特征的相似度计算,得到5组单演相似度函数,进而得到5幅单演相似度图,具体为:
第一步,计算riesz变换空间下各点的特征量,包括幅值、相位角、方向角,构造参考图像和失真图像在像素点(x,y)的单演信号相似度:
其中,
第二步,构造参数smi,作为单演信号相似度图:
其中,i=1,2,3,4,5。
步骤s4中利用csf曲线设置5幅单演相似度图的权重,然后进行加权求和,得到1幅单演相似度图sm,具体为:
利用csf曲线设置5幅单演相似度图的权重ki(i=1,2,3,4,5),然后进行加权求和,得到1幅单演相似度图sm:
sm=k1sm1+k2sm2+k3sm3+k4sm4+k5sm5(11)
其中k1=0.3370,k2=0.8962,k3=0.9809,k4=0.9753,k5=0.7411。
本发明设计的log-gabor滤波器组将对应csf曲线分为5段,以滤波器的半功率点作为带宽界限,得到log-gabor滤波器组的5个频段,取每个频段的最大值作为该段单演相似度图的权重,得到k1=0.3370,k2=0.8962,k3=0.9809,k4=0.9753,k5=0.7411。
本发明采用mannos等人提出的csf模型,所采用的视觉对比敏感度函数(contrastsensitivityfunction,csf)模型为:
a(fr)≈2.6[0.0192+0.114fr]exp[-(0.114fr)1.1](12)
其中,
所述步骤s5计算参考图像r和失真图像d的梯度幅度相似度图sg,实现步骤如下:
第一步,计算参考图像r和失真图像d的各像素点的梯度幅度g;
利用scharr算子来获得像素点的水平梯度gh和垂直梯度gv,然后计算参考图像和失真图像各像素点的梯度幅度g;以参考图像r为例,梯度幅度g计算公式如(13)、(14)所示:
第二步,计算参考图像r梯度幅度图gr和失真图像梯度幅度图gd的梯度幅度相似度图sg:
其中,gr(x,y)和gd(x,y)分别为参考图像和失真图像在像素点(x,y)的梯度幅值,c2和c3为小的正常数,c2=87497.64,c3=65025。sg(x,y)的数值范围为(0,1],数值越小,表明梯度幅度失真越严重。当sg(x,y)=1时,说明参考图像和失真图像在该像素点的梯度幅度无失真。常数c3可以防止公式(11)出现奇异现象,c2和c3在调节低梯度区域的对比响应中起到重要作用。
所述步骤s6中利用参考图像r得到的单演信号进行单演相位一致性计算,得到池化函数mpc,实现步骤如下:
根据公式(7),可算出局部能量信息的总和为:
其中,
局部幅度的总和为:
则单演相位一致性模型(也就是池化函数)表示为:
其中,符号
其中,c是过滤响应滤波器的截止值,在此之下相位一致性将会受到抑制;g是一个控制截止明锐度的增益因子;
扩展函数s(x)为:
其中,amax(x,y)是滤波器在(x,y)处的最大响应的振幅,ε是一个小的正常数,其中,g=1.8182,
为了更直观的展示权重函数,图5以live库中的monarch为例,绘制了w(x,y)和s(x,y)的三维曲面图。图5中,从左到右依次为monarch图、w(x,y)三维曲面图和w(x,y)三维曲面旋转图。
所述步骤s7中利用池化函数mpc对单演相似度图sm和梯度幅度相似度图sg的乘积进行卷积,得到最终的图像质量评价指数,实现步骤如下:
联合公式(11)(15),同时利用单演相位一致性评价模型作为池化函数mpc,得到最终的质量评价分数:
rvsim=sm·sg*mpc(21)
公式(21)即为提出的图像质量评价指数。
常数c1/c2/c3的确定如下:
常数c1、c2和c3除了可以规避sabi(i=1,2,3,4,5)和sg的定义式出现奇异现象外,还起到了微调rvsim指数的作用。为了确定常数c1、c2和c3的最优取值,本发明利用评价指标srocc在live数据库上进行正交实验。为了在实验的复杂度和参数确定之间取得平衡,我们进行了两轮正交实验。类似于ssim模型,首先令[c1,c2,c3]=[(k1*l)2,(k2*l)2,(k3*l)2],l=255。
第一轮:第一步固定k2=k2=1.0,当k1取不同值时将rvsim算法应用于live数据库,得到k1-srocc曲线,如图6a所示,可见k1=1.0时srocc取得最大值。第二步固定k1=1.0,k3=1.0,当k2取不同值时将rvsim算法应用于live数据库,得到k2-srocc曲线,如图6b所示,可见k2=1.2时srocc取得最大值。第三步固定k1=1.0,k2=1.2,当k3取不同值时将rvsim算法应用于live数据库,得到k3-srocc曲线,如图6c所示,可见k3=1.0时srocc取得最大值。至此,第一轮实验结束,得到的参数为:k1=1.0,k2=1.2,k3=1.0。
第二轮:在第一轮实验得到参数的基础上,重复第一轮的实验,得到如图7d/e/f所示的结果。至此,第两轮实验结束,最终确定的参数为:k1=1.09,k2=1.16,k3=1.00。
本发明的优点可以通过以下实验具体说明:
为了验证rvsim算法的性能,本发明首先通过四个图像数据库live、csiq、tid2008和tid2013进行指标测试。live数据库包含29张参考图像和779张失真图像,失真类型包括jpeg2000、jpeg压缩、白噪声、高斯模糊和瑞利衰落,数据库中提供了每个图像的差分平均意见分数(differentialmeanopinionscores,dmos),较小的dmos值表示较高的图像质量。csiq数据库包含30张参考图像,每张参考图像具有6种失真类型,每种失真有4-5种程度。tid2008数据库包括25幅彩色参考图像,失真图像包括17种不同的失真,每个失真含4个不同的等级,共计1700幅图像。库中的17种失真分别为:加性高斯噪声、颜色分量强于照明分量的加性噪声、空间位置相关噪声、掩模噪声、高频噪声、脉冲噪声、量化噪声、高斯模糊、图像噪声、jpeg压缩、jpeg2000压缩、jpeg传输错误、jpeg2000传输错误、非偏心式噪声、不同强度的局部块失真、强度均值偏移和对比度变化。tid2013包括25幅参考图像和3000幅失真图像,图像的分辨率为384×512。所有图像均为彩色bmp格式。该数据库中共有24种类型的失真,每幅图像的每种失真有5个不同的强度。除tid2008中的17种失真外,增加了7种新的失真类型:色彩饱和度改变、乘性高斯噪声、舒适噪声、含噪声的有损压缩、带抖动的图像颜色量差、色差和稀疏采样与重构。
本发明使用5参数非线性逻辑回归函数拟合数据,使用四个相应的指数来客观地比较算法的性能:pearson线性相关系数plcc,均方根误差rmse,kendall秩序相关系数krocc和spearman秩相关系数srocc。
其中,x是客观iqa分数,f(x)是iqa回归拟合分数,β1、β2、β3、β4和β5是回归函数参数。
自然界中的图像是丰富多彩的,仅对这四个数据库的实验还不足以说明问题,为了以更客观公平的方式测试算法性能,本发明还在waterlooexploration数据库中开展了竞争排名(groupmaximumdifferentiationcompetition,gmad)。waterlooexploration数据库包含4744张原始自然图像和94880张失真图像。并且gmad竞争可以自动地从数据库中选择图像对的子集,其提供最强的测试以使iqa模型彼此竞争。对所选子集的主观测试揭示了iqa模型的相对性能。
性能比较如下:
表2列出了rvsim算法在live、csiq、tid2008、tid2013四个数据库上的表现。从表2可以看到,rvsim算法的iqa指标在所有4个数据库上都有非常好的表现,体现了rvsim的优异性能。
表2不同数据库中不同算法的srocc/krocc/plcc/rmse性能比较
图7示出了live数据库上的psnr、ssim、ms-ssim、iw-ssim、fsim和rvsim指数相对于dmos的预测质量得分散点图。从图7可以发现,rvsim的散点图在整个坐标系统中均匀分布,并且与dmos具有很强的线性关系,这表明rvsim模型与人眼具有很强的一致性。
图8给出了17个算法的gmad竞争排名情况,除rvsim外,其余16个iqa模型数据是官方对已有的知名模型测试得到的结果。从图8可见,rvsim的综合性能排名第一,而且远超第二名ms-ssim和第三名cornia。
本发明针对全参考图像质量评价问题,提出了一种结合riesz变换和视觉对比敏感度的图像质量评价方法(rvsim)。rvsim充分利用单演信号理论和log-gabor滤波器,结合csf特性曲线来分配不同频段的权重,同时引入梯度幅度相似度,得到相似度矩阵,然后利用mpc图来构造池化函数,得到最终的质量评价指数。rgsim算法首次将riesz解析空间中的局部特征,包括幅度值,局部相位和局部方向,应用于图像质量评价,取得了良好的效果。5个流行数据库上的测试结果表明,rgsim算法整体性能优异。
实施例2
本实施例2是系统实施例,上述实施例1是方法实施例,需要说明的是,本实施例2与上述实施例1属于同一技术构思,在本实施例中未详尽描述的内容,请参见方法实施例1.
如图9所示,一种采用所述的图像质量评价方法的系统,包括:
分解单元,用于利用log-gabor滤波器分别对参考图像r和失真图像d进行5尺度分解,
riesz变换单元,用于基于单演信号理论对分解后的图像进行riesz变换,得到5对单演信号;
第一计算单元,用于对5对单演信号进行局部特征的相似度计算,得到5组单演相似度函数,进而得到5幅单演相似度图;
第二计算单元,用于利用csf曲线设置5幅单演相似度图的权重,然后进行加权求和,得到1幅单演相似度图sm;
第三计算单元,用于计算参考图像r和失真图像d的梯度幅度相似度图sg;
第四计算单元,用于利用riesz变换单元得到的参考图像r的单演信号进行单演相位一致性计算,得到池化函数;
图像质量评价指数获取单元,用于最后利用池化函数对单演相似度图sm和梯度幅度相似度图sg的乘积进行卷积,得到最终的图像质量评价指数。
本领域的技术人员容易理解,以上所述仅为本发明较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所做的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。