探地雷达成像算法研究

探地雷达成像算法研究

摘要

探地雷达(Ground Penetrating Radar,简称GPR)集无损检测、穿透能力强、分辨率高等众多优点而成为检测和识别地下目标的一种有效技术手段。性能优良的探地雷达成像方法有助于精确定位地下目标，同时提高对目标的检测和识别能力，从而推动探地雷达在城市质量监控、地质灾害、考古挖掘、高速公路无损检测、地雷探测等各个方面得到更广泛的应用。

本文以中国电波传播研究所的探地雷达LD-2000为实验设备，从中读取探测数据。以MATLAB为软件平台，实现了探地雷达数据的显示、处理、成像几个部分。其中数据显示方式包括数据的波形堆积图，剖面面色阶图以及带数据波形图；数据处理部分包括直达波的去除、背景噪声的去除、振幅增益等；雷达成像算法部分主要采用波前成像算法和投影层析成像算法。

Imaging Algorithm of Ground Penetrating Radar

ABSTRACT

GPR (Ground Penetrating Radar, referred GPR) set of non-destructive testing, penetration ability, many advantages of high resolution detection and identification of underground and become the target of an effective technical means. Excellent performance GPR imaging approach helps pinpoint underground targets, while increasing the target detection and identification capabilities, thereby promoting the quality of ground penetrating radar surveillance in the city, geological disasters, archaeological excavation, highway nondestructive testing, mine detection, etc. aspects to be more widely used.

In this paper, China Institute of Radiowave Propagation GPR LD-2000 for the experimental apparatus, reads probe data. MATLAB as the software platform to achieve a ground-penetrating radar data display, processing, imaging several parts. Wherein the data includes a data waveform display stacked, with a cross-sectional side view and a gradation data waveform; data processing section includes the removal of the direct wave, the background noise removal, the amplitude gain, etc.; radar imaging algorithm some of the major imaging algorithm and the wavefront projection tomography algorithms.

1 绪论

1．1 选题的背景及意义

雷达是利用电磁波探测目标的电子设备。发射电磁波对目标进行照射并接收其回波，由此获得目标至电磁波发射点的距离、距离变化率（径向速度）、方位、高度等信息。而探地雷达则是利用电磁波在地下媒质中的传播与反射特性进行地下不可见目标体或界面的探查与定位分辨。探地雷达实质是向地下发射电磁波，通过电磁波的回波来分析地下的情况[2]。由于地下的物质与空气相比反射情况会复杂的多，因此，与雷达相比，探地雷达的回波很大程度上会受到各种噪声的干扰。探地雷达在工作的过程中不断的接收回波的数据信息，数据量十分庞大，同时，接收的回波信息中，不可避免的会收到外界或者仪器本身产生的干扰[2]。所以在一般情况下，未经过处理的信号几乎无法判别，因此探地雷达处理准确的可视化成像图形分析越来越受到使用者的重视。再者，由于探地雷达现在的应用已经开拓了很大的应用领域，尤其是在工程地质领域的应用。探地雷达开拓应用领域的优点在于:

(1)探地雷达是一种非破坏性的探测技术，可以安全地用于城市和正在建设中的工程现场。工作场地条件宽松，适应性强 (对于轻便类的仪器)；

(2)抗电磁干扰能力强，可在城市各种噪声环境下工作，环境干扰影响小；

(3)具有工程上较满意的探测深度和分辨率。现场直接提供实时剖面记录图，图像清晰直观；

(4)便携微机控制数字采集、记录、存储和处理。轻便类仪器现场仅需3人或更少人员工作，工作效率高。

本文通过系统研究和阅读大量文献资料，着重的研究了从数据的提取到，到探地雷达回波信号的预处理，再到探地雷达成像算法的研究。

1．2 探地雷达成像算法的发展

探地雷达（Ground Penetrating Radar）是利用电磁波探测地下目标，通过分析电磁信号与地下目标的相互作用，提取目标的性质、形状等信息[1]。由于地球表层和人们生产生活的密切关系，电磁波在地下媒质中的传播与反射比空间更为复杂，有许多新的课题需要研究解决，所以探地雷达技术愈来愈受到人们的重视，得到了迅速的发展与提高，电磁波探测地下目标成为电磁波应用的一个重要领域。

国外发达国家自20世纪60年代以后，探地雷达技术得到迅速发展，国自20世纪90年代以来也开始重视和发展探地雷达技术的研究和应用，并开发出了实用的产品。但是国产品在分辨率、使用方便性、对雷达信号成像和图像解译技术等方面与国外产品存在差距，因此，国探地雷达的应用中绝大多数采用国外的产品。当前随着探地雷达技术的飞速发展，先进的高分辨数据处理和成像技术成为探地雷达技术发展的关键，成像方法也趋于多样化。

探地雷达合成孔径聚焦成像技术自上个世纪90年代初以来已得到逐步应用，同时基于雷达波和地震波在运动学上的相似性，反射地震学上的波动偏移成像技术也逐步应用于探地雷达数据处理和成像[1]。上个世纪70年代初Standford大学的J．Claerbout教授首先提出了用有限差分法解单程波动方程的近似式，用地面观测的地震数据重建地震波在地下传播过程中的波场，从这些传播过程的波场中提取地震界面剖面像的数据，组成地震偏移剖面，这种成像方法即为有限差分偏移成像。在上个世纪70年代后期，R．Stolt和J．Gazdag等人又先后提出了基于波动方程偏移的Stolt偏移成像方法和Phase—Shift偏移成像方法，由于此两种方法在计算中可以充分利用快速傅立叶交换，因此计算简单，效率高，很快得到推广。波动方程偏移成像在最近十年间迅速发展并不断完善，许多人对此做出了有益的贡献。其中，Loewenthal等人的爆炸反射面的概念，Hubral和Lamer等人的深度偏移的概念，A．J．Berkhout提出的偏移过程是一个空间卷积的概念，我国的马在田院士的高阶方程的分裂算法等都为波动方程偏移成像技术的发展作出了贡献，同时促进了探地雷达成像技术的进一步发展。当前仍有许多学者还在探索波动方程偏移成像，以期更加完善该方法，这也必将为探地雷达成像技术的发展注入新的活力。

最近几年，应用于随机不均匀介质中的时逆(Time Reverse)方法成为国际上研究的热点。时逆方法最早是由法国巴黎大学的Mathias Fink教授提出，并在超声医学中进行研究，其可以分为物理时逆和虚拟或计算时逆[3]。在物理时逆中，逆传播场是由发射天线向未知的真实介质辐射的；而在虚拟或计算时逆中，逆传播是在一个虚拟的参考介质中进行数值仿真。物理时逆在医学上已有许多应用，如肾结石的粉碎等，最近也出现了把物理时逆应用到如地雷引爆、通信等方面的探讨。而虚拟时逆主要应用于成像，最近两年，由Rice大学的Liliana Borcea教授、Standford 大学G．Papanicolaou教授、美国Lawrence Livermore国家实验室的J．Berryman

等学者组成的研究小组对随机介质中时逆成像方法的理论进行了大量深入的研究，并基于超声对此理论进行了大量的数值仿真，为时逆成像方法在各个领域的应用提供了坚实的基础。最近，Northeastern大学的Anthony J．Devaney教授在Standford 大学作报告时首先开始了时逆成像方法在探地雷达对地下目标成像中的探讨，为探地雷达成像技术开辟了新的方向[4]。

当前，探地雷达成像技术由于受到各自数学模型的影响，都有自己的不足和

优点。成孔径聚焦成像方法适用条件广，但计算复杂；Stolt偏移成像方法和Phase．Shift偏移成像方法计算简单，但无法适用于复杂的地下结构，复杂的速

度模型；时逆成像方法要求地下介质有随机变化，但变化幅度必须限制在一个很

小的围[5]。另外，探地雷达像基本上是一个二维截面图，探地雷达三维成像技术

还处于研究之中。探地雷达成像方法作为探地雷达的关键技术之一，其发展方向

是高分辨力、高效率、高精度成像。在此要求之下，在国外大量科研人员的不懈

努力下，探地雷达技术将不断提高，其应用围也将逐步拓展。在我国，探地雷达

除了用于各种建筑和公路质量监控外，还将在国防和国家安全部门有很大的应用

前景，现实的或潜在的应用还包括：地雷探测、地下掩体的探测、货物安全检查、打击毒品走私等方面[7]。

1．3 本文的研究容

本文的主要研究容如下：

第一章绪论

本章介绍了课题的选题背景和原因。介绍了探地雷达发展过程，探地雷达成像算法的国外发展状况以及今后的发展趋势。

第二章探地雷达基础

本章是从电磁波的基础理论出发，讲述了电磁波在岩石介质中的传播特性。同时分析了探地雷达的工作原理和探地雷达回波模型，研究了探地雷达的三种数据显示方式（A-scan,B-scan,C-scan），三种显示方式（Wiggle曲线，灰度图，彩色图）。第三章探地雷达数据预处理

探地雷达所接收到的信号十分复杂，脉冲在通过地下介质的过程中，波形和波幅将发生较大的变化，而脉冲余振、系统部干扰、地表不光滑或地下介质不均匀等引起的散射以及剖面旁侧的绕射等干扰，均使得实时记录图像多变和不易分辨。为

此我们需要先对我们所测量的数据进行一些处理。做奇异值分解以去除探地雷达最大的干扰成分直达波；做道增益以补偿电磁波在传播的过程的衰减；对探地雷达数据做中值滤波以去除背景噪声的影响。

第四章探地雷达成像算法的研究

本章开始介绍的两种探地雷达的扫描方式（合成孔径，实孔径）。然后主要介绍基于合成孔径扫描方式下两种成像算法（波前成像算法，投影层析成像）。

投影层析成像与CT类似，通过接收空间各方位的散射场，就可以通过层析技术实现目标的反演。只是在CT中接收的是前向透射场，而在雷达扫描中接收的是后向散射场。目标函数可以用逆拉顿变换得到[8]。

波前重建最早作为地震波迁徙技术出现在地震学中。基于地震波与电磁波传播的相似性，Cafforio和Soumekh提出了两个基本相同的SAR模型。在系统模型中Cafforio和Soumekh运用一阶波恩近似，假定目标电磁散射可以模型化为目标各个离散散射中心相应的矢量和[8]。

2 探地雷达基础

探地雷达操作规程

探地雷达操作规程 （文件编号：****-010） 共1页第1页版本/版次：D/ 0 生效日期：2016-01-01 1. 目的 为了使检测员更好地熟悉和掌握检测仪器的操作方法，保证检测数据的科学、公正和准确性，特制定本规程。 2. 适用范围 适用于探地雷达仪器 3 操作步骤 3.1测试前的安装准备 检查所有部件是否带齐，包括：电池、雷达主机、数据线、处理器电源线、信号线、工具箱、备件、固定用绑扎带、记录本； 3.2试验/检测的工作程序 （1）测试连接。将地质雷达天线通过支架安装。 （2）在扫描前调试主机并对主机进行参数设置。 （3）打开电源，控制天线移动的人员根据操作主机的人员口令，将天线紧贴待测界面上匀速移动。 （4）测试结束。按下stop结束测试点，保存文件并退出； （5）拆除信号线，拆除天线，支架。 3.3扫描之前的仪器调试和参数设置 （1）菜单系统—>设置—>调用，选择所用的天线。 （2）系统—>单位垂直刻度设为时间，单位为ns （3）测程：900M天线探测混凝土的量程约为15纳秒，为使所有有效信号完全显示，一般设置为20ns （4）采样点数：一般设为512或1024 采样点数越多，扫描曲线越光滑，垂直分辨率越好。但是采样点数增大，使得扫 描速率下降 （5）每秒扫描数：64 （6）增益点数：2 （7）垂向高通滤波器：225MHz

（8）垂向低通滤波器：2500MHz （9）数据位：16位 （10）发射率：100 KHz，发射功率越高，采集速度越快，但若采集过高，易损坏雷达系统 （11）信号位置设为手动 （12）表面设为0 （13）调出完整的直达波（首波），调整延时参数 若检测结构与上次相同，可不再次设置以上参数，系统默认上次检测参数。 （14）增益设置为自动，增益函数手动设置，可以改变增益点数多少、并且可以调整各增益点的函数大小，进而调整信号强度。增益函数调整过大，在探测资料中可能 人为造成假象。设置方法为先设为手动，再设为自动。 编制/日期：批准/日期：

探地雷达成像算法研究

探地雷达成像算法研究 摘要 探地雷达(Ground Penetrating Radar,简称GPR)集无损检测、穿透能力强、分辨率高等众多优点而成为检测和识别地下目标的一种有效技术手段。性能优良的探地雷达成像方法有助于精确定位地下目标，同时提高对目标的检测和识别能力，从而推动探地雷达在城市质量监控、地质灾害、考古挖掘、高速公路无损检测、地雷探测等各个方面得到更广泛的应用。 本文以中国电波传播研究所的探地雷达LD-2000为实验设备，从中读取探测数据。以MATLAB为软件平台，实现了探地雷达数据的显示、处理、成像几个部分。其中数据显示方式包括数据的波形堆积图，剖面面色阶图以及带数据波形图；数据处理部分包括直达波的去除、背景噪声的去除、振幅增益等；雷达成像算法部分主要采用波前成像算法和投影层析成像算法。

Imaging Algorithm of Ground Penetrating Radar ABSTRACT GPR (Ground Penetrating Radar, referred GPR) set of non-destructive testing, penetration ability, many advantages of high resolution detection and identification of underground and become the target of an effective technical means. Excellent performance GPR imaging approach helps pinpoint underground targets, while increasing the target detection and identification capabilities, thereby promoting the quality of ground penetrating radar surveillance in the city, geological disasters, archaeological excavation, highway nondestructive testing, mine detection, etc. aspects to be more widely used. In this paper, China Institute of Radiowave Propagation GPR LD-2000 for the experimental apparatus, reads probe data. MATLAB as the software platform to achieve a ground-penetrating radar data display, processing, imaging several parts. Wherein the data includes a data waveform display stacked, with a cross-sectional side view and a gradation data waveform; data processing section includes the removal of the direct wave, the background noise removal, the amplitude gain, etc.; radar imaging algorithm some of the major imaging algorithm and the wavefront projection tomography algorithms.

探地雷达的发展与现状

探地雷达的发展与现状 探地雷达的历史最早可追溯到20世纪初。1904年，德国人Hülsmeyer首次将电磁波信号应用于地下金属体的探测。1910年，Leimback和L?wy以专利形式提出将雷达原理用于探地，他们用埋设在一组钻孔中的偶极天线探测地下相对高导电性质的区域，正式提出了探地雷达的概念。1926年Hülsenbeck第一个提出应用脉冲技术确定地下结构的思路，他指出介电常数不同的介质交界面会产生电磁波反射。由于地下介质具有比空气强得多的电磁衰减特性，加之地下介质情况的多样性，电磁波在地下的传播比空气中复杂的多，之后二三十年尽管在美国出现过一些相关的专利，这项技术很少被运用到其它领域,直到50年代后期，探地雷达技术才慢慢重新被人们所重视。探地雷达在矿井(1960,、冰层厚度(1963,、地下粘土属性(1965,Barringer)、地下水位(1966,Lundien)的探测方面得到了应用。1967年，一个与stern最初用于冰川探测的仪器类似的系统被设计研制出来，1972年Procello将其于探测月球表面结构。同样在1972年，Rex Morcy和Art Drake开创了GSSI(Geophysical Survey Systems Inc.)公司，主要从事商业探地雷达的销售。随着电子技术的发展，数字磁带记录问世，加之现代数据处理技术的应用，特别是拟反射地震处理的应用，探地雷达的实际应用范围在70年代以后迅速扩大，其中有：石灰岩地区采石场的探测(1971,Takazi;1973,kithara;)、淡水和沙漠地区的探测(1974,、工程地质探测(1976,和、煤矿井探测(1975,、泥炭调查(1982,、放射性废弃物处理调查(1982,、以及地面和井中雷达用于地质构造填图(1997, )、水文地质调查(1996, ;1997,Chieh-Hou Yang )、地基和道路下空洞及裂缝调查、埋设物探测、水坝的缺陷检测、隧道及堤岸探测等。自70年代以来、许多商业化的通用数字探地雷达系统先后问世，其中有代表性的有：美国Geophysical Survey System Inc公司的SIR系统、Microwave Associates 的MK系列，加拿大Sensor & Software的Pulse Ekko系列，瑞典地质公司(SGAB)的RAMAC/GPR系列，日本应用地质株式会社OYO公司的GEORADAR系列及一些国内产品(电子工业部LTD系列，北京爱迪尔公司CR-20、CBS-900等)。这些雷达仪器的基本原理大同小异，主要功能有多通道采集、多维显示、实时处理、变频天线、多次叠加、多波形处理等，另外还有井中雷达系统，多态雷达系统，层析成像雷达系统等。国内探地雷达的研究始于70年代初。当时，地矿部物探所、煤炭部煤科院，以及一些高校和其他研究部门均做过探地雷达设备研制和野外试验工作，但由于种种原因，这些研究未能正式用于实际。90年代以来，由于大量国外仪器的引进，探地雷达得到了广泛的应用与研究。1990-1993年，中国地质大学(武汉)在国家自然科学基金

浅谈探地雷达的原理与特点

浅谈探地雷达的原理与特点 摘要：地下管线系统的建立是城市现代化建设的重要因素，但由于地下管线中的非金属管线的大量存在以及城市建设快速安全的需要，探地雷达探测技术的独特优势就显现出来，本文通过对探底雷达和地下管线的分析，为应用探底雷达在城市地下管线建设提供参考。 关键字：探地雷达；地下管线；探测技术 0 引言 随着城市现代化的发展，地下管线的密集程度也在不断地扩大。地下管线作为城市的重要基础设施之一，它一方面关系着城市居民生活及城市工业的发展，担负着巨大的社会责任，另一方面又由于它深埋于地下，具有不透明性，纵横交错、结构复杂。近年来，在许多大城市出现施工时挖断通信、电力电缆导致通讯中断、区域性停电、停产事故，这些事故给该地区经济和人们的生产生活带来了巨大的损失。因此，地下工程在施工时如何避免破坏这些地下管线就变得越来越重要，建立完整的城市地下管线系统成为现代城市快速建设的关键因素。 探地雷达(Ground Penetrating Radar，简称GPR)是一种对地下或物体内不可见的目标或界面进行定位的电磁法，并以其探测的高分辨率和高效率而成为地球物理勘探的一种有力工具。随着信号处理技术和电子技术的不断发展及工程实践的增多和经验的不断积累，探地雷达技术也得到极大提高，仪器也不断更新，探地雷达检测技术具有分辨率高、采集速度快、后期数据处理简便等特点。因此在铁路、公路、建筑、市政、考古等领域得到广泛的应用,并受到广大现场技术人员的认可和喜爱。 1 探地雷达的发展 国外探底雷达技术最早可追溯到二十世纪初，西方国家以专利形式提出将雷达原理用于探地，正式提出了探地雷达的概念。但是直到50年代后期探地雷达技术才被慢慢重视起来。探地雷达在矿井、冰层厚度、地下粘土属性、地下水位等方面的得到了应用。1967年，一个与Stern最初用于冰川探测的仪器类似的系统被设计研制出来，1972年，Rex Morcy和Art Drake开创了GSSI公司，主要从事商业探底雷达的销售。随着电子技术的发展，电子存储设备的问世，加之现代数据处理技术的应用，特别是拟反射地震处理的应用，探底雷达的应用领域迅速扩大，其中有：石灰岩地区采石场的探测、淡水和沙漠地区的探测、工程地质探测、煤矿井探测、泥灰调查、放射性废弃物处理调查、埋设物探测、水文地质调查、地基和道路下空洞及裂缝调查、水坝的缺陷检测、隧道及堤坝探测等。 自70年代以来，许多商业化的通用数字探地雷达系统先后问世，其中有代表性的有：美国GSSI公司的SIR和MK系列，加拿大Sensor&Software公司的Pulse Ekko系列。这些雷达的基本原理大同小异主要功能有多通道采集、多维显示、实时处理、变频天线、多次叠加、多波形处理等，另外还有井中雷达系统，

道路探地雷达检测方法

道路探地雷达检测方法 1 检测仪器及设备 1.1 探地雷达主机技术指标应符合下列规定: 1 系统增益应不低于120dB； 2 信噪比应不低于60dB； 3 模数转换应不低于12位； 4 信号迭加次数应可选择； 5 采样间隔宜不大于0.5ns； 6 实时滤波功能可选择； 7 应具有点测与连续测量功能； 8 应具有手动或自动位置标记功能； 9 应具有现场数据处理功能。 1.2 探地雷达天线可采用不同频率的天线组合，技术指标应符合下列规定: 1 应具有屏蔽功能； 2 最大探测深度应大于2m； 3 垂直分辨率应高于2cm。 1.3 探地雷达工作环境应符合下列规定: 1 工作环境温度-10℃～+40℃； 2 工作环境湿度<90%。 2 现场检测 2.1 检测前的准备应符合下列规定: 1 检测前应对被检工程进行现场调查，搜集设计、施工资料，了解工作条件及环境安全状况。 2 应调查施工过程中特殊施工段，记录结构物位置和影响检测工作障碍物和电磁干扰源的位置。 3 应调查已发病害，并记录其位置和类型。 4 检测前应正确连接雷达系统，并在检测前进行试运行。 5 检测前应准确标记检测里程桩号及测线位置。 6 测量轮连续采集时应保持测量轮随检测距离运转良好，计程准确。

7 在不间断通行道路检测时，检测仪器车后应跟保通警示车辆，检测车内应有专人负责安全。 2.2测线布置应符合下列规定： 1 测线布置应以纵向布线为主，横向布线为辅。 2 每车道应布设一条纵测线。应选取有代表特征部位布设部分横测线，对于重点病害异常区段宜进行加密测线，必要时应根据缺陷目标体形状布置横向测线。 3 测线每1km应有一个里程标记，标记应清晰。 2.3 介质参数标定应符合下列规定: 1 检测前应对道路结构层的介电常数或电磁波速做现场标定，每同类道路结构层宜不少于3处，取平均值为该类道路结构层的介电常数或电磁波速。当检测长度大于10km时应适当增加标定点数。 2 标定宜采用钻孔实测方法，标定记录中的界面反射信号应清晰、准确。 3 标定结果应按下列公式计算： （2.3-1） （2.3-2） 式中：ε—相对介电常数； v—电磁波速(m/s)； t—双程旅行时间(ns)； d—标定目标体厚度或距离(m)。 2.4现场参数应根据检测目标体特征而设定，检测参数应包括中心频率、时窗、采样率、测点点距等。现场检测参数设定应符合下列规定: 1 天线中心频率应按下列公式选定: （2.4-1） 式中：f —天线中心频率(MHz )； x—要求的空间分辨率； ε—相对介电常数。 2 时窗应按下列公式确定:

ReflexW读取意大利IDS探地雷达数据及常用信号处理方法 - Copy

ReflexW读取意大利IDS雷达数据DT和处理方法 Dr.Zhang/ 2020.2.23 1雷达数据 1.1雷达数据文件 .DT格式（意大利RIS系列探地雷达的数据格式） 国内其他常用的GPR数据还包括：，DT1格式（加拿大Pulse-Ekko 系列探地雷达的数据格式）以及RD3格式(瑞典MALA系列探地雷达的数据格式)。GSSI 公司的TerraSIRch SIR3000地质雷达系统(简称SIR-3000) ，格式dzt（美国SIR系列探地雷达的数据格式）。 1.2 雷达文件.DT数据 目IDS的探地雷达（GPR）数据格式。

2REFLEXW软件简介 2.1 REFLEXW软件简介 本文基于MATLAB开发的GUI界面，主要实现了REFLEXW软件的相应功能。 REFLEXW是地质雷达数据（类地震）数据处理及解释软件，应用于地质雷达的数据处理以及资料解释。Reflexw软件兼容了世界上大多数雷达的数据格式，在欧美地区，Reflexw已经成为了地质雷达数据处理的标准软件。随着地质雷达行业的发展，在国内也越来越多的人开始使用Reflexw软件。 软件特点：功能强大，可做多种滤波处理可对个雷达数据进行批量处理导入GPS数据，可绘制测线轨迹、修正地形可显示测线中的标记对不需要的雷达数据可进行删除可做2D剖面处理和3D时间切片处理。 2.2 REFLEXW软件功能分析 通过REFLEXW软件，可以实现对探地雷达数据的读入，一维滤波，校正，二维滤波，波形图观察等功能，REFLEXW提供了较为全面的滤波手段，可以将探地雷达图谱处理的更加容易观察 在REFLEXW软件中，在显示数据方面，也同样提供了大量的处理方式。其中，最主要的处理方式为Plot Options。在Plot Options选项中，包括了绘图模式，点模式比例，能量衰减和振幅比例。其中，点模式比例又分为了XY比例绘图，每样点像素，每道像素，能量衰减，振幅比例。以上的这些选项，在用户导入探地雷达图谱之后，可以对图谱进行一个基本的预处理，或者通过不同的方式来观察图谱。 3 读取.DT文件 3.1 建立预处理的工程文件 雷达信号分析时会产生一系列的文件，因此最好为每个雷达文件建立一个工程，并指定其存放位置。ReflexW会自动生成相应的文件。

基于MATLAB的二维探地雷达三维成像方法与设计方案

本技术公开了一种基于MATLAB的二维探地雷达三维成像方法，包括以下步骤：S1、对二维探地雷达进行组装和调试；S2、对采集路线中的各条车道进行数据采集；S3、三维成像模型建立；S4、对三维成像模型进行外观优化和尺寸调整；S5、三维数据导出；S6、数据格式转换；S7、三维反射电压数据模拟；S8、对最大及最小值的颜色进行设定；S9、将三维数据导入三维成像模型，得到三维成像；S10、对异常位置进行切片化显示，直观地显示出道路深层病害。本技术通过三维模型的建立与优化，形成能够容纳三维雷达数据的成像模型，在将模拟的三维数据导入后，实现三维雷达成像，并具备图像切割功能，以便观察道路内部异常图像。 权利要求书 1.一种基于MATLAB的二维探地雷达三维成像方法，其特征在于，该方法包括以下步骤： S1、准备工作：对二维探地雷达进行组装和调试，将其安装在雷达车上，确认进行数据采集的路段和采集时间，并制定数据采集计划； S2、数据采集：通过雷达车按照采集计划中的采集路线，对采集路线中的各条车道进行数据采集，为让雷达车按直线行驶，在采集左轮迹带时，将左轮轧在车道左标线上；在采集右轮迹带时，将右轮轧在车道右标线上；通过二维探地雷达配套数据处理软件Pavecheck对数据进行数据完整性检查； S3、三维成像模型建立：通过生成超立方体定义域中的数据点矩阵的方法，在MATLAB中建立用于存储三维数据的三维矩阵，作为三维成像模型；S4、模型优化：对三维成像模型进行外观优化和尺寸调整，使其符合道路外形； S5、三维数据导出：通过二维探地雷达配套数据处理软件Pavecheck对数据进行导出； S6、数据格式转换：对导出的数据进行格式转换，转换为MATLAB中三维成像模型能够识别的Excel格式数据； S7、三维反射电压数据模拟：目前采集到的单车道测线为3条，将两条测线间的数据进行线性扩充模拟，每两条测线间的数据扩充4组，提高其分辨率； S8、颜色设定：对最大及最小值的颜色进行设定；

探地雷达作业指导书

浙江公路技师学院试验检测中心 作业指导书 （版本：第一版第0次修订） 文件编号：GLJX/ZD-32-2017 试验类别：探地雷达测路面厚度检测作业指导书 编写： 审核： 批准： 受控状态： 分发号： 20XX年XX月XX日发布20XX年XX月XX日实施浙江公路技师学院试验检测中心

目录 目录 (2) 1、目的与适用范围 (3) 1.1 (3) 1.2 (3) 2、依据标准 (3) 3、仪具与材料技术要求 (3) 3.1 (3) 3.2 (3) 4、方法与步骤 (4) 4.1 准备工作 (4) 4.2 测试步骤 (8) 5、数据导出 (8) 5.1 (8) 5.2 (8) 6、数据处理 (11) 6.1 (11) 6.2 (11) 7、报告 (11)

GSSI Sir30数据采集作业指导书 1、目的与适用范围 1.1本方法适用于采用短脉冲雷达检测沥青路面面层厚度，可供道路施工过程质量控制、质量评定及旧路调查使用。 1.2本方法不适用于潮湿路面或用富含铁矿渣集料等介电常数较高的材料铺筑的路面。 2、依据标准 《路基路面现场测试规程》JTG E60-2008 3、仪具与材料技术要求 3.1雷达测试系统由承载车、发射天线、接收天线和控制单元等组成。 3.2测试系统技术要求和参数 (1)距离标定误差：≤0.1%； (2)最小分辨层厚：≤40mm； (3)系统测量精度要求如下表： 测量精度技术要求

(4)天线：采用空气耦合方式，带宽能适应所选择的发射脉冲频率。 4、方法与步骤 4.1 准备工作 (1) 检测前应收集设计图纸、施工配合比等资料，以便合理确定标定路段。 (2) 给电池及笔记本电脑充满电，保证野外数据采集能够正常进行。 (3) 把主机、天线、大缆及其他配件列一个详细清单，出野外之前由专人照着清 单清点一遍，保证需要带的设备全部齐全。 (4) 先连接系统各个部件，笔记本电脑+主机+电缆+天线+标记杆+测距轮； 注：设备连接要求无电操作不能带电连接、插拔各个部件，电缆连接防止虚接。联机时注意电缆接口方向，电缆接头应与面板垂直，拧紧，与主机端旋转至红线处。天线端旋转至三个小卡槽露出，同时注意固定电缆。电缆与天线应用环行扣连接。 (5)数据采集---参数设置 双击采集软件图标，进入如下系统设置界面：

探地雷达数据采集与解释

探地雷达数据采集以及解释 山东大学岩土中心

第1章．探地雷达简介 1．1工作基本原理 探地雷达（Ground Penetrating Radar，简称GPR）是利用频率介于106～109Hz的无线电波来确定地下介质的一种地球物理探测仪器。随着微电子技术和信号处理技术的不断发展，探地雷达技术被广泛应用于工程地质勘察、建筑结构调查、公路工程质量检测、地下管线探测等众多领域。 探地雷达的基本原理如图1所示。发射天线将高频短脉冲电磁波定向送入地下，电磁波在传播过程中遇到存在电性差异的地层或目标体就会发生反射和透射，接收天线收到反射波信号并将其数字化，然后由电脑以反射波波形的形式记录下来。对所采集的数据进行相应的处理后，可根据反射波的旅行时间、幅度和波形，判断地下目标体的空间位置、结构及其分布。探地雷达是在对反射波形特性分析的基础上来判断地下目标体的，所以其探测效果主要取决于地下目标体与周围介质的电性差异、电磁波的衰减程度、目标体的埋深以及外部干扰的强弱等。其中，目标体与介质间的电性差异越大，二者的界面就越清晰，表现在雷达剖面图上就是同相轴不连续。可以说，目标体与周围介质之间的电性差异是探地雷达探测的基本条件。 图1 探地雷达基本原理 1．2电磁波传播特征 探地雷达的电磁脉冲在介质中的传播速度为： v = 其中c为电磁波在空气中的传播速度，ε为介质的介电常数，常见介质的介电常数如表1所示。

材质相对介电常数材质相对介电常数 粉质粘土 6 水81 干砂3～5 灰岩4～8 湿砂20～30 花岗岩4～7 金属300 砂岩 6 PVC塑料 3.3 页岩5～15 混凝土 6.4 淤泥5～30 空气 1 海水80 粘土5～40 表1各种常见介质的介电常数 电磁波脉冲在地质界面上的反射系数为： k= 根据电磁脉冲的传播规律，在地质界面上如果反射系数为负，则相位与发射脉冲相反，若反射系数为正，则相位与反射脉冲一致。如图2和图3，可以清除看到反射波相位的变化规律。 图2

探地雷达阵列成像

摘要 有大量的记载探地雷达的相关书籍和论文，我们可以对探地雷达的发展概况有一个大致的了解，不难看出探地雷达未来的发展方向。 多输入多输出（Multiple-Input and Multiple-Output，MIMO）雷达是在数字阵列雷达、多基雷达以及现代通信技术基础上发展起来的一种新体制雷达，具有显著的技术优势和巨大的发展潜力，逐渐成为下一代雷达的主要发展方向之一。 由于采用了波形分集技术，与传统相控阵雷达相比，MIMO雷达具备多项优点，如突破阵元间距半波长限制、提高角度分辨率、提高系统自由度等。MIMO雷达在提高信号处理灵活性的同时，也由于其新体制特征及结构特点，带来了新的问题，如多通道积累检测问题，波形设计和阵列结构优化问题以及系统自由度过高导致的自适应算法性能下降等问题。 雷达成像方面有很多相关的软件，比如GprMax，一款基于FDTD（时域有限差分）的探地雷达仿真软件，其生成文件包括几何图和数据剖面图，这些生成文件是由MATLAB 读取的。 雷达成像最关键的部分就是成像的算法，通过对算法的不断研究和改进，提出了很多算法，有些算法可以相互结合，从而改善成像效果，基于不同的理论基础，这些算法又衍生出很多不同的算法，诸如RD-BP算法、TCC-BP算法等，针对不同的算法，借助MATLAB 对算法进行仿真，并得到仿真图像。 关键词：多输入多输出雷达；数字波束形成；反向投影算法；MATLAB

Abstract There are a large number of records of books and papers related Ground Penetrating Radar (GPR).W e can have a general understanding of the development of GPR, it is not difficult to see the development direction of the GPR in the future. Multiple-input and multiple-output (MIMO) radar is a new mode of radar system developed on the basis of digital array radar, multistatic radar and modern communication techniques. Since its obvious technical advantages and huge development potential, MIMO radar will become one of the major directions of the radar system development in future. By using waveform diversity technique, Through proper improvements traditional array signal processing techniques can be applied to MIMO radar. Moreover, this new type of radar system offers a new paradigm for radar signal processing. the strong demands for waveform and array structure optimization, and a noticeable performance decline of adaptive algorithms for the excessive degree of freedom, etc. There are a lot of software related to radar imaging, such as GprMax, a Ground Penetrating Radar (GPR) simulation software based on FDTD (finite difference time domain). The generated files including geometric figure and data section are read by the MA TLAB.. The most critical part of radar imaging is the imaging algorithm. Based on different theoretical basis, the algorithm derives a lot of different algorithms, such as RD-BP algorithm, the TCC-BP algorithm, etc. For different algorithms ,we can get the simulation images with the help of MA TLAB. Key words: MIMO radar, digital beamforming, back projection，MATLAB

探地雷达论

信息检索结课论文 题目：探地雷达的发展现状 学院：信息与通信学院 专业：电子与通信工程 学生姓名：甘睿 学号：1502303016 授课教师：蔡国永

探地雷达的发展现状 【中文摘要】(GPR)是一种利用电磁波对地下区域进行无损探测的装置，它根据电磁波在地下介质不连续处产生的反射和散射等现象来反演地下场景和目标的信息，实现对地下目标的检测和识别。通过GPR成像技术，可以对地下目标的信息进行直观显示，便于人们对地下目标的解译。由于GPR成像技术具有非破坏性、高分辨率、探测速度快、操作方便、安全性高等优点，其在军事和民用领域中得到广泛的应用，具有很好的发展前景。与传统雷达成像不同，GPR成像场景较为特殊，具体表现为GPR工作于近场条件，且工作场景和目标特性均比较复杂。 【英文摘要】GPR(GroundPenetratingRadar)uses electromagnetic waves to probe theunderground region nondestructively. GPR can be used to reconstruct the information ofunderground scenes and targets by the reflection and scattering phenomena stimulatedby the discontinuity of the underground medium, as well as underground targetdetection and recognition. The information of underground targets can be displayeddirectly by GPR imaging technique, which is also convenient for the interpretation ofunderground targets. Because of the advantages such as nondestructive, high resolution,fast probing and high safety, GPR has been widely used both in military and civilianapplications, and it possesses great developing potential. Unlike traditional radarimaging, the observing scene of GPR is relatively special because GPR works in nearfield conditions and the characteristics of both the imaging scene and targets are rathercomplicated. In order to solve these problems, an in-depth and detailed research hasbeen carried out in this thesis in the aspects of time domain imaging technique,frequency domain imaging technique, super resolutionimagingtechniqueandcompressive sensing imaging technique. 【关键词】探地雷达发展现状研究技术 【正文】 1.探地雷达的研究背景和意义 探地雷达(Ground Penetrating Radar, GPR)是用来对地下目标或场景进行探测的雷达系统。利用电磁波在地下介质电磁特性不连续之处产生的反射和散射象，GPR可以获取地下目标或场景的几何和物理信息。与传统的对地探测方式相比，GPR具有不破坏探测场景、可重复性强、采样密集、安全可靠、高分辨率、操作灵活、经济性强等优势，己在军事和民用的诸多领域得到广泛应用。 在系统层面，GPR存在着参数设计、信号形式选择、系统结构设计、天线设计、数据采集等方面的问题。这些问题是GPR系统性能和应用范围的决定性因素，在GPR相关研究中具有基础性的作用。在应用层面，GPR己在地雷或未爆炸武器(Unexploded Ordance, UXO)探测、无损探伤、冰川冰层探测、冻土探测、地质水文探测、市政工程、管线定位、考古勘测、海岸环境探测、河流沉积探测、风蚀沙丘探测等领域发挥了巨大作用。在信号处理层面，GPR 主要存在介电常数估计、直祸波与直达波去除、噪声抑制、射频干扰抑制、目标检测与识别、成像等方面的问题。信号处理方法的性能从一定程度上影响了GPR的性能，对降低数据处理和解译的难度具有重要的意义，其中目标检测和识别是信号处理最本质的目的，而其余方面的问题则是目标检测和识别的预处理过程或辅助手段。 由于GPR探测的根本目标是为了对地下目标或场景进行检测和识别，而最初人们直接使用一维或二维回波来完成检测或识别工作，这种工作方式既不方便直观，又对操作人员的专