CN105354804A - 一种基于最大化自相似性质的图像超分辨率重建方法 - Google Patents
一种基于最大化自相似性质的图像超分辨率重建方法 Download PDFInfo
- Publication number
- CN105354804A CN105354804A CN201510701726.2A CN201510701726A CN105354804A CN 105354804 A CN105354804 A CN 105354804A CN 201510701726 A CN201510701726 A CN 201510701726A CN 105354804 A CN105354804 A CN 105354804A
- Authority
- CN
- China
- Prior art keywords
- image
- sigma
- super
- resolution
- image sheet
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于最大化自相似性质的图像超分辨率重建方法,包括:对输入图像<i>Y</i>进行高斯低通滤波和双三次上采样,得到高斯低通滤波图像<i>Y0</i>和双三次上采样图像<i>X0</i>;从图像<i>Y</i>和图像<i>Y0</i>中抽取图像片集合{<i>yp</i>,<i>y0 p</i>},进而构建相应的训练样本数据集合D;采用混合高斯模型来描述数据集合<i>D</i>的概率密度函数,然后采用期望最大化算法进行参数训练,得到概率密度函数的参数;从图像<i>X0</i>中抽取图像片x0<i>p</i>,然后根据概率密度函数的参数对图像片x0<i>p</i>进行求解,得到在待求取的超分辨率图像X中以最大概率重现的图像片<i>zi</i>;根据图像片<i>zi</i>以及图像X需满足的退化处理约束方程求出最终的超分辨率图像X。本发明具有失真度小和噪声小的优点,可广泛应用于图像处理领域。
Description
技术领域
本发明涉及图像处理领域,尤其是一种基于最大化自相似性质的图像超分辨率重建方法。
背景技术
图像超分辨率重建技术是图像处理领域里一个重要且基础的课题,它指的是通过软件计算的方式处理输入的一幅或多幅同一场景的低分辨率图像来估算其对应的高分辨率图像。这种重建技术在数字电视高清显示、视频监控等领域有着广泛的应用,其实现方法主要可分为基于插值的方法、基于重建的方法和基于学习的方法这三种。近年来,随着机器学习、数据挖掘等技术的崛起,基于学习的方法逐渐为人们所重视,并成为图像超分辨率重建技术的主流。
D.Glasner等人通过统计实验表明:在自然图像中,绝大多数的小图像片(如5×5大小的图像片,7×7大小的图像片,…)在该幅图像本身或其不同尺度的图像中存在大量的重复,这一性质被称作自相似性。D.Glasner等人利用这一性质,将输入图像构造成高斯金字塔序列,对输入图像中的每一个小图像片,利用不同尺度间所能找到的相似图像片来构造学习样本,把在同一尺度内能找到的图像片看成是该小图像片的不同观察值,最终通过学习样本和观察值这两种方式的结合来构造图像的超分辨率结果。MariaZontak通过统计实验进一步指出:在自然图像中,在距离某图像片位置越近的区域内重现该图像片的可能性越大,同时,平滑区域内的图像片比边缘或纹理位置的图像片更易在整幅图像中重现。量化上述先验知识作为约束条件加入到前述的自相似性框架中,其结果较D.Glasner等人提出的方法有所改善。但这两种方法存在着以下共同的缺点:其估计出的高分辨率图像中会产生一些新的“奇异”的图像片,这些“奇异”的图像片在其自身或其低尺度图像内不会重复,使得图像失真较大。这些“奇异”的图像片是因为采用上述框架的算法只使得所抽取的图像片具有自相似性质,并不能保证在所估计的超分辨率重建结果中任意抽取的每个图像片都在自身或较低尺度内存在自相似性质(图像片与图像片的重叠部分因加权求和而改变了所求结果的值,或者抽取的图像片为两个相邻估计结果片内各自独立的一部分);同时,在图像细节丰富的位置附近(如边缘、纹理等位置)会产生一些噪声来影响图像的视觉质量。
发明内容
为了解决上述技术问题,本发明的目的是:提供一种失真度小和噪声小的,基于最大化自相似性质的图像超分辨率重建方法。
本发明解决其技术问题所采用的技术方案是:
一种基于最大化自相似性质的图像超分辨率重建方法,包括:
S1、对输入图像Y进行高斯低通滤波和双三次上采样,得到高斯低通滤波图像Y0和双三次上采样图像X0;
S2、从输入图像Y和高斯低通滤波图像Y0中抽取图像片集合{yp,y0 p},进而构建相应的训练样本数据集合其中,yp为Y中抽取的图像片,y0 p为Y0中抽取图像片,p表示图像片位置,T为转置;
S3、采用混合高斯模型来描述训练样本数据集合D的概率密度函数,然后采用期望最大化算法进行参数训练,得到概率密度函数的参数;
S4、从双三次上采样图像X0中抽取图像片x0 p,然后根据概率密度函数的参数对图像片x0 p所需满足的最大值方程进行求解,得到在待求取的超分辨率图像X中以最大概率重现的图像片zi;
S5、根据图像片zi以及超分辨率图像X需满足的退化处理约束方程求出最终的超分辨率图像X。
进一步,所述步骤S3,其包括:
S31、采用混合高斯模型来描述训练样本数据集合D的概率密度函数,所述训练样本数据集合D的概率密度函数p(d)的表达式为:
其中,d为D对应的多元随机变量,K为设定的高斯函数个数,αk、μk和Σk分别为第k个高斯函数的权值、均值向量和协方差矩阵,N为正态分布函数;
S32、采用期望最大化算法进行参数训练,得到概率密度函数p(d)的参数αk、μk和Σk。
进一步,所述步骤S4,其包括:
S41、从双三次上采样图像X0中抽取图像片x0 p;
S42、对图像片x0 p所需满足的最大值方程进行求解,得到图像片zi,所述所需满足的最大值方程为:
进一步,所述步骤S42,其具体为:
将μk改写为Σk改写为 然后代入x0 p所需满足的最大值方程中,从而得到图像片zi,所述图像片zi的表达式为:
其中,i=1,2,…,N;为双三次上采样图像X0中的第i个图像片。
进一步,所述步骤S5,其具体为:
固定图像片zi的集合{zi}不变,对最终的超分辨率图像X所需满足的退化处理约束方程进行求解,得到最终的超分辨率图像X,所述X所需满足的退化处理约束方程为:
其中,D0为设定的下采样矩阵,H为设定的低通滤波矩阵,β为设定的权重系数,Pj为抽取X中第j个图像片的概率,zj为集合{zi}中的第j个图像片;
所述最终的超分辨率图像X的表达式为:
进一步,所述输入图像Y为Ycbcr空间的图像。
本发明的有益效果是:包括根据概率密度函数的参数对图像片x0 p所需满足的最大值方程进行求解,得到在待求取的超分辨率图像X中以最大概率重现的图像片zi的步骤,根据图像的最大化自相似性质进行图像超分辨率重建,使得最终所求取的超分辨率图像X中的任一图像片zi在输入图像Y中以最大的概率存在,减少了不存在重复的“奇异”图像片的数量,失真度小;增设了退化处理约束,使得最终得到的超分辨率图像X经退化处理后的图像与输入图像Y尽可能接近,在图像细节丰富的位置附近对噪声有很强的抑制作用,使得重建后的图像在视觉上更加清晰自然。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1为本发明一种基于最大化自相似性质的图像超分辨率重建方法的整体流程图;
图2为本发明实施例一基于最大化自相似性质的图像超分辨率重建方法的框架图;
图3为采用双三次插值算法得到小孩的超分辨率重建图像;
图4为采用本发明的方法得到小孩的超分辨率重建图像。
具体实施方式
参照图1,一种基于最大化自相似性质的图像超分辨率重建方法,包括:
S1、对输入图像Y进行高斯低通滤波和双三次上采样,得到高斯低通滤波图像Y0和双三次上采样图像X0;
S2、从输入图像Y和高斯低通滤波图像Y0中抽取图像片集合{yp,y0 p},进而构建相应的训练样本数据集合其中,yp为Y中抽取的图像片,y0 p为Y0中抽取图像片,p表示图像片位置,T为转置;
S3、采用混合高斯模型来描述训练样本数据集合D的概率密度函数,然后采用期望最大化算法进行参数训练,得到概率密度函数的参数;
S4、从双三次上采样图像X0中抽取图像片x0 p,然后根据概率密度函数的参数对图像片x0 p所需满足的最大值方程进行求解,得到在待求取的超分辨率图像X中以最大概率重现的图像片zi;
S5、根据图像片zi以及超分辨率图像X需满足的退化处理约束方程求出最终的超分辨率图像X。
进一步作为优选的实施方式,所述步骤S3,其包括:
S31、采用混合高斯模型来描述训练样本数据集合D的概率密度函数,所述训练样本数据集合D的概率密度函数p(d)的表达式为:
其中,d为D对应的多元随机变量,K为设定的高斯函数个数,αk、μk和Σk分别为第k个高斯函数的权值、均值向量和协方差矩阵,N为正态分布函数;
S32、采用期望最大化算法进行参数训练,得到概率密度函数p(d)的参数αk、μk和Σk。
进一步作为优选的实施方式,所述步骤S4,其包括:
S41、从双三次上采样图像X0中抽取图像片x0 p;
S42、对图像片x0 p所需满足的最大值方程进行求解,得到图像片zi,所述所需满足的最大值方程为:
进一步作为优选的实施方式,所述步骤S42,其具体为:
将μk改写为Σk改写为 然后代入x0 p所需满足的最大值方程中,从而得到图像片zi,所述图像片zi的表达式为:
其中,i=1,2,…,N;为双三次上采样图像X0中的第i个图像片,和由均值向量μk拆分而来,和由协方差矩阵Σk拆分而来。
进一步作为优选的实施方式,所述步骤S5,其具体为:
固定图像片zi的集合{zi}不变,对最终的超分辨率图像X所需满足的退化处理约束方程进行求解,得到最终的超分辨率图像X,所述X所需满足的退化处理约束方程为:
其中,D0为设定的下采样矩阵,H为设定的低通滤波矩阵,β为设定的权重系数,Pj为抽取X中第j个图像片的概率,zj为集合{zi}中的第j个图像片;
所述最终的超分辨率图像X的表达式为:
进一步作为优选的实施方式,所述输入图像Y为Ycbcr空间的图像。
下面结合说明书附图和具体实施例对本发明作进一步详细说明。
实施例一
从前面背景技术部分的介绍可知,自然图像中绝大多数的图像片都存在着自相似性质,这对于估计出的超分辨率结果图像中的任意一个图像片而言,在输入图像中都至少有一个图像片和它的外观非常相似;同时,对超分辨率结果图像进行滤波以及下采样处理后得到的图像应和输入图像非常接近。
基于上述分析,本实施例提出了一种新的基于自相似性质的图像超分辨率重建方法,如图2所示,该方法先将输入的低分辨率图像Y进行高斯低通滤波处理得到图像Y0,将Y通过双三次插值方法放大到目标尺寸得到X0,以Y0和X0作为输入图像Y和高分辨率图像X的低频成分,然后通过Y,Y0,和X0重建Y的高分辨率图像X,具体的实现过程如下:
a.首先从Y和Y0中抽取图像片集合{yp,y0 p},p表示图像片位置,并以向量的形式连接Y和Y0中的图像片来构造训练样本数据集合接着,本发明用混合高斯模型来描述D的概率密度函数p(d),p(d)的表达式为:d为多元随机变量,K为指定的高斯函数个数,αk、μk和Σk分别为第k个高斯函数的权值、均值向量和协方差矩阵。最后,在训练概率密度函数p(d)的αk、μk和Σk这三个参数的过程中,本发明采用了期望最大化算法(即EM算法)。
b.在图像X0中抽取一定数量的图像片x0 p,其所抽取的数量为X0中的总像素数。对于图像X0中的每一个图像片而言,其所需满足的最大值方程为: 求解最大值方程得到的图像片zi在X中应以最大的概率重现。
在求解zi的过程中,本实施例将μk改写为Σk改写为 然后代入最大值方程,求得的zi为:
c.为保证最终的结果图像在进行滤波与下采样处理后与输入的低分辨率图像尽可能接近,在固定图像片zi的集合的{zi}不变前提下,还需要求解如下退化处理约束方程:
其中,D0为设定的下采样矩阵,H为设定的低通滤波矩阵,β为设定的权重系数,Pj为抽取X中第j个图像片的概率,zj为集合{zi}中的第j个图像片。通过对退化处理约束方程进行求解,可得X的表达式为:
d.为了保证得到的结果稳定,本发明还需将步骤b和c迭代运行4次或5次(在每次运行前将β增大,以保证其能逐步收敛),得到最终的结果图像。
此外,若输入图像Y为RGB彩色图像时,本发明还需要先将图像转换到Ycbcr空间。其中,Y通道采用本发明提出的方法,而cb、cr则采用简单的双三次插值方法,最后将得到的结果图像再转换回RGB空间。
本发明基于图像的最大化自相似性质提出了一种新的图像超分辨率重建算法,该方法使得所求解的超分辨率图像X中的任意一个小图像片zi在输入图像Y中以最大的概率存在,同时使得到的超分辨率图像X在经过退化处理后与输入的图像Y尽可能接近。对比图3和图4可以看出,本发明的方法在边缘等图像细节丰富的位置附近对噪声有很强的抑制作用,使得重建后的图像在视觉上更加清晰自然。
以上是对本发明的较佳实施进行了具体说明,但本发明创造并不限于所述实施例,熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替换,这些等同的变形或替换均包含在本申请权利要求所限定的范围内。
Claims (6)
1.一种基于最大化自相似性质的图像超分辨率重建方法,其特征在于:包括:
S1、对输入图像Y进行高斯低通滤波和双三次上采样,得到高斯低通滤波图像Y0和双三次上采样图像X0;
S2、从输入图像Y和高斯低通滤波图像Y0中抽取图像片集合{yp,y0 p},进而构建相应的训练样本数据集合其中,yp为Y中抽取的图像片,y0 p为Y0中抽取图像片,p表示图像片位置,T为转置;
S3、采用混合高斯模型来描述训练样本数据集合D的概率密度函数,然后采用期望最大化算法进行参数训练,得到概率密度函数的参数;
S4、从双三次上采样图像X0中抽取图像片x0 p,然后根据概率密度函数的参数对图像片x0 p所需满足的最大值方程进行求解,得到在待求取的超分辨率图像X中以最大概率重现的图像片zi;
S5、根据图像片zi以及超分辨率图像X需满足的退化处理约束方程求出最终的超分辨率图像X。
2.根据权利要求1所述的一种基于最大化自相似性质的图像超分辨率重建方法,其特征在于:所述步骤S3,其包括:
S31、采用混合高斯模型来描述训练样本数据集合D的概率密度函数,所述训练样本数据集合D的概率密度函数p(d)的表达式为:
其中,d为D对应的多元随机变量,K为设定的高斯函数个数,αk、μk和Σk分别为第k个高斯函数的权值、均值向量和协方差矩阵,N为正态分布函数;
S32、采用期望最大化算法进行参数训练,得到概率密度函数p(d)的参数αk、μk和Σk。
3.根据权利要求2所述的一种基于最大化自相似性质的图像超分辨率重建方法,其特征在于:所述步骤S4,其包括:
S41、从双三次上采样图像X0中抽取图像片x0 p;
S42、对图像片x0 p所需满足的最大值方程进行求解,得到图像片zi,所述所需满足的最大值方程为:
4.根据权利要求3所述的一种基于最大化自相似性质的图像超分辨率重建方法,其特征在于:所述步骤S42,其具体为:
将μk改写为 Σk改写为 然后代入x0 p所需满足的最大值方程中,从而得到图像片zi,所述图像片zi的表达式为:
其中,i=1,2,…,N;为双三次上采样图像X0中的第i个图像片。
5.根据权利要求4所述的一种基于最大化自相似性质的图像超分辨率重建方法,其特征在于:所述步骤S5,其具体为:
固定图像片zi的集合{zi}不变,对最终的超分辨率图像X所需满足的退化处理约束方程进行求解,得到最终的超分辨率图像X,所述X所需满足的退化处理约束方程为:
其中,D0为设定的下采样矩阵,H为设定的低通滤波矩阵,β为设定的权重系数,Pj为抽取X中第j个图像片的概率,zj为集合{zi}中的第j个图像片;
所述最终的超分辨率图像X的表达式为:
6.根据权利要求1-5任一项所述的一种基于最大化自相似性质的图像超分辨率重建方法,其特征在于:所述输入图像Y为Ycbcr空间的图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510701726.2A CN105354804B (zh) | 2015-10-23 | 2015-10-23 | 一种基于最大化自相似性质的图像超分辨率重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510701726.2A CN105354804B (zh) | 2015-10-23 | 2015-10-23 | 一种基于最大化自相似性质的图像超分辨率重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105354804A true CN105354804A (zh) | 2016-02-24 |
CN105354804B CN105354804B (zh) | 2018-07-31 |
Family
ID=55330771
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510701726.2A Active CN105354804B (zh) | 2015-10-23 | 2015-10-23 | 一种基于最大化自相似性质的图像超分辨率重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105354804B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110223229A (zh) * | 2019-05-23 | 2019-09-10 | 广东外语外贸大学 | 基于尺度间自相似性的图像放大方法、系统及存储介质 |
CN112907748A (zh) * | 2021-03-31 | 2021-06-04 | 山西大学 | 一种基于非降采样剪切波变换与深度图像纹理特征聚类的三维形貌重建方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102750677A (zh) * | 2012-06-12 | 2012-10-24 | 清华大学 | 基于同尺度结构自相似与压缩感知的单图像超分辨率方法 |
CN103093444A (zh) * | 2013-01-17 | 2013-05-08 | 西安电子科技大学 | 基于自相似性和结构信息约束的图像超分辨重建方法 |
CN103455988A (zh) * | 2013-09-17 | 2013-12-18 | 西安电子科技大学 | 基于结构自相似性与稀疏表示的超分辨率图像重构方法 |
US20140093185A1 (en) * | 2012-09-28 | 2014-04-03 | Luhong Liang | Apparatus, system, and method for multi-patch based super-resolution from an image |
CN104112286A (zh) * | 2014-08-01 | 2014-10-22 | 桂林电子科技大学 | 基于几何结构特征和自相似性的图像压缩感知重构方法 |
-
2015
- 2015-10-23 CN CN201510701726.2A patent/CN105354804B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102750677A (zh) * | 2012-06-12 | 2012-10-24 | 清华大学 | 基于同尺度结构自相似与压缩感知的单图像超分辨率方法 |
US20140093185A1 (en) * | 2012-09-28 | 2014-04-03 | Luhong Liang | Apparatus, system, and method for multi-patch based super-resolution from an image |
CN103093444A (zh) * | 2013-01-17 | 2013-05-08 | 西安电子科技大学 | 基于自相似性和结构信息约束的图像超分辨重建方法 |
CN103455988A (zh) * | 2013-09-17 | 2013-12-18 | 西安电子科技大学 | 基于结构自相似性与稀疏表示的超分辨率图像重构方法 |
CN104112286A (zh) * | 2014-08-01 | 2014-10-22 | 桂林电子科技大学 | 基于几何结构特征和自相似性的图像压缩感知重构方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110223229A (zh) * | 2019-05-23 | 2019-09-10 | 广东外语外贸大学 | 基于尺度间自相似性的图像放大方法、系统及存储介质 |
CN112907748A (zh) * | 2021-03-31 | 2021-06-04 | 山西大学 | 一种基于非降采样剪切波变换与深度图像纹理特征聚类的三维形貌重建方法 |
CN112907748B (zh) * | 2021-03-31 | 2022-07-19 | 山西大学 | 一种基于非降采样剪切波变换与深度图像纹理特征聚类的三维形貌重建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105354804B (zh) | 2018-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110119780A (zh) | 基于生成对抗网络的高光谱图像超分辨重建方法 | |
CN103871041B (zh) | 基于认知正则化参数构建的图像超分辨率重构方法 | |
CN103077511B (zh) | 基于字典学习和结构相似的图像超分辨率重建方法 | |
CN102156875B (zh) | 基于多任务ksvd字典学习的图像超分辨率重构方法 | |
CN103093444B (zh) | 基于自相似性和结构信息约束的图像超分辨重建方法 | |
CN106952228A (zh) | 基于图像非局部自相似性的单幅图像的超分辨率重建方法 | |
CN108428212A (zh) | 一种基于双拉普拉斯金字塔卷积神经网络的图像放大方法 | |
CN101794440B (zh) | 图像序列的加权自适应超分辨率重建方法 | |
CN105046672A (zh) | 一种图像超分辨率重建方法 | |
CN106709875A (zh) | 一种基于联合深度网络的压缩低分辨率图像复原方法 | |
CN109741256A (zh) | 基于稀疏表示和深度学习的图像超分辨率重建方法 | |
WO2021022929A1 (zh) | 一种单帧图像超分辨率重建方法 | |
CN109671023A (zh) | 一种人脸图像超分辨率二次重建方法 | |
CN106228512A (zh) | 基于学习率自适应的卷积神经网络图像超分辨率重建方法 | |
CN105069825A (zh) | 基于深度置信网络的图像超分辨率重建方法 | |
CN106920214A (zh) | 空间目标图像超分辨率重建方法 | |
CN104103052B (zh) | 一种基于稀疏表示的图像超分辨率重建方法 | |
CN107220957B (zh) | 一种利用滚动导向滤波的遥感图像融合方法 | |
CN104574328A (zh) | 一种基于直方图分割的彩色图像增强方法 | |
CN105631807A (zh) | 基于稀疏域选取的单帧图像超分辨重建方法 | |
CN105243670A (zh) | 一种稀疏和低秩联合表达的视频前景对象精准提取方法 | |
CN102842115A (zh) | 基于双重字典学习的压缩感知图像超分辨率重建方法 | |
CN105550989A (zh) | 基于非局部高斯过程回归的图像超分辨方法 | |
CN105427259A (zh) | 多方向加权tv和非局部自相似性正则化图像去模糊方法 | |
CN102063712A (zh) | 基于子带结构的多曝光图像融合方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |