CN101639932B - 一种数字图像分辨率增强的方法及系统 - Google Patents

Abstract

本发明公开了一种数字图像分辨率增强的方法及系统，以解决现有视频图像分辨率低，影响图像进一步处理和分析的问题。所述方法包括：图像插值，对低分辨率图像的亮度图进行最近邻插值；改进贝叶斯图像插值模型；插值校正，用改进的贝叶斯图像插值型校正高分辨率图像的亮度图。针对彩色图像，本发明提供了图像预处理和后处理，前者对输入的低分辨率彩色图像进行颜色空间转换，后者将校正后的高分辨率亮度图和色度分量图反向转换成彩色图像。本发明针对现有贝叶斯插值方法计算量庞大、执行速度慢的问题，通过快速确定下降方向和步长，加速了插值的执行，同时本发明具有边缘保持能力，能够有效恢复图像细节，提高图像清晰度。

Description

一种数字图像分辨率增强的方法及系统

技术领域

本发明涉及图像处理领域，特别涉及一种数字图像分辨率增强的方法及系统。

背景技术

在图像处理领域，图像插值是一个重要研究方向，在图像分辨率增强、图像校正、图像旋转和去马赛克等研究方向，图像插值都是一个必不可少的阶段。图像是由像素点组成的，当放大、缩小或者旋转图像时，原始像素就会重新排列，生成一些新的像素，图像插值就是确定这些已知像素中间的未知像素，通过图像插值能减少几何变换对原有图像的破坏。

现有的图像插值方法主要分为两类：自适应方法和非自适应方法。自适应插值方法只考虑了近邻像素的相关性，而没有考虑图像的全局相关性，没有较好地识别边缘而导致边缘模糊或锯齿现象，从而引起图像细节信息的质量下降。由于人眼对于边缘细节部分特别敏感，因此，插值后图像边缘对图像质量有重要的影响。非自适应插值中研究较多的是面向边缘的插值方法。图像插值往往会使图像损失细节，造成图像边缘模糊。

发明内容

本发明所要解决的技术问题是提供一种数字图像分辨率增强的方法及系统，以解决现在有技术条件下图像插值方法中存在的插值后图像模糊、图像边缘部分丢失的问题。

为解决上述技术问题，根据本发明提供的具体实施列，本发明公开了以下技术方案：

本发明提供了一种数字图像分辨率增强的方法，包括：

图像插值，对低分辨率图像的亮度图用最近邻插值方法进行插值，得到高分辨率图像的亮度图；

贝叶斯图像插值模型改进；

插值校正，利用改进的贝叶斯图像插值模型校正高分辨率图像的亮度图，得到校正后的高分辨率图像的亮度图；

优选的，在图像插值之前，对输入的低分辨率彩色图像进行颜色空间转换图像预处理，获得低分辨率彩色图像的亮度图和色度分量图。

优选的，在插值校正之后，对高分辨彩色图像进行图像后处理，将经过校正后的高分辨率亮度图和色度分量图，反向转换成彩色图像。

优选的，图像插值时，对低分辨率彩色图像的亮度图用最近邻插值方法进行插值，对色度分量图进行插值，得到高分辨率彩色图像的亮度图和色度分量图。

优选的，所述亮度图校正的步骤包括：计算代价函数，对高分辨率图像亮度图的平滑性度量值求和；搜索代价函数下降方向；将下降方向投影到可行方向；确定下降步长；像素亮度值校正，使代价函数值减小。

优选的，所述搜索代价函数下降方向的步骤包括：利用拉普拉斯导数算子，将高分辨率图像的亮度图转换为拉普拉斯图像；将高分辨率亮度图的像素的二阶方向差分模板和拉普拉斯图像像素亮度值的卷积图像排列成向量；对像素水平、垂直、对角线和反对角线四个方向的卷积图像排成的向量求和，得到代价函数的下降方向。

优选的，所述确定下降的步骤包括：设定初始步长，判断步长在下降方向上是否可行；如果不可行，减小步长为初始步长值的二分之一；重复执行上面步骤，直到步长在下降方向上可行。

优选的，所述判断步长在下降方向上是否可行，如果下降后新的代价函数值小于原代价函数值，步长在下降方向上可行。

一种数字图像分辨率增强系统，包括：

图像插值单元，用于将低分辨率图像重建成高分辨率图像；

模型改进单元，用于改进贝叶斯图像插值模型；

插值校正单元，用于校正高分辨率图像的亮度图；

优选的，所述插值校正单元进一步包括：

下降方向搜索子单元，用于代价函数值的下降方向；

空间映射子单元，用于将下降方向投影到可行方向上；

步长确定子单元，用于确定向量的下降步长；

校正值选取子单元，用于确定高分辨率未知像素的亮度值。

与现有技术相比，本发明具有以下优点：

本发明对贝叶斯图像插值模型进行了改进，贝叶斯图像插值问题实际上是一个约束最优化的求解问题，而搜索方向和步长的确定是最优化问题求解速度的关键。发明针对精确一维搜索确定步长计算量很大，且对加速收敛作用不大的问题，提出采用简单后退准则的非精确一维搜索确定步长。对于代价函数的梯度不能用解析式表示的情形，提出寻找一个可解析表达的下降方向，使代价函数沿该下降方向收敛到极小点。搜索方向和步长的快速确定，减少了算法的庞大计算量，加速了插值算法的执行。

附图说明

图1是数字图像分辨率增强方法的流程图；

图2是低分辨率和高分辨率图像的像素对应关系图；

图3是高分辨率亮度分量初值图校正子流程图；

图4a是向量Z_k，l水平方向二阶方向差分模板；

图4b是向量Z_k，l垂直方向二阶方向差分模板；

图4c是向量Z_k，l反对角线方向二阶方向差分模板；

图4d是向量Z_k，l正对角线方向二阶方向差分模板；

图5是下采样因子q取2时，下采样矩阵D的示意图；

图6a是输入计算机的低分辨率原始图像；

图6b是原始图像分辨率增强后得到的高分辨率图像；

图6c是采用本发明提供的方法得到高分辨率图像最终效果图；

图6d是采用双线性插值方法对彩色图像分辨率增强后效果图；

图7是数字图像分辨率增强系统结构图。

具体实施方式

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

下面详细说明本发明对贝叶斯图像插值模型的改进：

首先介绍现有Huber-Markov图像模型。Huber-Markov图像模型是贝叶斯图像插值模型所采用的一种用来衡量插值图像像素平滑度的模型，要使插值后的图像有良好的平滑度，就需要对Huber-Markov模型中的指数核最优化求解。

Huber-Markov图像模型描述如下：

设低分辨率图像的分辨率为N₁×N₂，将其按行排列成N₁N₂维列向量，记作y。则高分辨率图像的分辨率为qN₁×qN₂，将其按行排列成q²N₁N₂×1维列向量，记作z。其中，q为水平和垂直方向的下采样/插值因子。下采样模型通过对高分辨率图像取q×q像素块的均值得到，

$y(i,j)=\frac{1}{q^2}\left(\underset{x=\mathrm{qi}}{\overset{q(x+1)-1}{\mathrm{\Σ}}}\underset{y=\mathrm{qj}}{\overset{q(j+1)-1}{\mathrm{\Σ}}}z(x,y)\right)---\left(1\right)$

用矩阵-向量形式表示为：

y＝Dz    (2)

其中，D为尺寸为N₁N₂×q²N₁N₂下采样矩阵。它的每一行有q²个值为1/q²的元素，将高分辨率图像的一个q×q像素块映射为一个低分辨率像素。

Huber-Markov图像模型具有边缘保持能力，其指数核定义为

$\mathrm{\Ω}(z,T)=\underset{k=0}{\overset{q{N}_1-1}{\mathrm{\Σ}}}\underset{l=0}{\overset{q{N}_2-1}{\mathrm{\Σ}}}\underset{d=0}{\overset{3}{\mathrm{\Σ}}}\mathrm{\ρ}(d_{k,l,d}^{t}z,T)---\left(3\right)$

式中，d_k，l，d ^tz表示图像局部平滑度，在图像的平滑区域，值较小，而在边缘区域，值较大。高分辨率图像的局部平滑性用二阶方向导数 衡量，每个像素z_k，l的4个二阶方向导数由二阶差分近似为

$d_{k,l,d}^{t}z=m_{k,l,d}*z_{k,l},d=\mathrm{0,1,2,3}---\left(4\right)$

式中，m_d(d＝0，1，2，3)分别表示水平、垂直、对角线和反对角线方向的二阶方向差分模板，“*”表示空间卷积。

Huber边缘惩罚函数定义为

$\mathrm{\ρ}=(x,T)=\left\{\begin{array}{cc}{x}^2,& \left|x\right|\≤T\\ 2T\left|x\right|-T^2,& \left|x\right|>T\end{array}\right.---\left(5\right)$

其中，T为阈值，当|x|≤T时，ρ(x，T)定义为二次惩罚函数，当|x|＞α时，ρ(x，T)定义为线性惩罚函数。因此，噪声的差分较小，被二次惩罚，边缘的差分较大，被线性惩罚，所以Huber惩罚函数既可以平滑插值图像的噪声，又可以保持插值图像的边缘。

相应地，贝叶斯图像插值模型等价于解决如下的约束最优化问题：

$\underset{z}{\min }\mathrm{\Ω}(z,T)=\underset{k=0}{\overset{q{N}_1-1}{\mathrm{\Σ}}}\underset{l=0}{\overset{q{N}_2-1}{\mathrm{\Σ}}}\underset{d=0}{\overset{3}{\mathrm{\Σ}}}\mathrm{\ρ}(d_{k,l,d}^{t}z,T)---\left(6\right)$

s.t.y＝Dz

图像插值算法就是在约束空间

S＝{z:y＝Dz}    (7)

上将凸的、但非二次函数Ω(z，T)最小化，得到估计的插值图像 

现有的贝叶斯图像插值模型，通过下面方法搜索下降方向和确定步长，解决Huber-Markov图像模型的约束最优化问题。

搜索下降方向：

设z为d维向量，函数Ω(z)在z处的偏导数近似为

$\∂\mathrm{\Ω}\left(\stackrel{\‾}{z}\right)/{\∂z}_j=\mathrm{\Ω}({\stackrel{\‾}{z}}_0,L,{\stackrel{\‾}{z}}_j+1,L,{\stackrel{\‾}{z}}_{d-1})-\mathrm{\Ω}({\stackrel{\‾}{z}}_0,L,{\stackrel{\‾}{z}}_j,L,{\stackrel{\‾}{z}}_{d-1})---\left(8\right)$

式中，表示向量z的第i个分量。

现有贝叶斯图像插值模型利用差分近似偏导数的数值方法计算 梯度，梯度下降最快的方向确定为搜索方向，由于代价函数Ω(z_(n)，T)的梯度没有解析形式，梯度的每一个分量仅仅能用数值方法求得，运算效率非常低。

确定步长：

现有贝叶斯图像插值模型采用精确一维搜索确定步长，通过对第n次状态的代价函数进行二阶泰勒级数近似，得到如下步长：

${\mathrm{\α}}_n=-\frac{\▿\mathrm{\Ω}{(z_{\left(n\right)},T)}^T{d}_{\left(n\right)}}{d_{\left(n\right)}^T{\▿}^2\mathrm{\Ω}(z_{\left(n\right)},T)d_{\left(n\right)}}---\left(9\right)$

其中， 是代价函数Ω(z_(n)，T)的海塞矩阵。在每一步迭代，步长α_n是代价函数最大下降到局部极小点。

现有算法确定步长计算量庞大：第一，需要计算离散型变量的海塞矩阵，第二，需要计算大型矩阵和高维向量之间的乘法，因此计算复杂且耗时。

本发明对贝叶斯图像插值模型的搜索下降方向和确定步长方法进行了改进。

改进后搜索下降方向的方法如下：

用 

表示z_k，l的拉普拉斯图像，定义为

式中，k＝0，1，L，qN₁-1；l＝0，1，L，qN₂-1。

将 

的卷积图像按行排列成q²N₁N₂维列向量，记作v_d(z_(n))。式中m_d(d＝0，1，2，3)分别表示水平、垂直、对角线和反对角线方向的二阶方向差分模板，“*”表示空间卷积。

本发明找到了一个下降方向，用g_(n)表示为

$g_{\left(n\right)}=\underset{d=0}{\overset{3}{\mathrm{\Σ}}}{\mathrm{\υ}}_d\left(z_{\left(n\right)}\right).---\left(11\right)$

通过与现有方法的比较可以看出：Huber-Markov图像模型中代价函数Ω(z_(n)，T)的梯度没有解析形式，梯度的每一个分量仅仅能用数值方法求得。现有贝叶斯插值模型利用差分近似偏导数的数值方法计算梯度，运算效率非常低。在最优化理论中，最快下降方向未必是最好的搜索方向，通常情况下，一维搜索仅仅寻找一个使目标函数值减小的下降方向。为此，本发明寻找到了代价函数的一个可解析表示的下降方向。显然，这样很大程度上减少了计算量。

改进后确定步长的方法如下：

本发明提出一种非精确的一维搜索方法确定步长。设初始步长α＝1，若αd_(n)不是可行下降方向，则减小步长α为原值的1/2(称之为后退)，直至αd_(n)是可行下降方向。

判断αd_(n)是否为可行下降方向的准则：

1)若z_(n)+αd_(n)满足约束条件，则αd_(n)是可行方向；

2)若新代价函数值小于原代价函数值，即Ω(z_(n)+αd_(n))＜Ω(z_(n))，则αd_(n)是下降方向。

很显然，本方法克服了精确一维搜索方法做出大量不必要的计算的缺陷，同时又保证了每次迭代沿可行下降方向前进，本方法的使用在很大程度上节约了时间开销。

下面结合实例详细说明完整的数字图像插值方法。

如图1所示，公布了本发明一种数字图像分辨率增强的方法的具体实施例的流程图，具体可以包括以下步骤：

步骤S101，图像插值，对低分辨率图像的亮度图进行最近邻插值，对得到高分辨率图像的亮度图。

参见图2，是本实施例下采样因子q取2时，插值低分辨率和插值后高分辨率图像的像素对应关系图。

优选的，如果输入的图像是彩色图像而非亮度图像，在图像插值之前，对输入的低分辨率彩色图像进行图像预处理，将输入的低分辨 率彩色图像进行颜色空间转换，得到低分辨率彩色图像的亮度图和色度分量图。

对于彩色图像插值问题，通常有两种解决办法：一种方法是融合彩色信息提高运动估计的精度，然后独立地在每个通道上使用灰度超分辨率方法，再从分别处理的分量图像合成彩色图像；一种方法是转换颜色空间，将色度分量和亮度分量分开，仅对亮度分量进行图像插值。

在本实施例中，采用将彩色图像从RGB空间的图像转换到YCbCr(亮度分量、蓝色色度分量、红色色度分量)，分别插值后再还原成RGB彩色图像的方法。

将RGB空间的图像转换到YCbCr的公式如下：

Y＝0.299*R+0.587*G+0.114*B    (12)

Cr＝(R-Y)*0.713+128           (13)

Cb＝(B-Y)*0.564+128           (14)

优选的，图像插值时对低分辨率彩色图像的亮度图用最近邻插值方法进行插值，对两幅色度分量图进行插值，得到高分辨率彩色图像的亮度图和色度分量图；

在本实施例中，采用简单双三次插值的方法，对低分辨率彩色图像的蓝色色度分量和红色色度分量进行超分辨率插值重建，得到高分辨率彩色图像的蓝色色度分量图和红色色度分量图。

步骤S102，改进贝叶斯图像插值模型；

本发明对贝叶斯图像插值模型进行了改进，贝叶斯图像插值问题实际上是一个约束最优化的求解问题，而搜索方向和步长的确定是最优化问题求解速度的关键。发明针对精确一维搜索确定步长计算量很大，且对加速收敛作用不大的问题，提出采用简单后退准则的非精确一维搜索确定步长。对于代价函数的梯度不能用解析式表示的情形，提出寻找一个可解析表达的下降方向，使代价函数沿该下降方向收敛到极小点。搜索方向和步长的快速确定，减少了算法的庞大计算量， 加速了插值算法的执行。

步骤S103，插值校正，利用改进的贝叶斯图像插值模型校正高分辨率图像的亮度图，得到校正后的高分辨率图像的亮度图，其子步骤流程图见图3所示。

优选的，本步骤进一步包括如下子步骤：

子步骤A1，计算代价函数，即亮度图像平滑性度量值之和。代价函数的计算公式如下：

$\mathrm{\Ω}(z,T)=\underset{k=0}{\overset{q{N}_1-1}{\mathrm{\Σ}}}\underset{l=0}{\overset{q{N}_2-1}{\mathrm{\Σ}}}\underset{d=0}{\overset{3}{\mathrm{\Σ}}}\mathrm{\ρ}(d_{k,l,d}^{t}z,T)---\left(15\right)$

式中，d_k，l，d ^tz表示图像局部平滑度，在图像的平滑区域，值较小，而在边缘区域，值较大。高分辨率图像的局部平滑性用二阶方向导数衡量，每个像素z_k，l的4个二阶方向导数由二阶差分近似为

$d_{k,l,d}^{t}z=m_{k,l,d}*z_{k,l},d=\mathrm{0,1,2,3}---\left(16\right)$

式中，m_d(d＝0，1，2，3)分别表示水平、垂直、对角线和反对角线方向的二阶方向差分模板，“*”表示空间卷积。

参见图4a-4d是本实施列中一像素分量Z_k，l的四个方向的差分模板。图4a是向量Z_k，l水平方向二阶方向差分模板，图4b是向量Z_k，l 垂直方向二阶方向差分模板，图4c是向量Z_k，l反对角线方向二阶方向差分模板，图4d是向量Z_k，l正对角线方向二阶方向差分模板。

子步骤A2，搜索代价函数下降方向；

用 

表示高分辨率亮度分量初值图z_k，l的拉普拉斯图像，定义为

式中，k＝0，1，L，qN₁-1；l＝0，1，L，qN₂-1。

将 

的卷积图像按行排列成q²N₁N₂维列向量，记 作v_d(z_(n))。

本发明找到的下降方向，用g_(n)表示为

$g_{\left(n\right)}=\underset{d=0}{\overset{3}{\mathrm{\Σ}}}{\mathrm{\υ}}_d\left(z_{\left(n\right)}\right)\left(18\right)$

子步骤A3，将下降方向投影到可行方向，投影公式如下：

d_(n)＝-Pg_(n)    (19)

其中，投影矩阵P定义如下：

P＝I-q²D^TD    (20)

其中I为单位矩阵，q为水平和垂直方向的下采样/插值因子，D为尺寸为N₁N₂×q²N₁N₂的下采样矩阵。它的每一行有q²个值为1/q²的元素，将高分辨率图像的一个q×q像素块映射为一个低分辨率像素。参见图5，是本实施列下采样因子q取2时，下采样矩阵D的示意图。

子步骤A4，确定下降步长，本发明采用非精确一维搜索方法确定下降步长，设初始步长α＝1，若αd_(n)不是可行下降方向，则减小步长α为原值的1/2(称之为后退)，直至αd_(n)是可行下降方向。

判断αd_(n)是否为可行下降方向的准则：

1)满足约束条件，则αd_(n)是可行方向；

2)代价函数值小于原代价函数值，即Ω(z_(n)+αd_(n))＜Ω(z_(n))，则αd_(n)是下降方向。

子步骤A5，像素亮度值校正，使得代价函数值减小。

令：

z_(n+1)＝z_(n)+α_nd_(n) 。    (21)

设定阈值ε，

若连续两次的迭代结果满足

e_r＝|Ω(z_(n+1)，T)-Ω(z_(n)，T)|/Ω(z_(n)，T)≤ε    (22)

则迭代终止，估计的高分辨率亮度图结果为 $\hat{z}=z_{(n+1)};$

否则转向子步骤A1继续迭代。

优选的，如果输入的原始图像为彩色图像，在插值校正后，对高分辨彩色图像进行图像后处理，将经过校正后的高分辨率亮度图和色度分量图，反向转换成彩色图像。

实施例中，将步骤S103得到的高分辨率亮度图和优选步骤中得到的色度分量图，反向转换成RGB图像，公式如下：

R＝Y+1.403*(Cr-128)                   (23)

G＝Y-0.344*(Cr-128)-0.714*(Cb-128)    (24)

B＝Y+1.773*(Cb-128)                   (25)

参见图6a，是输入计算机的低分辨率原始图像；图6b，是分辨率增强后得到的高分辨率图像；图6c，是本实施例采用本发明提供的方案对彩色图像进行分辨率增强后的最终效果图；通过比较不难看出，经过本发明的处理校正，得到了较为清晰的高分辨率图像。图6d，是采用双线性插值方法对图像进行分辨率增强后的效果图，通过本发明处理后图像与双线性插值方法处理后图像的对比，本技术方案的效果更为明显。

针对上述数字图像插值方法，本发明还提供了一种数字图像分辨率增强系统，参照图7，是实施例所述数字图像分辨率增强系统结构图，所述系统主要包括图像插值单元U101、模型改进单元U102、插值校正单元U103。

其中，所述图像插值单元U101，用于对低分辨率图像的亮度图超分辨率重建，本实施例中采用最近邻插值方法对亮度图进行插值。

优选的，本系统还可以包括图像预处理单元，用于对输入的低分辨率彩色图像进行颜色空间转换。

优选的，本系统还可以包括图像后处理单元，用于将高分辨率彩色图像亮度图、蓝色色度分量图以及红色色度分量图转换成RGB彩色图像。

本实施例中图像预处理单元采用将彩色图像从RGB颜色空间转换到YcrCb颜色空间的方法进行转换，转换后得到彩色图像的亮度图和红色、蓝色色度分量图。

所述模型改进单元U102，用于改进贝叶斯图像插值模型，针对精确一维搜索确定步长计算量很大，提出一种非精确的一维搜索方法确定步长；对于代价函数的梯度不能用解析式表示的情形，提出寻找一个可解析表达的下降方向，使代价函数值沿该下降方向收敛到极小点。搜索方向和步长的快速确定，减少了算法的庞大计算量，加速了插值算法的执行。

所述插值校正单元U103，用于校正高分辨率图像的亮度图，得到校正后的高分辨率图像的亮度图。

优选的，所述插值校正单元U103进一步包括：

下降方向搜索子单元，用于搜索代价函数值的下降方向；

空间映射子单元，用于将下降方向映射到可行方向上；

步长确定子单元，用于确定向量的下降步长；

像素亮度值校正子单元，用于确定高分辨率未知像素的亮度值。

图7所示系统中未详述的部分可以参见图1-图6所示方法的相关部分，为了篇幅考虑，在此不再详述。

Claims (6)

1.一种数字图像分辨率增强的方法，其特征在于，包括如下步骤：

图像插值，对低分辨率图像的亮度图用最近邻插值方法进行插值，得到高分辨率图像的亮度图；

贝叶斯图像插值模型改进，包括：改进搜索代价函数的下降方向的方法和确定下降步长的方法；

插值校正，利用改进的贝叶斯图像插值模型校正高分辨率图像的亮度图，得到校正后的高分辨率图像的亮度图；

改进搜索代价函数的下降方向的方法包括：

利用拉普拉斯导数算子，将高分辨率图像的亮度图转换为拉普拉斯图像；

将高分辨率亮度图像素的二阶方向差分模板和拉普拉斯图像像素亮度值的卷积图像排列成向量；

对像素水平、垂直、对角线和反对角线四个方向的卷积图像排成的向量求和，得到代价函数的下降方向；

改进确定下降步长的方法包括：

设定初始步长，判断步长在下降方向上是否可行；

如果不可行，减小步长为初始步长值的二分之一；

重复执行上面步骤，直到步长在下降方向上可行；

所述插值校正具体包括：

计算代价函数，对高分辨率图像亮度图的平滑性度量值求和；

搜索代价函数下降方向；

将下降方向投影到可行方向上；

确定下降步长；

像素亮度值校正，使代价函数值减小。

2.根据权利要求1所述的方法，其特征在于，在所述图像插值之前，对输入的低分辨率彩色图像进行颜色空间转换的图像预处理，获得低分辨率彩色图像的亮度图和色度分量图。

3.根据权利要求2所述的方法，其特征在于，在插值校正之后，对高分辨彩色图像进行图像后处理，把经过校正后的高分辨率亮度图和色度分量图，反向转换成彩色图像。

4.根据权利要求2或3所述的方法，其特征在于，图像插值时对低分辨率彩色图像的亮度图用最近邻插值方法进行插值，对色度分量图进行插值，得到高分辨率彩色图像的亮度图和色度分量图。

5.根据权利要求1所述的方法，其特征在于，所述判断步长在下降方向上是否可行具体是指：如果下降后新的代价函数值小于原代价函数值，步长在下降方向上可行。

6.一种数字图像分辨率增强系统，其特征在于，包括：

图像插值单元，用于对低分辨率图像的亮度图进行插值；

模型改进单元，用于改进贝叶斯图像插值模型，所述模型改进单元包括：搜索方向改进子单元，用于改进搜索代价函数的下降方向的方法；步长确定方法改进子单元，用于改进确定下降步长的方法；

插值校正单元，用于校正高分辨率图像的亮度图；

其中，搜索方向改进子单元包括：

转换子单元，用于利用拉普拉斯导数算子，将高分辨率图像的亮度图转换为拉普拉斯图像；

排列子单元，用于将高分辨率亮度图像素的二阶方向差分模板和拉普拉斯图像像素亮度值的卷积图像排列成向量；

求和子单元，用于对像素水平、垂直、对角线和反对角线四个方向的卷积图像排成的向量求和，得到代价函数的下降方向；

步长确定方法改进子单元包括：

步长判断子单元，用于设定初始步长，判断步长在下降方向上是否可行；

步长后退子单元，用于在设定的初始步长不可行时，减小步长为初始步长值的二分之一；

重复子单元，用于重复执行判断步长和减小步长步骤，直到步长在下降方向上可行；

所述插值校正单元包括：

代价函数计算子单元，用于计算代价函数，对高分辨率图像亮度图的平滑性度量值求和；

下降方向搜索子单元，用于搜索代价函数的下降方向；

空间映射子单元，用于将下降方向映射到可行方向；

步长确定子单元，用于确定向量的下降步长；

像素亮度值校正子单元，用于确定高分辨率未知像素的亮度值。

Images