发明名称 基于DCT域的彩色图像可觉察失真度计算方法
摘要 本发明涉及视频图像通信技术领域,尤其涉及一种基于DCT域的彩色图像可觉察失真度计算方法。本发明将选定的彩色图像由RGB色彩空间转换到Lab色彩空间,对Lab色彩空间下的L、a、b三个分量分别进行8×8的DCT变换,将其由空域变换到DCT域,在DCT域,根据亮度基本阈值和调制因子的乘积计算获得亮度分量L的可觉察失真度JND值,在DCT域,根据色度基本阈值和调制因子的乘积计算获得色度分量a、b的可觉察失真度JND值。本发明彩色图像可觉察失真度计算方法整体上具有更低的峰值信噪比(PSNR)值,其能容忍更多的客观失真。
申请公布号 CN101621708B 申请公布日期 2011.01.26
申请号 CN200910063359.2 申请日期 2009.07.29
申请人 武汉大学 发明人 胡瑞敏;陈皓;胡金晖;毛丹;王师峥;汪欢;张岿;钟睿
分类号 H04N17/00(2006.01)I 主分类号 H04N17/00(2006.01)I
代理机构 武汉华旭知识产权事务所 42214 代理人 周宗贵
主权项 1.一种基于DCT域的彩色图像可觉察失真度计算方法,其特征在于,包括以下步骤:①将选定的彩色图像由RGB色彩空间转换到Lab色彩空间;②对Lab色彩空间下的L、a、b三个分量分别进行8×8的DCT变换,将其由空域变换到DCT域;③在DCT域,根据亮度基本阈值和调制因子的乘积计算获得亮度分量L的可觉察失真度JND值;④在DCT域,根据色度基本阈值和调制因子的乘积计算获得色度分量a、b的可觉察失真度JND值;步骤③中所述亮度分量L的可觉察失真度JND值采用如下公式计算获得:JND<sub>lum</sub>(n,i,j)=T<sub>lum_basic</sub>(n,i,j)·F<sub>lum_M</sub>(n,i,j)其中,T<sub>lum_basic</sub>(n,i,j)为亮度基本阈值,F<sub>lum_M</sub>(n,i,j)为调制因子,n为DCT块的索引,<img file="FSB00000312081900011.GIF" wi="310" he="103" />W为图像的宽度,H为图像的高度,i、j为DCT块中系数的索引,1≤i≤8,1≤j≤8;所述亮度基本阈值采用如下公式计算获得:T<sub>lum_basic</sub>(n,i,j)=s·T(n,i,j)其中,<maths num="0001"><![CDATA[<math><mrow><mi>T</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>&phi;</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>&phi;</mi><mi>j</mi></msub></mrow></mfrac><mo>&CenterDot;</mo><mfrac><mrow><mi>exp</mi><mrow><mo>(</mo><msub><mi>cw</mi><mi>ij</mi></msub><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><mi>a</mi><mo>+</mo><msub><mi>bw</mi><mi>ij</mi></msub><mo>)</mo></mrow></mrow><mrow><mi>r</mi><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>r</mi><mo>)</mo></mrow><msup><mi>cos</mi><mn>2</mn></msup><msub><mi>&upsi;</mi><mi>ij</mi></msub></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>w</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>N</mi></mrow></mfrac><msqrt><msup><mrow><mo>(</mo><mi>i</mi><mo>/</mo><msub><mi>&theta;</mi><mi>x</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>j</mi><mo>/</mo><msub><mi>&theta;</mi><mi>y</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>&phi;</mi><mi>m</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msqrt><mn>1</mn><mo>/</mo><mi>N</mi></msqrt></mtd><mtd><mi>m</mi><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msqrt><mn>2</mn><mo>/</mo><mi>N</mi></msqrt></mtd><mtd><mi>m</mi><mo>></mo><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>&upsi;</mi><mi>ij</mi></msub><mo>=</mo><mi>arcsin</mi><mrow><mo>(</mo><mn>2</mn><msub><mi>w</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><msub><mi>w</mi><mrow><mn>0</mn><mi>j</mi></mrow></msub><mo>/</mo><msubsup><mi>w</mi><mi>ij</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><img file="FSB00000312081900016.GIF" wi="654" he="140" />w<sub>i,j</sub>表示DCT子带的空间频率,R<sub>vd</sub>表示观察距离和图像高度的比例,3≤R<sub>vd</sub>≤6,P<sub>ich</sub>表示在图像高度上的像素点的个数,N为DCT的块大小,N=8,参数a=1.33,b=0.11,c=0.005,s=1.1,r=0.6;所述调制因子采用如下公式计算获得:F<sub>lum_M</sub>(n,i,j)=F<sub>lum</sub>(n,i,j)·F<sub>contrast</sub>(n,i,j)其中亮度自适应因子F<sub>lum</sub>的计算公式如下:<maths num="0005"><![CDATA[<math><mrow><msub><mi>F</mi><mi>lum</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>+</mo><mrow><mo>(</mo><mn>60</mn><mo>-</mo><mover><mi>I</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>/</mo><mn>150</mn></mtd><mtd><mover><mi>I</mi><mo>&OverBar;</mo></mover><mo>&le;</mo><mn>60</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>60</mn><mo>&lt;</mo><mover><mi>I</mi><mo>&OverBar;</mo></mover><mo>&lt;</mo><mn>170</mn></mtd></mtr><mtr><mtd><mn>1</mn><mo>+</mo><mrow><mo>(</mo><mover><mi>I</mi><mo>&OverBar;</mo></mover><mo>-</mo><mn>170</mn><mo>)</mo></mrow><mo>/</mo><mn>425</mn></mtd><mtd><mover><mi>I</mi><mo>&OverBar;</mo></mover><mo>&GreaterEqual;</mo><mn>170</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中,<img file="FSB00000312081900022.GIF" wi="28" he="55" />为DCT块中的平均亮度值;其中,亮度对比掩蔽因子F<sub>contrast</sub>的计算公式如下:对于块中的高频系数:<img file="FSB00000312081900023.GIF" wi="1655" he="238" />对于块中的低频系数:<img file="FSB00000312081900024.GIF" wi="1768" he="255" />ε=1.9,s<sub>1</sub>=0.9,s<sub>2</sub>=0.7;其中,<img file="FSB00000312081900025.GIF" wi="32" he="37" />的定义如下:将图像中的所有8×8大小DCT块进行分类,分成纹理、平坦和边缘三类;用canny算子来检测计算图像中的边缘信息,通过计算每个8×8块中的边缘像素密度ρ<sub>edge</sub>来对块进行分类,具体过程如下:<maths num="0006"><![CDATA[<math><mrow><msub><mi>&rho;</mi><mi>edge</mi></msub><mo>=</mo><mrow><mo>(</mo><munder><mi>&Sigma;</mi><mrow><mn>8</mn><mo>&times;</mo><mn>8</mn></mrow></munder><mi>edge</mi><mo>)</mo></mrow><mo>/</mo><msup><mi>N</mi><mn>2</mn></msup></mrow></math>]]></maths>其中,<img file="FSB00000312081900027.GIF" wi="154" he="95" />表示每个8×8块中通过canny算子得到的边缘像素的总数;N为DCT块大小,取N=8;<img file="FSB00000312081900028.GIF" wi="787" he="226" />由8×8的DCT块类型可知,块间掩蔽效应因子<img file="FSB00000312081900029.GIF" wi="30" he="38" />为:<img file="FSB00000312081900031.GIF" wi="847" he="212" />步骤④中色度分量a、b的可觉察失真度JND值采用如下公式计算获得:JND<sub>colar</sub>(n,i,j)=T<sub>color_basic</sub>(n,i,j)·F<sub>color_M</sub>(n,i,j)其中,T<sub>color_basic</sub>(n,i,j)为色度基本阈值,F<sub>color_M</sub>(n,i,j)为调制因子,n为DCT块的索引,i、j为DCT块中系数的索引。
地址 430072 湖北省武汉市武昌珞珈山