发明名称 基于GSSIM的立体图像客观质量评价算法
摘要 本发明属图像处理领域,涉及一种基于GSSIM的立体图像质量评价方法,包括(1)对于左图像和右图像,分别求取梯度结构相似度值,再求取两值的均值,得到立体图像质量评价值QE;(2)采用下列方法进行图像立体感觉客观评价:计算原始图像和处理后图像的绝对差值图像;求取出绝对差值图像的μ1,μ2;用Sobel算子求出绝对差值图像的梯度幅值图像;对原始图像的绝对差值图像进行滤波,计算原图双眼视差分布情况;求出双眼视差处的Dl(x,y),Dcg(x,y)与Dsg(x,y)值;计算出双眼视差处的DSSIM值;计算双眼视差处图像梯度结构相似度,即图像立体感觉客观评价值DE。本发明能很好地用于立体图像的质量评价,客观评价结果与主观评价结果相关性很强。
申请公布号 CN101833766A 申请公布日期 2010.09.15
申请号 CN201010169050.4 申请日期 2010.05.11
申请人 天津大学 发明人 杨嘉琛;王斌;韦娜;范超伟;武强一;李杰
分类号 G06T7/00(2006.01)I;G06T7/40(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 程毓英
主权项 1.一种基于GSSIM的立体图像质量评价方法,包括下列步骤:(1)对于左图像和右图像,分别利用梯度结构相似度GSSIM算法,求取梯度结构相似度值,再求取两值的均值,得到立体图像质量评价值QE;(2)采用下列方法进行图像立体感觉客观评价:第一步:把原始图像的视点对(L1,R1)和处理后图像的视点对(L2,R2)分别进行相减操作,得出绝对差值图像X1与X2;第二步,进行亮度测试,通过公式下列求出绝对差值图像X1与X2的μ<sub>1</sub>,μ<sub>2</sub>:<img file="FDA0000021233190000011.GIF" wi="562" he="198" /><img file="FDA0000021233190000012.GIF" wi="490" he="198" />式中,w为中心对称高斯加权窗<img file="FDA0000021233190000013.GIF" wi="501" he="135" />N取值为11;第三步,用Sobel算子求出绝对差值图像X1与X2的梯度幅值图像M1与M2,并求出图像M1与M2亮度均值μ<sub>1′</sub>,μ<sub>2′</sub>,再通过公式下列求出标准差σ<sub>1′</sub>,σ<sub>2′</sub>,以及协方差σ<sub>1′2′</sub>:<maths num="0001"><![CDATA[<math><mrow><msub><mi>&sigma;</mi><msubsup><mn>1</mn><mi>mn</mi><mo>&prime;</mo></msubsup></msub><mo>=</mo><msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>m</mi><mo>-</mo><mfrac><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow><mrow><mi>m</mi><mo>+</mo><mfrac><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>n</mi><mo>-</mo><mfrac><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow><mrow><mi>n</mi><mo>+</mo><mfrac><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow></munderover><msub><mi>w</mi><mi>ij</mi></msub><msup><mrow><mo>(</mo><msub><mi>M</mi><mrow><mn>1</mn><mi>ij</mi></mrow></msub><mo>-</mo><msub><mi>&mu;</mi><msubsup><mn>1</mn><mi>mn</mi><mo>&prime;</mo></msubsup></msub><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="0002"><![CDATA[<math><mrow><msub><mi>&sigma;</mi><msubsup><mn>2</mn><mi>mn</mi><mo>&prime;</mo></msubsup></msub><mo>=</mo><msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>m</mi><mo>-</mo><mfrac><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow><mrow><mi>m</mi><mo>+</mo><mfrac><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>n</mi><mo>-</mo><mfrac><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow><mrow><mi>n</mi><mo>+</mo><mfrac><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow></munderover><msub><mi>w</mi><mi>ij</mi></msub><msup><mrow><mo>(</mo><msub><mi>M</mi><mrow><mn>2</mn><mi>ij</mi></mrow></msub><mo>-</mo><msub><mi>&mu;</mi><msubsup><mn>2</mn><mi>mn</mi><mo>&prime;</mo></msubsup></msub><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="0003"><![CDATA[<math><mrow><msub><mi>&sigma;</mi><mrow><msup><mn>1</mn><mo>&prime;</mo></msup><msubsup><mn>2</mn><mi>mn</mi><mo>&prime;</mo></msubsup></mrow></msub><mo>=</mo><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>m</mi><mo>-</mo><mfrac><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow><mrow><mi>m</mi><mo>+</mo><mfrac><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>n</mi><mo>-</mo><mfrac><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow><mrow><mi>n</mi><mo>+</mo><mfrac><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow></munderover></mrow><msub><mi>w</mi><mi>ij</mi></msub><mrow><mo>(</mo><msub><mi>M</mi><mrow><mn>1</mn><mi>ij</mi></mrow></msub><mo>-</mo><msub><mi>&mu;</mi><msubsup><mn>1</mn><mi>mn</mi><mo>&prime;</mo></msubsup></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>M</mi><mrow><mn>2</mn><mi>ij</mi></mrow></msub><mo>-</mo><msub><mi>&mu;</mi><msubsup><mn>2</mn><mi>mn</mi><mo>&prime;</mo></msubsup></msub><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>第四步,对原始图像的绝对差值图像X1进行滤波、门限值判决操作,计算原图双眼视差分布情况;第五步,用μ<sub>1</sub>,μ<sub>2</sub>,σ<sub>1′</sub>,σ<sub>2′</sub>与σ<sub>1′2′</sub>,结合计算出的原图视差分布情况,求出双眼视差处的D1(x,y),Dcg(x,y)与Dsg(x,y)值,计算公式如下:<maths num="0004"><![CDATA[<math><mrow><mi>D</mi><mn>1</mn><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>2</mn><msub><mi>&mu;</mi><mn>1</mn></msub><msub><mi>&mu;</mi><mn>2</mn></msub><mo>+</mo><msub><mi>C</mi><mn>1</mn></msub></mrow><mrow><msubsup><mi>&mu;</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&mu;</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msub><mi>C</mi><mn>1</mn></msub></mrow></mfrac></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>Dc</mi><mi>g</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>2</mn><msub><mi>&sigma;</mi><msup><mn>1</mn><mo>&prime;</mo></msup></msub><msub><mi>&sigma;</mi><msup><mn>2</mn><mo>&prime;</mo></msup></msub><mo>+</mo><msub><mi>C</mi><mn>2</mn></msub></mrow><mrow><msubsup><mi>&sigma;</mi><msup><mn>1</mn><mo>&prime;</mo></msup><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><msup><mn>2</mn><mo>&prime;</mo></msup><mn>2</mn></msubsup><mo>+</mo><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>Ds</mi><mi>g</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>&sigma;</mi><mrow><msup><mn>1</mn><mo>&prime;</mo></msup><msup><mn>2</mn><mo>&prime;</mo></msup></mrow></msub><mo>+</mo><msub><mi>C</mi><mn>3</mn></msub></mrow><mrow><msub><mi>&sigma;</mi><msup><mn>1</mn><mo>&prime;</mo></msup></msub><msub><mi>&sigma;</mi><msup><mn>2</mn><mo>&prime;</mo></msup></msub><mo>+</mo><msub><mi>C</mi><mn>3</mn></msub></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>第六步,通过公式DSSIM(x,y)=[D1(x,y)]<sup>α</sup>[Dc<sub>g</sub>(x,y)]<sup>β</sup>[Ds<sub>g</sub>(x,y)]<sup>γ</sup>(3-34)计算出双眼视差处的DSSIM值,其中α,β,γ取值都为1;第七步,通过DSSIM值,利用公式下列计算双眼视差处图像梯度结构相似度DMSSIM,即图像立体感觉客观评价值DE:<img file="FDA0000021233190000023.GIF" wi="870" he="120" />第八步,对立体图像质量进行评价:QE值越大,图像质量越好,DE值越大,立体感越好。
地址 300072 天津市南开区卫津路92号天津大学