发明名称 基于压缩感知的图像解码方法
摘要 基于压缩感知的图像解码方法,是在解码端对反量化得到的图像信号进行压缩感知的重构,可以通过求解下式范数的优化问题对图像信号进行压缩感知的重构:然后将解出的列向量转化成矩阵实现图像的解码。为了提升CS重构的质量,在对图像信号进行压缩感知的重构前,先对图像进行块合并,即将p×p个图像块合并成一个图像块,合并后图像块的行数/列数是合并前图像块的行数/列数的p倍。本发明还进一步改进TV算子在块边缘的值,定义图像块矩阵In×n在i=n时的水平算子为Ii-1,j-Iij,在j=n时的垂直算子为Ii,j-1-Iij。本发明所有的改进都集中在解码端,编码端不需要作任何改动,相比现有的图像压缩标准,能够获得更好的效果。
申请公布号 CN101931814B 申请公布日期 2012.11.21
申请号 CN201010271764.6 申请日期 2010.09.03
申请人 北京工业大学 发明人 尹宝才;施云惠;张臻;李倩
分类号 H04N7/26(2006.01)I;H04N7/30(2006.01)I 主分类号 H04N7/26(2006.01)I
代理机构 北京中北知识产权代理有限公司 11253 代理人 卢业强
主权项 1.基于压缩感知的图像解码方法,其特征在于,在解码端对反量化得到的图像信号进行压缩感知的重构;具体通过求解下式范数的优化问题对图像信号进行压缩感知的重构:<maths num="0001"><![CDATA[<math><mrow><mover><mi>x</mi><mo>~</mo></mover><mo>=</mo><mi>arg</mi><mi>min</mi><msub><mrow><mo>|</mo><mo>|</mo><msup><mi>&Psi;</mi><mi>T</mi></msup><mi>x</mi><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub><mo>,</mo><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><msub><mrow><mo>|</mo><mo>|</mo><mi>&Phi;x</mi><mo>-</mo><mi>y</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub><mo>&le;</mo><mi>&epsiv;</mi><mo>,</mo></mrow></math>]]></maths>上式表示在满足||Φx-y||<sub>2</sub>≤ε条件下<img file="FSB00000828372900012.GIF" wi="27" he="36" />取使||ψ<sup>T</sup>x||<sub>1</sub>最小的x值;式中,x表示图像块矩阵I<sub>n×n</sub>经过列扫描后得到的N×1维向量,<img file="FSB00000828372900013.GIF" wi="26" he="37" />表示x的重构,ψ表示信号x的稀疏基,使||ψ<sup>T</sup>x||<sub>1</sub>稀疏,ε表示由量化噪声引起的误差,y表示观测值,Φ表示标准观测矩阵;然后将解出的列向量<img file="FSB00000828372900014.GIF" wi="26" he="38" />转化成矩阵<img file="FSB00000828372900015.GIF" wi="51" he="51" />实现图像的解码;其中图像信号x采用梯度稀疏基使TV(x)稀疏,通过求解<maths num="0002"><![CDATA[<math><mrow><mover><mi>x</mi><mo>~</mo></mover><mo>=</mo><mi>arg</mi><mi>min</mi><mi>TV</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>,</mo><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><msub><mrow><mo>|</mo><mo>|</mo><mi>Ax</mi><mo>-</mo><mi>y</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub><mo>&le;</mo><mi>&epsiv;</mi></mrow></math>]]></maths>实现图像信号x的重构;上式表示在满足||Ax-y||<sub>2</sub>≤ε条件下<img file="FSB00000828372900017.GIF" wi="25" he="37" />取使TV(x)最小的x值,其中,x表示图像块矩阵I<sub>n×n</sub>经过列扫描后得到的N×1维向量,<img file="FSB00000828372900018.GIF" wi="25" he="37" />表示x的重构,TV(x)表示x在梯度意义下稀疏,ε表示由量化噪声引起的误差,y表示观测值,A表示观测矩阵;式中,全变分:<maths num="0003"><![CDATA[<math><mrow><mi>TV</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>ij</mi></munder><msqrt><msup><mrow><mo>(</mo><msub><mi>D</mi><mrow><mi>h</mi><mo>;</mo><mi>ij</mi></mrow></msub><mi>x</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>D</mi><mrow><mi>v</mi><mo>;</mo><mi>ij</mi></mrow></msub><mi>x</mi><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>=</mo><munder><mi>&Sigma;</mi><mi>ij</mi></munder><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>D</mi><mi>ij</mi></msub><mi>x</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub></mrow></math>]]></maths>其中D<sub>h;ij</sub>x为水平算子,D<sub>v;ij</sub>x为垂直算子;以n×n的图像块矩阵I表示,全变分:<maths num="0004"><![CDATA[<math><mrow><mi>TV</mi><mrow><mo>(</mo><mi>I</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>ij</mi></munder><msqrt><msup><mrow><mo>(</mo><msub><mi>D</mi><mrow><mi>h</mi><mo>;</mo><mi>ij</mi></mrow></msub><mi>I</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>D</mi><mrow><mi>v</mi><mo>;</mo><mi>ij</mi></mrow></msub><mi>I</mi><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>=</mo><munder><mi>&Sigma;</mi><mi>ij</mi></munder><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>D</mi><mi>ij</mi></msub><mi>I</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub></mrow></math>]]></maths>定义i<n时水平算子D<sub>h;ij</sub>I=I<sub>i+1,j</sub>-I<sub>ij</sub>,j<n时垂直算子D<sub>v;ij</sub>I=I<sub>i,j+1</sub>-I<sub>ij</sub>,I<sub>ij</sub>表示图像块矩阵I<sub>n×n</sub>第i行、第j列的值;相应地,以N×1维的向量x表示:i<n时水平算子D<sub>h;ij</sub>x=x<sub>(j-1)n+i+1</sub>-x<sub>(j-1)n+i</sub>,j<n时垂直算子D<sub>v;ij</sub>x=x<sub>nj+i</sub>-x<sub>(j-1)n+i</sub>;在对图像信号进行压缩感知的重构前,先减少位于图像块边缘区域的像素,使重构图像块的尺寸变大。
地址 100124 北京市朝阳区平乐园100号