发明名称 基于统计的图像压缩感知低复杂度还原方法
摘要 本发明针对基于压缩感知的图像还原问题,给出了一种结合自然图像小波表示下统计信息的还原方法和还原系统,所述还原方法包括图像高速还原方法和图像高精度还原方法,可根据实际要求实现高速还原和高精度还原。其中,图像高速还原方法包括列维度一步快速还原,行维度一步快速还原;图像高精度还原方法包括列维度一步快速还原,行维度一步快速还原,行维度二次迭代还原。本发明的有益效果是:本发明充分利用图像的统计信息,在还原质量与还原速度上,较传统方法有很大提高,并可根据实际需求选择是否进行二次迭代,从而提升了整个图像压缩感知系统的性能和自由度。
申请公布号 CN104952046A 申请公布日期 2015.09.30
申请号 CN201510151657.2 申请日期 2015.04.01
申请人 哈尔滨工业大学深圳研究生院 发明人 吴绍华;杨竞然;陈大薇;张钦宇;刘云路
分类号 G06T5/00(2006.01)I;G06T9/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 深圳市科吉华烽知识产权事务所(普通合伙) 44248 代理人 刘显扬;胡玉
主权项 一种基于图像统计信息的压缩感知图像还原方法,其特征在于:所述还原方法利用自然图像统计信息,采用两个维度分别对图像进行高速还原,具体为针对一幅大小为n×n的图像X,经过投影后,获得大小为m×m的投影值Y,即Y=ΦXΦ<sup>T</sup>=ΦΨSΨ<sup>T</sup>Φ<sup>T</sup>=ASA<sup>T</sup>,其中A=ΦΨ且Φ为随机矩阵,Ψ为小波变换矩阵,S是原始图像X在小波域下的稀疏表示,即X=ΨSΨ<sup>T</sup>,对其进行还原,所述还原方法包括以下步骤:(1)列维度一步快速还原:利用表征图像小波域下衰减信息的权重矩阵K,对图像的列维度进行一步快速还原;令S<sub>1</sub>=SA<sup>T</sup>,Y=Y<sub>1</sub>则有Y<sub>1</sub>=A×S<sub>1</sub>,即可获得列维度还原结果<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>S</mi><mo>^</mo></mover><mn>1</mn></msub><mo>=</mo><msup><mi>K</mi><mn>2</mn></msup><mo>&times;</mo><msup><mi>A</mi><mi>T</mi></msup><mo>&times;</mo><msup><mrow><mo>(</mo><mi>A</mi><mo>&times;</mo><msup><mi>K</mi><mn>2</mn></msup><mo>&times;</mo><msup><mi>A</mi><mi>T</mi></msup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>&times;</mo><msub><mi>Y</mi><mn>1</mn></msub><mo>;</mo></mrow>]]></math><img file="FDA0000692484390000011.GIF" wi="706" he="93" /></maths>其中,所述权重矩阵K表征自然图像统计信息,其非主对角线元素均为0,主对角线元素一定程度上符合自然图像小波域表示下由上至下,由左至右衰减的特点;(2)转置:将求解结果<img file="FDA0000692484390000012.GIF" wi="60" he="86" />转置,<img file="FDA0000692484390000013.GIF" wi="190" he="86" />(3)行维度一步快速还原:利用表征图像小波域下衰减信息的权重矩阵K,对图像的行维度进行一步快速还原;因<img file="FDA0000692484390000014.GIF" wi="366" he="87" />令S<sub>2</sub>=S<sup>T</sup>,故可转化为对式Y<sub>2</sub>=A×S<sub>2</sub>求解,其求解结果为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>S</mi><mo>^</mo></mover><mn>2</mn></msub><mo>=</mo><msup><mi>K</mi><mn>2</mn></msup><mo>&times;</mo><msup><mi>A</mi><mi>T</mi></msup><mo>&times;</mo><msup><mrow><mo>(</mo><mi>A</mi><mo>&times;</mo><msup><mi>K</mi><mn>2</mn></msup><mo>&times;</mo><msup><mi>A</mi><mi>T</mi></msup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>&times;</mo><msub><mi>Y</mi><mn>2</mn></msub><mo>;</mo></mrow>]]></math><img file="FDA0000692484390000015.GIF" wi="718" he="86" /></maths>(4)转置:对上步求解结果<img file="FDA0000692484390000016.GIF" wi="68" he="84" />转置,<img file="FDA0000692484390000017.GIF" wi="145" he="86" />即获得图像的小波域稀疏表示(5)逆稀疏变换步骤:将获得图像的小波域表示转换为空域图像信息,<img file="FDA0000692484390000018.GIF" wi="238" he="76" />即还原出图像的原始信号。
地址 518000 广东省深圳市南山区西丽镇深圳大学城哈工大校区