发明名称 基于多尺度结构自相似与压缩感知的单图像超分辨率方法
摘要 基于多尺度结构自相似与压缩感知的单图像超分辨率方法,先设置高分辨率重构图像的初始估计值,设置迭代中止的误差,迭代最大的次数,根据图像的降质过程确定降采样矩阵和模糊矩阵,构建图像金字塔,并将其作为K-SVD方法的训练样本建立字典;按照Nonlocal方法在当前高分辨率重构图像中搜索具有相同尺度的相似图像块并确定权值矩阵;更新高分辨率重构图像的估计值,更新稀疏表示系数,更新高分辨率重构图像的估计值;然后进行下一次迭代,直到连续两步的高分辨率重构图像满足相应要求或达到最大的迭代次数,本发明通过压缩感知框架将蕴含在图像多尺度自相似结构中的附加信息加入到高分辨率重构图像中,具有更高的运算效率。
申请公布号 CN103020909B 申请公布日期 2015.02.18
申请号 CN201210519587.8 申请日期 2012.12.06
申请人 清华大学 发明人 潘宗序;禹晶;孙卫东
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 西安智大知识产权代理事务所 61215 代理人 贾玉健
主权项 基于多尺度结构自相似与压缩感知的单图像超分辨率方法,包括如下步骤:步骤1:设置高分辨率重构图像的初始估计值<img file="FDA0000610646110000011.GIF" wi="124" he="84" />k=0,设置迭代中止的误差∈,迭代最大的次数K<sub>max</sub>;步骤2:根据图像的降质过程确定降采样矩阵D和模糊矩阵H;步骤3:构建图像金字塔,并将其作为K‑SVD方法的训练样本建立字典Ψ;步骤4:按照Nonlocal方法在当前高分辨率重构图像中搜索具有相同尺度的相似图像块并确定权值矩阵B;步骤5:更新高分辨率重构图像的估计值<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mn>2</mn><mo>)</mo></mrow></msup><mo>=</mo><msup><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>+</mo><msup><mi>K</mi><mi>T</mi></msup><mrow><mo>(</mo><mover><mi>Y</mi><mo>~</mo></mover><mo>-</mo><mi>K</mi><msup><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>)</mo></mrow><mo>=</mo><msup><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>+</mo><mrow><mo>(</mo><msup><mrow><mo>(</mo><mi>DH</mi><mo>)</mo></mrow><mi>T</mi></msup><mi>Y</mi><mo>-</mo><mi>U</mi><msup><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>-</mo><mi>V</mi><msup><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000610646110000012.GIF" wi="1724" he="119" /></maths>其中,U=(DH)<sup>T</sup>DH,V=η<sup>2</sup>(I‑B)<sup>T</sup>(I‑B);步骤6:更新稀疏表示系数<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>&alpha;</mi><mi>i</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mn>2</mn><mo>)</mo></mrow></msubsup><mo>=</mo><msup><mi>&Psi;</mi><mi>T</mi></msup><msub><mi>R</mi><mi>i</mi></msub><msup><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mn>2</mn><mo>)</mo></mrow></msup><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>p</mi><mo>,</mo></mrow>]]></math><img file="FDA0000610646110000013.GIF" wi="905" he="95" /></maths><img file="FDA0000610646110000014.GIF" wi="608" he="112" />其中R<sub>i</sub>为抽取矩阵,p为图像块的个数,soft(x,τ)=sign(x)max(|x|‑τ,0)为含有阈值τ的软阈值函数,sign(x)表示符号函数;步骤7:更新高分辨率重构图像的估计值<img file="FDA0000610646110000015.GIF" wi="1300" he="116" />步骤8:k=k+1,进行下一次迭代,重复步骤4至步骤7,直到连续两步的高分辨率重构图像满足<img file="FDA0000610646110000016.GIF" wi="531" he="124" />或迭代次数k达到K<sub>max</sub>;其中:<img file="FDA0000610646110000017.GIF" wi="99" he="70" />表示第k次迭代后重构图像的估计值,是第k+1次迭代开始时重构图像的估计值;<img file="FDA0000610646110000021.GIF" wi="154" he="76" />表示在第k+1次迭代过程中对重构图像估计值的更新,即利用观测方程和非局部约束,通过误差反向投影进行更新从而获取的重构图像的估计值;<img file="FDA0000610646110000022.GIF" wi="131" he="77" />表示第k+1次迭代后重构图像的估计值,即通过收缩稀疏表示系数获得的重构图像的估计值;K为观测模型中的降采样矩阵、模糊矩阵与Nonlocal方法中的权值矩阵所组成的矩阵,对应观测模型与非局部约束对超分辨率重构的限制;<img file="FDA0000610646110000023.GIF" wi="59" he="75" />表示经过0延拓的低分辨率图像;Y表示低分辨率图像;η表示控制非局部约束项在代价函数中权重的参数;α<sub>i</sub><sup>(k+1/2)</sup>表示第k+1次迭代过程中图像块的稀疏表示系数,即利用<img file="FDA0000610646110000024.GIF" wi="158" he="70" />进行更新从而获得的图像块的稀疏表示系数;α<sub>i</sub><sup>(k+1)</sup>表示第k+1次迭代后图像块的稀疏表示系数,即将α<sub>i</sub><sup>(k+1/2)</sup>进行系数收缩从而获取的图像块的稀疏表示系数;α<sup>(k+1)</sup>表示重构图像的稀疏表示系数,由图像块的稀疏表示系数拼接而成。
地址 100084 北京市海淀区100084信箱82分箱清华大学专利办公室