发明名称 一种基于双层模型的单幅图像超分辨率重建方法
摘要 一种基于双层模型的单幅图像超分辨率重建方法,其包括(1)利用L<sub>0</sub>梯度最小化方法和HoG算子生成K个训练簇,然后训练它的对应字典对;(2)依据HoG算子,自适应的选取用于测试的低分辨率图像块对应的几何字典对,并求解出低分辨率图像对应的高分辨率纹理细节图像;(3)利用L<sub>0</sub>梯度最小化方法求解出用于测试的低分辨率图像对应的高分辨率边缘结构图像;(4)把求解出的高分辨率纹理细节图像加到高分辨率边缘结构图像上得到初始高分辨率图像;(5)对初始高分辨率图像进行全局和局部约束得到最终高分辨率图像。本发明可以使得重建后的图像轮廓更清晰,细节信息更丰富,提高重建后的图像质量。
申请公布号 CN103279933B 申请公布日期 2016.10.26
申请号 CN201310225721.8 申请日期 2013.06.07
申请人 重庆大学 发明人 龚卫国;李进明;李伟红;王立;潘飞宇;李正浩;杨利平
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 重庆华科专利事务所 50123 代理人 康海燕
主权项 一种基于双层模型的单幅图像超分辨率重建方法,所述方法包括以下步骤:(1)从高分辨图像库中随机选取n幅作为训练图像X<sub>1</sub>,X<sub>2</sub>,…X<sub>n</sub>,然后运用退化模型式(1)生成对应的低分辨率训练图像Y<sub>1</sub>,Y<sub>2</sub>,…Y<sub>n</sub>;Y=UBX  (1)其中,向量X表示高分辨率图像,向量Y表示对应的低分辨率图像,矩阵U表示下采样算子,矩阵B表示模糊算子;(2)对所有低分辨率训练图像Y<sub>i</sub>进行L<sub>0</sub>梯度最小化操作,生成高分辨率边缘结构图像<img file="FDA0001061487690000016.GIF" wi="99" he="63" />然后执行X<sub>i</sub>与<img file="FDA0001061487690000017.GIF" wi="74" he="63" />的相减操作,生成高分辨率纹理细节图像<img file="FDA0001061487690000018.GIF" wi="95" he="68" />(3)首先对所有的低分辨率训练图像Y<sub>i</sub>执行尺寸为3×3的分块操作,随机选取P个低分辨率图像块,并以向量y<sub>i</sub>来表示,获取对应的尺寸为3s×3s高分辨率纹理细节图像块,并以向量<img file="FDA0001061487690000019.GIF" wi="51" he="60" />来表示;最后聚集它们生成P个训练对<img file="FDA00010614876900000110.GIF" wi="899" he="70" />其中,P经验的选取为9万块,s表示重建倍数;(4)计算训练对<img file="FDA00010614876900000111.GIF" wi="880" he="70" />中y<sub>i</sub>的方向梯度直方图HoG,然后依据y<sub>i</sub>的HoG中的方向值,用K‑means算法对训练对<img file="FDA00010614876900000112.GIF" wi="704" he="71" /><img file="FDA00010614876900000113.GIF" wi="139" he="74" />进行分簇,得到K个簇C<sub>1</sub>,C<sub>2</sub>,...C<sub>K</sub>,其中C<sub>K</sub>表示第K个簇中具有相同方向值的训练对集合<img file="FDA00010614876900000114.GIF" wi="912" he="79" />k表示在第K簇中训练对的数目;(5)用式(2)训练K个簇C<sub>1</sub>,C<sub>2</sub>,...C<sub>K</sub>的字典对(D<sup>i</sup><sub>l</sub>,D<sup>i</sup><sub>T</sub>),得到K个字典对{(D<sup>1</sup><sub>l</sub>,D<sup>1</sup><sub>T</sub>),(D<sup>2</sup><sub>l</sub>,D<sup>2</sup><sub>T</sub>),...(D<sup>K</sup><sub>l</sub>,D<sup>K</sup><sub>T</sub>)};<maths num="0001"><math><![CDATA[<mrow><mi>arg</mi><munder><mi>min</mi><mrow><msubsup><mi>D</mi><mi>l</mi><mi>i</mi></msubsup><mo>,</mo><msubsup><mi>D</mi><mi>T</mi><mi>i</mi></msubsup><mo>,</mo><mi>S</mi></mrow></munder><mo>|</mo><mo>|</mo><msup><mi>X</mi><mi>i</mi></msup><mo>-</mo><msup><mi>D</mi><mi>i</mi></msup><mi>S</mi><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mi>&lambda;</mi><mo>|</mo><mo>|</mo><mi>S</mi><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001061487690000011.GIF" wi="1572" he="122" /></maths>其中,<img file="FDA0001061487690000012.GIF" wi="702" he="278" /><img file="FDA0001061487690000013.GIF" wi="46" he="63" />表示第i个簇中所有<img file="FDA0001061487690000014.GIF" wi="51" he="55" />集合,y<sup>i</sup>表示第i个簇中所有y<sub>i</sub>集合,λ为正则化参数,N和M为<img file="FDA0001061487690000015.GIF" wi="58" he="55" />和y<sub>k</sub>中像素的数目,S为稀疏系数;(6)对于一幅用于测试的尺寸为H×L(高×宽)的低分辨率图像Z,将Z分成尺寸为3×3的块,并以向量Z<sub>i</sub>来表示,块与块之间重叠2个像素;(7)对所有Z<sub>i</sub>进行以下遍历操作:首先,任取一个Z<sub>i</sub>,计算Z<sub>i</sub>与C<sub>1</sub>,C<sub>2</sub>,...C<sub>K</sub>的簇中心距离:d<sub>1</sub>,d<sub>2</sub>,...,d<sub>K</sub>,提取最小的d<sub>i</sub>对应簇的字典对(D<sup>i</sup><sub>l</sub>,D<sup>i</sup><sub>T</sub>),然后,利用式(2)求解出对应的高分辨率纹理细节图像块X<sub>Ti</sub>;(8)拼接由步骤(7)求出的所有高分辨率纹理细节图像块X<sub>Ti</sub>,组成初始的高分辨率纹理细节图像X<sub>T</sub>;(9)用L<sub>0</sub>梯度最小化方法求解测试图像Z的高分辨率边缘结构图像X<sub>E</sub>;(10)获取Z的初始高分辨率图像X<sup>0</sup>,X<sup>0</sup>=X<sub>E</sub>+X<sub>T</sub>,然后对X<sup>0</sup>执行全局和局部的约束,得到最终的高分辨率重建图像X<sup>*</sup>。
地址 400030 重庆市沙坪坝区沙正街174号