发明名称 一种多视觉敏感特征相似度的全参考图像质量评价方法
摘要 本发明涉及一种基于多视觉敏感特征相似度的全参考型图像质量评价方法。本发明的具体步骤是:(一)输入失真图像和参考图像;(二)对失真图像和参考图像提取多角度的梯度特征,并计算梯度相似度;(三)对失真图像和参考图像提取相位特征,并计算相位特征相似度;(四)对失真图像和参考图像提取锐度特征,并计算锐度特征相似度;(五)将梯度相似度、相位相似度和锐度相似度融合成最终的图像质量评价指标。本发明采用对图像提取梯度、相位和锐度等视觉敏感特征,对视觉敏感特征相似度组合得到客观评价参数,提高了全参考图像质量评价的精度。
申请公布号 CN104504676A 申请公布日期 2015.04.08
申请号 CN201410650125.9 申请日期 2014.11.07
申请人 嘉兴学院 发明人 汪斌
分类号 G06T7/00(2006.01)I;G06K9/62(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 代理人
主权项 一种基于多视觉敏感特征相似度的全参考图像质量评价方法,其特征在于包含以下步骤:(1)输入参考图像和失真图像,对失真图像和参考图像提取多角度下的梯度特征。具体步骤如下:对失真图像和参考图像分别采用0°,45°,90°,和135°方向的Sobel算子提取出梯度信息,采用如下公式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>S</mi><mn>0</mn></msub><mo>=</mo><mi>R</mi><mo>&CircleTimes;</mo><msub><mi>P</mi><mn>0</mn></msub><mo>;</mo><msub><mi>S</mi><mn>45</mn></msub><mo>=</mo><mi>R</mi><mo>&CircleTimes;</mo><msub><mi>P</mi><mn>45</mn></msub><mo>;</mo><msub><mi>S</mi><mn>90</mn></msub><mo>=</mo><mi>R</mi><mo>&CircleTimes;</mo><msub><mi>P</mi><mn>90</mn></msub><mo>;</mo><msub><mi>S</mi><mn>135</mn></msub><mo>=</mo><mi>R</mi><mo>&CircleTimes;</mo><msub><mi>P</mi><mn>135</mn></msub><mo>;</mo></mrow>]]></math><img file="FSA00001104032600000110.GIF" wi="1341" he="84" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>T</mi><mn>0</mn></msub><mo>=</mo><mi>D</mi><mo>&CircleTimes;</mo><msub><mi>P</mi><mn>0</mn></msub><mo>;</mo><msub><mi>T</mi><mn>45</mn></msub><mo>=</mo><mi>D</mi><mo>&CircleTimes;</mo><msub><mi>P</mi><mn>45</mn></msub><mo>;</mo><msub><mi>T</mi><mn>90</mn></msub><mo>=</mo><mi>D</mi><mo>&CircleTimes;</mo><msub><mi>P</mi><mn>90</mn></msub><mo>;</mo><msub><mi>T</mi><mn>135</mn></msub><mo>=</mo><mi>D</mi><mo>&CircleTimes;</mo><msub><mi>P</mi><mn>135</mn></msub><mo>;</mo></mrow>]]></math><img file="FSA00001104032600000111.GIF" wi="1344" he="80" /></maths>其中D为失真图像,S为参考图像;P<sub>0</sub>、P<sub>45</sub>、P<sub>90</sub>和P<sub>135</sub>为0°,45°,90°,和135°方向的sobel算子;<img file="FSA00001104032600000112.GIF" wi="46" he="61" />表示二维卷积运算;S<sub>0</sub>、S<sub>45</sub>、S<sub>90</sub>和S<sub>135</sub>分别是参考图像R经过0°,45°,90°,和135°方向的Sobel算子P<sub>0</sub>、P<sub>45</sub>、P<sub>90</sub>和P<sub>135</sub>卷积滤波后得到的图像;T<sub>0</sub>、T<sub>45</sub>、T<sub>90</sub>和T<sub>135</sub>分别是参考图像D经过0°,45°,90°,和135°方向的Sobel算子P<sub>0</sub>、P<sub>45</sub>、P<sub>90</sub>和P<sub>135</sub>卷积滤波后得到的图像;(2)计算失真图像和参考图像提取的多角度梯度特征之间的相似度,采用如下公式:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>G</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mrow><mn>2</mn><mi>S</mi></mrow><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><msub><mi>T</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mi>c</mi></mrow><mrow><msub><mi>S</mi><mn>0</mn></msub><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>T</mi><mn>0</mn></msub><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>c</mi></mrow></mfrac></mrow>]]></math><img file="FSA0000110403260000011.GIF" wi="679" he="160" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>G</mi><mn>45</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mrow><mn>2</mn><mi>S</mi></mrow><mn>45</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><msub><mi>T</mi><mn>45</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mi>c</mi></mrow><mrow><msub><mi>S</mi><mn>45</mn></msub><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>T</mi><mn>45</mn></msub><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>c</mi></mrow></mfrac></mrow>]]></math><img file="FSA0000110403260000012.GIF" wi="721" he="155" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>G</mi><mn>90</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mrow><mn>2</mn><mi>S</mi></mrow><mn>90</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><msub><mi>T</mi><mn>90</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mi>c</mi></mrow><mrow><msub><mi>S</mi><mn>90</mn></msub><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>T</mi><mn>90</mn></msub><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>c</mi></mrow></mfrac></mrow>]]></math><img file="FSA0000110403260000013.GIF" wi="724" he="160" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>G</mi><mn>135</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mrow><mn>2</mn><mi>S</mi></mrow><mn>135</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><msub><mi>T</mi><mn>135</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mi>c</mi></mrow><mrow><msub><mi>S</mi><mn>135</mn></msub><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>T</mi><mn>135</mn></msub><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>c</mi></mrow></mfrac></mrow>]]></math><img file="FSA0000110403260000014.GIF" wi="873" he="160" /></maths>(3)计算失真图像和参考图像提取的多角度梯度特征之间的相似度。采用如下公式:将0°,45°,90°,和135°方向的梯度相似度合成最终的合成梯度相似度,采用公式如下:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>[</mo><msub><mi>G</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo></mrow><msub><mi>&alpha;</mi><mn>1</mn></msub></msup><msup><mrow><mo>[</mo><msub><mi>G</mi><mn>45</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo></mrow><msub><mi>&alpha;</mi><mn>2</mn></msub></msup><msup><mrow><mo>[</mo><msub><mi>G</mi><mn>90</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo></mrow><msub><mi>&alpha;</mi><mn>3</mn></msub></msup><msup><mrow><mo>[</mo><msub><mi>G</mi><mn>135</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo></mrow><msub><mi>&alpha;</mi><mn>4</mn></msub></msup></mrow>]]></math><img file="FSA0000110403260000015.GIF" wi="1182" he="91" /></maths>(4)对输入失真图像和参考图像提取相位特征,具体采用步骤和公式如下:将失真图像D和参考图像S伸缩成像素为横向和纵向像素大小相等的图像D′和S′,对图像D′和S′进行二维离散傅里叶变换,并分别得到相位,具体采用公式如下:<img file="FSA0000110403260000016.GIF" wi="594" he="160" /><img file="FSA0000110403260000017.GIF" wi="615" he="160" />其中Y<sub>1</sub>(u,v)是图像D′的二维离散傅里叶变换系数;Y<sub>2</sub>(u,v)是图像S′的二维离散傅里叶变换系数;<img file="FSA0000110403260000018.GIF" wi="167" he="85" />和<img file="FSA0000110403260000019.GIF" wi="170" he="90" />是Y<sub>1</sub>(u,v)和Y<sub>2</sub>(u,v)的相位;re(·)表示取实部运算,im(·)表示取虚部运算;(5)计算图像D′和图像S′之间的相位相似度信息,计算公式如下:<img file="FSA0000110403260000021.GIF" wi="688" he="159" />最后采用二维傅立叶反变换将相位相似度信息转换到空域,并调整其大小到原始尺寸,得到空域相位相似度P(i,j),其中(i,j)为图像中某一像素的位置;(6)对输入失真图像和参考图像提取锐度特征,具体采用步骤和公式如下:首先,计算失真图像和参考图像(i,j)像素为中心的3×3邻域的平均灰度值,采用计算公式如下:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mover><mrow><mi>D</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mn>9</mn></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>u</mi><mo>=</mo><mo>-</mo><mn>1</mn></mrow><mn>1</mn></munderover><munderover><mi>&Sigma;</mi><mrow><mi>v</mi><mo>=</mo><mo>-</mo><mn>1</mn></mrow><mn>1</mn></munderover><mi>D</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mi>u</mi><mo>,</mo><mi>j</mi><mo>+</mo><mi>v</mi><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000110403260000022.GIF" wi="783" he="160" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mover><mrow><mi>R</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mn>9</mn></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>u</mi><mo>=</mo><mo>-</mo><mn>1</mn></mrow><mn>1</mn></munderover><munderover><mi>&Sigma;</mi><mrow><mi>v</mi><mo>=</mo><mo>-</mo><mn>1</mn></mrow><mn>1</mn></munderover><mi>R</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mi>u</mi><mo>,</mo><mi>j</mi><mo>+</mo><mi>v</mi><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000110403260000023.GIF" wi="783" he="170" /></maths>然后,计算失真图像和参考图像(i,j)像素处的最大邻域灰度差值S(i,j)和T(i,j),以最大领域灰度差值作为图像的锐度亮度,采用公式如下:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mi>S</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>max</mi><mo>|</mo><mi>D</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mi>u</mi><mo>,</mo><mi>j</mi><mo>+</mo><mi>v</mi><mo>)</mo></mrow><mo>-</mo><mover><mrow><mi>D</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mo>|</mo><mo>,</mo><mo>&ForAll;</mo><mo>-</mo><mn>1</mn><mo>&le;</mo><mi>u</mi><mo>&le;</mo><mn>1</mn><mo>,</mo><mo>-</mo><mn>1</mn><mo>&le;</mo><mi>v</mi><mo>&le;</mo><mn>1</mn></mrow>]]></math><img file="FSA0000110403260000024.GIF" wi="1277" he="85" /></maths><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mi>T</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>max</mi><mo>|</mo><mi>R</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mi>u</mi><mo>,</mo><mi>j</mi><mo>+</mo><mi>v</mi><mo>)</mo></mrow><mo>-</mo><mover><mrow><mi>R</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mo>|</mo><mo>,</mo><mo>&ForAll;</mo><mo>-</mo><mn>1</mn><mo>&le;</mo><mi>u</mi><mo>&le;</mo><mn>1</mn><mo>,</mo><mo>-</mo><mn>1</mn><mo>&le;</mo><mi>v</mi><mo>&le;</mo><mn>1</mn></mrow>]]></math><img file="FSA0000110403260000025.GIF" wi="1259" he="95" /></maths>(7)对输入失真图像和参考图像之间计算锐度相似度,具体采用步骤和公式如下:<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>2</mn><mi>S</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mi>T</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mi>c</mi></mrow><mrow><mi>S</mi><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>T</mi><msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>c</mi></mrow></mfrac></mrow>]]></math><img file="FSA0000110403260000026.GIF" wi="640" he="163" /></maths>(8)组合梯度相似度、相位相似度和锐度相似度组成联合相似度判据,求其均值得到最终的相似度数值,具体公式如下:GPS_SSIM(i,j)=[G(i,j)]<sup>α</sup>[P(i,j)]<sup>β</sup>[H(i,j)]<sup>γ</sup><maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><mi>GPS</mi><mo>_</mo><mi>SSIM</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mi>M</mi><mo>&times;</mo><mi>N</mi></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>GPS</mi><mo>_</mo><mi>SSIM</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000110403260000027.GIF" wi="979" he="234" /></maths>
地址 314001 浙江省嘉兴市越秀南路56号嘉兴学院