[发明专利]一种基于非线性预测稀疏编码的单幅图像超分辨率快速重建方法

摘要

本发明公开了一种基于非线性预测稀疏编码的图像超分辨率重建方法，在训练过程中，本发明在经典的基于稀疏编码算法的准则函数基础上叠加了非线性预测编码项，并且设计了自己的优化策略来最小化该目标函数，在重建过程中，本发明仅仅对输入的低分辨率图像块和预先训练得到的低分辨率字典采用一个非线性迭代步骤直接来逼近所要求的稀疏编码，避免了对每一个小图像块求解稀疏表示系数的问题。实验结果表明，与经典的基于稀疏编码的图像超分辨率算法相比，本发明在保证了重建结果具有充分竞争力的同时大大降低了实验所需的时间。

主权项

一种基于非线性预测稀疏编码的图像超分辨率重建方法，其特征在于：重建阶段不需要对每个图像块求解一个稀疏表示问题，而仅仅采用输入的低分辨率图像块和预先训练得到低分辨率字典直接预测出稀疏编码，采用如下步骤：1)如下给定采样的训练块对{X_h,X_l}，这里X_h是高分辨率图像块，X_l是经过对应的高分辨率图像块下采样得到的低分辨率图像块。本文考虑如下准则函数：$\begin{array}{c}\min f(D_h,D_l,A)\\ \begin{array}{cc}s.t.& \left|\right|D_{h,i}|{|}_2^2\≤1,|\left|D_{l,i}\right|{|}_2^2\≤1\end{array}\end{array}---\left(1\right)$上式中，i表示字典的第i列，A代码稀疏编码；其中：$\begin{array}{c}f(D_h,D_l,A)=\frac{1}{N}\left|\right|X_h-D_{h}A|{|}_2^2+\frac{1}{N}||\tanh ({\mathrm{WX}}_h+B)-A|{|}_2^2\\ +\frac{1}{M}\left|\right|X_l-D_{l}A|{|}_2^2+\frac{1}{M}||\tanh \left({\mathrm{WX}}_l+B\right)-A|{|}_2^2+\mathrm{\λ}\left(\frac{1}{N}+\frac{1}{M}\right)\left|\right|A|{|}_1\end{array}---\left(2\right)$上式中，N和M分别是高、低分辨率图像块以列向量形式表达的维数，W、B是非线性预测模型的参数；2)优化如下问题：$\begin{array}{c}\underset{\left\{D_h,D_l,A,W,B\right\}}{\min }=\frac{1}{N}\left|\right|X_h-D_{h}A|{|}_2^2+\frac{1}{N}||\tanh \left({\mathrm{WX}}_h+B\right)-A|{|}_2^2+\frac{1}{M}\left|\right|X_l-D_{l}A|{|}_2^2+\frac{1}{M}||\tanh \left({\mathrm{WX}}_l+B\right)-A|{|}_2^2+\mathrm{\λ}\left(\frac{1}{N}+\frac{1}{M}\right)\left|\right|A|{|}_1\\ \begin{array}{cc}s.t.& \left|\right|D_{h,i}|{|}_2^2\≤1,|\left|D_{l,i}\right|{|}_2^2\≤1\end{array}\end{array}---\left(3\right)$整个优化过程如下：2‑1)令t＝0，采用高斯随机矩阵对字典D_l和D_h初始化处理，并且将字典的每一个列都作单位归一化，随机初始化W和B；2‑2)固定采用ADMM方法，更新A^(t)：$A^{\left(t\right)}=\underset{A}{\mathrm{argmin}}f(B^{\left(t\right)},W^{\left(t\right)},D_h^{\left(t\right)},D_l^{\left(t\right)},A)$2‑3)固定和A^(t)采用梯度下降法，更新W和B：$W^{(t+1)},B^{(t+1)}=\underset{W,B}{\mathrm{argmin}}f(D_h^{\left(t\right)},D_l^{\left(t\right)},A^{\left(t\right)},W^{\left(t\right)},B^{\left(t\right)})$2‑4)固定A^(t)，W^(t),B^(t)，更新和$\begin{array}{c}{D}_h^{(t+1)},D_l^{(t+1)}=\underset{D_h,D_l}{\mathrm{argmin}}f(A^{\left(t\right)},W^{(t+1)},B^{(t+1)},D_h^{\left(t\right)},D_l^{\left(t\right)})\\ \begin{array}{cc}s.t.& \left|\right|D_{h,i}|{|}_2^2\≤1,|\left|D_{l,i}\right|{|}_2^2\≤1\end{array}\end{array}$这一步利用联合字典训练思路进行优化；2‑5)令t:＝t+1，迭代2‑2)到2‑4)，直到收敛；3)利用非线性预测稀疏编码进行图像超分辨率重建时，在重建阶段不必再为每一个图像块求解一个稀疏表示问题，设输入的图像为Y，每个输入的图像块为y，每个图像块的稀疏编码可以进行如下逼近：A＝tanh(Wy+B)   (4)高分辨率块重建表达式如下：x_h＝D_hA＝D_h(tanh(Wy+B))   (5)下面给出了非线性预测稀疏编码的超分辨率流程；步骤3)中基于非线性预测稀疏编码的图像超分辨率算法：3‑1)输入：预训练得到的字典D_h和D_l，一幅低分辨率图像Y，参数W和B；3‑2)For对于Y中的每个5*5小块y，从左到右，从上到下，在每个方向上保证4；个像素重叠：①计算每个块中的平均像素值m，每个块的像素值减去均值；②计算高分辨率块x_h＝D_h(tanh(Wy+B))；③计算x_h+m，填入高分辨率图像X₀；End3‑3)在解空间中，找到与X₀最匹配的图像：$X^{*}=\underset{X}{\mathrm{argmin}}\left|\right|SHX-Y|{|}_2^2+c||X-X_0|{|}_2^2$H是一个模糊算子，S是一个下采样算子，c是平衡参数；3‑4)输出：超分辨率图像X^*。