[发明专利]基于最大期望算法的自然图像超分辨方法

说明书

技术领域

本发明属于图像处理技术领域，更进一步涉及自然图像处理技术领域中的基于最大期望算法的自然图像超分辨方法。本发明是将低分辨的自然图像进行超分辨，以得到一幅清晰的高分辨自然图像，以便为图、像后续的解译、目标识别、目标检测提供更准确、全面的信息。

背景技术

图像超分辨技术是指从单幅或多幅低分辨率的图像中重建出一幅清晰的高分辨率的图像的过程。低分辨图像空间分辨率较低，影响了对景物更全面、清晰的描述。图像超分辨的目的是为了得到高分辨图像，增强和丰富了景物的细节内容，以便为后续的图像处理与解译提供更准确、全面的信息。要获得高分辨的图像，可通过以下两个途径实现：一是采用更新的传感器，这样能获取更高质量的图像，但新型传感器的成本较高；二是通过算法对低分辨图像进行超分辨，这样能在现有传感器技术的基础上获取更高质量的图像，且技术成本低。图像超分辨技术主要可以分为基于插值、基于重建和基于学习三大类。目前，单幅图像主要基于重构误差和先验信息约束的方法来实现图像超分辨。

Yang,Wright等人发表的论文“Image Super-Resolution Via Sparse Representation”(IEEE Trans.on Image Processing vol.19no.11pp.2861-2873.2010)中提出一种基于稀疏表示的图像超分辨方法。该方法将压缩感知的思想引入超分辨重构中，通过稀疏表示的方法来获得低分辨和高分辨的字典对。当低分辨图像的稀疏表示通过压缩感知的方法得到时，则高分辨图像的稀疏表示也相应的得到，从而可以有效地对低分辨图像进行重建。但是该方法存在的不足之处是，该方法在字典对的构造中，需要采集大量的外部训练数据，而这是不现实的，同时这种假设高低分辨图像之间表示是存在误差的，这样使得恢复得到的整体效果不是很好。

He和Siu二人发表的论文“Single Image Super-Resolution using Gaussian Process Regression”(IEEE Conference on Computer Vision and Pattern Recognition,pp.449–456.2011)中公开了一种使用高斯过程先验来解决图像超分辨问题的方法。该方法使用高斯过程来学习得到低分辨图像和高分辨图像之间的映射关系。高斯过程能自适应地通过高斯分布模型来寻找得到低分辨图像块和高分辨图像块之间的内在联系。该方法的优势在于在学习映射关系的时候，并不需要高分辨图像进行训练。但是，该方法仍然存在的不足之处是，该方法仅仅利用了图像自身的局部信息作为训练数据，同时该方法没有充分挖掘局部图像块的先验信息，这样使得当高斯过程能提供的信息不足的情况下，造成了重构得到的图像结果不稳定，细节和边缘恢复的不是很好，局部区域重建质量下降。

发明内容

本发明的目的在于针对上述已有技术的不足，提出一种基于最大期望算法的自然图像超分辨方法，充分地结合图像的先验信息，以在图像超分辨重建中，能够更好地去除振铃，极大地提高了图像超分辨的恢复效果。

为实现上述目的，本发明在基于最大期望算法的基础上实现自然图像超分辨，其技术方案是通过最大期望算法将图像超分辨过程拆成两个子过程，即隐图像期望最大过程和估计图像后验概率最大过程。在估计图像后验概率最大过程中，我们使用高斯过程方法学习得到估计图像块的均值和方差，然后训练得到隐图像块的字典，最后通过均方误差最小化的方法得到估计图像的最大后验值。通过对隐图像期望最大化过程和估计图像后验概率最大过程这两个过程进行循环迭代，当迭代满足终止条件，则跳出循环，最终得到最优的高分辨图像。

本发明的具体步骤如下：

(1)输入一幅待恢复的低分辨图像；

(2)插值图像：

使用matlab软件中的imresize函数，将待恢复的低分辨图像插值到待恢复的低分辨图像的3倍，得到插值后的低分辨图像；

(3)按照下式，获得隐图像：

Z＝L+λH^T(Y-HL)

其中，Z表示隐图像，L表示插值后的低分辨图像，λ表示迭代步长，λ＝0.8，H表示观测矩阵，T表示转置操作，Y表示待恢复的低分辨图像；

(4)将隐图像切成W个隐图像块：

将隐图像进行滑窗处理，其中隐图像块的大小设为6×6像素，滑窗步长设为1个像素，得到W个隐图像块集；

(5)获得每一个隐图像块的相似矩阵：