去除探地雷达信号中不相关噪声的方法.pdfVIP

第24卷第4期 华东交通大学学报 Vol24 No4 2加7年8月 J0urnal“Eds【China JiaoIollg UniversiIy A“g．，2007 文章编号：J005一0523(2007)04一0056—03 去除探地雷达信号中不相关噪声的方法 廖立坚，杨新安，黄小平 (同济大学城市轨道与铁道工程系，上海200092) 摘要：根据有艘信号在相邻道之间在波彤和能量上有相关性这一特点，作者改进了常规KL变换只能提取水平相关信号的不 足，使用小波域上的K工．变换对探地雷达剖面进行滤波，不仅使倾斜轴和弯由轴的到了加强，而且在低信噪比时，有效地区分 出有效信号和干扰信号，提高了雷达剖面的分辨率．本文运用铁路沪宁线路基检测的数据，使用该方法滤波．其处理后的剖面 更能有效地反映出路基病害的信息． 关键词：KL变换；滤波；探地雷达；小波 中图分类号：Tu45 文献标识码：A 1引言 方差阵为G∞(￡7膏)，有： 探地雷达勘探中，用数字仪记录雷达波时，为了 =工氇(互一E(_]|『))(互一E(x))屯 保持更多的波的特征，通常采用宽频带进行记录，因 =￡钆 (1) 此，在宽频带范围内记录了各种反射渡的同时．也记 实际计算协方差阵时一可用嬲7代替．当L为一的 录了各种干扰波(噪声)．根据相邻道反射波的有效 特征向量矩阵Ⅳ时(1)式取最大值． 成分(去除噪声的有用信号)在波形及能量上具有较 口=U啊 (2) 强的相关性，使用KL变换可以保留某一方向(如水 (2)式是x一口的变换式，这就是通常所说的KL变 平方向)的相关信号，从而压制不相关噪声和其他 换．重建信号能量的大小由特征值的大小决定，用一 渡_l J．但是在信噪比较低时，完全从时域的角度来考 的较大的几个特征值对应的特征向量组成的矩阵 虑信号的能量相关性就存在一定的缺陷，因此，作者 D7重构原始信号宣，写成矩阵形式为： 把具有很强时频分析能力的小波变换与KL变换结 童=D审 (3) 合起来，对探地雷达信号去噪．这种方法既能很好地 2．2倾斜KL变换的步骤 保持波的动力学特征，又能适用于低信噪比的雷达 使用常规KL变换去噪声时要求剖面中同相轴 资料去噪． 接近水平，否则不仅去噪效果较差，同相轴也会遭到 破坏．这里我们将常规的KL变换改进为倾斜KL变 2去噪原理及实现方法 换，它是设置一定大小的道间窗口，在窗口中将倾斜 轴和弯曲轴逼近水平轴再进行KL变换，其具体步 2．1常规KL变换的原理 骤如下12J： 设给定的探地雷达数据x为m×n的矩阵，第 对雷达数据剖面划分道时窗(如取道宽为21 l行表示第道数据，第，列表示第J个时刻的采样 道)，逐道滑动道时窗．对每个道取出前10道数据和 值．为了最大限度地解释变量变化的原因，在数学上 后10道数据构成时窗(前面不足10道从后面取，后 就是求￡，使得L啊的协方差阵最大，记L啊的协 面不足10道从前面取，取够2l道数据为止)做如下 收稿日期：2007—01—15 作者简介：廖立坚(J982一)，男，朔南郴州人，在读硕士，研究方向