CN110414182A - 引入天线方向图的探地雷达frtm算法 - Google Patents
引入天线方向图的探地雷达frtm算法 Download PDFInfo
- Publication number
- CN110414182A CN110414182A CN201910732391.9A CN201910732391A CN110414182A CN 110414182 A CN110414182 A CN 110414182A CN 201910732391 A CN201910732391 A CN 201910732391A CN 110414182 A CN110414182 A CN 110414182A
- Authority
- CN
- China
- Prior art keywords
- antenna
- frequency
- source
- imaging
- algorithm
- 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.)
- Granted
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 42
- 230000005855 radiation Effects 0.000 title claims abstract description 33
- 230000000149 penetrating effect Effects 0.000 title claims abstract description 28
- 238000003384 imaging method Methods 0.000 claims abstract description 65
- 230000005012 migration Effects 0.000 claims abstract description 33
- 238000013508 migration Methods 0.000 claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 21
- 238000010586 diagram Methods 0.000 claims description 43
- 238000004613 tight binding model Methods 0.000 claims description 31
- 238000005070 sampling Methods 0.000 claims description 21
- 238000001228 spectrum Methods 0.000 claims description 18
- 230000005684 electric field Effects 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000010287 polarization Effects 0.000 claims description 10
- 239000000523 sample Substances 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 6
- 230000005284 excitation Effects 0.000 claims description 6
- 238000000354 decomposition reaction Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000005672 electromagnetic field Effects 0.000 claims description 3
- 238000013213 extrapolation Methods 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 230000035939 shock Effects 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 6
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 230000004075 alteration 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
- 238000001514 detection method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Electromagnetism (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computing Systems (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了引入天线方向图的探地雷达FRTM算法,使得成像效果和精度获得进一步的改善。修正前的频域逆时偏移成像算法将辐射源视作完美点源,忽略了天线在目标介质中的能量辐射特性随辐射角度和介质参数的变化,本次方法在计算源波场和接收波场时分别点乘成像点到收发天线各自所成夹角的方向图函数进行修正,由修正后的源波场和接收波场得到最终的成像条件,并将所有炮的数据进行叠加进行成像。引入收发天线方向图函数修正后的频域逆时偏移成像算法使得成像效果和成像精度获得了进一步的提升。
Description
技术领域
本发明涉及探地雷达技术领域,具体为引入天线方向图的探地雷达FRTM算法。
背景技术
在中国的嫦娥五号月球探测任务中,MIMO探地雷达***工作在静止模式并会记录132道探地雷达数据。要求在3个小时内从探地雷达数据中解译出一个电大地下体积的高分辨率三维图像,作为重要的表土钻探任务的指导。为了大大提高计算效率,促进逆时偏移算法在月球探测以及其他工程领域的实际应用,近日相关学者提出了一种基于分层介质格林函数的频域逆时偏移算法。
在探地雷达进行数据采集的过程中,雷达发射天线向目标介质辐射高频电磁波脉冲信号,当电磁波遇到不同介质的分界面时会形成散射并最终由接收天线接收。天线向地下不同方向辐射的电磁波在不同方向上的能量分布一般是不均匀的,这种电磁波的能量辐射特性由天线的方向图来表征。在对雷达数据实施偏移成像时,辐射源通常用完美的理想点源辐射进行代替,然而实际上雷达***在目标介质中的辐射特性会随着天线型号和背景介质参数等因素发生变化。由于在一般的成像算法中,实际天线的辐射特性被忽略,从而对偏移成像算法的成像效果产生影响。
发明内容
本发明的目的在于提供引入天线方向图的探地雷达FRTM算法。
为实现上述目的,本发明提供如下技术方案:引入天线方向图的探地雷达FRTM算法,具体包括以下步骤:
步骤1、计算分层介质格林函数,格林函数是成像域中每个网格的每个频率的二阶张量;探地雷达通常在垂直宽边模式下工作,并使用线性极化天线进行发送和接收;这样,探地雷达仅记录电场的一个极化分量;因此,只需计算格林函数的一个分量,如gyy;而且分层介质格林函数具有对称性和水平平移不变性,因此有:
步骤2、对水平面上的格林函数进行切比雪夫插值,使得分层格林函数的计算进一步加速;
步骤3、利用分层介质格林函数的移不变性,完整成像空间各点的格林函数,并存入磁盘;
步骤4、计算源波场和接收波场,频域上空间的源波场和接收波场频谱可以通过与格林函数做简单的点乘运算获得:
其中,和分别为源波场频谱和接收波场频谱, 为激励源在发射天线位置和接收天线位置的并矢格林函数,和分别为接收天线激励信号源频谱和接收电磁场频谱的复共轭,e-jωT项由逆时外推产生;
步骤5、确定频域逆时偏移成像所需的最少频点数,在实施频率域逆时偏移成像时,频率点数往往需要人为决定,频率点数的多少一方面直接影响并矢格林函数的计算数量,另一方面也是我们在实际利用矢量网络分析仪采集数据的频率采样点数的设置依据;如果采样点数过多将增大并矢格林函数的计算时间和存储量,进一步的在引入天线方向图修正时会造成点源阵列综合方向图数量增多;如果采样点数过少将降低成像结果的精度并产生虚像;考虑到频率域与时间域的傅里叶变换关系,FRTM的频点数应至少满足采样定理要求,采样定理要求时间域的采样频率至少为信号最大频率的两倍才不会发生信号混叠从而保证信号不失真恢复,即:
fs≥2*fmax
其中,fs为采样频率,fmax为信号的最大频率。由时域和频域的对应关系可得:
其中,fmax和fmin为FRTM的最大频率和最小频率,N为频率逆时偏移使用的频点数,Tmax为成像的最大时窗。因此,频点数N需满足以下关系式:
N≥2*Tmax(fmax-fmin)+1
步骤6、获取实际雷达***中天线的方向图,真实天线的能量辐射与电磁波的频率和介质的性质有关,最后由接收天线接收的能量信号更是连续频点下能量辐射特性共同作用下的结果。
优选的,获取实际雷达***中天线的方向图的方法为三种,具体如下:
方法一、建立实际雷达***中天线的模型,利用电磁仿真工具进行模拟,得到天线在目标介质下的方向图;
方法二、将辐射源视作一个无限长的线源,在实际的探地雷达应用中,电磁波传播环境可以等效为分层均匀介质,天线近似位于空气和地下介质的交界处,无限长线源在两层介质交界处(即半空间环境)辐射电磁波;探地雷达常用线极化天线作为收发天线,包括偶极子天线和各种偶极子天线的变形如蝶形天线、Vivaldi天线等,这些天线的H面辐射方向图与线源方向图相近,因此线源半空间的方向图是实际天线半空间的方向图的一种良好的近似;利用无限长线源在半空间模型中方向图函数的远场解析解对偏移算法进行修正。
优选的,利用无限长线源在半空间模型中方向图函数的远场解析解对偏移算法进行修正,具体步骤为:
步骤1:考虑沿z轴放置的线源,令y=0为介质层分解面,空气层位于y>0半空间,介电常数为ε的介质层位于y<0半空间,线源可表示为
其中,为z方向的单位矢量,I为总电流,δ(x)和δ(y)为狄拉克冲击函数;
步骤2:为了表述的方便将坐标系转换到柱坐标系下,根据坐标转换公式:由对称结构可知,电场只有z分量磁场包含ρ,分量由麦克斯韦方程可得
由和可得,
其中,k2=ω2εμ=n2k0 2,k为传播常数;
步骤3:根据傅里叶积分变换关系:
将傅里叶积分变换关系带入上式得:
上式的解存在固定形式,由辐射边界条件和电场连续性并带入原方程可得上下半空间电场Ez1和Ez2分别为:
步骤4:在远场条件下,即k0ρ→∞,应用固定相位法可以求得积分解,由此而得到无限长线源上下半空间的方向图为:
其中为折射率,为临界角;
方法三、通过实测的方式得到雷达***中天线的方向图的方法。
优选的,通过实测的方式得到雷达***中天线的方向图的方法的具体步骤为:
步骤1:在目标介质内部埋设探头;
步骤2:使用雷达***中的天线在实际的工作高度对覆盖探头的一定区域进行足够密度的采样;
步骤3:记录下每个采样位置接收到的信号能量,拟合出天线的方向图。
步骤4:未经优化的频域逆时偏移成像算法忽略了实际情况下天线的能量辐射特性随目标与收发天线相对位置的影响,所以引入天线的方向图对频域逆时偏移成想算法进行修正;
步骤5:由修正后的源波场和接收波场,得到最终的成像条件为:
xT为雷达***的发射天线位置,xR为雷达***的接收天线位置,Us(x,z,ω)和Ur(x,z,ω)分别是源波场和接收波场的频谱。
步骤6:最后将所有的数据叠加即为最终逆时偏移所成的像。
优选的,引入天线的方向图对频域逆时偏移成想算法进行修正的具体做法是在计算源波场和接收波场时分别点乘散射点到收发天线各自所成夹角的方向图函数进行修正:
其中,fT(x,z,ω)和fR(x,z,ω)分别为发射天线到成像点的入射角对应的方向图幅度和成像点到接收天线的出射角对应的方向图函数,Gs(x,z,ω)和Gr(x,z,ω)分别是成像点相对于收发天线位置下的格林函数。
与现有技术相比,本发明的有益效果是:在偏移成像算法中引入天线辐射方向图对偏移成像算法进行修正,有助于提高地下目标,尤其是近地表复杂介质中小目标体的成像精度。使得成像效果和精度获得进一步的改善。修正前的频域逆时偏移成像算法将辐射源视作完美点源,忽略了天线在目标介质中的能量辐射特性随辐射角度和介质参数的变化,本次方法在计算源波场和接收波场时分别点乘成像点到收发天线各自所成夹角的方向图函数进行修正,由修正后的源波场和接收波场得到最终的成像条件,并将所有炮的数据进行叠加进行成像。引入收发天线方向图函数修正后的频域逆时偏移成像算法使得成像效果和成像精度获得了进一步的提升
附图说明
图1引入天线方向图修正后的频域逆时偏移成像算法流程图;
图2考虑实际天线辐射下频域逆时偏移算法原理示意图。
具体实施方式
本发明提供的一种实施例:
对频域逆时偏移成像算法,算法原理图如图2所示,xT为雷达***的发射天线位置,xR为雷达***的接收天线位置,fT(x,z,ω)和fR(x,z,ω)分别为发射天线和接收天线的方向图函数,Us(x,z,ω)和Ur(x,z,ω)分别是源波场和接收波场的频谱,Gs(x,z,ω)和Gr(x,z,ω)分别是成像点相对于收发天线位置下的格林函数。
基于分层介质格林函数的频域逆时偏移成像算法(FRTM)继承了时域逆时偏移成像算法高精度和高成像精度的特点的同时将计算量和成像时间降低了一个数量级,本次发明将天线方向图引入频域逆时偏移成像算法使得成像效果得到进一步的改善,具体包括以下步骤:
步骤1、计算分层介质格林函数,格林函数是成像域中每个网格的每个频率的二阶张量。探地雷达通常在垂直宽边模式下工作,并使用线性极化天线进行发送和接收。这样,探地雷达仅记录电场的一个极化分量。因此,只需计算格林函数的一个分量,如gyy。而且分层介质格林函数具有对称性和水平平移不变性,因此有,
步骤2、对水平面上的格林函数进行切比雪夫插值,使得分层格林函数的计算进一步加速。
步骤3、利用分层介质格林函数的移不变性,完整成像空间各点的格林函数,并存入磁盘。
步骤4、计算源波场和接收波场,频域上空间的源波场和接收波场频谱可以通过与格林函数做简单的点乘运算获得:
其中,和分别为源波场频谱和接收波场频谱, 为激励源在发射天线位置和接收天线位置的并矢格林函数,和分别为接收天线激励信号源频谱和接收电磁场频谱的复共轭,e-jωT项由逆时外推产生。
步骤5、确定频域逆时偏移成像所需的最少频点数,在实施频率域逆时偏移成像时,频率点数往往需要人为决定,频率点数的多少一方面直接影响并矢格林函数的计算数量,另一方面也是我们在实际利用矢量网络分析仪采集数据的频率采样点数的设置依据。如果采样点数过多将增大并矢格林函数的计算时间和存储量,进一步的在引入天线方向图修正时会造成点源阵列综合方向图数量增多;如果采样点数过少将降低成像结果的精度并产生虚像。考虑到频率域与时间域的傅里叶变换关系,FRTM的频点数应至少满足采样定理要求,采样定理要求时间域的采样频率至少为信号最大频率的两倍才不会发生信号混叠从而保证信号不失真恢复,即:
fs≥2*fmax
其中,fs为采样频率,fmax为信号的最大频率。由时域和频域的对应关系可得:
其中,fmax和fmin为FRTM的最大频率和最小频率,N为频率逆时偏移使用的频点数,Tmax为成像的最大时窗。因此,频点数N需满足以下关系式:
N≥2*Tmax(fmax-fmin)+1
步骤6、获取实际雷达***中天线的方向图,真实天线的能量辐射与电磁波的频率和介质的性质有关,最后由接收天线接收的能量信号更是连续频点下能量辐射特性共同作用下的结果。下面给出三种获取天线方向图的方法:
方法一、建立实际雷达***中天线的模型,利用电磁仿真工具进行模拟,得到天线在目标介质下的方向图;
方法二、将辐射源视作一个无限长的线源,在实际的探地雷达应用中,电磁波传播环境可以等效为分层均匀介质,天线近似位于空气和地下介质的交界处,无限长线源在两层介质交界处(即半空间环境)辐射电磁波。探地雷达常用线极化天线作为收发天线,包括偶极子天线和各种偶极子天线的变形如蝶形天线、Vivaldi天线等,这些天线的H面辐射方向图与线源方向图相近,因此线源半空间的方向图是实际天线半空间的方向图的一种良好的近似。利用无限长线源在半空间模型中方向图函数的远场解析解对偏移算法进行修正,具体步骤为:
步骤1:考虑沿z轴放置的线源,令y=0为介质层分解面,空气层位于y>0半空间,介电常数为ε的介质层位于y<0半空间,线源可表示为
其中,为z方向的单位矢量,I为总电流,δ(x)和δ(y)为狄拉克冲击函数;
步骤2:为了表述的方便将坐标系转换到柱坐标系下,根据坐标转换公式:由对称结构可知,电场只有z分量磁场包含ρ,分量由麦克斯韦方程可得
由和可得,
其中,k2=ω2εμ=n2k0 2,k为传播常数;
步骤3:根据傅里叶积分变换关系:
将傅里叶积分变换关系带入上式得:
上式的解存在固定形式,由辐射边界条件和电场连续性并带入原方程可得上下半空间电场Ez1和Ez2分别为:
步骤4:在远场条件下,即k0ρ→∞,应用固定相位法可以求得积分解,由此而得到无限长线源上下半空间的方向图为:
其中为折射率,为临界角;
方法三、通过实测的方式得到雷达***中天线的方向图的方法,具体步骤为:
步骤1:在目标介质内部埋设探头;
步骤2:使用雷达***中的天线在实际的工作高度对覆盖探头的一定区域进行足够密度的采样;
步骤3:记录下每个采样位置接收到的信号能量,拟合出天线的方向图。
步骤7、未经优化的频域逆时偏移成像算法忽略了实际情况下天线的能量辐射特性随目标与收发天线相对位置的影响,所以引入天线的方向图对频域逆时偏移成想算法进行修正。具体做法是在计算源波场和接收波场时分别点乘散射点到收发天线各自所成夹角的方向图函数进行修正:
其中,fT(x,z,ω)和fR(x,z,ω)分别为发射天线到成像点的入射角对应的方向图幅度和成像点到接收天线的出射角对应的方向图函数。
步骤8、由修正后的源波场和接收波场,得到最终的成像条件为:
步骤9、最后将所有炮的数据叠加即为最终逆时偏移所成的像。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (5)
1.引入天线方向图的探地雷达FRTM算法,其特征在于,具体包括以下步骤:
步骤1、计算分层介质格林函数,格林函数是成像域中每个网格的每个频率的二阶张量;探地雷达通常在垂直宽边模式下工作,并使用线性极化天线进行发送和接收;这样,探地雷达仅记录电场的一个极化分量;因此,只需计算格林函数的一个分量,gyy;而且分层介质格林函数具有对称性和水平平移不变性,因此有:
步骤2、对水平面上的格林函数进行切比雪夫插值,使得分层格林函数的计算进一步加速;
步骤3、利用分层介质格林函数的移不变性,完整成像空间各点的格林函数,并存入磁盘;
步骤4、计算源波场和接收波场,频域上空间的源波场和接收波场频谱可以通过与格林函数做简单的点乘运算获得:
其中,和分别为源波场频谱和接收波场频谱, 为激励源在发射天线位置和接收天线位置的并矢格林函数,和分别为接收天线激励信号源频谱和接收电磁场频谱的复共轭,e-jωT项由逆时外推产生;
步骤5、确定频域逆时偏移成像所需的最少频点数,在实施频率域逆时偏移成像时,频率点数往往需要人为决定,频率点数的多少一方面直接影响并矢格林函数的计算数量,另一方面也是我们在实际利用矢量网络分析仪采集数据的频率采样点数的设置依据;如果采样点数过多将增大并矢格林函数的计算时间和存储量,进一步的在引入天线方向图修正时会造成点源阵列综合方向图数量增多;如果采样点数过少将降低成像结果的精度并产生虚像;考虑到频率域与时间域的傅里叶变换关系,FRTM的频点数应至少满足采样定理要求,采样定理要求时间域的采样频率至少为信号最大频率的两倍才不会发生信号混叠从而保证信号不失真恢复,即:
fs≥2*fmax
其中,fs为采样频率,fmax为信号的最大频率。由时域和频域的对应关系可得:
其中,fmax和fmin为FRTM的最大频率和最小频率,N为频率逆时偏移使用的频点数,Tmax为成像的最大时窗。因此,频点数N需满足以下关系式:
N≥2*Tmax(fmax-fmin)+1
步骤6、获取实际雷达***中天线的方向图,真实天线的能量辐射与电磁波的频率和介质的性质有关,最后由接收天线接收的能量信号更是连续频点下能量辐射特性共同作用下的结果。
2.根据权利要求1所述的引入天线方向图的探地雷达FRTM算法,其特征在于,获取实际雷达***中天线的方向图的方法为三种,具体如下:
方法一、建立实际雷达***中天线的模型,利用电磁仿真工具进行模拟,得到天线在目标介质下的方向图;
方法二、将辐射源视作一个无限长的线源,在实际的探地雷达应用中,电磁波传播环境可以等效为分层均匀介质,天线近似位于空气和地下介质的交界处,无限长线源在两层介质交界处(即半空间环境)辐射电磁波;探地雷达常用线极化天线作为收发天线,包括偶极子天线和各种偶极子天线的变形如蝶形天线、Vivaldi天线等,这些天线的H面辐射方向图与线源方向图相近,因此线源半空间的方向图是实际天线半空间的方向图的一种良好的近似;利用无限长线源在半空间模型中方向图函数的远场解析解对偏移算法进行修正。
3.根据权利要求2所述的引入天线方向图的探地雷达FRTM算法,其特征在于,利用无限长线源在半空间模型中方向图函数的远场解析解对偏移算法进行修正,具体步骤为:
步骤1:考虑沿z轴放置的线源,令y=0为介质层分解面,空气层位于y>0半空间,介电常数为ε的介质层位于y<0半空间,线源可表示为
其中,为z方向的单位矢量,I为总电流,δ(x)和δ(y)为狄拉克冲击函数;
步骤2:为了表述的方便将坐标系转换到柱坐标系下,根据坐标转换公式:由对称结构可知,电场只有z分量磁场包含ρ,分量由麦克斯韦方程可得
由和可得,
其中,k2=ω2εμ=n2k0 2,k为传播常数;
步骤3:根据傅里叶积分变换关系:
将傅里叶积分变换关系带入上式得:
上式的解存在固定形式,由辐射边界条件和电场连续性并带入原方程可得上下半空间电场Ez1和Ez2分别为:
步骤4:在远场条件下,即k0ρ→∞,应用固定相位法可以求得积分解,由此而得到无限长线源上下半空间的方向图为:
其中为折射率,为临界角;
方法三、通过实测的方式得到雷达***中天线的方向图的方法。
4.根据权利要求3所述的引入天线方向图的探地雷达FRTM算法,其特征在于,通过实测的方式得到雷达***中天线的方向图的方法的具体步骤为:
步骤1:在目标介质内部埋设探头;
步骤2:使用雷达***中的天线在实际的工作高度对覆盖探头的一定区域进行足够密度的采样;
步骤3:记录下每个采样位置接收到的信号能量,拟合出天线的方向图。
步骤4:未经优化的频域逆时偏移成像算法忽略了实际情况下天线的能量辐射特性随目标与收发天线相对位置的影响,所以引入天线的方向图对频域逆时偏移成想算法进行修正;
步骤5:由修正后的源波场和接收波场,得到最终的成像条件为:
Us(x,z,ω)和Ur(x,z,ω)分别是源波场和接收波场的频谱。
步骤6:最后将所有的数据叠加即为最终逆时偏移所成的像。
5.根据权利要求4所述的引入天线方向图的探地雷达FRTM算法,其特征在于,引入天线的方向图对频域逆时偏移成想算法进行修正的具体做法是在计算源波场和接收波场时分别点乘散射点到收发天线各自所成夹角的方向图函数进行修正:
其中,fT(x,z,ω)和fR(x,z,ω)分别为发射天线到成像点的入射角对应的方向图幅度和成像点到接收天线的出射角对应的方向图函数,Gs(x,z,ω)和Gr(x,z,ω)分别是成像点相对于收发天线位置下的格林函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910732391.9A CN110414182B (zh) | 2019-08-09 | 2019-08-09 | 引入天线方向图的探地雷达frtm算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910732391.9A CN110414182B (zh) | 2019-08-09 | 2019-08-09 | 引入天线方向图的探地雷达frtm算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110414182A true CN110414182A (zh) | 2019-11-05 |
CN110414182B CN110414182B (zh) | 2021-10-22 |
Family
ID=68366580
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910732391.9A Expired - Fee Related CN110414182B (zh) | 2019-08-09 | 2019-08-09 | 引入天线方向图的探地雷达frtm算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110414182B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112946373A (zh) * | 2021-02-01 | 2021-06-11 | 北京邮电大学 | 一种基于紧缩场***的无相位测量方法及装置 |
CN113868852A (zh) * | 2021-09-22 | 2021-12-31 | 中国人民解放军63660部队 | 一种开阔场地电磁波近地面传播特性快速获取方法 |
CN113887103A (zh) * | 2021-09-30 | 2022-01-04 | 北京智芯仿真科技有限公司 | 基于不同介质特性的集成电路全波电磁仿真方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7511654B1 (en) * | 2006-01-12 | 2009-03-31 | The United States Of America As Represented By The Secretary Of The Army | Systems and methods for mine detection |
CN104535999A (zh) * | 2015-01-02 | 2015-04-22 | 中国人民解放军国防科学技术大学 | 一种校正天线方向图影响的雷达成像数据预处理方法 |
CN105137408A (zh) * | 2015-07-23 | 2015-12-09 | 电子科技大学 | 一种最优天线方向图选取的雷达角超分辨方法 |
CN107390181A (zh) * | 2017-07-26 | 2017-11-24 | 电子科技大学 | 一种基于多波束扫描的雷达高分辨成像方法 |
CN108845317A (zh) * | 2018-06-07 | 2018-11-20 | 广州大学 | 一种基于分层介质格林函数的频域逆时偏移算法 |
-
2019
- 2019-08-09 CN CN201910732391.9A patent/CN110414182B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7511654B1 (en) * | 2006-01-12 | 2009-03-31 | The United States Of America As Represented By The Secretary Of The Army | Systems and methods for mine detection |
CN104535999A (zh) * | 2015-01-02 | 2015-04-22 | 中国人民解放军国防科学技术大学 | 一种校正天线方向图影响的雷达成像数据预处理方法 |
CN105137408A (zh) * | 2015-07-23 | 2015-12-09 | 电子科技大学 | 一种最优天线方向图选取的雷达角超分辨方法 |
CN107390181A (zh) * | 2017-07-26 | 2017-11-24 | 电子科技大学 | 一种基于多波束扫描的雷达高分辨成像方法 |
CN108845317A (zh) * | 2018-06-07 | 2018-11-20 | 广州大学 | 一种基于分层介质格林函数的频域逆时偏移算法 |
Non-Patent Citations (4)
Title |
---|
EDEMSKIY F D 等: "Qualitative model of spatio-temporal radiation pattern of GPR antenna", 《GROUND PENETRATING RADAR (GPR), 2012 14TH INTERNATIONAL CONFERENCE》 * |
HAI LIU等: "Frequency Domain Reverse Time Migration of GPR with Attenenuation Compensation", 《2018 17TH INTERNATIONAL CONFERENCE ON GROUND PENETRATING RADAR (GPR)》 * |
孙延鹏 等: "基于OFDM的探地雷达压缩感知成像方法", 《沈阳航空航天大学学报》 * |
龙志军: "探地雷达逆时偏移成像及地下目标智能识别算法", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112946373A (zh) * | 2021-02-01 | 2021-06-11 | 北京邮电大学 | 一种基于紧缩场***的无相位测量方法及装置 |
CN112946373B (zh) * | 2021-02-01 | 2024-02-09 | 北京邮电大学 | 一种基于紧缩场***的无相位测量方法及装置 |
CN113868852A (zh) * | 2021-09-22 | 2021-12-31 | 中国人民解放军63660部队 | 一种开阔场地电磁波近地面传播特性快速获取方法 |
CN113868852B (zh) * | 2021-09-22 | 2024-04-19 | 中国人民解放军63660部队 | 一种开阔场地电磁波近地面传播特性快速获取方法 |
CN113887103A (zh) * | 2021-09-30 | 2022-01-04 | 北京智芯仿真科技有限公司 | 基于不同介质特性的集成电路全波电磁仿真方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN110414182B (zh) | 2021-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Annan | Electromagnetic principles of ground penetrating radar | |
CN110414182B (zh) | 引入天线方向图的探地雷达frtm算法 | |
CN107783092A (zh) | 基于链条关系式的近场后向rcs测量***及方法 | |
Alvarez et al. | An inverse fast multipole method for geometry reconstruction using scattered field information | |
Liu et al. | Migration of ground penetrating radar with antenna radiation pattern correction | |
Maruddani et al. | The development of ground penetrating radar (GPR) data processing | |
CN113534140B (zh) | 基于波场互相关的探地雷达三维成像方法 | |
Valle et al. | Ground penetrating radar antennas: theoretical and experimental directivity functions | |
CN110297237B (zh) | 考虑天线方向图的探地雷达绕射叠加成像方法及*** | |
Wang et al. | GPR imaging using the generalized Radon transform | |
CN117406190A (zh) | 基于雷达信号的无开挖拉棒腐蚀检测方法、装置及设备 | |
Liu et al. | Buried targets detection from synthetic anc measured B-scan ground penetrating radar data | |
Van der Kruk et al. | Improved three-dimensional image reconstruction technique for multi-component ground penetrating radar data | |
Góes et al. | Refraction Effect in SAR Processing for Focused Subsurface Tomography | |
Ding et al. | Near-field phase cross correlation focusing imaging and parameter estimation for penetrating radar | |
Bai et al. | Weakly convex discontinuity adaptive regularization for microwave imaging | |
Biancheri-Astier et al. | Modeling the configuration of HF electrical antennas for deep bistatic subsurface sounding | |
Serrano et al. | Alternative cost function for full waveform inversion of GPR data | |
Grinev et al. | Multi-channel Ground Penetrating Radar based on ultra-wideband short-pulse signal: Hardware and software | |
Zheglova et al. | 3D radio-frequency electromagnetic wave modelling with an optimized frequency-domain finite-difference scheme | |
Guzel | A Fast Matched Filtered Method for Ground Penetrating Radar Tomographic Imaging | |
Lambot et al. | Near-field modeling of radar antennas for wave propagation in layered media: when models represent reality | |
Bevacqua et al. | Virtual experiments and compressive sensing for subsurface microwave tomography | |
Sharma et al. | FFT implementation of Kirchhoff's migration for ground penetrating radar image focussing | |
Estatico et al. | Variable Exponent Lebesgue Space Inversion for Cross-Borehole Subsurface Imaging |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20211022 |