发明名称 一种基于局部回归模型的图像超分辨率重建方法
摘要 本发明公开一种基于局部回归模型的图像超分辨率重建方法,首先,对输入低分辨率图像,高斯低通滤波得到其低频带图像,双立方插值得到输出高分辨率图像的近似低频带图像。其次,重建时对高分辨率图像的低频带图像中每个图像块应用一阶回归模型,回归模型中高/低图像间的映射函数可对输入图像通过机器学的方法得到,即利用输入低分辨率图像及其低频带图像对应位置采样得到对应位样本图像块对进行字典训练。最后,重建图像块的非局部自相似块分别应用一阶回归模型,加权综合得到重建的高分辨率图像块。本发明所提出方法无需利用外部图像模型,而是利用输入图像自身学得到先验模型,依此模型重建出的高分辨图像取得了较好的主客观重建效果。
申请公布号 CN105488759A 申请公布日期 2016.04.13
申请号 CN201510902909.0 申请日期 2015.12.09
申请人 南京邮电大学 发明人 李欣;崔子冠;干宗良;唐贵进;朱秀昌
分类号 G06T3/40(2006.01)I 主分类号 G06T3/40(2006.01)I
代理机构 南京知识律师事务所 32207 代理人 汪旭东
主权项 一种基于局部回归模型的图像超分辨率重建方法,其特征在于,包括如下步骤,步骤1:读入待重建低分辨率图像X<sub>0</sub>,放大因子s;步骤2:对X<sub>0</sub>高斯低通滤波得到其低频带图像Y<sub>0</sub>,对X<sub>0</sub>双立方插值近似输出高分辨率图像的低频带图像Y;步骤3:对Y划分成互相重叠的大小为a×a的图像块y;步骤4:根据步骤3得到的图像块y的位置(i,j)找到其在Y<sub>0</sub>中的对应位(i<sub>r</sub>,j<sub>r</sub>),其中<img file="FDA0000872036480000011.GIF" wi="790" he="87" />得到y在Y<sub>0</sub>和X<sub>0</sub>中的对应位相似样本块对{y<sub>0</sub>,x<sub>0</sub>};步骤5:由一阶局部回归模型计算步骤3得到的图像块y的高分辨率图像块x,公式如下:<math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>x</mi><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mrow><msub><mi>y</mi><mn>0</mn></msub><mo>+</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>+</mo><mo>&dtri;</mo><msup><mi>f</mi><mi>T</mi></msup><mrow><mo>(</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub></mrow><mo>)</mo></mrow><mo>+</mo><mi>O</mi><mrow><mo>{</mo><mrow><mo>|</mo><mo>|</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup></mrow><mo>}</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>&ap;</mo><msub><mi>x</mi><mn>0</mn></msub><mo>+</mo><mo>&dtri;</mo><msup><mi>f</mi><mi>T</mi></msup><mrow><mo>(</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000872036480000012.GIF" wi="1486" he="359" /></maths>步骤6:计算步骤3得到的图像块y的非局部自相似块的对应位相似样本块对组<img file="FDA0000872036480000013.GIF" wi="254" he="111" />分别执行步骤5,将计算结果加权求和,公式如下:<math><![CDATA[<mrow><mi>x</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><mrow><mo>(</mo><msubsup><mi>x</mi><mn>0</mn><mi>p</mi></msubsup><mo>+</mo><mo>&dtri;</mo><msup><mi>f</mi><mi>T</mi></msup><mo>(</mo><msubsup><mi>y</mi><mn>0</mn><mi>p</mi></msubsup><mo>)</mo><mo>(</mo><mrow><mi>y</mi><mo>-</mo><msubsup><mi>y</mi><mn>0</mn><mi>p</mi></msubsup></mrow><mo>)</mo><mo>)</mo></mrow><msub><mi>w</mi><mi>p</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000872036480000014.GIF" wi="1397" he="181" /></maths>权值w<sub>p</sub>可通过共轭梯度法求解下列公式得到:<img file="FDA0000872036480000015.GIF" wi="1349" he="127" />其中,<img file="FDA0000872036480000016.GIF" wi="517" he="110" />w=[w<sub>1</sub>,w<sub>2</sub>,...w<sub>p</sub>,...w<sub>J</sub>];步骤7:对所有像素重叠的高分辨率图像块x叠加平均,得到高分辨率图像X。
地址 210003 江苏省南京市鼓楼区新模范马路66号