一种基于权重小波变换的图像多分辨率重建方法技术

本发明专利技术提供一种基于权重小波变换的图像多分辨率重建方法，利用小波变换的多分辨率特性，将低分辨率图像分解成低频子图和高频子图，并根据高、低频子图各自的图像特征分别利用非均匀插值和改进的快速插值法进行图像重建，最后利用基于平均梯度的权重小波逆变换对高、低频子图的重建图像进行小波逆变换，得到最终的高分辨率图像；本发明专利技术利用改进的快速插值方法进行高频子图的超分辨率，利用高斯卷积核对图像的“空洞”信息进行估计填充，对不同方向的高频信息具有较好的处理效果；同时使用基于平均梯度的权重小波逆变换对高、低频子图的重建图像进行小波逆变换，可以自适应的对图像的细节信息进行增强，使重建图像信息丰富度高、细节突出。

A Multiresolution Image Reconstruction Method Based on Weighted Wavelet Transform

【技术实现步骤摘要】
一种基于权重小波变换的图像多分辨率重建方法
本专利技术涉及图像处理
，具体而言，尤其涉及一种基于权重小波变换的图像多分辨率重建方法。
技术介绍
近年来，空间图像超分辨率重建技术一直是学界的研究热点，其中多帧图像超分辨率重建技术是通过对多幅具有互补信息的低分辨图像进行处理，重建得到一幅高分辨率图像，在安全、监控、计算机视觉、军事侦查、医学成像等领域有重要的应用价值。小波变换具有多分辨率分析功能和逐步分解等特性，能够在图像重建过程中全面利用图像的信息，使得超分辨率重建后的图像信息丰富且细节清晰，在图像超分辨率重建中有广泛的应用。陈等.使用小波变换与FPGA实现了每秒30帧1024×1024图像的处理速度。张等.提出一种基于凸集投影法和复数小波包域的遥感图像上采样算法，实验证明，该方法对在遥感图像的高频细节部分重建和去模糊降噪方面，都有较好效果。孙等.将小波变换能够提取图像高频信息的特点与深度网络优异的图像重建性能相结合，提出一种结合小波变换与深度网络的图像超分辨率方法,能够良好地恢复图像细节。唐等.提出一种结合DWT和SWT的小波方法(D-S)，利用由SWT得到的高频频带来修正DWT得到的高频频带，并通过逆离散小波变换(IDWT)组合修正的高频频带和输入图像，得到高分辨率图像。张等.提出了在小波域内实现图像的超分辨率复原的方法，利用低分辨率图像序列、超分辨率图像和观测模型在不同尺度和方向上的特性，使超分辨率复原结果能够达到自适应边缘保持的目的。Shin等.提出一种结合离散小波变换(DWT)与基于实例的超分辨率技术，该方法通过搜索与输入的低分辨率图像的子带中的片段匹配的高频片来估计高频子带中的系数，实验表明该方法具有去模糊和细节增强作用。
技术实现思路
由于小波变换具有多分辨率分析功能和逐步分解等特性，本文提出一种基于小波变换的图像多分辨率重建的超分辨率方法，利用小波变换的多分辨率特性，将低分辨率图像分解成低频子图和高频子图，并根据高、低频子图各自的图像特征分别利用非均匀插值和改进的快速插值法进行图像重建，最后利用基于平均梯度的权重小波逆变换对高、低频子图的重建图像进行小波逆变换，得到最终的高分辨率图像。本专利技术采用的技术手段如下：一种基于权重小波变换的图像多分辨率重建方法，包括以下步骤：步骤S1：获取多张源图像，对源图像进行小波变换，得到低频子图LA和高频子图HV、HH、HD；步骤S2：利用非均匀插值方法对步骤S1中的低频子图LA进行多帧超分辨率重建，重建得到高分辨率图像HLA；步骤S3：利用改进的快速插值方法对步骤S1中的高频子图HV、HH、HD进行多帧超分辨率重建，重建得到对应的高分辨率图像HHV，HHH，HHD；步骤S4：利用基于平均梯度的权重小波逆变换对高、低频子图的重建图像HLA，HHV，HHH，HHD进行小波逆变换，得到最终的高分辨率图像。进一步地，所述步骤S3中基于改进的快速插值方法具体步骤如下：步骤S31：分别将对步骤S1中的高频子图HV、HH、HD进行图像配准，得到三组对应的配准参数W1、W2、W3；步骤S32：利用步骤S31中的配准参数W1、W2、W3对分别对高频子图进行映射插值，得到对应的初始估计图H01、H02、H03；步骤S33：使用高斯卷积核作为模板对步骤S32中的初始估计图H01、H02、H03进行迭代卷积，得到高分辨率图像HHV，HHH，HHD。进一步地，所述步骤S4中基于平均梯度的权重小波逆变换具体步骤如下：步骤S41：分别计算步骤S3中的高频子图重建图像HHV，HHH，HHD的平均梯度▽GHH、▽GHV和▽GHD，图像平均梯度▽G公式如下：其中，f表示图像的像素值，(i，j)是像素的坐标，M和N代表图像的行数和列数；步骤S42：利用步骤S41得到的图像平均梯度▽GHH、▽GHV和▽GHD，计算图像HHV，HHH，HHD的权重因子αHH、αHV和αHD，权重因子的计算公式如下：步骤S43：利用步骤S42中的权重因子αLH、αHL、αHH和高、低频子图的重建图像HLA，HHV，HHH，HHD进行基于平均梯度的权重小波逆变换，得到最后的高分辨率图像H，基于平均梯度的权重小波逆变换公式如下：H＝ISWT(HLL,αHHHHH,αHVHHV,αHDHHD)其中，ISWT代表平稳小波逆变换。较现有技术相比，本专利技术具有以下优点：1、为了解决均匀分布卷积模板在图像重构过程中，使对角线高频信息模糊的现象，使用适用性更强的高斯卷积模板，从此对各方向的高频信息都有良好的重建效果。2、在传统小波逆变换基础上提出使用基于图像平均梯度的权重小波逆变换，可以根据图像不同方向细节信息的强弱自适用的对图像高频信息进行增强，使重建图像信息丰富、细节凸显。基于上述理由本专利技术可在图像处理等领域广泛推广。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做以简单地介绍，显而易见地，下面描述中的附图是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术的流程示意图。图2为本专利技术与其他算法针对仿真图像的超分辨率重建效果对比图。图3为本专利技术与其他算法针对实拍图像的超分辨率重建效果对比图。具体实施方式为了使本
的人员更好地理解本专利技术方案，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分的实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本专利技术保护的范围。需要说明的是，本专利技术的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本专利技术的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。实施例为了验证本专利技术超分辨率重建的效果，选取仿真图像和实拍图像作为测试集，同时与Bicubic、IBP、POCS、MAP、L1-SAR、D-S算法的实验结果从主观角度和客观角度对算法的重建效果对比分析。具体步骤和原理如下：如图1所示，本专利技术提供了一种基于权重小波变换的图像多分辨率重建方法，包括以下步骤：步骤S1：获取多张源图像，对源图像进行小波变换，得到低频子图LA和高频子图HV、HH、HD；步骤S2：利用非均匀插值方法对步骤S1中的低频子图LA进行多帧超分辨率重建，重建得到高分辨率图像HLA；步骤S3：利用改进的快速插值方法对步骤S1中的高频子图HV、HH、HD进行多帧超分辨率重建，重建得到对应的高分辨率图像HHV，HHH，HHD；步骤S3中基于改进的快速插值方法重建过程，包括以下步骤：步骤S31：分别将对步骤S1中的高频子图HL、LH、HH进行图像配准，得到三组对应的配准参数W1、W2本文档来自技高网...

【技术保护点】
1.一种基于权重小波变换的图像多分辨率重建方法，其特征在于，包括以下步骤：步骤S1：获取多张源图像，对源图像进行小波变换，得到低频子图LA和高频子图HV、HH、HD；步骤S2：利用非均匀插值方法对步骤S1中的低频子图LA进行多帧超分辨率重建，重建得到高分辨率图像HLA；步骤S3：利用改进的快速插值方法对步骤S1中的高频子图HV、HH、HD进行多帧超分辨率重建，重建得到对应的高分辨率图像HHV，HHH，HHD；步骤S4：利用基于平均梯度的权重小波逆变换对高、低频子图的重建图像HLA，HHV，HHH，HHD进行小波逆变换，得到最终的高分辨率图像。

【技术特征摘要】
1.一种基于权重小波变换的图像多分辨率重建方法，其特征在于，包括以下步骤：步骤S1：获取多张源图像，对源图像进行小波变换，得到低频子图LA和高频子图HV、HH、HD；步骤S2：利用非均匀插值方法对步骤S1中的低频子图LA进行多帧超分辨率重建，重建得到高分辨率图像HLA；步骤S3：利用改进的快速插值方法对步骤S1中的高频子图HV、HH、HD进行多帧超分辨率重建，重建得到对应的高分辨率图像HHV，HHH，HHD；步骤S4：利用基于平均梯度的权重小波逆变换对高、低频子图的重建图像HLA，HHV，HHH，HHD进行小波逆变换，得到最终的高分辨率图像。2.根据权利要求1所述的基于权重小波变换的图像多分辨率重建方法，其特征在于，所述步骤S3中基于改进的快速插值方法具体步骤如下：步骤S31：分别将对步骤S1中的高频子图HV、HH、HD进行图像配准，得到三组对应的配准参数W1、W2、W3；步骤S32：利用步骤S31中的配准参数W1、W2、W3对分别对高频子图进行映射插值，得到对应的初始估计图H01、H02...

【专利技术属性】
技术研发人员：董丽丽，晋杰，姜宇航，张萌，许文海，
申请(专利权)人：大连海事大学，
类型：发明
国别省市：辽宁,21