CN113534140B - 基于波场互相关的探地雷达三维成像方法 - Google Patents
基于波场互相关的探地雷达三维成像方法 Download PDFInfo
- Publication number
- CN113534140B CN113534140B CN202110745399.6A CN202110745399A CN113534140B CN 113534140 B CN113534140 B CN 113534140B CN 202110745399 A CN202110745399 A CN 202110745399A CN 113534140 B CN113534140 B CN 113534140B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- penetrating radar
- ground penetrating
- fourier transform
- imaging
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 79
- 230000000149 penetrating effect Effects 0.000 title claims abstract description 62
- 238000004613 tight binding model Methods 0.000 claims abstract description 26
- 230000009466 transformation Effects 0.000 claims abstract description 17
- 238000000034 method Methods 0.000 claims description 33
- 238000001228 spectrum Methods 0.000 claims description 11
- 238000009825 accumulation Methods 0.000 claims description 6
- 230000001427 coherent effect Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000021615 conjugation Effects 0.000 claims 1
- 235000013399 edible fruits Nutrition 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 45
- 238000012545 processing Methods 0.000 abstract description 11
- 238000013519 translation Methods 0.000 abstract description 2
- 238000001514 detection method Methods 0.000 description 13
- 238000013508 migration Methods 0.000 description 12
- 230000005012 migration Effects 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 9
- 239000011295 pitch Substances 0.000 description 9
- 238000004088 simulation Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000005284 excitation Effects 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000005290 field theory Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提出了一种基于波场互相关的探地雷达三维成像方法，能够实现对地下场景进行快速高分辨三维成像处理，计算复杂度低，满足工程实际需求。本发明探地雷达***采用共中心模式采集回波数据，通过沿俯仰方位步进移动收发天线对得到成像场景C扫描回波数据省去了后续在时间维作傅里叶变换，降低了算法复杂度。考虑到分层介质格林函数的平移不变性，得到收发天线位置分层介质格林函数共轭相乘结果，然后进行二维傅里叶变换，降低了算法复杂度。
Description
技术领域
本发明涉及探地雷达技术领域，具体涉及一种基于波场互相关的探地雷达三维成像方法。
背景技术
探地雷达(Ground Penetrating Radar)技术是一种较新的地球物理探测方法，该技术根据高频电磁脉冲波在地下不同物性介质之间得反射、透射和绕射等来探测地下介质分布规律，在国内外已得到广泛的应用。相比于电阻率法、低频电磁感应法等常规无损探测方法，探地雷达技术具有探测效率高，可连续测量，分辨率高，操作方便灵活，探测费用低等优越性，已经成为目前浅表层目标探测的热门技术之一。目前，探地雷达技术的应用领域已经远远超出了“探地”的范畴，作为一种重要的无损探测方法，在隐藏危险物探查、工程质量检测、环境监测、地质勘查、考古以及军事反恐等领域中发挥着重要的作用。
由于GPR探测的根本目标是为了对地下目标或场景进行检测和识别，而最初人们直接使用一维或二维回波来完成检测或识别工作，这种工作方式既不方便直观，又对操作人员的专业技能和先验知识提出了较高层次的要求。为了降低目标检测和识别的难度，GPR成像技术应运而生并不断发展。通过改变探地雷达回波中能量的分布情况，GPR成像技术可以更为直观地表现地下场景及其中目标的几何或物理特征，便于目标解译工作的开展。由于GPR成像技术可使目标解译变得更加容易和准确，已成为当今非破坏性对地探测领域中应用最广泛、最具发展前景的技术之一。
根据成像过程中处理回波数据时所在域的不同，可将GPR成像算法分为时域算法和频域算法两类，其中最典型和最具代表性的算法有时域的后向投影(Back Projection,BP)成像算法和频域的距离偏移(Range Migration,RM)成像算法。BP算法起源于计算机辅助层析技术，是SAR成像中原理最直观简单的一种成像算法。由于原理类似，BP算法可直接应用到GPR成像之中，而且可以克服GPR分层介质带来的不利影响。比如[Halman J I,Shubert K A,Ruck G T.SAR processing of ground-penetrating radar data forburied UXO detection:results from a surface-based system[J].IEEE Transactionson Antennas and Propagation,1998,46(7):1023～1027.]使用BP成像算法对地下未***武器进行了三维成像，结果证明了该算法在对地下目标定位和识别时的有效性。Tanyer使用时域有限差分法对非均匀介质中的BP成像问题，通过在能量求和中引入加权步骤，改善了BP算法的成像精度。[雷文太,刘立业,粟毅.平面分层媒质中目标的TAM-BP成像算法[J].信号处理,2007,23(5):680～685.]在考虑折射和衰减效应的基础上，增加了对近场条件下球面波的传播衰减的考虑，使得BP算法可以适应有耗媒质场景。从已有算法可以发现：BP算法在采样孔径不均匀、运动误差较大时鲁棒性较强，可以获得比频域算法更优的效果，但是其自身具有运算速度慢、旁瓣和杂波能量高等缺点。RM算法起源于地震信号处理，[Stolt RH.Migration by transform[J].Geophysics,1978,43(1):23～48.]首先提出了RM算法。由于原理类似，基于波方程的RM算法同样也适合于GPR成像场景为均匀介质时的情况。比如[Mast J E.Microwave pulse-echo radar imaging for the nondestructiveevaluation of civil structures[D].Urbana,Illinois:University of Illinois atUrbana-Champaign,1993.]将SAR成像中的RM算法引入GPR成像中，由标量波动方程计算出电磁波在各项均匀介质中的传播函数，提出了适合于GPR成像的RM算法；[Gu K,Wang G,LiJ.Migration based SAR imaging for ground penetrating radar systems[J].IEEProceedings of Radar Sonar Navigation,2004,151(5):317～325.]于2004年讨论了使用RM算法对下视和前视GPR进行成像的方法，通过空间补零的方法克服了假频问题。但是，以上的研究都是基于波速不变模型进行推导的，即总是假设成像场景中只存在一种介质，而在实际GPR应用时，受人员安全或设备工作环境的考虑，天线总是距离地面一定的高度，从而使得上述基于单一场景模型的RM算法失效。
逆时偏移(Reverse-time Migration,RTM)最早在地震勘探领域提出，由于其成像精度高等优势逐渐被应用于探地雷达成像[Leuschen C J,Plumb R G.A matched-filter-based reverse-time migration algorithm for ground-penetrating radar data[J].IEEE Trans.Geo.Remote Sens,2001,39(5):929-936.]。探地雷达逆时偏移的基本原理是将地表记录到的接收点电磁波场在时间轴上进行逆向传播，当电磁波场逆推至零时刻，则所有反射波和绕射波的能量都回到最初被反射和绕射的空间位置，然后应用成像条件可获得最终的偏移剖面。逆时偏移基于双程波动方程，且不存在对波动方程的近似，允许电磁波在任意方向进行传播，能够对任意类型的电磁波场进行成像，因此逆时偏移方法相比其他偏移方法具有无法比拟的优越性。但是传统逆时偏移算法需要多次的正向和反向波场外推，然后应用成像条件将每个炮点的成像结果进行累加。因此传统时域逆时偏移需要很大的计算与存储资源，导致其并不能在实际场景中对大数据量进行实时成像处理，故并未得到广泛应用。
BP和RM成像算法往往只适用于远场环境成像，在近场条件下其性能将大幅恶化。传统时域RTM算法计算复杂度过高。在地探地雷达工程实践中，使用场景非常复杂，采集数据往往不均匀。
发明内容
有鉴于此，本发明提出了一种基于波场互相关的探地雷达三维成像方法，能够实现对地下场景进行快速高分辨三维成像处理，计算复杂度低，满足工程实际需求。
为实现上述目的，本发明的技术方案为：
本发明的一种基于波场互相关的探地雷达三维成像方法，包括如下步骤：
步骤1，将探地雷达***C扫描回波数据在俯仰方位时间维作三维傅里叶变换，得到三维傅里叶变换后的探地雷达扫描回波；
根据成像场景先验信息，计算探地雷达数据采集平面中心点处收发天线分层介质格林函数并共轭相乘，得到共轭相乘结果，然后对共轭相乘结果在俯仰、方位维作傅里叶变换，得到二维傅里叶变换后的共轭相乘结果；
步骤2，将探地雷达发射信号频谱、所述三维傅里叶变换后的探地雷达扫描回波以及所述二维傅里叶变换后的共轭相乘结果相乘，得到三者相乘结果；
将所述三者相乘结果作俯仰以及方位逆傅里叶变换得到各单频点的三维复图像；
将各单频点的三维复图像相参累加得到最终成像结果。
其中，所述步骤1中，所述探地雷达***采用共中心模式采集回波数据，通过沿俯仰方位步进移动收发天线对得到成像场景C扫描回波数据R([x,y,z＝0],t)，其中[x,y,z＝0]表示收发天线对中点坐标，t表示时间；当探地雷达***为频率步进体制，则接收信号为R([x,y,z＝0],ω)，ω表示角频率。
其中，所述步骤1中，对采集得到的C扫描回波数据R([x,y,z＝0],t)在俯仰方位时间维作三维傅里叶变换的具体方式为：
其中为采集C扫描回波数据俯仰方位时间维作三维傅里叶变换的结果。
其中，所述步骤1中，所述分层介质格林函数为G(r,r0,ω)，其中r＝(x,y,z)为三维成像场景的坐标，r0＝(x0,y0,z＝0)为探地雷达数据采集平面中心点坐标；
利用分层介质格林函数的平移不变性，rs＝(xs,ys,z＝0)位置处的格林函数分布为G(r,rs,ω)＝G(r-rs,r0,ω)；将收发天线位置分层介质格林函数共轭相乘的表达式为：M(r,r0,ω)＝G(r,r0,ω)G*(r+rd,r0,ω)
其中*表示复共轭，rd表示收发天线对间距；
将M(r,r0,ω)作俯仰方位二维傅里叶变换的表达式为：
其中为收发天线位置分层介质格林函数共轭相乘后作二维傅里叶变换结果。
其中，所述步骤2中，探地雷达数据采集中心位置处分层介质格林函数共轭相乘结果与作三维傅里叶变换后的探地雷达扫描回波、发射波形频谱三者相乘的表达式为：
其中S(ω)为探地雷达发射信号s(t)的频谱，B([kx,ky,z],ω)为三维复图像作俯仰方位傅里叶变化结果。
其中，对B([kx,ky,z],ω)作俯仰方位逆傅里叶变换，得到地下场景各频点三维复图像，所述表达式为：
其中为地下场景各频点的三维复图像。
其中，各频点复图像相参累加得到最终成像结果的表达式为：
其中，I(x,y,z)为最终地下场景三维成像结果。
其中，所述步骤1中，采用非均匀傅里叶变换将回波数据变换到均匀波数域，然后再进行三维傅里叶变换。
其中，所述步骤1中，所述成像场景先验信息包括被探测地下介质的相对介电常数以及天线口面离地面高度。
其中，所述步骤1中，在俯仰、方位维作傅里叶变换的计算根据成像场景先验信息提前计算预存在磁盘中。
有益效果：
本发明考虑到分层介质格林函数的平移不变性，得到收发天线位置分层介质格林函数共轭相乘结果，然后进行二维傅里叶变换，降低了算法复杂度。本发明成像分辨率高，算法稳健性强；适用于远近场场景以及分层介质场景成像；算法复杂度低，收敛性稳定，能对探地雷达实测数据实时处理。
本发明探地雷达***采用共中心模式采集回波数据，通过沿俯仰方位步进移动收发天线对得到成像场景C扫描回波数据省去了后续在时间维作傅里叶变换，进一步降低了算法复杂度。
本发明采用非均匀傅里叶变换将回波数据变换到均匀波数域，只需精确记录收发天线位置的坐标即可，进一步降低了算法复杂度。
本发明中，二维傅里叶变换后的共轭相乘结果的计算根据成像场景先验信息提前计算预存在磁盘中，大幅提升所提算法数据处理效率，满足工程实际需求。
附图说明
图1是本发明的探地雷达***数据采集方式示意图
图2是本发明的算法流程图；
图3是本发明探地雷达仿真回波数据的三维场景示意图；
图4是本发明图3三维模型剖面图；
图5是本发明对仿真场景回波数据三维成像结果。
具体实施方式
下面结合附图并举实施例，对本发明进行详细描述。
本发明提出了一种基于波场互相关的探地雷达三维成像方法，技术原理为：
本发明使用共偏移距模式采集GPR数据，即收发天线之间距离恒定。基于波场零延时互相关成像条件，探地雷达成像结果可以表示为：
其中r＝(x,y,z)为成像网格点坐标，对于二维垂直截面和水平切片成像则网格坐标分别取r＝(x,z)，r＝(x,y)；rd＝(xd,yd,zd＝0)是收发天线之间的坐标差值，并且天线仅在xoy平面以rd步进移动；ur(·)与us(·)分别是收发天线位置的后向与前向递推波场。对每个天线收发位置的数据进行时域逆时偏移并进行累加便得到最终成像结果I(r)。
根据卷积定理：
其中Us(·)与Ur(·)分别为前向与后向时域波场值ur(·)、us(·)的频域形式。根据电磁场理论，波场分布是分层介质格林函数与相应收发位置处激励源频谱的乘积，即：
Us(r,rs,ω)＝G(r,rs,ω)S(ω)
其中G(·)频域分层介质格林函数，上标*表示复共轭，S(ω)为探地雷达发射信号s(t)的频谱，为接收信号时间反转后的傅里叶换。
考虑到分层介质格林函数的平移不变性，即：
G(r,rs,ω)＝G(r-rs,r0,ω)
其中r0＝(0,0,0)为数据采集平面的中心坐标，G(r,r0,ω)为对应坐标处的频域分层介质格林函数。则：
其中，M(r,r0,ω)＝G(r,r0,ω)G*(r+rd,r0,ω)
在某些远场探测场景中，收发天线间距很小，式(12)可简化为M(r,r0,ω)＝G(r,r0,ω)G*(r,r0,ω)
根据卷积定理：
其中，分别为M(r,r0,ω)、/>对x、y方向作二维傅里叶变换的结果。I(r)为最终三维成像结果。
本发明算法流程如图2所示，本发明包括如下步骤：
步骤1，将探地雷达***C扫描回波数据在俯仰方位时间维作三维傅里叶变换，得到三维傅里叶变换后的探地雷达扫描回波；
所述探地雷达***采用共中心模式采集回波数据，数据采集方式示意图如图1所示。通过沿俯仰方位步进移动收发天线对得到成像场景C扫描回波数据R([x,y,z＝0],t)，其中[x,y,z＝0]表示收发天线对中点坐标，t表示时间；当探地雷达***为频率步进体制，则接收信号为R([x,y,z＝0],ω)，ω表示角频率，省去了后续在时间维作傅里叶变换；
对采集得到的C扫描回波数据R([x,y,z＝0],t)在俯仰方位时间维作三维傅里叶变换的表达式为：
其中为采集C扫描回波数据俯仰方位时间维作三维傅里叶变换的结果，在实际数据采集过程中，数据采集点往往非均匀分布，可以采用非均匀傅里叶变换将回波数据变换到均匀波数域，只需精确记录收发天线位置的坐标即可。
根据成像场景先验信息，计算探地雷达数据采集平面中心点处收发天线分层介质格林函数并共轭相乘，得到共轭相乘结果，然后对共轭相乘结果在俯仰、方位维作傅里叶变换，得到二维傅里叶变换后的共轭相乘结果。
本实施例中，根据被探测地下介质的相对介电常数，以及天线口面离地面高度等成像场景的先验信息，计算探地雷达数据采集平面中心点处分层介质格林函数G(r,r0,ω)，其中r＝(x,y,z)为三维成像场景的坐标，r0＝(x0,y0,z＝0)为探地雷达数据采集平面中心点坐标；
利用分层介质格林函数的平移不变性，rs＝(xs,ys,z＝0)位置处的格林函数分布为G(r,rs,ω)＝G(r-rs,r0,ω)；将收发天线位置分层介质格林函数共轭相乘的表达式为：
M(r,r0,ω)＝G(r,r0,ω)G*(r+rd,r0,ω)
其中*表示复共轭，rd表示收发天线对间距；
将M(r,r0,ω)作俯仰方位二维傅里叶变换的表达式为：
其中为收发天线位置分层介质格林函数共轭相乘后作二维傅里叶变换结果，该步骤/>的计算可以根据成像场景先验信息提前计算预存在磁盘中，大幅提升所提算法数据处理效率，满足工程实际需求。
步骤2，将探地雷达发射信号频谱、所述三维傅里叶变换后的探地雷达扫描回波以及所述二维傅里叶变换后的共轭相乘结果相乘，得到三者相乘结果；具体地，探地雷达数据采集中心位置处分层介质格林函数共轭相乘结果与作三维傅里叶变换后的探地雷达扫描回波、发射波形频谱三者相乘的表达式为：
其中S(ω)为探地雷达发射信号s(t)的频谱，B([kx,ky,z],ω)为三维复图像作俯仰方位傅里叶变化结果；
将所述三者相乘结果作俯仰以及方位逆傅里叶变换得到各单频点的三维复图像；对B([kx,ky,z],ω)作俯仰方位逆傅里叶变换，得到地下场景各频点三维复图像，所述表达式为：
其中为地下场景各频点的三维复图像。
将各单频点的三维复图像相参累加得到最终成像结果，各频点复图像相参累加得到最终成像结果的表达式为：
其中，I(x,y,z)为最终地下场景三维成像结果。
本发明的效果可以通过以下仿真结果进一步说明：
仿真条件：采用探地雷达正演软件gprMax来模拟真实场景回波数据，激励信号为中心频率为1.5GHz，-3dB带宽约为1.3GHz的Ricker脉，仿真场景为2m*2m*0.7m的三维半空间模型，地下介质的相对介电常数为6，相对磁导率为1，电导率为0.001S/m。
仿真实验一：图3给出了仿真三维场景示意图，目标为间距不同的3对半径为3cm的金属圆柱，其具体参数如图4所示。
图5为所提算法对仿真场景回波数据进行处理得到的三维成像结果，结果证明所提方法具有突出的地下目标重构能力，很好还原了地下目标的几何形态与物理尺寸，能够为地下目标探测提供技术支撑。此外由于所提算法基于波场传播特性进行推导，能够很好适用于近场环境成像。通过与BP三维成像算法进行比较，所提算法具有更低的算法复杂度，适用于工程实际中探地雷达数据的实时成像处理。
综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于波场互相关的探地雷达三维成像方法，其特征在于，包括如下步骤：
步骤1，将探地雷达***C扫描回波数据在俯仰方位时间维作三维傅里叶变换，得到三维傅里叶变换后的探地雷达扫描回波；
根据成像场景先验信息，计算探地雷达数据采集平面中心点处收发天线分层介质格林函数并共轭相乘，得到共轭相乘结果，然后对共轭相乘结果在俯仰、方位维作傅里叶变换，得到二维傅里叶变换后的共轭相乘结果；
步骤2，将探地雷达发射信号频谱、所述三维傅里叶变换后的探地雷达扫描回波以及所述二维傅里叶变换后的共轭相乘结果相乘，得到三者相乘结果；
将所述三者相乘结果作俯仰以及方位逆傅里叶变换得到各单频点的三维复图像；
将各单频点的三维复图像相参累加得到最终成像结果。
2.如权利要求1所述的方法，其特征在于，所述步骤1中，所述探地雷达***采用共中心模式采集回波数据，通过沿俯仰方位步进移动收发天线对得到成像场景C扫描回波数据R([x,y,z＝0],t)，其中[x,y,z＝0]表示收发天线对中点坐标，t表示时间；当探地雷达***为频率步进体制，则接收信号为R([x,y,z＝0],ω)，ω表示角频率。
3.如权利要求2所述的方法，其特征在于，所述步骤1中，对采集得到的C扫描回波数据R([x,y,z＝0],t)在俯仰方位时间维作三维傅里叶变换的具体方式为：
其中为采集C扫描回波数据俯仰方位时间维作三维傅里叶变换的结果。
4.如权利要求2或3所述的方法，其特征在于，所述步骤1中，所述分层介质格林函数为G(r,r0,ω)，其中r＝(x,y,z)为三维成像场景的坐标，r0＝(x0,y0,z＝0)为探地雷达数据采集平面中心点坐标；
利用分层介质格林函数的平移不变性，rs＝(xs,ys,z＝0)位置处的格林函数分布为G(r,rs,ω)＝G(r-rs,r0,ω)；将收发天线位置分层介质格林函数共轭相乘的表达式为：
M(r,r0,ω)＝G(r,r0,ω)G*(r+rd,r0,ω)
其中*表示复共轭，rd表示收发天线对间距；
将M(r,r0,ω)作俯仰方位二维傅里叶变换的表达式为：
其中为收发天线位置分层介质格林函数共轭相乘后作二维傅里叶变换结果。
5.如权利要求4所述的方法，其特征在于，所述步骤2中，探地雷达数据采集中心位置处分层介质格林函数共轭相乘结果与作三维傅里叶变换后的探地雷达扫描回波、发射波形频谱三者相乘的表达式为：
其中S(ω)为探地雷达发射信号s(t)的频谱，B([kx,ky,z],ω)为三维复图像作俯仰方位傅里叶变化结果。
6.如权利要求5所述的方法，其特征在于，对B([kx,ky,z],ω)作俯仰方位逆傅里叶变换，得到地下场景各频点三维复图像，所述表达式为：
其中为地下场景各频点的三维复图像。
7.如权利要求6所述的方法，其特征在于，各频点复图像相参累加得到最终成像结果的表达式为：
其中，I(x,y,z)为最终地下场景三维成像结果。
8.如权利要求1、2、3、5、6和7任意一项所述的方法，其特征在于，所述步骤1中，采用非均匀傅里叶变换将回波数据变换到均匀波数域，然后再进行三维傅里叶变换。
9.如权利要求1、2、3、5、6和7任意一项所述的方法，其特征在于，所述步骤1中，所述成像场景先验信息包括被探测地下介质的相对介电常数以及天线口面离地面高度。
10.如权利要求1、2、3、5、6和7任意一项所述的方法，其特征在于，所述步骤1中，所述二维傅里叶变换后的共轭相乘结果的计算根据成像场景先验信息提前计算预存在磁盘中。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110745399.6A CN113534140B (zh) | 2021-07-01 | 2021-07-01 | 基于波场互相关的探地雷达三维成像方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110745399.6A CN113534140B (zh) | 2021-07-01 | 2021-07-01 | 基于波场互相关的探地雷达三维成像方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113534140A CN113534140A (zh) | 2021-10-22 |
| CN113534140B true CN113534140B (zh) | 2024-02-20 |
Family
ID=78097628
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110745399.6A Active CN113534140B (zh) | 2021-07-01 | 2021-07-01 | 基于波场互相关的探地雷达三维成像方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113534140B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117392311B (zh) * | 2023-09-28 | 2024-06-25 | 北京化工大学 | 毁伤场景sar图像仿真方法及装置 |
| CN117554921B (zh) * | 2024-01-12 | 2024-03-29 | 西安中创云图科技有限公司 | 一种探地雷达三维场景正演方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5835054A (en) * | 1996-03-01 | 1998-11-10 | The Regents Of The University Of California | Ultra wideband ground penetrating radar imaging of heterogeneous solids |
| WO2002086542A2 (en) * | 2001-04-20 | 2002-10-31 | Witten Alan J | A method and apparatus for identifying buried objects using ground penetrating radar |
| CN102540184A (zh) * | 2011-12-30 | 2012-07-04 | 陕西延长石油(集团)有限责任公司研究院 | 探地雷达频域成像方法 |
| CN106291540A (zh) * | 2016-09-14 | 2017-01-04 | 河北省电力勘测设计研究院 | 一种基于doa估计的多输入多输出探地雷达逆向投影目标成像方法 |
| CN108845317A (zh) * | 2018-06-07 | 2018-11-20 | 广州大学 | 一种基于分层介质格林函数的频域逆时偏移算法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7841982B2 (en) * | 1995-06-22 | 2010-11-30 | Techniscan, Inc. | Apparatus and method for imaging objects with wavefields |
-
2021
- 2021-07-01 CN CN202110745399.6A patent/CN113534140B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5835054A (en) * | 1996-03-01 | 1998-11-10 | The Regents Of The University Of California | Ultra wideband ground penetrating radar imaging of heterogeneous solids |
| WO2002086542A2 (en) * | 2001-04-20 | 2002-10-31 | Witten Alan J | A method and apparatus for identifying buried objects using ground penetrating radar |
| CN102540184A (zh) * | 2011-12-30 | 2012-07-04 | 陕西延长石油(集团)有限责任公司研究院 | 探地雷达频域成像方法 |
| CN106291540A (zh) * | 2016-09-14 | 2017-01-04 | 河北省电力勘测设计研究院 | 一种基于doa估计的多输入多输出探地雷达逆向投影目标成像方法 |
| CN108845317A (zh) * | 2018-06-07 | 2018-11-20 | 广州大学 | 一种基于分层介质格林函数的频域逆时偏移算法 |
Non-Patent Citations (3)
| Title |
|---|
| Improved F-K Migration Based on Interpolation Method for GPR Imaging;zhongwei Yang 等;2019 6th Asia-Pacific Conference on Synthetic Aperture Radar (APSAR);1-6 * |
| 一种基于半空间扫描测量模式的隧道坍腔地质雷达三维成像技术;邓海明 等;现代隧道技术;第58卷(第3期);52-59 * |
| 基于波场互相关的探地雷达快速自聚焦成像;杨忠委 等;信号处理;第37卷(第9期);1663-1668 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113534140A (zh) | 2021-10-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8193967B2 (en) | Method and system for forming very low noise imagery using pixel classification | |
| Zhuge et al. | Modified Kirchhoff migration for UWB MIMO array-based radar imaging | |
| US7796829B2 (en) | Method and system for forming an image with enhanced contrast and/or reduced noise | |
| Özdemir et al. | A review on migration methods in B‐scan ground penetrating radar imaging | |
| González-Huici et al. | A comparative study of GPR reconstruction approaches for landmine detection | |
| CN113534140B (zh) | 基于波场互相关的探地雷达三维成像方法 | |
| Zhang | Two-dimensional microwave tomographic algorithm for radar imaging through multilayered media | |
| Pereira et al. | 3d tomographic image reconstruction for multistatic ground penetrating radar | |
| CN108761449B (zh) | 一种钢筋回波干扰下的灾害目标成像方法 | |
| Zhang et al. | Back-projection algorithm based on self-correlation for ground-penetrating radar imaging | |
| Meschino et al. | Spot-gpr: A freeware toolfor target detection and localizationin gpr data developedwithin the cost action tu1208 | |
| Osumi et al. | Detection of buried plant | |
| Liu et al. | Buried targets detection from synthetic anc measured B-scan ground penetrating radar data | |
| Hintz | SNR improvements in NIITEK ground-penetrating radar | |
| Yang et al. | Improved FK migration based on interpolation method for GPR imaging | |
| Almutiry | UAV tomographic synthetic aperture radar for landmine detection | |
| Agrawal et al. | GPR data analysis of weak signals using modified S-transform | |
| Meschino et al. | A practical guide on using SPOT-GPR, a freeware tool implementing a SAP-DoA technique | |
| Li et al. | MIMO borehole radar imaging based on high degree of freedom for efficient subsurface sensing | |
| CN110161500B (zh) | 一种基于Radon-Clean的改进圆周SAR三维成像方法 | |
| Ghaderi Aram et al. | Buried target imaging: A comparative study | |
| Rohman et al. | Basic Shape Classification of Buried Object Using Pattern Matching in Ultrawideband Radar Image | |
| Harkat et al. | Ground penetrating radar imaging for buried cavities in a dispersive medium: Profile reconstruction using a modified hough transform approach and a time-frequency analysis | |
| Song et al. | Ground-penetrating radar land mine imaging: Two-dimensional seismic migration and three-dimensional inverse scattering in layered media | |
| Qu et al. | Efficient back projection imaging approach for airborne GPR using NUFFT technique |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |