一种基于高斯过程回归的图像超分辨率重建方法与流程

本发明属于图像处理技术领域，特别涉及高斯过程回归的图像超分辨率方法和改进模糊c均值聚类方法以及两种方法的结合。

背景技术：

图像超分辨率重建简单的来说就是由一幅或多幅低分辨率图像经过一系列处理之后生成高分辨率图像的过程。图像超分辨率技术目前已经广泛的应用于社会生活生产的众多领域，例如医学成像、遥感雷达成像、车辆监控、刑事侦查以及高清数字电视等。

目前图像超分辨率技术可以大致分为三类：基于插值的方法、基于重建的方法、基于学习的方法。基于插值的方法能够简单快速的重建出高分辨率图像，但是由于该方法重建过程细节信息丢失严重，所以重建出的图像质量有限。基于重建的方法相较基于插值的方法重建的图像质量有一定的提高，但是图像重建工作量大难，重建图像高频信息不完整。基于学习的方法能够利用先验知识详细描述高分辨率图像与低分辨图像间的映射关系，可以重建出质量较高、细节信息完整高分辨图像。

高斯过程作为一种强大的统计学习工具常用解决非线性问题。近年来，研究人员发高斯过程回归在解决超分辨率问题上有非常良好的应用。何河等人在文献he,siu.singleimagesuper-resolutionusinggaussianprocessregression[c]//proceedingsoftheieeeconferenceoncomputervision&patternrecognition.ieee,2011:499-456.中利用图像自身特点，提出一种自学习框架，来学习图像内部关系，进而重建出高分辨率图像。王建军等人在文献wangh,gaox,zhangk,etal.imagesuper-resolutionusingnon-localgaussianprocessregression[j].neurocomputing,2016,194:95-106.提出借用网格采样的方法学习图像块上的非局部相似性的高斯过程回归方法，处理高分辨率重建问题。曲延云等人在文献quyy,liaomj,zhouyw,etal.imagesuper-resolutionbasedondata-drivengaussianprocessregression[c]//internationalconferenceonintelligentscience&bigdataengineering.2013:513-520.中提出了基于锚点的高斯过程回归超分辨方法。

上述这些方法虽然说能恢复出低分辨率图像缺失的细节信息，但由没有考虑图像特征块的相似性，以至于构建出高斯过程回归模型预测的准确性不高，进而导致最终超分辨重建图像质量不佳。

技术实现要素：

本发明的目的在于解决超分辨率重建技术目前存在的问题。提出一种结合高斯过程回归和改进模糊c均值聚类的图像超分辨率重建方法，该方法能够有效提高图像重建质量。

一种基于高斯过程回归的图像超分辨率重建方法，其特征在于：包括训练阶段和测试阶段两部分，其中所述训练阶段的步骤为：

s1、将高分辨率样本图像集h＝(h1,h2,...,hn)进行处理得到高频低分辨率图像集其中hi表示第i个高分辨率样本图像，表示第i个高频低分辨率图像；具体为，将高分辨率样本图像集h＝(h1,h2,...,hn)利用模糊核大小为7×7、标准方差为1.1的高斯模糊核函数进行处理，然后进行3倍下采样处理，得到低分辨率图像集l＝(l1,l2,...,ln)，其中li表示第i个低分辨率图像，对低分辨图像集进行双三次插值算法放大处理至高分辨率图像大小，得到高频低分辨率图像集

s2、分别对步骤s1中的高分辨率图像集h、高频低分辨率图像集lf进行分块操作，得到样本训练集其中pl^h表示高分辨率图像块，表示低分辨率图像块，m表示图像块的个数；

s3、将步骤s2得到的样本训练集聚类为c个簇类，用m＝{m1,m2,...,mc}表示，其中mi表示一个簇集合；

(s3a)设置加权指数n，终止阈值ε，最大迭代次数rmax，当前迭代次数r＝1；

(s3b)设置类簇个数c，以及初始聚类中心v⁽⁰⁾；

(s3c)设置di＝1,i＝1,2,...,c，di表示第i个类簇的分散度值；

(s3d)计算隶属度矩阵其中xk表示第k个样本，vt表示第t个聚类中心；

(s3e)计算各类的分散度值其中ci表示第i个类簇，vi表示第i个聚类中心，n表示样本的数量；

(s3f)计算各类簇的聚类中心其中表示第k个样本到t个聚类中心的隶属度值；

(s3g)若||vi-v(i-1)||＜ε或r＞rmax，则停止迭代，输出聚类结果；否则r＝r+1，返回步骤(s3d)。

s4、基于聚类集合m中的每个子集合分别训练高斯过程回归模型g；

(s4a)对于输入样本x与测试样本x^*，定义的核函数k(x,x^*)为：

其中，||·||表示约束范式，σn表示高斯核的标准差；

(s4b)σn由下式自适应确定，即：

其中，ρ为比例系数；

(s4c)假设训练样本集d＝{xi,yi|i＝[1,n]}得到观察值列向量y＝{y1,y2,...yn}，那么可以计算k(x,x^*)的协方差矩阵：

其中xi、yi分别表示第i个低分辨率样本和第i个高分辨率样本；

(s4d)高斯过程回归模型采用的似然函数定义如下：

其中y、y^t分别表示输入样本和y的转置；

(s4d)极小化l是一个非凸的优化问题，利用共轭梯度法优化算法，求得超参数最优解。

所述测试阶段的步骤为：

s5、将待测试低分辨率图像it插值放大至需要大小，记为低分辨率插值图像ih；

s6、将低分辨率插值图像ih进行分块x＝{x1,x2,...,xn}，将每个图像块xi作为输入，利用高斯过程回归模型预测xi中丢失的像素值y^*；

s7、将y^*插入图像块xi；

s8、根据图像块对应的位置将各个块组合起来重建出高分辨率图像。

本发明与现有的超分辨技术相比,本发明具有以下显著优点：

第一，本发明分析样本数据间的相似性与差异性，通过将相似度高的数据分为一类进行高斯过程回归模型训练,提升了高斯过程回归模型预测的精度。

第二，本发明与现有的其它图像超分辨方法相比,实现了在视觉体验上更接近原始高分辨图像，呈现出了非常优质的超分辨率重建效果。

第三，通过模糊c均值聚类将事先准备的训练样本按照图像块之间结构相似关系分为c类，使得类内的样本具有较高的相似性，类与类之间的样本具有较高的差异性。通过提高训练样本的紧密关系进而提高构建的高斯过程回归模型预测值的精确度，使得重建出的图像质量更好、更符合人们的视觉体验。

附图说明

为了使本发明的目的、技术方案和有益效果更加清楚，本发明提供如下附图进行说明：

图1为本发明提供实施方法的流程图；

图2为本发明利用baby图像进行超分辨率重建和其他算法进行超分辨率重建的对比图；

图3为本发明利用bird图像进行超分辨率重建和其他算法进行超分辨率重建的对比图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、详细地描述。所描述的实施例，仅仅是本发明的一部分实施例。

本发明的具体方法参见图1，

本发明的训练阶段具体方法步骤如下所示：

s1、将高分辨率样本图像集h＝(h1,h2,...,hn)，hi表示第i个高分辨率样本图像利用模糊核大小为7×7、标准方差为1.1的高斯模糊核函数进行处理，然后进行3倍下采样处理，得到低分辨率图像集l＝(l1,l2,...,ln)，li表示第i个低分辨率样本图像，对低分辨图像集进行双三次插值算法放大处理至高分辨率图像大小，得到高频低分辨率图像集表示第i个高频低分辨率样本图像；

s2、分别对步骤s1中得到的高分辨率图像集h、高频低分辨率图像集lf集进行分块操作，得到样本训练集其中m分别表示高分辨率图像块、低分辨率图像块以及图像块的个数；；

s3、将步骤s2得到的样本训练集聚类为c个簇类用m＝{m1,m2,...,mc}表示，其中mi表示一个簇集合；

(s3a)设置加权指数n，终止阈值ε，最大迭代次数rmax，当前迭代次数r＝1；

(s3b)设置类簇个数c，以及初始聚类中心v⁽⁰⁾；

(s3c)设置di＝1,i＝1,2,...,c，di表示第i个类簇的分散度值；

(s3d)计算隶属度矩阵其中xk、vt分别表示第k个样本、第t个聚类中心；；

(s3e)计算各类的分散度值其中ci表示第i个类簇，vi表示第i个聚类中心，n表示样本的数量；；

(s3f)计算各类簇的聚类中心其中表示第k个样本到t个聚类中心的隶属度值；；

(s3g)如果||vi-v(i-1)||＜ε或r＞rmax，则停止迭代，输出聚类结果；否则r＝r+1，返回步骤(s3d)；

s4、基于m聚类集合中的每个子集合训练高斯过程回归模型g。

(s4a)对于输入样本x与测试样本x^*，所定义的核函数k(x,x^*)为：

其中，||·||表示约束范式，σn表示高斯核的标准差；

(s4b)σn由下式自适应确定，即：

其中，ρ为比例系数；

(s4c)假设训练样本集d＝{xi,yi|i＝[1,n]}得到观察值列向量y＝{y1,y2,...yn}，那么可以计算k(x,x^*)的协方差矩阵：

其中xi、yi分别表示第i个低分辨率样本和第i个高分辨率样本；

(s4d)高斯过程回归模型采用的似然函数定义如下：

其中y、y^t分别表示输入样本和y的转置；

(s4d)极小化l是一个非凸的优化问题，利用共轭梯度法优化算法，求得超参数最优解。

本发明的测试阶段具体方法步骤如下所示：

s5、将待测试低分辨率图像it使用双三次插值算法插值放大至需要大小，记为低分辨率插值图像ih；

s6、将低分辨率插值图像ih进行分块x＝{x1,x2,...,xn}，其中xi表示第i个图像块；

s7、根据测试阶段的步骤s3计算每个图像块xi隶属于mi，测试阶段s4步骤学习到的该类样本高斯过程回归模型预测xi中丢失的像素值y^*；

s8、将y^*插入图像xi；

s9、根据图像块的位置将其组合起来重建出高分辨率图像。

为了评估本算法的有效性，本发明通过采用自然图像bird、人物图像baby两幅不同种类的图像进行超分辨率重建。本发明通过与bicubic算法、srgpr算法、scsr算法进行仿真测试结果的对比。

实验一，用本发明和上述现有的3种方法对人物图像baby，进行实验仿真测试，图2是实验结果。其中图2(a)是bicubic算法重建出的图像，图2(b)是srgpr算法重建出的图像，图2(c)是scsr算法重建出的图像，图2(d)是本发明算法重建出的图像，图2(e)是原始高分辨率图像。通过分析个算法重建结果图发现本发明处理baby图像上边缘细节丰富的眼部区域有着非常好的效果，相比其它算法更能接近原始高分辨率图像。

实验二，用本发明和上述现有的3种方法对人物图像baby，进行实验仿真测试，图3是实验结果。其中图3(a)是bicubic算法重建出的图像，图3(b)是srgpr算法重建出的图像，图2(c)是scsr算法重建出的图像，图3(d)是本发明算法重建出的图像，图3(e)是原始高分辨率图像。可以看出对比其它算法本发明重建出的图像质量更好

以上这些实施例应理解为仅用于说明本发明而不用于限制本发明的保护范围。在阅读了本发明的记载的内容之后，技术人员可以对本发明作各种改动或修改，这些等效变化和修饰同样落入本发明权利要求所限定的范围。