一种抗脉冲干扰的数学形态滤波方法.docxVIP

一种抗脉冲干扰的数学形态滤波方法 0 数学形态滤波 地磁观测容易受到外部环境的干扰。地磁站的日变观测数据通常会受到车辆、闪电和未知铁磁性材料的极大干扰。在某些情况下，噪声干扰可能大于真实的地磁信号本身。地磁观测仪器产生的仪器噪声也会干扰地磁数据。对地磁数据的噪声处理可以减少干扰对地磁质数据质量的影响，成为地磁数据预处理研究的重要内容之一。 数字滤波器作为最常用的信号处理手段可对地磁信号进行预处理,常用的数字滤波器有IIR滤波器、FIR滤波器等.但在数字滤波器设计中存在时滞、相移等缺陷,并且当信号与噪声差异较小且难以区分,或者干扰信号的频带分布较宽等情况下,数字滤波方法的应用受到限制.因此,我们尝试从干扰形态上寻找一种新的方法—数学形态滤波,来对受干扰地磁信号进行预处理. 数学形态学是从图像处理理论处理演变而来的新方法.作为一种重要的非线性滤波器,与传统的数字滤波器相比,数字形态滤波通过数学形态变换将复杂信号分解为具有物理意义的各个部分,在处理伴随有较强噪声干扰,甚至发生严重畸变的情况下,其基本形态依旧能别识别、重构及增强.Maragos等人在数学形态框架下系统总结形态滤波器集合理论,并阐明形态滤波器与线性平移不变、排序统计和层叠滤波器的关系,为数学形态学的发展奠定了坚实基础.目前数学形态学主要的应用领域包括图像处理、模式识别、电力系统监测、心(脑)电信号干扰抑制、机械振动信号特征提取等方面[9,10,11,12,13,14,15,16].但是,数学形态滤波应用在地磁信号预处理方面的工作开展较少. 地磁台站每天所产出的地磁数据包含平静变化信号和叠加在平静变化上的不规则扰动变化信号.平静期变化的频率周期一般以24、12、8、6小时的周期为主;而不规则扰动主要有磁扰和脉动两大类,其频率分布较广,成分较为复杂,频率分布从10-3～102Hz的数量级. 通过仪器采集的地磁数据中除了地磁信号本身以外,还会记录下外界环境干扰所产生的电磁场,以及仪器本身噪声.外界环境干扰主要包括车辆等铁磁性物质在距离仪器探头较近处运动产生的脉冲状干扰;雷暴雨天气所引起的连续的脉冲状干扰;部分大城市轨道交通距离地磁观测台站较近时,由于轨道交通漏电电流所产生的轨道交通干扰;以其一些未知的人工干扰等.其中由车辆等铁磁性物质运动所产生的脉冲状干扰最为常见,其干扰幅度一般较大,干扰持续时间从几秒到几十秒不等,而且干扰发生较为随机.由于是脉冲状干扰,其频谱分布较宽,采用通常的数字滤波器进行抑制效果并不理想. 本文根据数学形态滤波原理,提出一种地磁日变数据干扰抑制的新方法.讨论了不同结构元素对仪器干扰和暂态脉冲干扰的抑制效果.通过数值模拟,对地磁日变数据中添加随机噪声和不同形态的脉冲干扰噪声进行模拟抑制.最后对实际地磁观测仪器的仪器噪声和由车辆开行所产生的脉冲干扰进行处理,获得较好的应用效果.鉴于数学形态学用于地磁数据处理是一个新的应用领域,这方面的研究工作更需要进一步深入. 1 数学滤波的基本原理 1.1 腐蚀运算定义 数学形态学是一种以集合理论与逻辑运算为基础的非线性图像处理理论.其基本思想是用具有一定形态的结构元素作为“探针”元素对目标信号进行变换和匹配,保持信号基本形态,抑制不相干信号,以达到提取信号、保持细节和抑制噪声的目的. 假设地磁信号f(n)与结构元素g(n)分别定义在实数上的一维离散函数.f(n)的定义域D[f]∈{l,2,3,…,N},g(n)的定义域D[g]∈{0,1,2,3,…,P},且N?P.g(n)关于原点对称函数gs(n)=g(-n): 则腐蚀运算定义: 膨胀运算定义 式中,符号inf和sup分别表示函数的极小和极大值运算. 开运算定义: 闭运算定义 这里符号?和·分别表示开运算和闭运算. 开运算采用先膨胀后腐蚀运算,其作用是分离信号中比结构元素小的孤立部分,用于抑制信号中的正脉冲噪声;闭合运算采用先腐蚀后膨胀运算,其作用是补缺与内部连通,用于抑制信号的负脉冲噪声. 1.2 结构元素的选取 数学形态学方法是基于结构元素填充探测的思想,结构元素的选取是影响滤波器性能的关键因素,常见的结构元素为扁平型、矩形、三角形、半圆形、正(余)弦型等.大量仿真经验表明,较长的结构元素能够获得较好的低通效果,滤波的效果也较为明显,但结构元素的结构越长,复杂性越高,计算量也会相应增大. 地磁信号易受电磁瞬态脉冲干扰,因此在结构元素的选取上,结构元素应大于瞬态脉冲长度.而对于仪器背景噪声等随机噪声干扰,结构元素长度为地磁信号采样周期的5至10倍左右可获得较好的效果.在下一章节中将会具体讨论不同结构元素对于这些干扰的抑制效果.. 1.3 混合滤波器与交替混合滤波器 级联开、闭运算可构建交替滤波器foc(n)和fco(n). 同时,还能通过开、闭基本运算构建混合滤波