发明名称 一种基于图像边缘保持自适应分解的超分辨率重建方法
摘要 本发明公开了一种图像超分辨率重建方法。该方法包括两部分:1、图像的自适应分解。这里所谓的自适应是指图像数据自适应的驱动分解过程。2、在图像自适应分解的基础上进行超分辨率放大。把图像自适应分解为低频信息和高频的信息的基础上,在高频部分应用高斯混合模型作为先验模型进行超分辨率放大,低频部分采用样条插值,最后进行叠加。本发明可以实现图像的高品质放大。<pb pnum="1" />
申请公布号 CN106508045B 申请公布日期 2014.08.27
申请号 CN201110014854.1 申请日期 2011.12.05
申请人 中国科学院自动化研究所 发明人 谢启伟;韩华;陈曦;沈丽君
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 中国科学院西安专利中心 61001 代理人 任越
主权项 一种图像超分辨率重建方法,其特征在于,该方法包括以下步骤:步骤1,输入初始图像;步骤2,对初始图像进行图像边缘保持自适应分解,得到初始图像的卡通部分图像与振动部分图像;步骤3,对分解得到的卡通部分图像与振动部分图像分别进行放大;步骤4,对放大后的卡通部分图像和振动部分图像进行叠加得到最终的超分辨率重建图像;所述步骤2进一步包括以下步骤:步骤2.1,根据初始图像I确定邻域尺度;步骤2.2,根据初始图像I的灰度信息计算得到局部邻域的极值;步骤2.3,根据极值点信息构造上下包络图像;步骤2.4,取上下包络图像的平均值得到平均包络图像,即卡通部分图像;步骤2.5,求取初始图像与卡通部分图像的差,得到振动部分图像;所述步骤2.3中,使用下式根据极值点信息构造上下包络图像:<maths id="cmaths0001" num="0001"><math><![CDATA[<mrow><mi>J</mi><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow><mo>=</mo><mi>m</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><munder><mo>&Sigma;</mo><mrow><mi>r</mi><mo>&NotElement;</mo><mi>Q</mi></mrow></munder><msup><mrow><mo>(</mo><mrow><mi>E</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>-</mo><munder><mo>&Sigma;</mo><mi>s</mi></munder><msub><mi>w</mi><mrow><mi>r</mi><mi>s</mi></mrow></msub><mi>E</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FWW0000000072780000011.GIF" wi="718" he="101" /></maths>s.t.E(r)=I(r),r∈Q,其中,E表示待求的包络图像,J(E)为关于包络图像E的目标函数,r和s为初始图像中的任一点的位置,E(r)、E(s)分别为对应点r和s位置上的包络图像像素点,I(r)为对应点r位置上的初始图像像素点,Q表示局部极值位置,w<sub>rs</sub>表示所采用的局部邻域权重系数,w<sub>rs</sub>的加权和为1,如果Q表示局部极大值位置,则求解出的E(r)为上包络,如果Q表示局部极小值位置,则求解出的E(r)为下包络;所述步骤2.4中得到的卡通部分图像表示为:C=(E<sub>u</sub>+E<sub>d</sub>)/2,其中,C表示卡通部分图像,E<sub>u</sub>表示上包络图像,E<sub>d</sub>表示下包络图像;所述步骤2.5中得到的振动部分图像表示为:<math><![CDATA[<mrow><mover><msubsup><mi>I</mi><mi>T</mi><mi>m</mi></msubsup><mo>^</mo></mover><mo>=</mo><msub><mi>argmax</mi><mi>I</mi></msub><mfrac><mn>1</mn><mrow><msqrt><mrow><mn>2</mn><mi>&pi;</mi></mrow></msqrt><mi>&sigma;</mi></mrow></mfrac><mi>exp</mi><mrow><mo>(</mo><mfrac><mrow><mo>|</mo><mo>|</mo><mi>H</mi><mrow><mo>(</mo><msub><mi>I</mi><mi>s</mi></msub><mo>+</mo><msubsup><mi>I</mi><mi>T</mi><mi>m</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mi>L</mi><mi>o</mi><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup></mrow><mrow><mn>2</mn><msup><mi>&sigma;</mi><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><munderover><mo>&Pi;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>w</mi><mi>j</mi></msub><mi>g</mi><mrow><mo>(</mo><msubsup><mi>I</mi><mrow><mi>T</mi><mo>,</mo><mi>i</mi></mrow><mi>m</mi></msubsup><mo>|</mo><msub><mi>&mu;</mi><mi>j</mi></msub><mo>,</mo><msub><mo>&Sigma;</mo><mi>j</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FWW0000000072780000024.GIF" wi="1369" he="142" /></maths>其中,I为原始图像,<img file="FWW0000000072780000025.GIF" wi="52" he="53" />为每次迭代时想要得到的超分辨率图像的振动部分,<img file="FWW0000000072780000026.GIF" wi="349" he="63" />w<sub>j</sub>为权重系数,<img file="FWW0000000072780000027.GIF" wi="267" he="60" />为高斯混合模型部分,i为原始图像中任一点的位置,k为混合高斯模型中不同的高斯模型的个数,N为图像的像素数,H为模糊核矩阵,L<sub>O</sub>=I,<img file="FWW0000000072780000028.GIF" wi="60" he="58" />为<img file="FWW0000000072780000029.GIF" wi="51" he="54" />的对应位置的灰度值,m为迭代次数;所述步骤3.4中的条件为:<maths id="cmaths0003" num="0003"><math><![CDATA[<mrow><mo>|</mo><mo>|</mo><msubsup><mi>I</mi><mi>T</mi><mi>m</mi></msubsup><mo>-</mo><msubsup><mi>I</mi><mi>T</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub><mo>&lt;</mo><mi>&epsiv;</mi><mo>,</mo></mrow>]]></math><img file="FWW00000000727800000210.GIF" wi="315" he="85" /></maths>其中,m为迭代次数,||·||为2‑范数运算,ε为一个较小的值;所述权重系数w<sub>j</sub>的加权和为1;最终得到的超分辨率重建图像表示为:<maths id="cmaths0004" num="0004"><math><![CDATA[<mrow><msub><mi>I</mi><mi>L</mi></msub><mo>=</mo><msub><mi>I</mi><mi>s</mi></msub><mo>+</mo><msubsup><mi>I</mi><mi>T</mi><mi>m</mi></msubsup><mo>,</mo></mrow>]]></math><img file="FWW00000000727800000211.GIF" wi="245" he="56" /></maths>其中,I<sub>L</sub>为最终得到的超分辨率重建图像,I<sub>s</sub>为放大的卡通部分图像,<img file="FWW00000000727800000212.GIF" wi="50" he="55" />为放大的振动部分图像。pb pnum="2" />
地址 100190 北京市海淀区中关村东路95号