基于CS的脉冲GPR成像技术研究（20111）

分类号                                  密级

  U D C                                   编号

 

 

 

CENTRAL SOUTH  UNIVERSITY 

 

 

 

 硕士学位论文

 

 

 

 

 

论文题目基于CS的脉冲GPR成像技术研究

 

2011年10月

 

 

 

 

原创性声明

 

本人声明，所呈交的学位论文是本人在导师指导下进行的研究工作及取得的研究成果。尽我所知，除了论文中特别加以标注和致谢的地方外，论文中不包含其他人已经发表或撰写过的研究成果，也不包含为获得中南大学或其他单位的学位或证书而使用过的材料。与我共同工作的同志对本研究所作的贡献均已在论文中作了明确的说明。

 

 

 

       作者签名：             日期：     年  月  日

 

 

 

学位论文版权使用授权书

 

本人了解中南大学有关保留、使用学位论文的规定，即：学校有权保留学位论文并根据国家或湖南省有关部门规定送交学位论文，允许学位论文被查阅和借阅；学校可以公布学位论文的全部或部分内容，可以采用复印、缩印或其它手段保存学位论文。同时授权中国科学技术信息研究所将本学位论文收录到《中国学位论文全文数据库》，并通过网络向社会公众提供信息服务。

 

 

作者签名：           导师签名         日期：     年  月

 

摘 要

 

压缩传感(CS) 理论是在已知信号具有稀疏性或可压缩性的条件下对信号数据进行采集、编解码的新理论。压缩传感采用非自适应线性投影来保持信号的原始结构，能通过数值最优化问题准确重构原始信号. 压缩传感以远低于奈奎斯特频率进行采样，在高分辨压缩成像系统、视频图像采集系统、雷达成像以及MRI医疗成像等领域有着广阔的应用前景。压缩传感理论（CS），可通过求解一个凸最小ℓ1范数问题的非自适应线性测量重建稀疏信号。本文介绍了一种新型宽带合成孔径成像数据采集系统，它是基于CS理论利用图像空间中的点状目标具体稀疏性的特性来处理数据的技术。它代替了传统的Nyquist速率采样测量传感器，使用线性投影测量返回的信号和随机向量一起被用来作为测量数据。此外，沿合成孔径扫描点随机抽样可以合并在数据采集系统中。 CS测量所必需数据可以比均匀采样时数据少一个数量级。在地面探地雷达（GPR）地下成像中的应用中，典型的图像只包含几个目标，因此目标图像满足信号具有稀疏性的CS理论必要条件。我们通过使用CS测量仿真和探地雷达实验数据表明，与标准的成像效果相比，通过CS测量随机采样获取的稀疏目标空间图像也只有很少的杂乱的回波。

    利用电磁波在非金属媒质中的穿透性能，探地雷达可实现多种媒质中目标的检测、定位和识别。该技术具有非破坏性、穿透能力强、分辨率高、操作方便和费用低廉等优点，广泛应用于军事和民用的多种探测场合。本论文针对探地雷达应用中的目标二维和三维重构问题研究了脉冲探地雷达高分辨成像相关技术。

关键词：压缩传感，合成孔径，探地雷达，最小ℓ1范数法，随机采样，csuGPR数据处理专业软件

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Abstract

Compressed Sensing(CS) theory is a novel data collection and coding theory under the condition that signal is sparse or compressible. It first employs nonadaptive linear projections that preserve the structure of the signal, and then the signal reconstruction is conducted using an optimization process from these projections. Different from the traditional signal acquisition process, compressive sensing, which is a new theory that captures and represents compressible signals at a sampling rate significantly below the Nyquist rate. It has broad applications such as high resolution compressive imaging, image and video processing systems, Radar imaging, MRI imaging, etc. The theory of compressive sensing (CS) enables the reconstruction of sparse signals from a small set of non-adaptive linear measurements by solving a convexℓ1 minimization problem. This paper presents a novel data acquisition system for wideband synthetic aperture imaging based on CS by exploiting sparseness of point-like targets in the image space. Instead of measuring sensor returns by sampling at the Nyquist rate, linear projections of the returned signals with random vectors are used as measurements. Furthermore, random sampling along the synthetic aperture scan points can be incorporated into the data acquisition scheme. The required number of CS measurements can be an order of magnitude less than uniform sampling of the space–time data. For the application of underground imaging with ground penetrating radars (GPR), typical images contain only a few targets. Thus we show, using simulated and experimental GPR data, that sparser target space images are obtained which are also less cluttered when compared to standard imaging results.

Based on electromagnetic wave penetrating property through non-metal medium,ground penetrating radar(GPR)can implement the detection,location and shape reconstruction of targets buried under the earth's surface．It has many advantages such as nondestructive detection,strong penetrability,high resolution performance operational convenience low cost and so on. Aimed at subsurface targets two dimension and three dimension shape reconstruction,high resolution GPR imaging techniques are studied in this thesis．

Keywords: Compressive sensing，Synthetic aperture，Ground penetrating radar (GPR)，ℓ₁ Minimization，random sampling，csuGPR Data processing professional software

 

 

 

 

 

 

目 录

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

                        

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

第一章 绪论

探地雷达技术的出现和发展很大程度地扩展了人类认识和改造客观世界的能力，而人类对表层下隐蔽目标的感知需求又不断地促进着这一技术的发展。本章简要介绍了探地雷达技术的发展概况，分析了探地雷达系统研究的技术路线，并对冲激脉冲体制探地雷达和高分辨雷达成像技术做了较为详细的阐述，最后介绍了本文的内容安排。

 

1.l探地雷达简介

 

雷达最早出规在三十年代后期。早期的雷达只能发现目标和测量目标的距离[l5,16]。人们把它叫做“无线电探测和测距”(Radar:Radio Detection And Ranging)。经过几十年的发展，雷达的功能早已超出了探测和测距的范围，在宇宙探测、地球遥感、地下目标探测等多方面都显示出无可比拟的优越性，成为人们获取空间目标信息的重要手段[16-25]。随着人类对自然界认识的逐步深化，人们对地表下世界的探知要求变得越来越深入和迫切，从探测当今世界危及人生安全的地雷、考古发现到地质勘探、工程建设，要求的探测深度从几十厘米到数千米不等。在现代国民经济建设的众多部门和国防建设中，浅层地下目标的探测和识别己成为一个重要的研究方向[26-31]。在不破坏地表结构的情况下获取地下未知区域的信息显得尤为重要[32]，尤其是当把成像技术应用其中时。表层穿透雷达又称为探地雷达(GPR:Ground Penetrating Radar)是当今各种浅层地下目标非破坏性探测技术中最具有应用前景和发展前途的方法之一[33-43]。

 

1.1.1 探地雷达特点与应用

 

GPR又分为表层穿透雷达(SPR:Surface Penetrating Radar)和表层下雷达

(SSR:Subsurface Radar)，是指利用电磁波在媒质电磁特性不连续处产生的反射

和散射实现非金属覆盖区域中目标的成像、定位进而定性或者定量地辨识探测区

域中的电磁特性变化，实现对探测区域中目标的探测[44-57]。简单地说，GPR的任务就是描述目标的几何和物理性质。GPR所具有的优于其它遥感技术特点包括:快速、高空间分辨率、对目标的三维电磁特征敏感，可以测量前方区域甚至机载测量。其运动平台可分为飞机、地表车载和手持、地下凿洞以及多洞间探测[48-52]。本文主要讨论地表工作的GPR。

作为一种地下浅层探测工具，GPR在实际应用中与其它地下探测技术相比，具有以下优点:

.GPR是一种非破坏性探测技术，可对各种非金属覆盖区域内部进行无损探测。

.GPR可远距探测未知区域中的目标分布，避免了潜在的危险。

.GPR可反映出媒质的电磁参数变化，对金属和非金属目标都可以进行探测和识别。

.GPR可实现较高的空间分辨率，能满足多种探测需求。

配置合适的天线和收发装置，GPR几乎可以用于任何非金属体内的目标探测，

以下是目前己经成功的应用[29,54-61]:

.市政工程，如:坝、桥、高速公路的评估

.埋地管线的定位

.地下设施的确定，如:隧道和地下室

.考古探测

.地质与水文地质探测

.冰层与冻土探测

.地雷、UXO探测

.伪装的爆炸物和毒品的探测与识别

.无损评估/探测(NDE加DT)

 

1.1.2 探地雷达分类

 

按照数据采集方式的不同，GPR可分为地表式GPR、钻孔式GPR等;按照其装载方式，可分为手持式、车载式、机载式等。不同体制的GPR有不同的探测能力和应用环境。GPR的体制主要是由其调制方式决定的，调制方式则是根据穿透深度、分辨率、电磁干扰程度以及体积和成本等因素的特殊要求进行选择的[20，62-65]。理论上，只要工作频带相同，系统的分辨能力和穿透性能是相同的，获取目标的散射信息也是相同的。然而不同的信号形式有其不同的系统实现特点。目前最为常见的商业或实验性GPR系统信号形式是冲激脉冲，其次是线性调频，然后是步进变频和噪声信号。下面对这四种体制的探地雷达做一简单的介绍。

(一)冲激脉冲。为兼顾探测深度和分辨率，要求系统具备宽带特性。冲激

脉冲是一种典型的宽带信号，这种体制的探地雷达系统体积小，测量速度快，广

泛应用于浅表层探测中[66-68]。系统发射纳秒级脉冲，通过等效采样技术对回波脉冲进行采样，目标回波包含了宽带电磁散射特征，可进行成像和识别。由于冲激脉冲的辐射能量有限，因此这种体制的探地雷达系统不适用于深层探测。

(二)线性调频。系统发射线性调频信号，根据预知的地下介质的频率衰减

特性以及可能的地下目标的频率响应特性预先设定工作频带。其优点为分辨率高，发射频谱易于控制，具有很宽的动态范围;缺点是体积大、成本高，且系统比较复杂[69,70]。，目前大多数线性调频GPR仅限于对几十厘米以内表层的探测，如机场跑道和高速公路等表层中的结构异常或孔穴的探测。

(三)步进频率。发射信号可以是点频，也可以是一些特定间隔的频率，接

收端采用孔径天线在地表接收来自地下区域的后向散射信号并对信号的幅度和相位进行测量[71,71-74]。发射信号的窄带特性使天线的设计较之宽带系统相对容易但其工作频率须根据地下介质的频率衰减特性仔细选择。

(四)噪声信号。系统发射噪声信号，接收信号与经过时延的发射信号相关

处理实现对地下目标的检测和定位[75，76]。系统测量功率与目标反射系数幅度的平方成正比。

下表给出了四种体制的探地雷达性能比较[44，56]。

 

 

 

 

 

 

 

表1.1四种体制探地雷达的性能比较

 

性能   信号 形式	辐射 能量	发射 波形可控	动态 范围	硬件 复杂度	电磁 兼容	测量 速度
冲激脉冲	小	不可控	小	简单	差	快
线性调频	大	精确可控	较大	较复杂	较差	快
步进频率	大	精确可控	大	较简单	好	慢
噪声	较小	可控	较小	复杂	差	较慢

 

 

各种体制的GPR具有各自的特点，有些体制之间还存在着应用上的互补性。

对浅层勘探而言，冲激脉冲体制探地雷达可以通过时间窗截取技术减小直达波和

地面反射波的干扰，得到高分辨的地下区域成像结果。另外，冲激脉冲体制测量

速度快。因此目前绝大多数GPR产品是冲激脉冲体制系统。线性调频体制基于其

大动态范围、高辐射功率可有效的应用于深层有耗媒质中的目标探测。步进频率

体制探地雷达的实现原理简单，但受制于高性能的元器件，且后续的信号处理工

作较为繁琐。噪声信号体制探地雷达系统实现较为复杂[24,25,44,77]。当前国外已有适用于各种类型的军用和商用GPR，工作频带一般在10MHz到3GHz之间，探测深度从几厘米到数十米，探测分辨率最高可以达到厘米量级。

 

1.1.3 冲激脉冲探地雷达系统构成

 

在不同的应用背景中，GPR的实现形式和测量方式可能千差万别，但其基本

的系统构成大体类似。典型的冲激脉冲体制GPR系统构成和探测原理如图1.1所示:

图1.1冲激脉冲体制GPR系统构成

在计算机的控制下，时序控制电路(包括比较器、可变电平、快斜坡信号、

慢斜坡信号和控制电路等电路模块)输出同步脉冲和取样脉冲。同步脉冲触发脉

冲源发射纳秒级宽频带窄脉冲信号，经由位于地面上的宽带发射天线藕合到地下。当发射的脉冲波在地下传播过程中遇到电磁特性不同的介质界面、目标或局域介质不均匀体时，一部分脉冲波能量被反射回到地面，由地面上的宽带接收天线所接收。取样电路在取样脉冲的控制下，按等效采样原理将接收到高速重复的脉冲信号变换成低频信号。该信号送往数据采集卡，经过放大、滤波，再进行A\D变换，通过CompactPCI总线传输给计算机模组。计算机模组的应用软件对数据进行信号处理和成像，并在显示器上显示出来。

 

1.2国内外研究现状

 

1.2.1探地雷达系统

 

GPR技术起源于德国科学家在研究埋地特性时的专利技术。1910年，Leimbach和Lowy在德国专利中提出了利用电磁法探测埋藏物体的方法，他们将偶极子天线埋设在孔洞中进行发射和接收电磁波。由于高导电率的媒质对电磁波的衰减作用，通过比较不同孔洞之间接收信号的幅度差别，可以对媒质中电导率高的部分进行定位。1926年，Hulsenbeek在其研究工作中采用了脉冲技术来确定地下埋设物的结构特征。他注意到，不仅仅是电导率，任何介电特性的变化都将造成电磁波的反射。他的发现首次确定了地下目标电磁波回波信号与地下介质及目标间的本质联系，从而为GPR在技术和方法上提供了优于其它地球物理探测方法的理论依据[41，43，47，78-80]。直到上个世纪六十年代末、七十年代初，等效采样技术和亚纳秒脉冲产生技术的发展，从技术角度加速了探地雷达的发展。同时美国阿波罗月球表面探测实验从应用角度促进了探地雷达技术的广泛使用。1972年，第一个探地雷达设备制造公司(GcofogicalsurveysystemInc.)宣告成立，这直接导致了探地雷达研究投入的增加，反过来又促进了探地雷达产业的发展。20世纪70年代以后，探地雷达技术开始在市政工程、考古、地质、探雷等方面的研究和应用变得频繁起来[81]。早期的探地雷达系统属粗距离分辨和方位分辨系统。随着硬件技术的提高和近场合成孔径雷达技术的发展，高分辨率二维成像作为一项探地雷达基本功能在很多系统采用，阵列天线技术的采用，使得三维实时成像成为可能。二维或三维成像对于探地雷达数据解释是一个质的飞跃，因为它良好的可视化降低了对判读员的要求，使得原本不直观的数据易于从形态上和局部散射强弱上来理解。

自20世纪60年代以来，随着高速脉冲形成技术、采样接收技术及计算机技术的迅速发展，探地雷达技术也得到了迅猛的发展，地下浅层目标探测得以实现。国外开始出现适用于各种类型的军用和商用探地雷达，探测深度从几厘米到数百米，探测分辨率可以达到厘米量级。美、日、加拿大等国家的公司相继推出了自己的产品，民用GPR代表性的系统如美国GSSI公司SIR系列产品如图1.2(a)所示、加拿大SSI公司pulseEKKO系列产品如图1.2(b)、日本OYO公司的YLRZ系列产品等。这些系统均为通用系统，根据不同探测目的和探测深度采用不同频段选用装备不同的天线，用于各种地下目标成像探测。

               图1.2(a) 美国GSSI-SIR系列GPR及数据处理软件

图1.2(b)加拿大SSI-pulseEKKO系列GPR及数据处理软件

国内对探地雷达的研究起步较晚，国内产品在分辨率、使用方便性、对雷达信号成像和图像解释技术等方面与国外产品存在差距。从总体情况分析，硬件系统性能与国外产品有阶段性差异；并且没有自主知识产权的解译软件；技术指标不能满足应用需要等。从应用情况分析，由于系统性能原因，国内研发的系统基本没有使用价值，尤其是在公路探测（定量）方面。

但是由于不断引进和借鉴国外的先进技术，近年来在该领域内也取得了较为显著的研究成果。航天部25所从20纪80年代中期开始进行用于地雷检测的探地雷达技术的研究和实验。20世纪90年代以来，我国开始重视探地雷达技术的研究和应用，并且开发出新的实用产品，如北京爱尔迪公司的CBS-9000和CR-20 GPR系统如图1.2(c)所示，中国电波传播研究所研制的LTD-3探地雷达、煤炭科学总院重庆分院开发的 KDL-3、4 等[53]。

图1.2(c)国内爱尔迪公司CR-20 GPR系统及实际的数据处理图像

 

目前国际上从事GPR技术研究的公司及产品主要有:美国GSSl公司的SIR

系列;加拿大SSI公司的pulseEKKO系列;瑞典MalaGcoscience公司的RAMAC

系列等，其典型的GPR产品及其信号体制和工作频段见下表[82-85]。

 

表1.2典型的商用GPR产品

产品型号	制造商	信号体制	频段/MHZ
SIR-2000	GSSI,USA	冲激脉冲	16-1500
PulseEKKO 1000	SSI,Canada	冲激脉冲	110-1200
RAMAC	MalaGeoscience,Sweden	冲激脉冲	25-1000
KSD-21	Koden,Japan	冲激脉冲	50-2000
Zond 12c	Radar Systems,Latvia	冲激脉冲	38-2000
GeoRadar	GeoRadar,USA	步进频率	100-1000
SPR-Scan	ERA Technology,UK	冲激脉冲	500-1000

 

这些设备的共同特点是均配有不同主频的天线以满足不同探测深度的需求。

以GSSI公司的SIR-ZOOO系统为例，它可以装配GSSl公司研发的各种天线并完成单通道、多通道车载或手持式探测。从表中也可以看出，目前GPR技术发展的主流便是冲激体制，本文也正是对这种体制GPR的高分辨成像技术进行研究的。

国内对探地雷达的研究起步较晚，近年来在该领域内也取得了一定的技术进步，不少研究单位也推出了自己的探地雷达样机。如中国电波传播研究所研制的LTD一3探地雷达;煤炭科学总院重庆分院开发的KDL-3、4型矿井防爆探地雷达;艾迪尔公司的CBS-900和CBS-2000探地雷达[86-91]。在国家“863”计划项目的资助下，国防科学技术大学电子科学与工程学院研制了一套冲激脉冲体制的探地雷达系统--RadarEye[46，92，93]，成功应用于公路质量评估和高分辨地下目标成像场合。

 

1.2.2 探地雷达信号处理

 

GPR是否可以有效应用，不仅取决于硬件系统的性能，同时取决于信号处理

和数据解译技术。相对于硬件技术的发展，探地雷达信号处理技术更加滞后。反

过来也表明信号处理技术的生命周期比硬件技术更长。探地雷达信号处理技术主

要包括目标检测、成像、识别和介质分层[94-96]。具体实现时，一般要在前端进行信号预处理工作，包括随机噪声抑制、直达波和射频干扰抑制[97-100]。为增强信号处理的直观性和人机交互性能，还要进行原始数据和处理结果的可视化处理。典型的探地雷达全系统设计路线如图1.3所示。

 

 

 

 

 

 

 

      

图1.3 GPR系统研究技术路线

目标检测主要包括能量检测、相关检测、特征检测(基于模型的特征方法)

和变换域检测算法(小波域检测算法、双谱域检测算法)[95,101-106]。检测算法要求有高检测率和低虚警率。对于探雷系统而言，则是要在保证漏报率为O的前提下尽可能降低虚警率[107]。检测结果要直观快速的显示在终端设备上以供操作员进行判决。

介质分层是GPR进行公路探测时最具价值的一个应用[73,108-110]。基于一维逆散射的层状介质电磁参数估计和层厚度估计技术己经发展的比较完善。但实际的公路各层并不是理想的均匀媒质，有必要引入随机媒质中波的解析理论和强起伏理论进行高分辨介质分层。

目标识别是GPR信号处理中最具挑战性的一个课题。探测地域电磁环境的复

杂性和埋地目标的多样性使得目标特征不变量的确定非常困难[111-113]。现有的识别算法一般是对确定己知探测区域内有限的多个己知目标的检测识别或是对某一类目标(地雷)的识别。可采用的方法有基于回波相位的模板匹配和基于高阶谱的特征不变量提取。当目标散射信号的空间采样足够多时，就可以借助于探地雷达成像技术进行目标的分类识别。

雷达成像技术则是GPR应用中最实用的处理方法，但同时它对目标散射信号

的采样要求也最高，包括大采样区域和高采样密度[22,114-117]。首先要求探地雷达天线要具有大主瓣宽度，这样才能有效地进行合成孔径处理。同时还要求探地雷达在天线有效照射范围内进行密集采样。典型的收发一体探地雷达一维合成孔径和二维合成孔径扫描示意图如图1.4所示。

图1.4 GPR一维和二维扫描示意图

 

与对空雷达不同，GPR的探测目标从数厘米的地雷、管道到横向扩展数米的

地下掩体、矿藏，探测深度从几十厘米到数百米不等。其应用需求也从目标的三

维形状和电磁特性精细重构到大尺度目标的低分辨反演而各不相同。这些因素造

成了GPR成像技术方面很大的差异。但各种成像技术都是基于目标散射信号的空

间采样，建立GPR成像基本理论框架如下图所示:

图1.5 GPR成像基本理论框架

 

探地雷达成像技术从电磁理论角度考虑是典型的逆散射问题，而从雷达技术

领域考虑则和SAR成像有相似之处。这就决定了这一技术沿两个方向向前发展。

求解逆散射问题时，目标通常被视作具有特定电磁特性的散射体，成像的目的便

是恢复出目标体的形状和电磁参数。而采用合成孔径成像技术时，得到的通常是

目标散射中心的重构结果，包括散射中心的位置和散射强度。理论上讲，在对探

地雷达全系统进行精细建模的基础上，逆散射成像技术可对探测区域的电磁参数

进行真实反演，有效应用于目标识别。但考虑到探地雷达系统的多样性、天线形

式的多样性、探测环境的复杂性和噪声的干扰，全系统的精细建模极难实现，这

就造成了基于逆散射理论的成像技术很难直接应用于工程探测。而合成孔径成像

技术通过适当的近似处理，在对探测场景合理建模的基础上可实现目标散射中心

的重构，这在一般的工程探测场合已足以满足成像探测需求，因此应用面较广。

下面对这两种成像技术做一简单的介绍。

逆散射成像技术将目标视为具有特定电磁特性的散射体，在这种建模方式下，目标散射场和目标函数之间通过LS(LS:LiPPmann-Schwinger)积分方程联系起来，可通过求解这一方程实现目标体的电磁特性反演。常用的求解方法有解析法和数值求解法。当目标和背景媒质的电磁特性相差不大时，通过一阶Bom近似可将LS方程进行一定的简化并导出目标体电磁特性的谱域形式和目标散射信号谱域形式的对应关系，这即为衍射层析(DT:Diffraction Tomography)成像算法。Devaney首先研究了均匀无耗媒质中二维目标的DT成像算法实现问题[l18]，继而采用了Tikhhonov正则伪逆算子实现了有耗媒质中目标的DT成像[l19，120]。Meincke通过引入天线辐射场的平面波谱分解函数实现了不同天线形式下有耗媒质中目标的高效DT成像算法[121]。为放宽DT成像算法对目标弱电磁散射特性的限制，cuiTiejun采用矩阵分析的方法实现了二维强散射体的DT成像[l22，123]。当目标的电磁散射特性不满足一阶Bom近似的适用条件时，还可以采用数值计算的方式进行LS方程的求解从而得到目标函数的重构结果。常用的数值计算方法有Bom迭代法(BIM:BomlterationMethod)和变型BIM(DBIM:DistortBIM)。euiTiejun运用DBIM对有耗媒质中介质目标的电磁逆散射问题进行了深入的研究，采用共扼梯度算法进行散射积分方程的求解，得到了二维目标的介电常数剖面和电导率剖面的重构结果[124]。NiezaiPing将多重网格技术和BIM法相结合提出了扩展BIM(vBIM:variationalBIM)法，算法比DBIM更高效，可应用于混合迭代反演中[l25,126]。基于小波分析的多分辨特性，Erer和YangYang等人分别将小波分析运用于DT成像和数值迭代成像中，利用传递矩阵的稀疏性质有效降低了成像算法的运算量[127-129]。上述成像算法通过对探测环境的精确建模，都在一定程度上解决了目标的电磁参数重构问题。但正因为其精确建模的要求，致使这种成像技术应用面受限，目前还很难直接应用于GPR实测数据的成像处理中。

在散射中心模型下，可采用合成孔径成像技术进行目标散射中心重构。此时

GPR成像和SAR成像有类似之处，但并不是所有的SAR成像算法都可以应用于

GPR成像。SAR成像中，目标回波信号的多普勒频移经常被用来进行高分辨成像，

而在GPR成像中这一点并不成立。GPR成像中，可通过标量波动方程建立目标散

射场和目标函数之间的关系进而进行成像处理，常用的成像算法有距离偏移(RM:

Range Migration)算法[130，131]、逆时偏移(RTM:Reverse Time Migration)算法[132]和后向投影(BP:BackProjection)算法[133-136]。JeffreyE.Mast从标量波动方程出发建立了探地雷达反向传播成像算法[l37]。Yingwei和xuxiaoyin通过分析波场连续性的谱域表现形式，导出了频率波数域偏移成像算法[138-140]。这两种成像算法都可归结为cafforio提出的RM成像算法[141]。RM成像算法将目标散射场数据转化到谱域并进行插值处理，再经逆傅里叶变换转化回空域即可得到目标的散射中心重构结果，算法流程简单，运算速度快，可有效应用于均匀无耗媒质中的目标成像问题。RTM成像算法是由Fisher首先提出的[142]，Leuschen运用匹配滤波理论也导出了算法的时域实现形式[l42]。与FDTD正向计算类似，该成像算法可处理复杂背景媒质中的目标成像问题，但运算量较大。BP成像算法的“延迟一求和”运算物理概念清晰，处理流程简单，广泛应用于SAR成像中。对探地雷达而言，媒质不均匀性对成像结果影响很大。WuRenbiao深入研究了这一问题，指出为避免速度误差引起的成像结果失真，成像处理需要一个较大的深度区间内进行[143]。Tanyer运用FDTD法分析了媒质的非均匀性对目标散射信号和BP成像结果的影响，在此基础上提出了“延迟一加权求和”的成像算法，有效改善了BP算法的性能[144]。在理想情况下，上述三种成像算法都可以得到高分辨的目标散射中心重构结果。但探地雷达对地下目标进行探测时，收发天线的间距、距地面的高度和背景媒质的均匀性都会对成像结果产生影响。考虑到这些因素，则三种成像算法各有优劣。

综上所述，目标建模方法的不同导致了成像算法和处理结果的不同。体散射模型下的成像结果真实反映了目标的材质构成和轮廓信息，可直接用于目标识别。但目标的高分辨重构依赖于雷达系统和探测环境的精确建模，包括天线形式、扫描方式、媒质电磁特性等，因此目前还很难直接应用于GPR成像中。而基于散射中心模型的建模方法简单直观，成像算法物理意义明确、概念清晰，成像结果可满足大多数工程应用的需求。

由于探地雷达(GPR)信号具有很好的稀疏性,可以很好的满足压缩传感(CS)的条件,因此将CS技术运用在GPR信号数据的获取上,可以解决海量数据采集和存储问题,能显著降低图像处理的计算量。此外还可以省去雷达接收端的脉冲压缩匹配滤波器,降低接收端对模数转换器件带宽的要求。设计重点由传统的设计昂贵的接收端硬件转化为设计新颖的信号恢复算法,从而简化雷达成像系统。

对GPR数据进行成像处理后，可采用二维和三维可视化技术对成像结果进行显示。二维成像结果的显示方式有三种:平面图、立体图和等值线图，前两种是较常用的显示方式。典型二维成像结果的三种显示方式如下图所示:

 

图1.6二维成像结果显示方式

 

对三维成像而言，图像显示变得复杂。设三维成像结果为 ，若按照二维成像结果的显示方法则只能呈现出各个表面的目标函数值，区域内部的目标函数值无法直观地显示出来。这时可以采用切片显示的方式，即任一x值或y值都对应一个二维成像结果，因此可以用一幅图像来表示。通过等间隔的在x方向或y方向进行数个二维成像剖面的抽取，并将抽取的成像剖面进行组合显示即可较为直观地显示出探测区域的目标走向及分布。也可以采用等值面图的方法进行三维显示:即将相等的各点用平滑曲面连接起来就形成了一个等值面。一般以等值面上的目标函数值和整个区域目标函数的极大值之比来描述等值面。记目标函数的极大值，则目标函数值为所在的等值面可记为dB等值面图，其中:。典型三维成像结果的各种显示方式如下图所示: