基于变分和小波变换图像放大算法.docVIP

基于变分和小波变换的图像放大算法 冯象初1) 姜东焕 2) 徐光宝2) 1) (西安电子科技大学 理学院数学系，陕西，西安，710071) 2) (山东科技大学 信息科学与工程学院，山东，青岛，266510 ) 摘 要 为了更好地放大图像，利用小波变换的思想，提出了一种变分和小波变换相结合的图像放大算法。该算法的思想是先构造一个用Besov范数估计图像正则性的变分泛函，然后在小波域中最小化变分泛函得到放大图像。小波变换后的高频分量具有丰富的细节边缘信息，因而能够重构出高质量的图像，而且小波的引入使得文中新算法具有运行时间短、速度快的特点。与传统的插值放大图像不同，该算法是用变分的思想进行图像放大。理论分析和实验仿真表明，该算法能达到和样条插值同样的放大效果。 关键词 图像放大；变分泛函；小波变换； 中图分类号 TN911.73 引言 图像放大是从一幅低分辨率图像获得其高分辨率版本的一种图像处理技术。图像的放大和缩小在数字图像处理中起着很重要的作用，为适应特殊的应用场合或者得到一个较好的视觉效果，例如要突出某些细节，常常需要一种可以有效改变已有图像大小的方法，使图像放大和缩小后仍有较高的质量。 各种插值技术是实现图像放大的常用方法。现有的图像处理理论及常用的处理软件都是采用插值方法对图像进行放大和缩小，如平移重复插值、双线性插值和样条插值等。这些插值方法有很多的优越性，用于图像插值也取得了较好的效果。但这些方法都是用一些已知的光滑函数根据一定的光滑性要求逼近源图像。由于这种固定方式的局限性，势必会在图像放大倍数较高时形成斑点以及在明暗区域出现偏移现象，而且放大倍数越大，这种现象越明显。为克服该问题，朱宁等[1]提出了图像放大的偏微分方程方法。Guichard和Maltouyres[2,3]提出了一种基于全变差的图像插值方法。其基本思想是通过给插值图像的傅里叶系数加某种约束条件使得插值方法是可逆的，即通过对插值图像下采样可以得到原始图像。满足这样约束条件的插值方法很多，他们选用全变差估计图像的正则性并选择正则性最好的插值方法。2004年Chambolle[4]给出了基于全变差的图像放大的实现算法，得到比较理想的放大效果，但是该算法运行时间长、速度慢。 本文在Chambolle图像放大算法的基础上，用Besov范数代替图像的全变差估计图像的正则性，得到一个新的变分泛函，然后通过最小化变分泛函得到放大图像。因为Besov范数可以用小波系数序列的范数等价描述，所以可以将该泛函的最小化问题转化到小波域中用迭代的方法求解，并给出了迭代方法的收敛性证明。Chambolle基于全变差的图像放大算法运行时间长、速度慢，而文中新算法是在小波域中进行求解，小波变换后的高频分量具有丰富的细节边缘信息，因而能够重构出高质量的图像，而且小波的引入使得该算法具有运行时间短、速度快的特点。所以新算法能在比较短的时间内获得质量更好的放大图像。与传统的基于插值的放大算法不同，新算法是用变分的思想放大图像。变分模型对放大后的图像强加了几何正则性要求，这就保证放大后图像边缘的光滑性。 基础知识 2.1 小波和Besov范数 假设是定义在区域上的一副灰度图像，用Besov范数刻画图像的正则性是比较合适的[5,6]。关于Besov空间的具体定义，参考文献[5]。本文只需要Besov半范数和小波系数范数之间的等价关系。选用上的紧支撑周期正交小波(如Daubechies正交小波[7])分解图像。令表示周期化的尺度函数，那么为1。令表示由一维正交小波的张量积构成的两维正交小波基。定义 那么函数的小波表达式为： 其中，。此时，可以用小波系数来表示范数： (1) 如果小波在空间中()，那么有等价关系： 此处仅考虑，时的特殊情况，即空间。为了方便起见，记，则 (2) 其中。 2.2 Chambolle图像放大模型 假设讨论的图像是的二维矩阵，表示空间，是的一个子空间，表示一幅粗糙的图像。例如，当放大因子为2时，用数学式子表示为： 。Chambolle[4]的图像放大模型为： 其中，是放大后的图像，是在空间上的正交投影。显然有及 从而，Chambolle图像放大模型变为： Chambolle给出了求解上述最小值问题的一种迭代算法。虽然该方法能得到比较好的图像放大效果，但是它的计算时间长、速度慢，所以寻求相对运行时间短与实现速度快的算法具有重要的实际和现实意义。 基于变分和小波变换的图像放大算法及其收敛性证明 基于变分和小波变换的图像放大算法 引理[8] 设是Hilbert空间中的一组正