发明名称 基于结构失真与空间频率指标的图像质量评价方法
摘要 本发明公开了一种基于结构失真与空间频率指标的图像质量评价方法,该方法在基于结构失真的图像质量评价方法的基础上,提出结构失真与空间频率指标结合的图像质量评价方法,从而实现了从影响图像质量评价的主要因素即亮度、清晰度以及相关度对图像质量进行评价。本发明基于结构失真与空间频率指标的图像质量评价方法在客观评价图像时更能与人眼主观感觉相一致,并且由于所选用的计算关系式较为简单,只涉及到原始图像与失真图像像素之间的计算,所以该方法能够快速有效的对图像质量进行评价。
申请公布号 CN102663764A 申请公布日期 2012.09.12
申请号 CN201210125400.6 申请日期 2012.04.25
申请人 武汉大学 发明人 易尧华;刘菊华;苏海;李帅
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 薛玲
主权项 1.一种基于结构失真与空间频率指标的图像质量评价方法,其特征是在于,包括以下步骤:步骤1,计算原始图像与失真图像的均值<img file="FDA0000157303510000011.GIF" wi="99" he="44" />利用关系式<img file="FDA0000157303510000012.GIF" wi="158" he="117" />计算原始图像与失真图像的平均亮度逼近度,其中,x,y分别代表原始图像与失真图像,图像中横向和纵向的像素数目分别为m和n,<maths num="0001"><![CDATA[<math><mrow><mover><mi>x</mi><mo>&OverBar;</mo></mover><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>x</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mover><mi>y</mi><mo>&OverBar;</mo></mover><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>y</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></math>]]></maths>x(i,j)和y(i,j)分别为原始图像和失真图像在第i行第j列的灰度值;步骤2,计算原始图像与失真图像的均方误差σ<sub>x</sub>,σ<sub>y</sub>及协方差σ<sub>xy</sub>,利用关系式<img file="FDA0000157303510000015.GIF" wi="94" he="101" />计算原始图像与失真图像的线性相关度,其中,x,y分别代表原始图像与失真图像,图像在横向和纵向的像素数目分别为m和n,<maths num="0003"><![CDATA[<math><mrow><msub><mi>&sigma;</mi><mi>x</mi></msub><mo>=</mo><msup><mrow><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><msup><mrow><mo>(</mo><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>&sigma;</mi><mi>y</mi></msub><mo>=</mo><msup><mrow><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><msup><mrow><mo>(</mo><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>&sigma;</mi><mi>xy</mi></msub><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><mrow><mo>(</mo><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mrow><mo>(</mo><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow></math>]]></maths>x(i,j)和y(i,j)分别为原始图像和失真图像在第i行第j列的灰度值;步骤3,根据步骤2所得原始图像与失真图像的均方误差σ<sub>x</sub>,σ<sub>y</sub>,利用关系式<img file="FDA0000157303510000019.GIF" wi="181" he="111" />计算原始图像与失真图像的相似度;步骤4,计算失真图像的空间行频率RF以及空间列频率CF,通过空间行频率RF及空间列频率CF计算失真图像空间频率指标<img file="FDA00001573035100000110.GIF" wi="404" he="62" />计算空间行频率RF和空间列频率CF的公式如下<maths num="0006"><![CDATA[<math><mrow><mi>RF</mi><mo>=</mo><msqrt><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>2</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><mi>CF</mi><mo>=</mo><msqrt><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>2</mn></mrow><mi>m</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths>其中,m,n分别为图像在横向和纵向的像素数目,y(i,j)为失真图像在第i行第j列的灰度值;步骤5,综合步骤1所得平均亮度逼近度、步骤2所得线性相关度、步骤3所得相似度及步骤4所得失真图像空间频率指标,计算基于结构失真与空间频率指标的图像质量评价指标<img file="FDA0000157303510000023.GIF" wi="763" he="114" />根据计算所得图像质量评价指标数值从图像的亮度、清晰度以及相关性三方面对图像质量进行综合评价。
地址 430072 湖北省武汉市武昌珞珈山武汉大学