专利名称：探地雷达数据处理方法
技术领域：
本发明涉及雷达探测及其应用技术领域，尤其涉及探地雷达数据处理方法。
背景技术：
二十世纪八十年代，探地雷达被引入到铁路路基检测中，迅速在世界各国开展了 相关的试验和研究，并取得了一定进展。虽然随着探地雷达硬件不断更新发展，其采样速 度、适时显示技术、探测深度和精度方面都有不同程度地提高。 然而几乎所有的探地雷达配套系统软件只能提供数据采集及简单的数据处理功 能，这一方面既限制了探地雷达的使用范围(只能局限于专业人员)，另一方面其最终成果 的表现也无法让现场的施工人员轻易理解。其成果解释始终处于传统的物探资料解释的束 缚之下，即先打印雷达时间剖面图一由专业人员对雷达剖面进行解释一从图上量测界面 分层和病害位置一绘制解释成果图。这对解释人员的地球物理专业知识及计算机水平提出 了更高的要求，在很大程度上限制了探地雷达在铁路路基检测领域内的应用与推广。
因此，有必要开发一种适用于铁路路基和轨道交通路基检测的探地雷达后处理方 法。
发明内容
本发明的目的就是为了克服上述现有技术存在的不足而提供一种具有多种功能 的探地雷达数据处理方法。 本发明的目的可以通过以下技术方案来实现 探地雷达数据处理方法，其特征在于，该方法RAILGPRl.O的界面开发环境选择的 是Visual 0++2005，计算核心语言采用的是區11^8 7. l，通过前台的VC调用MATLAB引擎， 使用MATLAB在后台计算；该方法包括
(1)程序操作一般步骤
①数据单项处理 探地雷达数据处理方法RAILGPR1. 0数据处理的主要作用是去除杂波和提高雷达 剖面分辨率；杂波可以分为两种，一种是由探地雷达硬件本身及其工作方式带来的，是固定 不变的，另一种是由检测时介质性质和周边环境造成的，是随外界条件不同而不同的；因 此，可以将探地雷达的数据处理归为三大类，一类是预处理，用来去除固有噪声；一类是噪 声和干扰的去除，根据外界杂波种类而定；一类是分辨率的提高，用来增强有效信号；
②自动化批处理 数据处理是一个相互穿叉的过程，对一组数据的处理顺序和往返次数应视具体情 况而定，对于大批量的数据，先进行样组处理过程试验，确定处理方法和流程后再做自动化 批处理； ③路基病害反演 判别路基结构是否发生病害，反演出路基层面图和路基病害位置和类型；
(2)程序设计 利用Visual 0++2005新建一个项目，项目名称为1^0^，项目类型为多文档；在整
个项目中，包含28个类，各类均有分工。 所述的28个类，各类的分工如下 1) CBXTFrame :"波形堆积图"的子框架； 2)CHDTFrame :"灰度堆积图"的子框架； 3) CCSTFrame :"彩色堆积图"的子框架； 4)CMainFrame :程序的主框架； 5) CJMSBFrame :"层面识别"的子框架； 6)CBHSBFrame :"病害识别"的子框架； 7) CAboutDlg :显示"关于"的对话框； 8) CBCSZDlg :显示"标尺设置"的对话框； 9) CCYDZLDlg :显示"采样点整理"的对话框； 10) CGTLBDlg :显示"高通滤波"的对话框； 11) CJMSZDlg :显示"层面设置"的对话框； 12) CKLBHDlg :显示"KL变换"的对话框； 13) CRMCHFDlg :显示"RMC雷达头文件"的对话框； 14)CSIRHFDlg :显示"SIR雷达头文件"的对话框； 15) CSJWLXLDlg :显示"神经网络训练"的对话框； 16) CXCGZYXDlg :显示"消除轨枕影响"的对话框； 17) CXCZOBDlg :显示"消除直耦波"的对话框； 18) CYWXBQZDlg :显示"小波去噪"的对话框； 19) CDataProDlg :显示"数据批处理"的对话框； 20)CBXTView :"波形堆积图"的视图； 21)CHDTView :"灰度堆积图"的视图； 22)CCSTView :"彩色堆积图"的视图； 23)CJMSBView :"层面识别"的视图； 24)CBHSBView :"病害识别"的视图； 25)CYView :"Y坐标"的视图； 26)CMSDITemplate :文档视图关联类； 27)CRWGPRA卯应用程序类； 28) CRWGPRDoc :文档类。 与现有技术相比，本发明具有以下优点 (1).读取显示功能 兼容SIR和REMAC两种雷达数据格式，能同时读入多个文件，能以三种方式(波形 堆积图、灰度堆积图、彩色堆积图)显示雷达数据，能任意放大和縮小图形。 (2).数据处理功能 能根据不同的雷达原始波形选择不同的参数进行波形处理，也可任意设置数据处 理的顺序和参数并自动保存，以实现数据的批量化处理。
(3).自动解释功能 根据不同的地质环境选择不同的病害类库，得到不同的判决规则，对数据进行解 释，能够自动跟踪层面位置，并画出路基病害分布图。
图1为本发明数据处理的分类示意图；
图2为本发明路基病害的反演示意图。
具体实施例方式
下面结合附图对发明作进一步说明。 如图1、图2所示，一种探地雷达数据处理方法，该方法RAILGPR1.0的界面开发 环境选择的是Visual 0++2005，计算核心语言采用的是區11^8 7. 1，通过前台的VC调用 MATLAB引擎，使用MATLAB在后台计算。
1.程序操作一般步骤
a.数据单项处理 探地雷达数据处理方法RAILGPR1. 0数据处理的主要作用是去除杂波和提高雷达 剖面分辨率；杂波可以分为两种，一种是由探地雷达硬件本身及其工作方式带来的，是固定 不变的，另一种是由检测时介质性质和周边环境造成的，是随外界条件不同而不同的；因 此，可以将探地雷达的数据处理归为三大类，如图1所示， 一类是预处理，用来去除固有噪 声；一类是噪声和干扰的去除，根据外界杂波种类而定；一类是分辨率的提高，用来增强有 效信号； b.自动化批处理 数据处理是一个相互穿叉的过程，对一组数据的处理顺序和往返次数应视具体情 况而定，对于大批量的数据，先进行样组处理过程试验，确定处理方法和流程后再做自动化 批处理； c.路基病害反演 判别路基结构是否发生病害，反演出路基层面图和路基病害位置和类型，反演过
程如图2所示； 2.程序设计 利用Visual C++2005新建一个项目，项目名称为RWGPR，项目类型为多文档；在整
个项目中，包含以下28个类，如图3所示，各类的分工如下 1) CBXTFrame :"波形堆积图"的子框架； 2)CHDTFrame :"灰度堆积图"的子框架； 3) CCSTFrame :"彩色堆积图"的子框架； 4)CMainFrame :程序的主框架； 5) CJMSBFrame :"层面识别"的子框架；6)CBHSBFrame :"病害识别"的子框架； 7) CAboutDlg :显示"关于"的对话框； 8) CBCSZDlg :显示"标尺设置"的对话框；
9) CCYDZLDlg :显示"采样点整理"的对话框； 10) CGTLBDlg :显示"高通滤波"的对话框； ll)CJMSZDlg :显示"层面设置"的对话框； 12) CKLBHDlg :显示"KL变换"的对话框； 13) CRMCHFDlg :显示"RMC雷达头文件"的对话框； 14)CSIRHFDlg :显示"SIR雷达头文件"的对话框； 15) CSJWLXLDlg :显示"神经网络训练"的对话框； 16) CXCGZYXDlg :显示"消除轨枕影响"的对话框； 17) CXCZOBDlg :显示"消除直耦波"的对话框； 18) CYWXBQZDlg :显示"小波去噪"的对话框； 19) CDataProDlg :显示"数据批处理"的对话框； 20)CBXTView :"波形堆积图"的视图； 21)CHDTView :"灰度堆积图"的视图； 22)CCSTView :"彩色堆积图"的视图； 23)CJMSBView :"层面识别"的视图； 24)CBHSBView :"病害识别"的视图； 25)CYView :"Y坐标"的视图； 26)CMSDITemplate :文档视图关联类； 27)CRWGPRA卯应用程序类； 28) CRWGPRDoc :文档类。
权利要求
探地雷达数据处理方法，其特征在于，该方法RAILGPR1.0的界面开发环境选择的是Visual C++2005，计算核心语言采用的是MATLAB 7.1，通过前台的VC调用MATLAB引擎，使用MATLAB在后台计算；该方法包括(1)程序操作一般步骤①数据单项处理探地雷达数据处理方法RAILGPR1.0数据处理的主要作用是去除杂波和提高雷达剖面分辨率；杂波可以分为两种，一种是由探地雷达硬件本身及其工作方式带来的，是固定不变的，另一种是由检测时介质性质和周边环境造成的，是随外界条件不同而不同的；因此，可以将探地雷达的数据处理归为三大类，一类是预处理，用来去除固有噪声；一类是噪声和干扰的去除，根据外界杂波种类而定；一类是分辨率的提高，用来增强有效信号；②自动化批处理数据处理是一个相互穿叉的过程，对一组数据的处理顺序和往返次数应视具体情况而定，对于大批量的数据，先进行样组处理过程试验，确定处理方法和流程后再做自动化批处理；③路基病害反演判别路基结构是否发生病害，反演出路基层面图和路基病害位置和类型；(2)程序设计利用Visual C++2005新建一个项目，项目名称为RWGPR，项目类型为多文档；在整个项目中，包含28个类，各类均有分工。
2. 根据权利要求1所述的探地雷达数据处理方法，其特征在于，所述的28个类，各类的 分工如下1) CBXTFrame :"波形堆积图"的子框架；2) CHDTFrame :"灰度堆积图"的子框架；3) CCSTFrame :"彩色堆积图"的子框架；4) CMainFrame :程序的主框架；5) CJMSBFrame :"层面识别"的子框架；6) CBHSBFrame :"病害识别"的子框架；7) CAboutDlg :显示"关于"的对话框；8) CBCSZDlg :显示"标尺设置"的对话框；9) CCYDZLDlg :显示"采样点整理"的对话框；10) CGTLBDlg :显示"高通滤波"的对话框；11) CJMSZDlg :显示"层面设置"的对话框；12) CKLBHDlg :显示"KL变换"的对话框；13) CRMCHFDlg :显示"RMC雷达头文件"的对话框；14) CSIRHFDlg :显示"SIR雷达头文件"的对话框；15) CSJWLXLDlg :显示"神经网络训练"的对话框；16) CXCGZYXDlg :显示"消除轨枕影响"的对话框；17) CXCZ0BDlg :显示"消除直耦波"的对话框；18) CYWXBQZDlg :显示"小波去噪"的对话框；19) CDataProDlg :显示"数据批处理"的对话框；20) CBXTView :"波形堆积图"的视图；21) CHDTView :"灰度堆积图"的视图；22) CCSTView :"彩色堆积图"的视图；23) CJMSBView :"层面识别"的视图；24) CBHSBView :"病害识别"的视图；25) CYView :"Y坐标"的视图;26) CMSDITemplate :文档视图关联类；27) CRWGPRApp :应用程序类；28) CRWGPRDoc :文档类。
全文摘要
本发明涉及探地雷达数据处理方法，该方法RAILGPR1.0的界面开发环境选择的是Visual C++2005，计算核心语言采用的是MATLAB 7.1，通过前台的VC调用MATLAB引擎，使用MATLAB在后台计算。与现有技术相比，本发明具有读取显示功能、数据处理功能、自动解释功能等多种功能。
文档编号G01V3/38GK101738642SQ20081020245
公开日2010年6月16日 申请日期2008年11月10日 优先权日2008年11月10日
发明者凌保林, 叶培旭, 廖立坚, 杨新安 申请人:同济大学