发明名称 一种基于硬件加速的单帧红外图像超分辨方法
摘要 本发明属于红外成像图像处理领域,具体涉及一种基于硬件加速的单帧红外图像超分辨方法。它包括:步骤一:双线性插值;步骤二:滤波;滤波得到四幅图像,步骤三:计算;计算得到最终的高分辨率图像。本发明的有益效果是:用本申请的方法可以在低分辨率图像的基础上,经过图像处理得到足够满足要求的高分辨率图像。该方法不仅适用于普通可见光图像,而且适用于红外图像。同时该方法可以很方便的用硬件实现,满足系统实时性的要求。
申请公布号 CN103310431A 申请公布日期 2013.09.18
申请号 CN201310240454.1 申请日期 2013.06.17
申请人 中国航天科工集团第三研究院第八三五八研究所 发明人 彭志勇;任亚芬;胡耀金;张瑞
分类号 G06T5/50(2006.01)I 主分类号 G06T5/50(2006.01)I
代理机构 核工业专利中心 11007 代理人 高尚梅;刘昕宇
主权项 1.一种基于硬件加速的单帧红外图像超分辨方法,其特征在于:包括下述步骤:步骤一:双线性插值针对输入的低分辨率图像,采用下述公式进行双线性插值:d<sub>p</sub>(i,j)=(1-w<sub>q</sub>)d(p,q)+w<sub>q</sub>d(p,q+1)d<sub>p+1</sub>(i,j)=(1-w<sub>q</sub>)d(p+1,q)+w<sub>q</sub>d(p+1,q+1)    (1)d(i,j)=(1-w<sub>p</sub>)d(i,j)+w<sub>p</sub>d<sub>p+1</sub>(i,j)其中,(i,j)表示新图像的像素点,d(i,j)表示点(i,j)的灰度值,(p,q),(p,q+1),(p+1,q),(p+1,q+1)表示原图像中的点,w<sub>p</sub>和w<sub>q</sub>为(p,q)点的权重,并且p,qw<sub>p</sub>和w<sub>q</sub>按照下述公式计算得到<img file="FDA00003355216700011.GIF" wi="342" he="151" /><img file="FDA00003355216700012.GIF" wi="334" he="148" /><img file="FDA00003355216700013.GIF" wi="1370" he="198" /><img file="FDA00003355216700014.GIF" wi="918" he="156" />其中<img file="FDA00003355216700015.GIF" wi="86" he="88" />表示取整,mod表示模除本步骤得到双插值后的灰度图像,图像中象元的灰度值用d(i,j)或者用d<sub>ij</sub>表示,步骤二:滤波分别用下述四个滤波模板对步骤一中得到的图像滤波,滤波得到四幅图像,滤波模板分别为hf<sub>1</sub>=[-1,0,1];vf<sub>1</sub>=[-1,0,1]'                    (4)hf<sub>2</sub>=[1,0,-2,0,1];vf<sub>2</sub>=[1,0,-2,0,1]'四幅图像分别为d<sub>ij</sub><sup>1</sup>=d<sub>ij</sub>*hf<sub>1</sub>d<sub>ij</sub><sup>2</sup>=d<sub>ij</sub>*vf<sub>1</sub>                            (5)d<sub>ij</sub><sup>3</sup>=d<sub>ij</sub>*hf<sub>2</sub>d<sub>ij</sub><sup>4</sup>=d<sub>ij</sub>*vf<sub>2</sub>其中d<sub>ij</sub><sup>1</sup>表示用第一个滤波模板进行滤波后得到图像象元的灰度值,类似的d<sub>ij</sub><sup>2</sup>、d<sub>ij</sub><sup>3</sup>、d<sub>ij</sub><sup>4</sup>就是用第二、三、四个滤波模板进行滤波后得到图像象元的灰度值,d<sub>ij</sub>为步骤一得到图像象元的灰度值,步骤三:计算在进行计算之前首先建立高分辨率图像字典D<sub>h</sub>和低分辨率图像字典D<sub>l</sub>,字典的建立过程为:首选选取高分辨率(640×512)红外图像,然后通过过采样的手段降低图像分辨率,得到低分辨率(320×256)红外图像,则高分辨率(640×512)红外图像被加入高分辨率图像字典D<sub>h</sub>,低分辨率(320×256)红外图像被加入低分辨率图像字典D<sub>l</sub>,同时建立高分辨率图像字典中特定图像与低分辨率图像字典中特定图像的一一对应关系,字典中图像的个数一般以25个为最优,当然也可以选取其他个数的字典,一般字典中图像个数选取的数量范围为15~30个,将所有高分辨率图像字典D<sub>h</sub>和低分辨率图像字典D<sub>l</sub>中的全部图像都抽取出来,然后将高分辨率图像字典D<sub>h</sub>的所有图像排列成三维矩阵,该三维矩阵中的横坐标和纵坐标为高分辨率图像字典D<sub>h</sub>中图像的位置坐标,三维矩阵中的高度坐标为图像序号,类似的也建立由低分辨率图像字典D<sub>l</sub>的所有图像排列成三维矩阵,此时低分辨率图像字典D<sub>l</sub>所有图像排列成三维矩阵的排列顺序与高分辨率图像字典D<sub>h</sub>所有图像排列成三维矩阵的排列顺序保持一致,最后,将步骤二中得到的四幅图像,也按照上述规则,构建三维矩阵,上述构建的三个三维矩阵分别记录为:d<sup>*</sup>表示用步骤二中图像构建的三维矩阵,D<sup>*</sup><sub>l</sub>表示用低分辨率图像字典D<sub>l</sub>所有图像排列成的三维矩阵,D<sup>*</sup><sub>h</sub>表示用高分辨率图像字典D<sub>h</sub>所有图像排列成的三维矩阵,用下述公式计算系数矩阵α<sup>*</sup><maths num="0001"><![CDATA[<math><mrow><msup><mi>&alpha;</mi><mo>*</mo></msup><mo>=</mo><mi>arg</mi><munder><mi>min</mi><mi>&alpha;</mi></munder><mi>f</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>)</mo></mrow><mo>=</mo><mi>arg</mi><mi>min</mi><mfrac><mn>1</mn><mn>2</mn></mfrac><msubsup><mrow><mo>|</mo><mo>|</mo><msup><mi>d</mi><mo>*</mo></msup><mo>-</mo><msub><msup><mi>D</mi><mo>*</mo></msup><mi>l</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mi>&lambda;</mi><msub><mrow><mo>|</mo><mo>|</mo><mi>&alpha;</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>上述公式中||||<sub>2</sub>表示运算符号中的所有元素平方之和再开方,D<sup>*</sup><sub>l</sub>和d<sup>*</sup>为前面建立的三维矩阵,λ为预先设定的常数,argmin表示使公式值最小的变量,最后用下述公式计算高分辨图像D=D<sup>*</sup><sub>h</sub>α<sup>*</sup>    (7)上述公式中α<sup>*</sup>为公式(6)计算得到的结果,D<sup>*</sup><sub>h</sub>表示用高分辨率图像字典D<sub>h</sub>所有图像排列成的三维矩阵,D为最终得到的高分辨率图像。
地址 300192 天津市南开区航天道58号