浅谈 中值滤波.docx

《浅谈 中值滤波.docx》由会员分享，可在线阅读，更多相关《浅谈 中值滤波.docx（14页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、浅谈 中值滤波 浅谈中值滤波 1.中值滤波的现状 在数字信号处理和数字图像处理的早期探讨中,线性滤波是主要的处理手段。线性滤波简洁的数学表达式以及某些志向特性使其很简单设计和实现。然而，当信号中含有非叠加性噪声时,例如非线性引起的噪声或非高斯噪声等,线性滤波的处理效果就很难令人满足。在处理图像时,线性滤波将破坏边缘,而且不能有效滤除脉冲噪声。为了克服线性滤波方法的局限性,探讨非线性滤波的方法为数字信号处理重要课题之一。 非线性滤波基于对输入信号序列的一种非线性映射关系,常可把某一特定的噪声近似映射为零而保留信号的重要特征,因而可以在肯定程度上克服线性滤波的不足。 1971年闻名学者J.W.Tu

2、key在他的开拓性论文中提出了中值滤波的概念并用作时间序列平滑。中值滤波一出现就因其具有对尖脉冲的良好抑制实力，在平滑加性噪声时能保持信号的边缘特征等优点而备受瞩目。 常用的中值滤波是非线性滤波的代表，由于经典的中值滤波算法在滤除噪声的同时会使信号中重要的细微环节信息受损,因此,很多改进的中值滤波算法相继被提出。 2.中值滤波 在数据处理中我们常常运用的是滑动中值滤波，即取定中值滤波的跨度N(一般 N 为奇数)，在数据序列中顺次取得 N 个数据，然后将该数据列的中值作为中心位置的值输出以形成新的数据序列，在滤波中应将原数据序列的两个边界各补充(N-1)/2(N为奇数时)个等于边界的点以使滤波后

3、的新数据序列长度与原始的数据序列长度一样。 2.1一般中值滤波 2.11一般中值滤波的基本原理 设有一个序列：x1,x2,x3,x4,x5 ，将它们根据肯定值大小重新排列此序列 x3, x5,x2, x4,x1 重排以后的中值是x2，此值就作为滤波的输出。明显，x2不能表示成输入数据和滤波系数的褶积的线性组合。其主要特点有： (1)一般中值滤波肯定阻挡噪声峰值，因为中值滤波只取中位数，肯定不会取异样数。例如有一组数 (x1,x2,x3,x4,x5) 正常数 axna， n=2,3,4,5 异样数 x1 a a 表示一个数，将以上数组自小到大排列后为 (x3, x5,x2, x4,x1 )取中位

4、数x2，决不会取异样数x1。 (2)一般中值滤波是低通滤波器，中值滤波取中值为序列的输出，可以看作是对数据序列进行局部平滑，这种局部平滑实质就是低通滤波。 (3)一般中值滤波不变更阶越函数在空间、时间上的位置，这一性质对于信号处理中的爱护边缘有着重要的作用。 (4)当中值滤波的滤波窗口足够长时，有限宽度的三角波和矩形波可以被完全平滑。 (5)中值滤波由于没有统计效应，对随机出现的小的振幅值有时不能完全平滑，所以通常信号在中值滤波处理以后须要再进行带通滤波。 2.1.2一般中值滤波(MF)的数学基础 中值滤波对数字序列有平滑作用，平滑也就是数据靠近，这样则存在误差，如何利用误差最小来确定平滑参数

5、，一般常见的有两种准则：(1)使误差的平方和达到最小；(2)使误差的肯定值和达到最小。平均值平滑的数学原理应用准则(1)，即符合误差的平方和最小。中值滤波则是利用准则(2)来实现对数据序列的平滑。 设x是 n 个数据序列的中位数，xi 表示一组序列。x与xi之差的肯定值和为： Q=x-xi （3.1.1） i=1n要使Q最小，则 Q=0 （3.1.2） xnnx-xiQn2即 =(x-xi)=Sign(x-xi)=0 xxi=1i=1x-xii=1式中：Sign符号函数。 当xi x 时，Sign 为负； 当xi 这样在选择x时，使得在 n 个数中，有 n/2 个xi大于x,同样有 n/2 个

6、xi小于x，中间的xi即为x；假如 n 为偶数，则取中间的两个xi的平均值为x。 2.2 加权中值滤波 2.2.1 加权中值滤波(WM)的基本原理 由上可以看出通过变更加权系数，完全可以变更中值滤波的性质，来达到我们的要求。 2.3 一维中值滤波对信号作用的结果分析 由于中值滤波是一种特别的非线性滤波手段，它对脉冲的响应为零，（在一个输入上施加一个脉冲函数引起的时间响应。）所以在傅氏域没有“真正的”振幅谱和相位谱。我们只能通过它对已知信号及其频谱特征的响应来分析其各种滤波特性。虽然中值滤波的理论比较完善，但是由于多数状况处理的信号是对称信号，所以并没有人留意到中值滤波对信号相位的影响。 2.3

7、.1一般中值滤波对对称信号相位的影响 （1）在频谱图中，一般中值滤波引入了假高频成分，并且在子波的频带范围内，滤波后子波的主要频带向低频方向移动，此特点在数据处理时应当着重留意，要依据数据处理时的详细要求来推断，同时也成为选择滤波长度的一个条件。 （2）经过一般中值滤波后对称信号的相位不发生移动，这使得我们在处理由对称信号(例如雷克子波、奥姆斯比子波等)作为子波的合成地震记录时，不须要考虑相移问题。但由此就得出结论说中值滤波处理后的全部类型的信号的相位都不发生移动则是片面甚至错误的。 2.3.2 一般中值滤波对非对称信号相位的影响 一般中值滤波对非对称信号的处理效果不同于处理对称信号，假如 用

8、处理对称信号的规律来对待非对称信号则往往不能达到预期的效果。对应滤波后的频谱同样向低频方向移动，但假高频现象却并不如对称信号滤波后明显。处理非对称信号的同时必需留意选择的滤波点数是否使相位的变更在要求的范围内。从滤波后和滤波前的最大振幅平方比来看滤波前后的能量改变，发觉在同等状况下，一般中值滤波对非对称信号的衰减实力大于对对称信号的衰减。 2.3.3 加权中值滤波对对称信号相位的影响 （1）在频谱图中，加权中值滤波也引入了假高频成分；并且滤波后的子波的主要频带向低频方向移动，说明白加权中值滤波的低通滤波特性。 (2) 同样加权中值滤波对对称信号的相位不产生影响。 2.3.4 加权中值滤波对非对

9、称信号相位的影响 当滤波长度大于肯定的子波宽度时，波形已经失去了原有的形态，但是在波形失去原有形态之前，经过加权中值滤波后的子波表现出较好的辨别率特性；在频谱上加权中值滤波仍旧表现出低通特性；信号的相位也因滤波而产生了畸变。 2.3.5 一般、加权中值滤波对不同信号作用的比较 一般中值滤波和加权中值滤波对于同一种信号表现出相像的特性：二者在处理对称信号时，都起到了衰减的作用，并且对信号的相位都不产生影响，同时使信号的频谱中掺入了假高频成分，还表现出了中值滤波的低通特性；在处理非对称信号时，除了对信号产生衰减作用外，还使信号的相位发生了畸变。 尽管一般和加权中值滤波有相像之处，但是它们还是存在着

10、较大的差异： 在处理对称信号时，一般和加权中值滤波分别对同一信号进行滤波以后，信号的峰值很接近，但是加权中值滤波比一般中值滤波更有利于提高信号的辨别率，在频谱上加权中值滤波比一般中值滤波表现出更加严峻的假高频现象。 在处理非对称信号时，加权中值滤波比一般中值滤波表现出的更好的提高辨别率的性质。而在相位谱分析中，尽管加权和一般中值滤波都使信号的相位发生变更，但是在滤波长度较小的状况下，经过一般中值滤波得到的信号的相位曲线虽然已经发生变更，但是仍旧与原始信号的相位曲线有相同的趋势，并没有偏离太多；而此时即使在滤波长度较小的状况下，经过加权中值滤波得到的信号的相位曲线已经变得不行分辨了。 经过以上的

11、探讨，我们可以相识到在实际地震资料的处理中，应用中值滤波除了应当考虑信噪比和辨别率以外，更加不容忽视的就是信号的相移问题，这对于资料的可信度起着至关重要的作用。 由于在实际中，常常用到的一维中值滤波是不加权的，所以经常把一维不加权中值滤波简说成一般中值滤波，但是随着对处理手段的进一步要求，加权中值滤波的地位日益突出，并且终归不加权的状况只是一种特别的加权中值滤波，所以一般中值滤波的概念也应当扩充为加权中值滤波。 通过对两种不同加权中值滤波（一般、加实数权）的探讨，总结出了一些关于一维中值滤波方面的阅历： (1) 通过不同权系数的选取，中值滤波表现出不同的特性，我们可以依据对实际状况的分析来选取

12、不同的权系数以适应各自的须要。 (2) 本次只是选取了两种特别的权系数来分析，而在实际中存在着更多的权系数的选取方法，但是不管权系数的形式如何，都可以仿照本文的方法加以探讨。 (3) 虽然中值滤波可以满意肯定的要求，但是我们同时也应当留意到它们存在的问题：中值滤波会引起信号形态上的畸变，而且畸变程度和滤波长度有关；中值滤波会引入假高频，因此信号在经过中值滤波后可以依据状况做一次低通滤波；中值滤波对非对称信号进行处理时，会引起相位畸变，因此在运用中值滤波之前应当试验相位畸变是否在处理的允许的范围内；虽然选取适当的权系数后，加权中值滤波可以使信号提高辨别率，但是同时带来“小台阶”效应，因此经过加权

13、中值滤波处理后的信号举荐做一次平滑处理。 由于中值滤波是一类特别的滤波方法，因此我们利用它进行信号处理时应当特别留意。为了得到预期的效果，处理之前做一下试验以确定最佳的滤波长度是特别必要的。 3中值滤波在地震资料处理中的应用 3.1中值滤波在井间地震资料处理中的应用 中值滤波是一种简便有效且信号失真较小的信号处理方法。在不同的道集域下，井间地震资料中的直达波、一次反射波和多次反射波在相邻道间的时差具有不同的表现形式，利用这一特点，应用中值滤波在不同道集域内对井间地震资料进行滤波处理，可以得到很好的效果。 对于井间地震资料，我们所须要的有效反射波是来自于激发点与接收点下方的一次反射波(上行反射)

14、和来自于激发点与接收点上方的一次反射波(下行反射)，其它波均视为相关干扰或无效信息。中值滤波是以正常时差不同为基础的多道滤波技术，在井间地震特别的观测系统中，中值滤波可以发挥其自身的优点。 通过对井间地震不同道集域下道间时差的分析可以知道，仅运用中值滤波即可达到较好的波场分别效果。为了验证不同道集域下中值滤波对数据处理的效果，进行直达波与多次波的衰减、一次反射波的增加以及上下行波场的分别。 3.1.1直达波和多次波的衰减 首先对数据进行带通滤波，消退有效频带之外的噪音干扰，将共炮点道集重排为共偏移距道集。在共偏移距道集下，依据(1)式和(3)式可知，直达波和多次波除了受速度影响外，其相邻道间时

15、差为0，通过共偏移距道集对初至时间拉平排齐，在肯定程度上消退了速度的影响，然后选择适当的时窗参数，采纳中值滤波消退相邻道时差为0的波组记录，使直达波和多次反射波得到衰减。 3.1.2反射波的增加 把衰减了直达波和多次波的数据体重新抽道组成共中心点道集，对于共中心点道集，由于s=-g，依据(2)式可知，一次反射波在不受速度影响的状况下其相邻道间的时差为0，通过共中心点道集对反射波时间拉平排齐，消退速度的影响，再次做中值滤波处理，本次中值滤波是为了保留相邻道时差为0的波组记录，而相邻道时差不为0的波组记录将被减弱，因而一次反射波同相轴得以增加，而其它 波场(如直达波和多次波)再次得到衰减(图4，虚

16、线圈)。 3.1.3上、下行波场的分别 首先对上行反射波进行拉平(图5中的)，然后通过中值滤波使上行反射增加而下行反射减弱(图5中的)，最终返回原始时间剖面得到上行反射波场(图5中的)；反之，得到下行反射波场(图5中的)。图6为通过中值滤波最终获得的上行和下行反射波场。 在波场分别中，对于资料相对较好的地震数据，仅应用中值滤波即可达到较好的波场分别效果；对于信噪比较低的资料，可以用中值滤波技术使资料的有效波场加强，并得到上、下行反射信息，然后再利用中值滤波进行波场分别。 井间观测系统所具有的特别性，使得同一种地震波在不同道集域下的时差表现形式不同，因此可以在井间地震处理中利用中值滤波技术实现直

17、达波和多次波的衰减，反射波的增加以及上、下行反射波的分别。同时，由于中值滤波处理对地震原始资料畸变程度较小，因此可以提高井间地震资料反射成像的质量。 3.2多道中值滤波在分别VSP波场中的应用？ 多道中值滤波处理流程图 4.中值滤波特点 中值滤波是一个非线性过程，最大优点是算法简洁且去噪效果明显。中值滤波具有 如下特点: (1)中值滤波肯定阻挡噪声峰值。 (2)中值滤波不变更阶跃函数在空间、时间上的位置。 （3）消退尖峰波以及增加部分有效波； （4）对野外原始地震资料信息的畸变和负面影响较小等优点。 (5)中值滤波平滑三角波,其平滑作用随着中值滤波长度N的增加而增加,当其达到肯定长度时,可将三

18、角波平滑为具有相等幅度的志向的直流重量。 (4)中值滤波平滑矩形波,若中值滤波足够长时,矩形波被完全平滑。 (6)中值滤波由于沒有像计算均值那样的统计效应,对随机出现的小的振幅值有时不能完全平滑。中值滤波事实上是一个平滑滤波,经过其处理之后,主频向低频移动,高频成分受到损害,正是由于其平滑作用,使处理后的地震数据波形过于一样而显呆板,有一些信息不行能客观反映出来。 5.总结 中值滤波器是一种特别的非线性滤波器，与线性滤波器不同。线性滤波器的振幅、相位谱完全确定了滤波器在频率域和时间域的特征。与之相比，中值滤波器对脉冲的响应是零，在频率域没有“真正”的振幅谱和相位谱。所以，尽管中值滤波器原理很简

19、洁，但了解其特性比了解线性滤波器特性要困难得多。 虽然中值滤波技术在信号处理领域得到很大的重视，特殊在非平稳信号的处理中取得了较大的胜利，然而中值滤波的一个严峻不足是引起相对滤波窗口而言较为“细小”的信号细微环节结构的破坏和丢失，在图象处理中，中值滤波的这一缺憾要比在一维信号的处理中更加显著。缘由主要来自两个方面： 第一，二维信号几乎没有根信号，也就是说几乎全部的二维信号经中值滤波以后都要受到不同程度的破坏； 其次，图象中的某些诸如细线，拐角等细微环节结构往往包含重要的信息，这些结构的破坏或丢失往往比噪声本身更为不行接受。 所以，爱护细微环节的中值类滤波器的探讨成为非线性滤波器探讨的一个重要方

20、面。多级中值滤波则正是人们在努力寻求的兼有细微环节爱护和噪声抑制的优良特性的滤波器结构。探讨表明，长度较小的窗口能够较好地爱护信号的细微环节信息，但却不能有效地滤除随机噪声；而长度较大的滑动窗口能更好地抑制噪声，同时却严峻地损失重要信息。依据噪声性质自动变更滤波窗口长度的滑动加权中值滤波器，更好地适应去噪的需求。 浅谈 中值滤波 数字图像处理_领域平均滤波_中值滤波 滤波技术 图像滤波总结 滤波电路教案 放大及滤波电路 滤波电路教案（举荐） 滤波和算子总结 Matlab频域滤波练习 逆滤波试验报告 本文来源：网络收集与整理，如有侵权，请联系作者删除，谢谢！第14页 共14页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页