滤波方法与流程

1.本发明涉及视频图像滤波技术领域，尤其涉及一种在已知噪声强度条件下的滤波方法。


背景技术：

2.现有的滤波技术，多数不考虑噪声的强度，只是通过设计合适的滤波器达到滤波的效果。但是，不同视频、图像中的噪声的强度不尽相同，通用的滤波方法缺乏针对性，对具体视频及图像的滤波效果还有提升的空间。
3.近年来，基于模板匹配的滤波方法取得了较好的效果，其核心思想是在像素周围取模板，模板包含若干模板像素，通过比较两个像素的模板来判断这两个像素是否属于同类像素，也就是真实值相同或非常相近的像素，同类像素的值之间的差别主要是由于噪声引起的，因此可以用待滤波像素的多个同类像素和待滤波像素自身的平均值来代替待滤波像素的值，达到消除噪声的目的。
4.基于模板匹配的滤波方法的一个关键步骤是判断同类像素。判断的准确率越高，滤波效果越好。现有的判断条件的设置多数是根据经验，造成滤波效果不够稳定，时好时坏。


技术实现要素：

5.本发明的目的在于提供一种滤波方法，能够提高同类像素判断的准确率，进而使滤波效果始终保持在较好水平。
6.为了达到上述目的，本发明提供了一种滤波方法，根据噪声的强度进行滤波，所述滤波方法包括以下步骤：
7.步骤s1：以待滤波像素为中心选取第一模板，所述第一模板中包含所有到所述待滤波像素的距离小于或等于第一半径的像素；
8.步骤s2：根据所述待滤波像素附近的噪声强度设定门限值；
9.步骤s3：以所述待滤波像素为中心选取搜索区，所述搜索区中包含所有到所述待滤波像素的距离小于等于第二半径的待判定像素，
10.步骤s4：对所述搜索区中的每一个待判定像素执行以下步骤：
11.步骤s41：以所述待判定像素为中心选取第二模板，所述第二模板中包含所有到所述待判定像素的距离小于等于第一半径的像素；
12.步骤s42：计算所述待判定像素与所述待滤波像素之间的模板距离；
13.步骤s43：若所述模板距离小于或等于所述门限值，则所述待判定像素为所述待滤波像素的同类像素，并将所述待判定像素加入同类像素集；
14.步骤s5：计算所述同类像素集中所有像素的像素值的平均值，并将所述平均值作为所述待滤波像素滤波后的值。
15.可选的，所述步骤s1、所述步骤s3及所述步骤s41中的两个像素之间的距离按照以
下公式进行计算：
16.d(p1,p2)＝|x1‑
x2|+|y1‑
y2|
17.其中，所述(x1，y1)和(x2，y2)分别为像素p1和像素p2的坐标。
18.可选的，所述步骤s2中，根据所述待滤波像素附近的噪声的标准差计算门限值。
19.可选的，所述门限值按照以下计算公式进行计算：
20.threshold＝σ
a
21.其中，threshold为门限值，σ为已知的标准差，a是一个正数。
22.可选的，所述a的取值介于0.5
‑
3之间。
23.可选的，所述第一模板中包含的像素数与所述第二模板中包含的像素数相同。
24.可选的，所述步骤s42中，所述待判定像素与所述待滤波像素之间的模板距离按照以下公式进行计算：
[0025][0026]
其中，m为所述待判定像素与所述待滤波像素之间的模板距离，n为所述第一模板中包含的像素数或所述第二模板中包含的像素数，p
i
和q
i
分别为把所述第一模板及所述第二模板中的像素按照相同的方法排序后的第i个像素的值，e
i
为所述第i个像素的权重。
[0027]
可选的，所述第i个像素的权重按照以下公式进行计算：
[0028][0029]
其中，d(p
i
,p)为第i个像素到所述待滤波像素的距离。
[0030]
可选的，所述第一半径的取值介于1
‑
4之间。
[0031]
可选的，所述第二半径的取值大于或等于4。
[0032]
本发明提供了一种滤波方法，通过根据噪声强度设定门限值，可以更加准确的将真实值相同的像素判定为同类像素，将真实值不同的像素判断为非同类像素，并取所有同类像素的平均值作为待滤波像素滤波后的像素值，能够达到更好的滤波效果。
附图说明
[0033]
本领域的普通技术人员应当理解，提供的附图用于更好地理解本发明，而不对本发明的范围构成任何限定。其中：
[0034]
图1为本发明实施例提供的滤波方法的流程图；
[0035]
图2为本发明实施例提供的第一模板的示意图；
[0036]
图3为本发明实施例提供的搜索区的示意图；
[0037]
其中，附图标记为：
[0038]
p0，p1，p2，p3，p4，p5，p6，p7，p8，p9，p
10
，p
11
，p
12
‑
第一模板中包含到待滤波像素的距离小于或等于第一半径的像素。
具体实施方式
[0039]
为使本发明的目的、优点和特征更加清楚，以下结合附图和具体实施例对本发明作进一步详细说明。需说明的是，附图均采用非常简化的形式且未按比例绘制，仅用以方便、明晰地辅助说明本发明实施例的目的。此外，附图所展示的结构往往是实际结构的一部分。特别的，各附图需要展示的侧重点不同，有时会采用不同的比例。还应当理解的是，除非特别说明或者指出，否则说明书中的术语“第一”、“第二”、“第三”等描述仅仅用于区分说明书中的各个组件、元素、步骤等，而不是用于表示各个组件、元素、步骤之间的逻辑关系或者顺序关系等。
[0040]
请参照图1，并结合图2
‑
图3，本实施例提供了一种滤波方法，根据噪声的强度进行滤波，所述滤波方法包括以下步骤：
[0041]
步骤s1：以待滤波像素为中心选取第一模板，所述第一模板中包含所有到所述待滤波像素的距离小于或等于第一半径的像素；
[0042]
步骤s2：根据所述待滤波像素附近的噪声强度设定门限值；
[0043]
步骤s3：以所述待滤波像素为中心选取搜索区，所述搜索区中包含所有到所述待滤波像素的距离小于等于第二半径的待判定像素，
[0044]
步骤s4：对所述搜索区中的每一个待判定像素执行以下步骤：
[0045]
步骤s41：以所述待判定像素为中心选取第二模板，所述第二模板中包含所有到所述待判定像素的距离小于等于第一半径的像素；
[0046]
步骤s42：计算所述待判定像素与所述待滤波像素之间的模板距离；
[0047]
步骤s43：若所述模板距离小于或等于所述门限值，则所述待判定像素为所述待滤波像素的同类像素，并将所述待判定像素加入同类像素集；
[0048]
步骤s5：计算所述同类像素集中所有像素的像素值的平均值，并将所述平均值作为所述待滤波像素滤波后的值。
[0049]
具体的，本实施例中提供了一种在已知噪声强度的条件下进行视频图像滤波的方法。所述视频图像可以看成一个二维矩阵，记为i[x][y]，x＝0，1，...，w
‑
1，y＝0，1，...，h
‑
1，其中w和h分别是视频图像的宽度和高度，以像素为单位。然后对所述视频图像中的每一个待滤波像素执行以下步骤。
[0050]
首先执行步骤s1，以待滤波像素p为中心选取第一模板ω(p)，所述第一模板ω(p)中包含所有到所述待滤波像素p的距离小于或等于第一半径r的像素，所述第一模板ω(p)中包含的像素数记为n。本实施例中，如图2所示，所述第一半径r＝2，所述第一模板ω(p)中包含的像素数n＝13，分别为p0，p1，p2，p3，p4，p5，p6，p7，p8，p9，p
10
，p
11
及p
12
，所述p0即所述待滤波像素p。
[0051]
本实施例中，所述第一模板ω(p)中的任一像素到所述待滤波像素p的距离的计算公式如下：
[0052]
d(p
i
,p)＝|x
i
‑
x0|+|y
i
‑
y0|
[0053]
其中，所述(x
i
，y
i
)和(x0，y0)分别为像素p
i
和待滤波像素p的坐标。
[0054]
然后执行步骤s2，根据所述待滤波像素p附近的噪声强度设定门限值。本实施例中，根据所述待滤波像素附近的噪声的标准差计算门限值，即所述噪声的标准差用于表征噪声强度。应当理解的是，所述待滤波像素附近的噪声的标准差为已知值，例如可通所述第
一模板内的像素的值计算获得，也可以通过整个视频图像的像素的值计算获得，本申请对此不作限制。
[0055]
所述门限值按照以下计算公式进行计算：
[0056]
threshold＝σ
a
[0057]
其中，threshold为门限值，σ为已知的标准差，a是一个正数。
[0058]
本实施例中，所述a的取值介于0.5
‑
3之间。
[0059]
然后执行步骤s3：以所述待滤波像素p为中心选取搜索区φ(p)，所述搜索区φ(p)中包含所有到所述待滤波像素p的距离小于等于第二半径r的待判定像素q，计算公式同上。本实施例中，如图3所示，所述第二半径r的取值大于或等于4。
[0060]
接着执行步骤s4，对所述搜索区φ(p)中的每一个待判定像素q执行以下步骤：
[0061]
步骤s41：以所述待判定像素q为中心选取第二模板ω(q)，所述第二模板ω(q)中包含所有到所述待判定像素q的距离小于等于第一半径r的像素，计算公式同上。本实施例中，所述第二模板ω(q)中包含的像素数n＝13，即所述第一模板ω(p)中包含的像素数与所述第二模板ω(q)中包含的像素数相同。
[0062]
步骤s42：计算所述待判定像素与所述待滤波像素之间的模板距离。本实施例中，所述待判定像素与所述待滤波像素之间的模板距离按照以下公式进行计算：
[0063][0064]
其中，m为所述待判定像素q与所述待滤波像素p之间的模板距离，n为所述第一模板ω(p)中包含的像素数或所述第二模板ω(q)中包含的像素数，p
i
和q
i
分别为把所述第一模板ω((p)及所述第二模板ω((q)中的像素按照相同的方法排序后的第i个像素的值，e
i
为所述第i个像素的权重。
[0065]
进一步的，所述第i个像素的权重按照以下公式进行计算：
[0066][0067]
其中，d(p
i
,p)为第i个像素到所述待滤波像素p的距离。
[0068]
本实施例中，e
i
＞0且
[0069]
步骤s43：若所述模板距离小于或等于所述门限值threshold，则所述待判定像素q为所述待滤波像素p的同类像素，并将所述待判定像素加入同类像素集。
[0070]
最后执行步骤s5，计算所述同类像素集中所有像素的像素值的平均值，并将所述平均值作为所述待滤波像素p滤波后的值。
[0071]
综上，本发明提供了一种滤波方法，通过根据噪声强度设定门限值，可以更加准确的将真实值相同的像素判定为同类像素，将真实值不同的像素判断为非同类像素，并取所有同类像素的平均值作为待滤波像素滤波后的像素值，能够达到更好的滤波效果。
[0072]
此外还应该认识到，虽然本发明已以较佳实施例披露如上，然而上述实施例并非用以限定本发明。对于任何熟悉本领域的技术人员而言，在不脱离本发明技术方案范围情
况下，都可利用上述揭示的技术内容对本发明技术方案作出许多可能的变动和修饰，或修改为等同变化的等效实施例。因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均仍属于本发明技术方案保护的范围。