基于先验信息指导的子空间低秩正则图像复原方法技术

本发明专利技术涉及一种基于先验信息指导的子空间低秩正则图像复原方法，包括以下步骤：步骤S1:获取自然图像，并利用混合高斯模型学习图像块的多个子空间结构；步骤S2:将待复原图像进行重叠分块，利用学习得到的多个子空间结构对这些重叠块进行聚类，得到多个相似块组；步骤S3:针对每一个相似块组,利用该块组的子空间结构指导线性变换建立子空间；步骤S4:在每个子空间进行低秩矩阵近似；步骤S5:将所有的块组进行聚合，再结合ADMM得到复原后的图像，步骤S5:循环步骤S2

【技术实现步骤摘要】
基于先验信息指导的子空间低秩正则图像复原方法


[0001]本专利技术涉及图像处理
，具体涉及一种基于先验信息指导的子空间低秩正则图像复原方法。

技术介绍

[0002]图像复原是典型的逆问题。通常需要借助先验信息实现图像复原，大致可以分为两类：一类是利用外部清楚图像的先验(例如，全变分、稀疏)来指导图像复原；另一类是利用图像内部先验(例如，图像自相似性)来复原图像。对于破损严重的图像，通常同时结合外部和内部先验实现更好地图像复原。近年来，有大量关于这方面的研究工作。比如(PELL方法利用混合高斯模型学习外部图像的块子空间结构并结合维纳滤波实现图像复原，WNNM利用块匹配寻找相似块并利用加权低秩实现图像复原，PCLR利用外部图像块先验信息指导块聚类并结合低秩方法实现图像复原。随着即插即用框架的快速发展，IDBP方法利用去噪模型(例如BM3D、IRCNN)作为先验结合ADMM算法实现图像复原。然后，这些方法在数据丢失或噪声污染严重的情况下，仅依赖图像内部信息不能够提供足够的信息用于图像复原，强制局部相似或非局部相似容易导致复原后的图像具有震荡边缘，产生虚假纹理，从而使复原的图像不够自然。

技术实现思路

[0003]有鉴于此，本专利技术的目的在于提供一种基于先验信息指导的子空间低秩正则图像复原方法，能够提升信息丢失或噪声污染严重图像的复原质量。
[0004]为实现上述目的，本专利技术采用如下技术方案：
[0005]一种基于先验信息指导的子空间低秩正则图像复原方法，包括以下步骤：
[0006]步骤S1:获取自然图像，并利用混合高斯模型学习图像块的多个子空间结构；
[0007]步骤S2:将待复原图像进行重叠分块，利用学习得到的多个子空间结构对这些重叠块进行聚类，得到多个相似块组；
[0008]步骤S3:针对每一个相似块组,利用该块组的子空间结构指导线性变换建立子空间；
[0009]步骤S4:在每个子空间进行低秩矩阵近似；
[0010]步骤S5:将所有的块组进行聚合，再结合ADMM得到复原后的图像，
[0011]步骤S5:循环步骤S2
‑
S5，迭代至满足预设复原要求的结果，即最终的复原图像。
[0012]进一步的，所述步骤S1具体为：
[0013]设外部自然图像块形成K个子空间，每个图像块x
i
用混合高斯模型表示：
[0014][0015]其中Θ＝(w1，
…
，w
K
，θ1，
…
，θ
K
)是相关参数集合，并满足参数θ
k
刻画了一个子空间表示高斯密度函数p
k
的均值为u
k
协方差矩阵为Σ
k
；
[0016]给定外部图像训练集，从中随机图像块，通过训练混合高斯模型，得到K个子空间
[0017]进一步的，所述步骤S2具体为：对于每个图像块x
i
＝P
i
x计算其条件概率得到所在类别,其中P
i
为一个d
×
n的矩阵，表示从图像x中取第i图像块
[0018][0019]对所有的类别都计算其条件概率，取最大概率所在的类别为该图像块所在的类别；
[0020]将所有属于同一类别的图像块进行合并构成一个块组，第k类块组记为
[0021]进一步的，所述步骤S3具体为：
[0022]对子空间协方差矩阵Σ
k
进行(SVD)奇异值分解
[0023][0024]其中Ψ
k
＝diag(ψ
k，1
，
…
，ψ
k，d
)是一个对角矩阵，是是中的第i列；
[0025]引入一个正交变换矩阵用于寻找一个子空间更好实现低秩矩阵复原，其中F
k
为固定子矩阵由Σ
k
最大前b个特征值对应的特征向量构成
[0026]为学习子矩阵，由块组自适应学习得到；
[0027]构建如下自然图像先验信息指导的子空间低秩正则项S(Γ
k
(x))，在固定项F
k
的指导下，通过子空间变换矩阵Q
k
寻找最优的子空间，通过低秩矩阵近似，实现块组复原
[0028][0029]其中，，Q
k
是正交矩阵，其中F
k
是由外部自然图像块先验得到，而L
k
是从内部块组数据计算得到。
[0030]进一步的，所述步骤S4具体为：
[0031]给定待复原的图像希望复原真实的图像为对相似块组Γ
k
(x)，实现子空间低秩正则复原，图像复原模型如下：
[0032][0033][0034]其中表示退化矩阵。
[0035]进一步的，所述退化矩阵，对于图像去噪问题，A为一个单位矩阵；对于图像修补问题，A为二值掩码矩阵。
[0036]进一步的，所述步骤S5具体为：
[0037]步骤S51：给定x
t
和有求解k＝1，
…
，K；
[0038][0039]令σ
j
(Z
k
)为Z
k
的第j个单值，通过奇异值(SVD)分解得到最优解
[0040][0041]其中β
j
＝1/(σ
j
(Z
k
)+∈)，∈为一个大于0的常数，且(x)
+
＝max{x，0}
[0042]步骤S52：给定x
t
和求解k＝1，
…
，K
[0043][0044]令其中该方程与如下方程存在相同解，
[0045][0046]做奇异值(SVD)分解,其中和为正交矩阵，上面方程的唯一解表示为
[0047]步骤S53：给定和有求解x
t+1
[0048][0049]由于
[0050][0051]其中表示从图像x
t
中抽取第k类第j个图像块，Z
k
(:，j)表示Z
k
中第j列。
[0052]本专利技术与现有技术相比具有以下有益效果：
[0053]1、本专利技术提出了外部自然图像先验信息指导图像块聚集，可以聚集结构相似的图
像块，同时外部先验信息指导寻找最优图像块子空间实现低秩正则，可以更好保持图像的结构，使得复原的图像更加自然；
[0054]2、本专利技术复原框架，可以根据不同复原问题(如图像去噪、图像修补)，设置不同矩阵A，从而实现图像复原；
[0055]3、本专利技术提出了一种交替迭代求解算法，寻找最优子空间实现低秩矩阵近似，可以通过两次奇异值分解求解，实现起来简单、快速。
附图说明
[0056]图1是本专利技术方法流程图；
[00本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于先验信息指导的子空间低秩正则图像复原方法，其特征在于，包括以下步骤：步骤S1：获取自然图像，并利用混合高斯模型学习图像块的多个子空间结构；步骤S2：将待复原图像进行重叠分块，利用学习得到的多个子空间结构对这些重叠块进行聚类，得到多个相似块组；步骤S3：针对每一个相似块组，利用该块组的子空间结构指导线性变换建立子空间；步骤S4：在每个子空间进行低秩矩阵近似；步骤S5：将所有的块组进行聚合，再结合ADMM得到复原后的图像，步骤S5：循环步骤S2
‑
S5，迭代至满足预设复原要求的结果，即最终的复原图像。2.根据权利要求1所述的基于先验信息指导的子空间低秩正则图像复原方法，其特征在于，所述步骤S1具体为：设外部自然图像块形成K个子空间，每个图像块x
i
用混合高斯模型表示：其中Θ＝(w1，
…
，w
K
，θ1，
…
，θ
K
)是相关参数集合，并满足参数θ
k
刻画了一个子空间表示高斯密度函数p
k
的均值为u
k
协方差矩阵为∑
k
；给定外部图像训练集，从中随机图像块，通过训练混合高斯模型，得到K个子空间3.根据权利要求1所述的基于先验信息指导的子空间低秩正则图像复原方法，其特征在于，所述步骤S2具体为：对于每个图像块x
i
＝P
i
x计算其条件概率得到所在类别，其中P
i
为一个d
×
n的矩阵，表示从图像x中取第i图像块对所有的类别都计算其条件概率，取最大概率所在的类别为该图像块所在的类别；将所有属于同一类别的图像块进行合并构成一个块组，第k类块组记为4.根据权利要求1所述的基于先验信息指导的子空间低秩正则图像复原方法，其特征在于，所述步骤S3具体为：对子空间协方差矩阵∑
k
进行(SVD)奇异值分解其中Ψ
k
＝diag(ψ
k，1
，
…
，ψ
k，d
)是一个对角矩阵，是是中的第i列；引入一个正交变换矩阵用于寻找一个子...

【专利技术属性】
技术研发人员：陈飞，曾勋勋，
申请(专利权)人：福州大学，
类型：发明
国别省市：