一种非局部哈尔变换图像去噪方法与流程

本发明涉及图像去噪技术领域，更具体地说，特别涉及一种非局部哈尔变换图像去噪方法。

背景技术：

图像去噪一直以来都是一个重要的基础研究课题，尽管图像去噪已经被研究了几十年，但直到今天图像去噪仍然是一个比较活跃的研究方向。这是因为无论在什么环境下，无论用多么高级的设备获取的图像总会不可避免地引入各种各样的噪声，噪声的存在会对计算机视觉、模式识别、机器学习甚至人工智能等方面的研究带来或多或少的影响，所以人们总在追求更加完美的图像去噪结果。

图像去噪研究总体分为空域方法与变换域方法，这两种方法在整个研究过程中总相互补充，更迭交替地向前发展。最典型的空域图像去噪方法是高斯滤波及维纳滤波，最典型的变换域方法是基于小波变换的图像去噪方法。而近十几年发展起来的非局部方法使图像去噪的发展有了很大的进展，其中最典型的nl-means(非局部平均)属于空域方法，而bm3d(块匹配三维变换)属于非局部变换域方法，尤其是bm3d方法为当前的最先进的图像去噪方法之一，尽管近几年又有多种优秀的算法被提出，但至今没有一种算法能全面超过bm3d方法。

技术实现要素：

本发明的目的在于解决现有技术中的不足，从而提供一种不是对相似像素进行简单加权平均，从而能很好保留图像细节；提出的方法也不像bm3d那样进行了大量的复杂的小波变换与离散余弦变换致使在去噪的同时引入大量假信号，本发明仅仅是对聚类的相似像素进行了简单的哈尔变换就实现了比bm3d好得多的图像去噪结果。

为了达到上述目的，本发明采用的技术方案如下：

一种非局部哈尔变换图像去噪方法，包括以下步骤：

s1、对输入原始含噪图像执行块匹配与列扫描后的行匹配；

s2、初步降噪：对图像行匹配结果执行哈尔变换，并对变换系数采用系数硬阈值与结构硬阈值的双硬阈值进行收缩，将结果执行哈尔逆变换后对图像块进行聚合获得初步降噪图像；

s3、精细降噪：将步骤s1的原始含噪图像和步骤s2中初步降噪图像同步执行块匹配与行匹配，执行哈尔变换，得到两组变换系数，再将两组变换系数执行维纳滤波去噪，将去噪后的变换系数执行逆哈尔变换，执行图像块聚合获得最终去噪后的图像。

进一步地，所述步骤s2的具体执行方法为：

s21、进行图像块匹配和行匹配，然后对相似像素进行聚类；

s22、进行哈尔变换；

s23、进行双硬阈值系数收缩；

s24、进行哈尔逆变换与像素聚合形成初步降噪图像；

所述步骤s3的具体执行方法为：

s31、将原始含噪图像和初步降噪图像进行块匹配和行匹配，然后对相似像素进行聚类，得到两个群组结果；

s32、将步骤s31中的两个群组结果分别进行哈尔变换；

s33、进行维纳滤波去噪；

s34、进行哈尔逆变换与像素聚合形成最终降噪图像。

进一步地，所述步骤s21的具体执行方法为：首先执行原始含噪图像的块匹配,按一个指定的步长n_step提取大小为n1×n1的图像块作为参考块，然后在以该参考块为中心的一个大小为ns×ns的邻域内进行块匹配获得数量为n2的相似图像块，将匹配的所有图像块进行列扫描拼接成一个大小为(n1×n1)×n2的矩阵m，在m上用每一行作为参考行与其余的所有行计算欧式距离d，将距离最小n3行群组获得最大程度的图像中相似像素的聚类。

进一步地，所述步骤s22的具体执行方法为：步骤s21获得的行群组执行可分的二维哈尔变换，即纵向与横向各自执行哈尔变换。

进一步地，所述步骤s23的具体执行方法为：将最相似的行之间执行哈尔变换后的结果系数进行结构硬阈值，即将这些系数全部设置为零；对其余哈尔变换的系数采用传统的系数硬阈值进行系数收缩。

进一步地，所述步骤s23中系数硬阈值参数为thr，

当输入图片噪声标准偏差小于等于25时，thr＝1.65，迭代执行步骤s21-s24，2-3次；

当输入图片噪声标准偏差大于25且小于等于75时，thr＝1.45，迭代执行步骤s21-s24，4-6次；

当输入图片噪声标准偏差大于75时，thr＝1.45，迭代执行步骤s21-s24，7-12次。

进一步地，所述步骤s24的具体执行方法为：对双硬阈值收缩后的哈尔变换系数执行逆哈尔变换，获得去噪后的二维矩阵，再将此二维矩阵扫描回原来的图像块，最终将这些图像块加权放回到图像中的原始位置获得去噪后的图像。

进一步地，所述步骤s31的具体执行方法为：用步骤s24的结果图像作为参考,同步执行步骤s21结果图像与原始含噪图像的块匹配操作及行匹配操作，按一个指定的步长n_step在第一阶段结果图像中提取大小为n1×n1的图像块作为参考块，然后在以该参考块为中心的一个大小为ns×ns的邻域内进行块匹配获得数量为n2的相似图像块，将匹配的所有图像块进行列扫描拼接成一个大小为(n1×n1)×n2的矩阵m，在m上用每一行作为参考行与其余的所有行计算欧式距离d，将距离最小n3行群组获得最大程度的图像中相似像素的聚类，原始含噪图像进行同步的同参数的操作，经过块匹配与行匹配后获得两个同样大小的行群组结果。

进一步地，所述步骤s32中，将步骤s31中的两个行群组各自执行二维哈尔变换获得两组变换系数,所述步骤s33中，用步骤s2中初步降噪的结果的变换系数θ1对原始含噪图像的变换系数θ2按下述公式执行维纳滤波：

其中σ为噪声标准偏差。

进一步地，对维纳滤波后的含噪图像的系数执行逆哈尔变换，获得去噪后的二维矩阵，再将此二维矩阵扫描回原来的图像块，将这些图像块加权放回到图像中的原始位置获得最终去噪后的图像。

本发明首先执行图像块匹配获得一定数目的图像块，将获得的图像按行或列扫描拼接成一个二维矩阵，再执行此二维矩阵上的行匹配或列匹配获得最相似的一定数目像素行。通过这种方法最大程度地将图像中的自相似信息进行了聚类。对聚类的像素执行简单的哈尔变换，由于哈尔变换存在简单的提升方法，即只是对这些行之间进行加权平均，所以很容易硬件实现。整个算法不需要任何矩阵乘法，与bm3d方法相比大大降低了运算复杂度。

在去噪过程中本发明提出一种双硬阈值策略，即不仅对变换系数进行一般的硬阈值实现系数收缩，而且将最相似的两行像素的哈尔变换的系数全部置为零(本发明称这种阈值为结构硬阈值)，从而大大提高了图像去噪能力。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明的非局部哈尔变换图像去噪方法的方法原理图；

图2是本发明的非局部哈尔变换图像去噪方法的方法流程图；

图3是本发明实施例中原始图像和含噪图像；

图4是本发明实施例中bm3d去噪结果图；

图5是本发明实施例中非局部哈尔变换图像去噪方法去噪结果图。

具体实施方式

下面结合附图对本发明的优选实施例进行详细阐述，以使本发明的优点和特征能更易于被本领域技术人员理解，从而对本发明的保护范围做出更为清楚明确的界定。

实施例一

参阅图1所示，一种非局部哈尔变换图像去噪方法，包括以下步骤：

s1、对输入原始含噪图像执行块匹配与列扫描后的行匹配；

s2、初步降噪：对图像行匹配结果执行哈尔变换，并对变换系数采用系数硬阈值与结构硬阈值的双硬阈值进行收缩，将结果执行哈尔逆变换后对图像块进行聚合获得初步降噪图像；

s3、精细降噪：将步骤s1的原始含噪图像和步骤s2中初步降噪图像同步执行块匹配与行匹配，执行哈尔变换，得到两组变换系数，再将两组变换系数执行维纳滤波去噪，将去噪后的变换系数执行逆哈尔变换，执行图像块聚合获得最终去噪后的图像。

所述步骤s2的具体执行方法为：

s21、进行图像块匹配和行匹配，然后对相似像素进行聚类；

s22、进行哈尔变换；

s23、进行双硬阈值系数收缩；

s24、进行哈尔逆变换与像素聚合形成初步降噪图像；

所述步骤s3的具体执行方法为：

s31、将原始含噪图像和初步降噪图像进行块匹配和行匹配，然后对相似像素进行聚类，得到两个群组结果；

s32、将步骤s31中的两个群组结果分别进行哈尔变换；

s33、进行维纳滤波去噪；

s34、进行哈尔逆变换与像素聚合形成最终降噪图像。

本发明首先执行图像块匹配获得一定数目的图像块，将获得的图像按行或列扫描拼接成一个二维矩阵，再执行此二维矩阵上的行匹配或列匹配获得最相似的一定数目像素行。通过这种方法最大程度地将图像中的自相似信息进行了聚类。对聚类的像素执行简单的哈尔变换，由于哈尔变换存在简单的提升方法，即只是对这些行之间进行加权平均，所以很容易硬件实现。整个算法不需要任何矩阵乘法，与bm3d方法相比大大降低了运算复杂度。

在去噪过程中本发明提出一种双硬阈值策略，即不仅对变换系数进行一般的硬阈值实现系数收缩，而且将最相似的两行像素的哈尔变换的系数全部置为零(本发明称这种阈值为结构硬阈值)，从而大大提高了图像去噪能力。

实施例二

参阅图2所示，本实施例的图像去噪方法整体分为两个阶段，其中第一阶段为初步去噪阶段，第二阶段为精细去噪阶段。

第一阶段：初步去噪

第一步：图像块匹配与行匹配，对应实施例一中步骤s21，

首先执行原始含噪图像的块匹配,按一个指定的步长n_step提取大小为n1×n1的图像块作为参考块，然后在以该参考块为中心的一个大小为ns×ns的邻域内进行块匹配获得数量为n2的相似图像块，将匹配的所有图像块进行列扫描拼接成一个大小为(n1×n1)×n2的矩阵m，在m上用每一行作为参考行与其余的所有行计算欧式距离d，将距离最小n3行群组获得最大程度的图像中相似像素的聚类。

本发明中，示例性的设置参数值为：n1＝9，n2＝16，ns＝43，nstep＝8，n3＝4。也可以根据需要设置其他参数值。

第二步：可分哈尔变换，对应实施例一中步骤s22，

将第一步获得的行群组执行可分的二维哈尔变换，即纵向与横向各自执行哈尔变换，由于哈尔变换矩阵的独特的规则结构，在执行变换过程中本发明应用了提升变换方法，即只对行或列之间执行加权的平均，在整个变换过程中不用矩阵乘法，从而大大降低了运算复杂度。

第三步：双硬阈值系数收缩，对应实施例一中步骤s23，

由于最相似的两行的哈尔变换的结果相当于对两行的加权求差，所以本发明发现这些求差后的结果系数基本全部来自于噪声，所以对这样的系数采用了结构硬阈值，即将这些系数全部设置为零；对其余哈尔变换的系数采用传统的系数硬阈值进行系数收缩，系数硬阈值参数为thr。由于采用了双硬阈值策略，本发明的方法大大提高了比传统方法的去噪能力。

第四步：哈尔逆变换与像素聚合，对应实施例一中步骤s24,

对双硬阈值收缩后的哈尔变换系数执行逆哈尔变换，获得去噪后的二维矩阵，再将此二维矩阵扫描回原来的图像块，最终将这些图像块加权放回到图像中的原始位置获得去噪后的图像。

由于每次参与去噪的最相似像素数目是有限的，为了最好保留图像细节，硬阈值参数又不能选择的太大，所以仅执行一次上述步骤很难将图像噪声全部去除，所以本发明采用了迭代策略，用一个恰好能最大程度保留图像细节的硬阈值参数对上述四步操作实行多次迭代以获得最好的图像去噪结果。对于弱噪声情形只迭代两次，对强噪情形迭代的次数随着噪声强度加强而增加。

优选的，

当输入图片噪声标准偏差小于等于25时，thr＝1.65，迭代执行步骤s21-s24，2-3次；

当输入图片噪声标准偏差大于25且小于等于75时，thr＝1.45，迭代执行步骤s21-s24，4-6次；

当输入图片噪声标准偏差大于75时，thr＝1.45，迭代执行步骤s21-s24，7-12次。

第二阶段：精细去噪

第一步：图像块匹配与行匹配，对应实施例一中步骤s31，

由于第一阶段的图像的噪声已经几乎全部去除，所以能获得更好的图像块匹配结果。用第一阶段的结果图像作为参考同步执行第一阶段记过图像与原始含噪图像的块匹配操作及行匹配操作。按一个指定的步长n_step在第一阶段结果图像中提取大小为n1×n1的图像块作为参考块，然后在以该参考块为中心的一个大小为ns×ns的邻域内进行块匹配获得数量为n2的相似图像块，将匹配的所有图像块进行列扫描拼接成一个大小为(n1×n1)×n2的矩阵m，在m上用每一行作为参考行与其余的所有行计算欧式距离d，将距离最小n3行群组获得最大程度的图像中相似像素的聚类。原始含噪图像进行同步的同参数的操作，经过块匹配与行匹配后获得两个同样大小的行群组结果。

在本实施例中，根据不同参数强度对所用的参数值分为3组：

低噪：n1＝8，n2＝64，ns＝129，nstep＝5，n3＝8；

中噪：n1＝16，n2＝64，ns＝129，nstep＝10，n3＝8；

强噪：n1＝20，n2＝64，ns＝129，nstep＝13，n3＝8。

也可以根据需要设置其他参数值。

第二步：可分哈尔变换，对应实施例一中步骤s32，

进一步地，所述步骤s32中，将第一步获得的两个行群组分别执行同样的可分二维哈尔变换。

第三步：经验维纳滤波，对应实施例一中步骤s33，

所述步骤s33中，用步骤s2中初步降噪的结果的变换系数θ1对原始含噪图像的变换系数θ2按下述公式执行维纳滤波：

其中σ为噪声标准偏差。

第四步：哈尔逆变换与像素聚合，对应实施例一中步骤s34，

对维纳滤波后的含噪图像的系数执行逆哈尔变换，获得去噪后的二维矩阵，再将此二维矩阵扫描回原来的图像块，最终将这些图像块加权放回到图像中的原始位置获得最终去噪后的图像。

为了最大程度地提高去噪性能，本发明在第二阶段同样采用了迭代策略，无论噪声强弱，本阶段都采用两步迭代。

上述第一阶段是为了对图像初步去噪，第二阶段是为了将第一阶段去噪的过程中损失的真实信号能量进行恢复以及进一步取出第一阶段残留的噪声。

本发明通过简单哈尔变换对最相似的像素群组实现了有效的图像去噪。通过本发明的方法不仅能较好保留图像细节而且极少引入假信号，获得的去噪结果图像不仅在客观评价指标方面比现有算法有了较大提升，主观视觉质量也明显优于现有方法。

在实施例二中例举了采用行匹配的方法，图像还可以按列扫描的替代方法为按行扫描，两种扫描方法对本发明的结果是同样的。

参阅图3所示，其是原始图像与含噪图像；

本发明对图像去噪研究中通用的多幅标准图像用matlab软件向图像进行加噪与去噪实验，每次向图像中添加不同强度的噪声，形成图3所示的原始图像和含噪图像，然后对加噪的图像按本发明方法去噪。图4为房屋图像在不同噪声强度情形下采用bm3d方法去噪结果，图5为房屋图像在不同噪声强度情形下采用本发明方法去燥结果，从图中可以看出，本发明方法的去噪结果不仅更好保留了图像细节，而且几乎不引入假信号。本发明用房屋图像给出实验结果的原因是此图像含有很小的图像细节，更能体现本发明的去噪保留细节信息的优势。从图中数据可以看出，本发明的去噪结果比bm3d方法有大幅度提高，psnr值最多的比bm3d方法提高1db以上。

虽然结合附图描述了本发明的实施方式，但是专利所有者可以在所附权利要求的范围之内做出各种变形或修改，只要不超过本发明的权利要求所描述的保护范围，都应当在本发明的保护范围之内。