基于幅度加权的探地雷达后向投影成像方法
【专利摘要】基于幅度加权的探地雷达后向投影成像方法，属于雷达成像领域。传统后向投影成像方法存在对杂波、噪声抑制能力不足的问题。一种基于幅度加权的探地雷达后向投影成像方法，对于任意成像点A，利用时延计算成像点在各通道的散射响应幅值xA,k＝rk(t)，t＝τA,k；通过边缘提取获得双曲线回波上的点对应的响应幅值构造窗函数；对各通道的散射响应幅值xA,k进行加权，然后相干叠加完成对A点的成像；遍历所有成像点完成整个成像过程。本发明利用窗函数对散射幅度进行加权处理，改进了目标检测的精度，提高了对噪声、杂波干扰的抑制能力。
【专利说明】
基于幅度加权的探地雷达后向投影成像方法
技术领域
[0001] 本发明属于雷达成像领域，具体涉及一种基于幅度加权的探地雷达后向投影成像 方法。
【背景技术】
[0002] 探地雷达(Ground Penetrating Radar,GPR)通过向地下媒质发射高频电磁波，利 用电磁波在媒质电磁特性不连续处产生的反射、散射来实现地下目标检测与识别。随着信 号处理技术的发展和电子器件水平的提高，高分辨率数据处理和目标成像技术成为了探地 雷达最有应用前景的处理方法。
[0003] 探地雷达关于理想点目标的回波模型呈双曲线形状。为了提高图像的方位分辨 率，一般采用合成孔径成像技术(SAR)来对双曲线进行聚焦处理，实现对目标的精确定位。 探地雷达在地面沿观测线方向移动时向地面发射高频电磁波，在每条孔径处的回波信号中 都记录了目标的位置信息，扫过整个区域后即完成了合成孔径扫描。
[0004] 目前最具有代表性的成像方法分别为时域上的后向投影(Back Projection，BP) 成像方法和频域上的距离偏移(Range Migration，RM)成像方法。
[0005] ( - )后向投影成像方法。这类方法通过求出成像点后向散射回波的相干叠加结果 后，再求其幅度即可得到该点的后向散射强度，其基本思想可以用"延迟-求和"概括。后向 投影成像方法最早应用于计算机层析成像，因原理类似，在Munson，D. C.，J. D. 0 'Brien，and ff.Jenkins.tomographic formulation of spotlight-mode synthetic aperture radar. 〃P:roceedings ofthe IEEE 71.8(1983) :917-925(聚束模式合成孔径雷达的层析成 像）.，中使用后向投影成像方法进行合成孔径成像；Halman，J. I .，K.Shubert，and G·T·Ruck·〃SAR processing of ground-penetrating radar data for buried UX0 detection:results from a surface-based system,IEEE Transactions on Antennas& Propagation 46.7(1998): 1023-1027(地下未爆炸武器探地雷达合成孔径成像）.中利用后 向投影方法对地下未爆炸武器进行了三维成像，证明了该方法对地下目标成像的准确性； Cui , Guo long , L. Kong , and J.Yang,A Back-Projection Algorithm to Stepped-Frequency Synthetic Aperture Through-the-ffal1 Radar Imaging."Synthetic Aperture Radar,2007.APSAR 2007.1st Asian and Pacific Conference on 2007:123-126( -种频率步进孔径雷达的反投影算法）.中将后向投影方法应用于穿墙雷达成像，证明 了该方法可以进行穿墙高质量成像。这类方法对地下各类随机噪声和杂波干扰的抑制能力 有限，导致回波信号的信噪比达不到要求，严重影响成像结果。
[0006] (二)距离偏移成像方法。这类方法的信号处理过程在频域中进行，其关键在于频 域上的Stolt插值处理。最早在Mast J E.Microwave pulse-echo radar imaging for the nondestructive evaluation of civil structures[D].Urbana, Illinois :University of Illinois at Urbana-Champaign，1993(微波脉冲回波成像对土木结构无损检测的研 究）·中提出应用于探地雷达的距离偏移成像方法;Gu，K.，G.Wang，and J.Li · "Migration based SAR imaging for ground penetrating radar systems·"Radar，Sonar and Navigation，IEE Proceedings-151.5(2004) :317-325(基于偏移成像的探地雷达系统）·中 指出相比于射线追踪合成孔径方法，距离偏移方法由于其在空气-地面对斜入射波处理简 单、直观的优点更适合对地下目标进行成像；Lopera, Olga, et al · "Filtering Soil Surface andAntennaEffects From GPRData to Enhance Landmine Detection·〃IEEE Transactions on Geoscience&Remote Sensing 45.3(2007) :707-717(提高地雷检测的方 法：消除土壤表面和GPR数据的天线效应）.中通过线性变换和格林函数滤波对探地雷达回 波数据进行了更好的预处理，提高了图像的成像质量。这类方法的问题在于当存在分层介 质时并不符合该方法的初始条件，从而使Stolt插值不能很好地找到频域的回波数值，导致 成像结果很差，双曲线无法聚焦，严重影响了目标的检测与定位。
【发明内容】
[0007] 本发明的目的是为了有效提高后向投影方法对杂波、噪声干扰的抑制能力，针对 地下目标成像，提出一种基于幅度加权的探地雷达后向投影成像方法。
[0008] -种基于幅度加权的探地雷达后向投影成像方法，所述方法通过以下步骤实现：
[0009] 步骤一、设定探地雷达的总通道数为M，计算成像区域内任意一个成像点A到各个 通道的双程时延，那么，成像点A到第k个通道的双程时延i A,k表示为：
[0010]
[0011] 其中，成像点A坐标为A(XA，zA);
[0012] 式中，k表示通道序号，且k=l，2,......，M;iA,k表示成像点A到第k个通道的双程 时延;V表示电磁波在介质中的传播速度，v = e /:$，C表示电磁波在真空中传播的速度，大 小为3.0 X 108m/s，ε为介质的介电常数；
[0013] 步骤二、通过公式XA,k = rk(t)计算成像点Α在各通道的散射响应幅值XA,k，其中t = TA，k;
[0014]式中，rk(t)表示探地雷达在各通道的回波数据；
[0015] 步骤三、构造窗函数W;
[0016]步骤四、利用步骤三中构建的窗函数扫描成像点A，进行成像点A的成像过程；
[0017]步骤五、重复步骤一至步骤四的过程，遍历成像区域所有成像点，则完成了成像过 程。
[0018] 本发明的有益效果为：
[0019] 本发明的基于幅度加权的探地雷达后向投影成像方法为探地雷达合成孔径成像 开辟了一条新途径，在经典的后向投影成像方法的基础上对探地雷达回波信号进行加权， 获得的信噪比和成像分辨率与现有经典后投影成像技术相比能提高40-60%;改善了经典 后向投影成像算法对噪声、杂波干扰抑制能力较差的缺点，与实物相比所成像的成像质量 达98%。本发明可应用于需要高分辨率成像的探地雷达成像领域，利用本文中所提出的窗 函数对散射幅度进行加权处理，改进了目标检测的精度，对噪声、杂波干扰的抑制能力提高 30%左右。
【附图说明】
[0020] 图1为本发明涉及的预处理后的探地雷达回波B-scan图像；
[0021] 图2为本发明涉及的经典后向投影成像方法的成像结果；
[0022] 图3为本发明涉及的基于幅度加权的后向投影方法（Delay-Weight-Sum algorithm，DWS algorithm)的成像结果，
[0023] 图4为本发明涉及的两种成像方法在X轴的剖面图对比；
[0024] 图5为本发明涉及的两种成像方法在Y轴的剖面图对比；
[0025]图6为本发明的流程图。
【具体实施方式】
[0026]【具体实施方式】一：
[0027] 本实施方式的基于幅度加权的探地雷达后向投影成像方法，结合图6的流程图，所 述方法通过以下步骤实现：
[0028] 步骤一、设定探地雷达的总通道数为M，计算成像区域内任意一个成像点A到各个 通道的双程时延，那么，成像点A到第k个通道的双程时延i A,k表示为：
[0029]
[0030] 其中，成像点A坐标为A(XA，zA);
[0031] 式中，k表示通道序号，且k=l，2,......，M;iA,k表示成像点A到第k个通道的双程 时延;V表不电磁波在介质中的传播速度，.V .= C· /士 ,c表不电磁波在真空中传播的速度，大 小为3.0 X 108m/s，ε为介质的介电常数；
[0032] 步骤二、通过公式XA,k = rk(t)计算成像点Α在各通道的散射响应幅值XA,k，其中t = TA，k;
[0033]式中，rk(t)表示探地雷达在各通道的回波数据；
[0034] 步骤三、构造窗函数W;
[0035] 步骤四、利用步骤三中构建的窗函数扫描成像点A，进行成像点A的成像过程；
[0036] 步骤五、重复步骤一至步骤四的过程，遍历成像区域所有成像点，则完成了成像过 程。
[0037]【具体实施方式】二：
[0038] 与【具体实施方式】一不同的是，本实施方式的基于幅度加权的探地雷达后向投影成 像方法，步骤二所述计算成像点A在各通道的散射响应幅值的过程为，首先，设探地雷达各 通道的回波用rk⑴表示，然后得成像点A在各通道的散射响应幅值^^分别为^^=^^)4 = TA,k;式中，k表示通道序号，且k=l，2,......，M。
【具体实施方式】 [0039] 三：
[0040] 与【具体实施方式】一或二不同的是，本实施方式的基于幅度加权的探地雷达后向投 影成像方法，步骤三所述构造窗函数W的过程为，
[0041] 步骤三一、设双曲线的水平区域对应m个通道，分别为lu，k2, ...，km;对获得的探地 雷达B-scan图像进行边缘提取，得到双曲线上的点对应的响应幅值分别为m，u2,......， um，并构造一个1 XM维的矢量W'，且：
[0042]
[0043] 其中，Μ表示总通道数，m表示探地雷达回波双曲线水平区域占据的通道数；
[0044] 步骤三二、设总通道数Μ和双曲线对应的通道数m都是偶数，
[0045]步骤三三、对获得的矢量W'进行归一化处理，则得到了基于幅度加权的探地雷达 后向投影成像方法中使用的窗函数W，窗函数W表示为：
[0046]
[0047]【具体实施方式】四：
[0048]与【具体实施方式】三不同的是，本实施方式的基于幅度加权的探地雷达后向投影成 像方法，步骤四所述利用步骤三中构建的窗函数扫描成像点A，进行成像点A的成像过程具 体为，
[0049] 利用步骤三中构建的窗函数W对步骤二中成像点A在各通道的散射响应幅值XA,k构 成的Xa=(xa,i，xa,2,......，χα,μ)进行加权;然后将加权后的各通道回波相干叠加，获得成 像点Α的散射幅度信息，完成对成像点Α的成像过程。
【具体实施方式】 [0050] 五：
[0051] 与【具体实施方式】四不同的是，本实施方式的基于幅度加权的探地雷达后向投影成 像方法，所述的完成对成像点A的成像过程后，得到成像结果表示为：
[0052]
【具体实施方式】 [0053] 六：
[0054] 与【具体实施方式】四或五不同的是，本实施方式的基于幅度加权的探地雷达后向投 影成像方法，所述的进行对成像点A的成像过程中，
[0055] 当窗函数扫过非目标成像点时，对该成像点幅度加权值改变不大；
[0056] 当窗函数扫过目标成像点时，窗函数与目标成像点形成的幅度曲线重合，基于波 束形成和匹配滤波器的原理可知，此时对幅值的改变最大，从而达到了抑制旁瓣与改进分 辨力的目的，提高了信噪比。
[0057]仿真实验：
[0058] 计算成像点对应于每个合成孔径位置的延时;找到成像点对应于每个合成孔径位 置的散射强度;利用窗函数对散射强度进行加权后相干叠加获得该点的成像结果;遍历所 有成像区域，完成对整个区域的成像。
[0059] 其中，图2为经典后向投影成像方法对图1处理后的聚焦结果，图3为本发明涉及的 后向投影成像方法对图1处理后的聚焦结果。
[0060] 仿真结果证明：
[0061] 图3与图2相比可知。本发明的基于幅度加权的探地雷达后向投影成像方法，在经 典的后向投影成像方法的基础上对探地雷达回波信号进行加权，能够获得更好的信噪比和 成像分辨率，改善了经典后向投影成像算法对噪声、杂波干扰抑制能力较差的缺点，提高了 成像质量。本发明可应用于需要高分辨率成像的探地雷达成像领域，利用本文中所提出的 窗函数对散射幅度进行加权处理，改进了目标检测的精度，提高了对噪声、杂波干扰的抑制 能力。
[0062]本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域 技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于 本发明所附的权利要求的保护范围。
【主权项】
1. 一种基于幅度加权的探地雷达后向投影成像方法，其特征在于:所述方法通过W下 步骤实现： 步骤一、设定探地雷达的总通道数为M，计算成像区域内任意一个成像点A到各个通道 的双程时延，那么，成像点A到第k个通道的双程时延TA,k表示为：其中，成像点A坐标为A(xa，za); 式中，k表示通道序号，且k=l，2,......，M;T、k表示成像点A到第k个通道的双程时延； V表示电磁波在介质中的传播速度，.V - e / ，C表示电磁波在真空中传播的速度，大小为 3.0X108m/s，e为介质的介电常数； 步骤二、通过公式XA,k = rk(t)计算成像点A在各通道的散射响应幅值XA,k，其中t = TA,k; 式中，rk (t)表示探地雷达在各通道的回波数据； 步骤Ξ、构造窗函数W; 步骤四、利用步骤Ξ中构建的窗函数扫描成像点A，进行成像点A的成像过程； 步骤五、重复步骤一至步骤四的过程，遍历成像区域所有成像点，则完成了成像过程。2. 根据权利要求1所述基于幅度加权的探地雷达后向投影成像方法，其特征在于:步骤 二所述计算成像点A在各通道的散射响应幅值的过程为，首先，设探地雷达各通道的回波用 rk(t)表示，然后得成像点A在各通道的散射响应幅值XA,k分别为別,1< ='如），* =叫山式中^ 表示通道序号，且k=l，2,......，M。3. 根据权利要求1或2所述基于幅度加权的探地雷达后向投影成像方法，其特征在于： 步骤Ξ所述构造窗函数W的过程为， 步骤Ξ -、设双曲线的水平区域对应m个通道，分别为kl，k2, ...，km;对获得的探地雷达 B-scan图像进行边缘提取，得到双曲线上的点对应的响应幅值分别为m，U2,......，加，并 构造一个1XM维的矢量：其中，Μ表示总通道数，m表示探地雷达回波双曲线水平区域占据的通道数； 步骤Ξ二、设总通道数Μ和双曲线对应的通道数m都是偶数， 步骤ΞΞ、对获得的矢量W'进行归一化处理，则得到了基于幅度加权的探地雷达后向 投影成像方法中使用的窗函数W，窗函数W表示为：4. 根据权利要求3所述基于幅度加权的探地雷达后向投影成像方法，其特征在于:步骤 四所述利用步骤Ξ中构建的窗函数扫描成像点A，进行成像点A的成像过程具体为，利用步 骤^中构建的窗函数W对步骤二中成像点A在各通道的散射响应幅值XA,k构成的Xa=(xa,i， XA.2,......，XA,Μ)进行加权;然后将加权后的各通道回波相干叠加，获得成像点A的散射幅 度信息，完成对成像点A的成像过程。5. 根据权利要求4所述基于幅度加权的探地雷达后向投影成像方法，其特征在于:所述 的完成对成像点A的成像过程后，得到成像结果表示为：6. 根据权利要求4所述基于幅度加权的探地雷达后向投影成像方法，其特征在于:所述 的进行对成像点A的成像过程中， 当窗函数扫过非目标成像点时，对该成像点幅度加权值改变不大； 当窗函数扫过目标成像点时，窗函数与目标成像点形成的幅度曲线重合，对幅值的改 变最大。
【文档编号】G01S13/90GK105974405SQ201610289347
【公开日】2016年9月28日
【申请日】2016年5月4日
【发明人】李高鹏, 张明昊, 谢青青, 田文龙, 赵彬
【申请人】哈尔滨工业大学