发明名称 基于对偶约束的联合学超分辨方法
摘要 本发明公开了一种基于对偶约束的联合学超分辨方法,主要解决现有邻域嵌入超分辨方法导致图像模糊的问题,其步骤是:(1)使用高分辨图像生成训练集,建立低分辨和高分辨图像块的成组块对;(2)提取低分辨图像块的特征向量,在训练集中查找最匹配的成组块对;(3)对成组块对中低分辨和高分辨块的特征向量,进行联合学,构造投影矩阵,生成联合特征子空间;(4)在生成的联合特征子空间中使用邻域嵌入估计高分辨图像块;(5)将所有估计的高分辨图像块合成高分辨图像;(6)对合成的高分辨图像,利用全局重构约束和反向投影算法改善图像质量。实验结果表明,本发明具有更强的超分辨能力,可用于低分辨自然图像的放大。
申请公布号 CN101976435B 申请公布日期 2012.10.24
申请号 CN201010298564.X 申请日期 2010.10.07
申请人 西安电子科技大学 发明人 高新波;张凯兵;李洁;邓成;田春娜;路文;王茜;沐广武
分类号 G06T5/50(2006.01)I 主分类号 G06T5/50(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种基于对偶约束的联合学习超分辨方法,包括如下步骤:(1)采集高分辨灰度自然图像,模拟图像降质过程,生成低分辨图像,对低分辨图像使用双立方插值进行2倍放大,得到插值图像,对插值图像和原始高分辨图像,分别随机选择大小为6×6的低分辨图像块和对应大小为9×9的高分辨图像块,并提取水平和垂直方向上的一阶和二阶梯度特征,构造低分辨图像块训练集<img file="FDA00001778911800011.GIF" wi="229" he="84" />和高分辨图像块训练集<img file="FDA00001778911800012.GIF" wi="243" he="84" />以及高分辨图像块的高频细节训练集<img file="FDA00001778911800013.GIF" wi="238" he="84" />其中<img file="FDA00001778911800014.GIF" wi="40" he="58" />表示第i个低分辨图像块特征向量,<img file="FDA00001778911800015.GIF" wi="39" he="58" />表示第i个高分辨图像块特征向量,<img file="FDA00001778911800016.GIF" wi="46" he="58" />表示第i个高分辨图像块的高频细节特征向量,N表示训练集中样本总个数;(2)将低分辨图像块训练集<img file="FDA00001778911800017.GIF" wi="229" he="84" />和高分辨图像块训练集<img file="FDA00001778911800018.GIF" wi="219" he="84" />的每个特征对<img file="FDA00001778911800019.GIF" wi="140" he="76" />进行组合,经归一化后,得到对偶特征c<sub>i</sub>,将N个对偶特征构成对偶特征集<img file="FDA000017789118000110.GIF" wi="219" he="72" />对每个对偶特征c<sub>i</sub>,选择K个最近邻,构成成组块对集合<img file="FDA000017789118000111.GIF" wi="156" he="84" />其中G<sup>i</sup>表示对偶集C中第i个对偶特征c<sub>i</sub>的K个最近邻组成的成组块对,K=128;(3)输入待处理的低分辨彩色图像F,并将图像F从红、绿、蓝的RGB颜色空间转换到YCbCr颜色空间,其中Y表示亮度分量,两个色差分量Cb和Cr表示颜色信息;(4)对YCbCr颜色空间中的亮度分量Y使用双立方插值进行2倍放大,得到插值图像I<sub>int</sub>;(5)将插值图像I<sub>int</sub>自上而下、从左到右划分成大小为6×6的图像块,使相邻块之间重叠2个像素,提取该图像块在水平和垂直方向的一阶和二阶梯度特征构成低分辨图像块测试集<img file="FDA000017789118000112.GIF" wi="260" he="91" />其中<img file="FDA000017789118000113.GIF" wi="46" he="57" />表示低分辨图像块测试集X<sub>t</sub>中的第j个特征向量,M表示该测试集X<sub>t</sub>中样本总个数;(6)对低分辨图像块测试集<img file="FDA000017789118000114.GIF" wi="237" he="91" />中每个特征向量<img file="FDA000017789118000115.GIF" wi="69" he="57" />在低分辨图像块训练集<img file="FDA000017789118000116.GIF" wi="230" he="84" />中查找与<img file="FDA000017789118000117.GIF" wi="46" he="58" />最匹配的样本<img file="FDA000017789118000118.GIF" wi="63" he="58" />得到该样本所对应的成组块对G<sup>i</sup>,从该成组块对G<sup>i</sup>中分离出K个低分辨图像块和高分辨图像块特征向量,分别构成低分辨图像块成组矩阵J<sub>l</sub>和高分辨图像块成组矩阵J<sub>h</sub>,其中<img file="FDA00001778911800021.GIF" wi="209" he="59" /><img file="FDA00001778911800022.GIF" wi="221" he="58" />d和m分别是低分辨图像块和对应的高分辨图像块特征向量的维数,K表示每个成组块对中低、高分辨率特征向量组成的对偶个数,d=144,m=324;(7)使用联合学习方法分别构造低分辨图像块成组矩阵J<sub>1</sub>的投影矩阵P<sub>1</sub>和高分辨图像块成组矩阵J<sub>h</sub>的投影矩阵P<sub>h</sub>,将低分辨图像块测试集X<sub>t</sub>中的特征向量<img file="FDA00001778911800023.GIF" wi="47" he="58" />通过投影矩阵P<sub>l</sub>变换成一个<img file="FDA00001778911800024.GIF" wi="79" he="45" />(r+1≤d)维的变换特征<img file="FDA00001778911800025.GIF" wi="122" he="59" />将高分辨图像块成组矩阵J<sub>h</sub>中的K个特征向量通过投影矩阵P<sub>h</sub>变换成<img file="FDA00001778911800026.GIF" wi="97" he="45" />维的变换特征<img file="FDA00001778911800027.GIF" wi="128" he="59" />将低分辨图像块测试集中的特征向量<img file="FDA00001778911800028.GIF" wi="46" he="58" />的变换特征<img file="FDA00001778911800029.GIF" wi="102" he="59" />与高分辨图像块成组矩阵J<sub>h</sub>的K个变换特征<img file="FDA000017789118000210.GIF" wi="103" he="59" />构成一个<img file="FDA000017789118000211.GIF" wi="145" he="46" />维的联合特征子空间S;(8)在联合特征子空间S中选择特征向量<img file="FDA000017789118000212.GIF" wi="46" he="57" />的变换特征<img file="FDA000017789118000213.GIF" wi="103" he="58" />的q个邻域,使用邻域嵌入算法,估计最优权值,将变换特征<img file="FDA000017789118000214.GIF" wi="103" he="58" />的q个邻域在高频细节训练集<img file="FDA000017789118000215.GIF" wi="214" he="85" />中对应的q个高频细节特征向量进行线性组合,并加上输入的低分辨图像块的灰度均值<img file="FDA000017789118000216.GIF" wi="54" he="54" />得到低分辨图像块的特征向量<img file="FDA000017789118000217.GIF" wi="47" he="58" />的高分辨图像块的估计<img file="FDA000017789118000218.GIF" wi="74" he="58" />低分辨图像块测试集<img file="FDA000017789118000219.GIF" wi="237" he="92" />中所有特征向量的高分辨图像块估计构成了预测结果集<img file="FDA000017789118000220.GIF" wi="246" he="91" />q=5;(9)按插值图像I<sub>int</sub>中图像块的划分顺序合并所有预测结果集<img file="FDA000017789118000221.GIF" wi="246" he="91" />对重叠区域的像素取平均值进行融合,得到高分辨图像的初始化估计X<sub>0</sub>;(10)利用全局重构约束F=DBX和反向投影算法,通过如下迭代公式进一步改善超分辨图像的质量:X<sub>n+1</sub>=X<sub>n</sub>+α[B<sup>T</sup>D<sup>T</sup>(F-DBX<sub>n</sub>)+μ(X<sub>n</sub>-X<sub>0</sub>)],其中n,n≥0表示迭代次数,X<sub>0</sub>表示高分辨图像的初始化估计,X<sub>n</sub>表示第n次迭代后的高分辨图像,X<sub>n+1</sub>表示第n+1次迭代的高分辨图像,F表示输入的低分辨图像,α表示梯度下降的迭代步长,B和D分别表示模糊和下采样操作,T表示矩阵的转置操作,B<sup>T</sup>和D<sup>T</sup>分别表示去模糊和上采样操作,μ表示反向投影约束项与先验惩罚项的平衡参数;(11)对输入的低分辨彩色图像F在YCbCr颜色空间中的两个色差分量Cb和Cr直接采用双立方插值进行3倍放大,并组合亮度分量Y的超分辨估计X<sub>n</sub>,将YCbCr颜色空间的超分辨结果转换到RGB颜色空间,得到输入的低分辨图像F被放大了3倍的RGB图像。
地址 710071 陕西省西安市太白南路2号