一种基于改进图像块匹配的非局部均值图像去噪方法
【专利摘要】本发明的一种基于改进图像块匹配的非局部均值图像去噪方法，属于图像处理技术领域。有变差系数值图的建立、变差系数值图的分类、平坦区域上去噪处理、结构区域上去噪处理四个步骤。本发明具有运算时间短、计算权重准确、去噪效果好等优点。
【专利说明】
-种基于改进图像块匹配的非局部均值图像去噪方法
技术领域
[0001] 本发明属于图像处理技术领域，具体设及一种改进图像块匹配的非局部均值图像 去噪方法。
【背景技术】
[0002] 人类从外界所获取的大部分信息是通过视觉上对事物的反映，进而在大脑中形成 图像。由于图像具有直观、易懂、形象和信息量大等特点，因此应用领域非常广泛，并且发挥 着越来越重要的作用。因此，图像已经成为人类获取信息、感知世界进而改变世界的一种重 要手段。但是，由于成像系统、传输介质和记录仪器的不完善，在获取、传输和接收的过程 中，图像会被噪声污染，影响人类对图像观察和正确的理解，也会影响机器对图像的识别 等，所W对图像进行去噪处理就变得尤为重要。图像去噪的目的就是采用一定的技术手段， 对含有噪声图像利用相应的算法进行处理，同时也要求满足一定条件的噪声平滑程度，从 而能尽量把图像的细节等重要信息保存下来，近似恢复原图像，提高含噪声图像的清晰度 和质量。
[0003] 2005年，Buades等人首次提出了非局部均值(Non-local means，NLM)去噪方法，充 分利用了加性高斯白噪声模型和图像中像素存在自相似的特性，通过对各个像素配W相应 的权重并且相乘做和来近似估计目标像素灰度值。该方法W其特有的优势和理论基础受到 了学者的广泛关注，从而对NLM改进的一系列方法被提出来，不同程度的改善了图像去噪效 果。在经典NLM去噪方法中权重大小是由当前像素点所在的图像块与邻域内其他像素点所 在的图像块之间的相似性来度量的，图像块能把更加丰富的信息包含进来，运样就能描述 图像的特征，也能更好的估计当前目标像素的灰度值，从而达到有效的保护图像边缘、细节 等重要特征。但是在像素变化较大的区域上，两个图像块之间的相似性程度与图像块尺寸 有关系，进而影响了权重和最终的去噪效果。
【发明内容】
[0004] 本发明要解决的问题是:克服非局部均值在自适应匹配相似图像块尺寸方面的不 足，基于两个相似图像块之间的相关系数，提供一种在结构区域上能自适应匹配相似图像 块尺寸的非局部均值图像去噪方法。
[0005] 上述的技术问题通过W下的技术方案实现：
[0006] -种基于改进图像块匹配的非局部均值图像去噪方法，有变差系数值图的建立、 变差系数值图的分类、平坦区域上去噪处理、结构区域上去噪处理四个步骤；
[0007] 所述的变差系数值图的建立是，将噪声图像利用变差系数进行重新构图，得到变 差系数值图，所述的变差系数是图像上任一点的标准差相对于平均数大小的一个相对量， 由图像上所有点的变差系数值组成一个新的图像就得到变差系数值图；
[000引所述的变差系数值图的分类是，利用聚类算法对所述的变差系数值图进行分析， 将变差系数值图中的区域分成两大类，并记为平坦区域和结构区域；
[0009] 所述的平坦区域上去噪处理是，在变差系数值图的分类之后得到的平坦区域上， 用经典非局部均值图像去噪方法进行处理，得到平坦区域的去噪结果图；
[0010] 所述的结构区域上去噪处理是，在变差系数值图的分类之后得到的结构区域上， 先利用两个图像块之间相关系数的值来选取块尺寸，然后再利用经典非局部均值图像去噪 方法进行处理。
[0011] 本发明的一种基于改进图像块匹配的非局部均值图像去噪方法中，所述的变差系 数值图的建立，具体是，利用变差系数CV来衡量图像的匀质性，变差系数CV的计算公式为：
[0012]
[001引其中Fi是W像素 i为中屯、的正方形邻域，vs(j)是正方形邻域内任意一点像素 j的 灰度值，Ave化)为其正方形邻域内所有像素灰度值的均值，Var(Fi)是其方形邻域内所有像 素灰度值的方差;对整幅噪声图像进行计算后，建立变差系数值图。
[0014] 本发明的一种基于改进图像块匹配的非局部均值图像去噪方法中，所述的变差系
数值图的分类，具体是，根据公5 量数据与 聚类中屯、之间的距离，并根据公 中屯、值；其 中w(A,)表示同类中数据的个数;利用K-means聚类算法对CV值分成平坦区域和结构区域； 所述的K-means聚类算法，具体可参考文献：王千，王成，冯振元，叶金凤.K-means聚类算法 研究综述[J].电子设计工程，2012，07:21 -24。
[0015] 本发明的一种基于改进图像块匹配的非局部均值图像去噪方法中，所述的平坦区 域上去噪处理，具体是利用公式
[0016]
[0017]
[001 引
[0019] ^<'，7)=|，(^,)-1'的,'|对图像的平坦区域进行处理得到图像平坦区域的去噪结果 图；其中d(i，j)表示两个像素之间的欧式距离，w(i，j)表示的是像素点j所赋予的相似权 重，Z(i)表示欧氏距离的归一化常数，韦i表示的是目标像素点i的灰度估计值。
[0020] 本发明的一种基于改进图像块匹配的非局部均值图像去噪方法中，所述的结构区 域上去噪处理，具体是利用公^5
计 算图像块之间的相似系数Similarity,其中CR(a，b)表示两个图像块的相关系数;然后设定 阔值Tl = O. 25,T2 = 0.5，当SimiIarityCTl时,选择3 X 3的图像块；当Similarity >T2时， 选择7X7的图像块；当Tl《Similarity《T2时，选择5X5的图像块，由此来自适应匹配相似 图像块的尺寸。
[0021] 本发明的技术方案有W下有益效果：
[0022] 1、本发明在图像像素变化较大的区域上，能快速匹配相似图像块的尺寸，W减少 方法的运算时间。
[0023] 2、本发明可W实现相似图像块尺寸的自适应选择，使得在计算权重时能更加准 确，去噪效果更好。
【附图说明】
[0024] 图1是实施例2中所述的噪声图像。
[0025] 图2是实施例2中所述的CV值图。
[0026] 图3是实施例3中所述的CV值分类图。
[0027] 图4是实施例5中所述的最终去噪后的整幅图像。
[0028] 图5是本发明方法的流程图。
【具体实施方式】
[0029] 实施例1本发明的总体步骤
[0030] 本发明是一种改进图像块匹配的非局部均值图像去噪方法，有变差系数值图的建 立、变差系数值图的分类、平坦区域去噪处理、结构区域去噪处理四个步骤，首先利用变差 系数公式将要处理的带噪声的图片转换成变差系数值图，然后将变差系数值图划分成平坦 区域和结构区域，在平坦区域上，用经典非局部均值图像去噪方法进行处理，在结构区域 上，先利用两个图像块之间相关系数的值来选取块尺寸，然后再利用经典非局部均值图像 去噪方法进行处理，最终得到去噪的图像。
[0031] 实施例2变差系数值图的建立
[0032] 在本发明中，利用变差系数(Coefficient of Variation,CV)来衡量图像的匀质 性，所谓的变差系数是标准差相对于平均数大小的相对量，其公式可W表示为：
[0033]
(1)
[0034] 其中Fi是W像素 i为中屯、的正方形邻域，vs(j)是正方形邻域内任意一点像素 j的 灰度值，Ave化)为其正方形邻域内所有像素灰度值的均值，Var(Fi)是其方形邻域内所有像 素灰度值的方差。将图1所示的噪声图利用公式(1)处理后得到如图2所示的变差系数值图。
[0035] 实施例3变差系数值图的分类
[0036] 在实施例2中的变差系数值图中，CV值会有大小之分，其中较小的CV值为平坦区 域，较大的CV值为结构区域。利用K-means聚类算法对CV值分成两类并显示变差系数值图的 分类结果图。下面就简单的介绍K-means算法，该算法接受参数k，即k表示的是要进行分类 的种数。为了方便描述算法，定义W表示为样本包含n个数据，即
[0037]
(2)
[0038] 其中，Xa = Xal, Xa2...Xad是一个d维向量，表示的是第a个数据的d个属性，n是样本 容量。然后还要对将数据样本W划分为k类，则聚类中屯、可W表示为：
[0039] C=kb|cb=kbi,cb2. . .cbd),b=l,2. . .k} (3)
[0040] 其中Cb=(Cbl，Cb2. . .Cbd)表示的是第b个聚类中屯、，每个中屯、都有d种属性，k表示 的是k种聚类。要求满足同一个类的距离较小，若是不同的类则距离较大。其衡量距离采用 的公式如下：
[0041 ]
(4)
[0042] 然后把属于相同类的数据做相应均值处理，得到新的聚类中屯、的值，其计算公式 如下所示：
[0043]
日）
[0044] 其中W知)表示该类中数据的个数。
[0045] 该聚类方法基本思想是:先在图像空间中选择k个点做为质屯、然后进行聚类，先把 在位置上最靠近聚类中屯、的对象进行类别划分，归属到属于它们的那一类上。然后通过均 值迭加的方法，依次更新各不同聚类的中屯、值，直到得到最好的聚类结果。具体算法描述如 下：
[0046] (1)选择k个类的初始中屯、，通常情况下可W设置初始中屯、或者算法自动选择；
[0047] (2)在进行迭代处理中，对一个任意已知数据的样本，求各个数据到k个聚类中屯、 的距离，将该样本数据划分到距离最短的那一类；
[0048] (3)利用均值叠加方法更新该类的中屯、值；
[0049] (4)对于所有的k个聚类中屯、，当利用第(2)、（3)步的迭加更新后，其聚类中屯、值最 终保持不变或者满足一定要求条件情况下，则迭代结束，否则继续进行迭代。
[0050] CV值图经过K-means算法聚类之后，就把图像分成平坦区域和结构区域两类。如图 3所示，其中，图像中白色部分表示的是结构区域，黑色部分表示的是平坦区域，最终得到变 差系数值分类图。
[0051] 实施例4平坦区域上采用经典方法进行去噪处理
[0052] 在实施例3所得到的图像平坦区域上，采用经典非局部均值图像去噪方法(NLM)进 行去噪处理，NLM方法对图像去噪的公式如下所示：
[005；3]
(6)
[0化4]式中而]表示的是目标像素点i的灰度估计值，Q (i)表示的是W像素点i为中屯、， 尺寸为固定不变的捜索邻域，像素点j表示的是在捜索邻域Q (i)内所有像素中，任意选择 的一个像素，而v(j)表示像素点j的灰度值，w(i，j)表示的是像素点j所赋予的相似权重，其 权重公式表不为：
[0055]
(7)
[0056] Z(i)表示欧氏距离的归一化常数，其公式可W表示为：
[0057]
(8)
[005引 d(i，j)表示两个像素之间的欧式距离，用来衡量像素之间的相似性。h表示平滑参 数，对权重函数的衰减速度进行控制，去噪后图像平滑程度正是受到该值的控制。求欧氏距 离的公式如下：
[0化9]
(9)
[0060]其中Ni表示W目标像素 i为中屯、所在固定尺寸的相似图像块。V(Ni)表示目标像素 i 的灰度值，V(叫)表示的是其邻域中任意一像素点j的灰度值。另外0《w(i，j)《l，并且满足
-般情况下，在平坦区域相似图像块尺寸为7 X 7的正方形，平滑参数h = 5 X O。 然后利用公式(6)对整幅图像进行计算处理，就可W得到噪声图像在平坦区域的去噪结果 图。
[0061 ]实施例5结构区域上采用改进的方法进行去噪处理
[0062] 在实施例3所得的图像结构区域上，采用改进的非局部均值图像去噪方法进行去 噪处理。在图像结构区域上，相似信息比较复杂，所W相似图像块尺寸不同会造成它们之间 的相似性不同。目标像素所在的图像块为目标块，将捜索域内其他像素所在的图像块定义 为参与块。在变化较慢的结构区域上，大尺寸的目标快与参与块之间的相似度比较大;而在 变化较为激烈的结构域上，小尺寸的目标块与参与块之间有更高的相似度。相似图像块之 间的相似度随着两个图像块尺寸的变化而变化，由此说明含有最多相似信息的图像块尺寸 是不固定的。
[0063] 因此，本发明在非局部均值去噪方法的整体框架下，针对图像结构区域上的特点， 在计算相似图像块权重之前先自适应的匹配相似图像块的尺寸。其具体方法是利用两个图 像块的相似系数Similarity值的大小来衡量其之间相似性的大小，其公式如下：
[0064]
(10)
[0065] 其中CR(a，b)表示两个图像块的相关系数，可W很好地评价图像块之间的相似性， 其定义公式如下：
[0066]
(II
[0067] 其中，志表示目标块像素灰度值的平均值，^：表示参与块像素灰度值的平均值，并 且CR(a，b)的取值范围为[-1,1]。又根据公式（10)可得，Similarity的取值范围是[0，1]。当 Similarity的值越大时，说明目标块和参与块之间的相似度越低；当Similarity值越小时， 说明目标块和参与块之间的相似度越高，另外SimiIarity = O时说明目标块和参与块完全 相同。
[006引设置阔值Tl和T2，T1 < T2，并且在结构区域上选择3 X3、5X5、7X7的相似图像块 尺寸。有若Similarity<Tl，则说明相似度较高，选择3X3的图像块;若 说明相似度较低，选择7 X 7的图像块;若Tl <Simi Iarity <T2，则选择5 X 5的图像块。本发 明中，捜索域的尺寸为21 X21，T1 = 0.25，T2 = 0.5.当目标块和参与块选择相应的尺寸后， 再利用公式(6)对目标像素进行去噪，进而通过该方法对图像在结构区域的去噪结果图。最 终通过把图像在平坦区域和结构区域处理后结果图进行重新组合，得到最终去噪后的图 像，如图4所示。
[0069] 实施例2~5的方法流程图如图5所示，其步骤如下：
[0070] 步骤1:首先利用公式1计算整幅图像得变差系数值并得到CV值图，然后用Kmeans 算法根据CV值图将噪声图像分为平坦区域和结构区域，得到CV值分类图。
[0071] 步骤2:对噪声图像的平坦区域上的每点像素，采用固定尺寸的图像块，利用公式 (9)计算图像块间的相似性，然后用公式(8)计算权重，最终估计所有像素的灰度值得到平 坦区域上的去噪图像。
[0072] 步骤3:对噪声图像的结构区域上的每点像素，利用公式（10)计算相似系数值，根 据值的大小来自适应选取最相似的图像块尺寸，然后在经典非局部均值的框架下进行去噪 处理，最终得到整幅图像去噪后的结果图。
[0073] 实施例6本发明实施例2~5的相关程序
[0074]
【主权项】
1. 一种基于改进图像块匹配的非局部均值图像去噪方法，有变差系数值图的建立、变 差系数值图的分类、平坦区域上去噪处理、结构区域上去噪处理四个步骤； 所述的变差系数值图的建立是，将噪声图像利用变差系数进行重新构图，得到变差系 数值图，所述的变差系数是图像上任一点的标准差相对于平均数大小的一个相对量，由图 像上所有点的变差系数值组成一个新的图像就得到变差系数值图； 所述的变差系数值图的分类是，利用聚类算法对所述的变差系数值图进行分析，将变 差系数值图中的区域分成两大类，并记为平坦区域和结构区域。 所述的平坦区域上去噪处理是，在变差系数值图的分类之后得到的平坦区域上，用经 典非局部均值图像去噪方法进行处理，得到平坦区域的去噪结果图； 所述的结构区域上去噪处理是，在变差系数值图的分类之后得到的结构区域上，先利 用两个图像块之间相关系数的值来选取块尺寸，然后再利用经典非局部均值图像去噪方法 进行处理。2. 根据权利要求1所述的一种基于改进图像块匹配的非局部均值图像去噪方法，其特 征在于，所述的变差系数值图的建立，具体是，利用变差系数CV来衡量图像的匀质性，变差 系数CV的计算公式为：其中？1是以像素i为中心的正方形邻域，vs(j)是正方形邻域内任意一点像素j的灰度 值，Ave(F1)为其正方形邻域内所有像素灰度值的均值，Var(F1)是其方形邻域内所有像素灰 度值的方差;对整幅噪声图像进行计算后，建立变差系数值图。3. 根据权利要求1所述的一种基于改进图像块匹配的非局部均值图像去噪方法，其特 征在于，所述的变差系数值图的分类，具体是，根据公式2. . .n;b=l，2. . .k衡量数据与聚类中心之间的距离，并根据么2. · .d;b = l ,2. · .k更新聚类中心值;其中及(匕)表示同类中数据的个数;利用Κ-means聚类 算法对CV值分成平坦区域和结构区域。4. 根据权利要求1所述的一种基于改进图像块匹配的非局部均值图像去噪方法，其特 征在于，所述的平坦区域上去噪处理，具体是利用公式外,7)=|，(义)1('：|对图像的平坦区域进行处理得到图像平坦区域的去噪结果图；其 中d(i，j)表示两个像素之间的欧式距离，w(i，j)表示的是像素点j所赋予的相似权重，Z(i) 表示欧氏距离的归一化常数，项表示的是目标像素点i的灰度估计值。5.根据权利要求1~4任一所述的一种基于改进图像块匹配的非局部均倌图像去噪方 法，其特征在于，所述的结构区域上去噪处理，具体是利用公：和计算图像块之间的相似系数Simi Iarity，其中CR(a，b)表示 两个图像块的相关系数;然后设定阈值1'1 = 0.25，了2 = 0.5，当3^1&4七7<1'1时，选择3\3 的图像块；当Similarity>T2时，选择7\7的图像块；当1'1<3丨1^1&4七7<了2时，选择5\5 的图像块。
【文档编号】G06T5/00GK106023093SQ201610298113
【公开日】2016年10月12日
【申请日】2016年5月9日
【发明人】于银辉, 周恒 , 刘生成, 张羽丰, 陈兆强
【申请人】吉林大学