CN109283524B - 一种用于提高地质雷达信号分辨率的方法 - Google Patents
一种用于提高地质雷达信号分辨率的方法 Download PDFInfo
- Publication number
- CN109283524B CN109283524B CN201810868475.0A CN201810868475A CN109283524B CN 109283524 B CN109283524 B CN 109283524B CN 201810868475 A CN201810868475 A CN 201810868475A CN 109283524 B CN109283524 B CN 109283524B
- Authority
- CN
- China
- Prior art keywords
- frequency
- band
- scan signal
- scan
- harmonic
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000009466 transformation Effects 0.000 claims abstract description 14
- 238000003672 processing method Methods 0.000 claims abstract description 7
- 238000006243 chemical reaction Methods 0.000 claims abstract description 6
- 238000001228 spectrum Methods 0.000 claims description 42
- 238000004364 calculation method Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 11
- 238000010183 spectrum analysis Methods 0.000 abstract description 7
- 230000000694 effects Effects 0.000 description 6
- 230000003044 adaptive effect Effects 0.000 description 5
- 229910001294 Reinforcing steel Inorganic materials 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000007480 spreading Effects 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000002087 whitening effect Effects 0.000 description 1
Images
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本公开揭示了一种用于提高地质雷达信号分辨率的信号处理方法,所述方法在连续小波变换域内自适应的拓展地质雷达信号的频率,以得到高分辨率的拓频信号。首先采用基于地质雷达信号的频谱分析,确定各拓频频段和相应的加权系数,然后对地质雷达信号进行连续小波变换,得到相应的时间‑尺度(频率)分布,根据频率和小波尺度的换算关系,对相应尺度范围内的小波变换系数进行加权,最后基于加权的连续小波变换,重构对应的拓频信号,达到提高地质雷达信号分辨率的目的。所述方法能够为隧道衬砌和铁路路基检测提供高分辨率地质雷达数据,从而增加后续检测的准确性。
Description
技术领域
本发明属于地质雷达领域,特别涉及一种用于提高地质雷达信号分辨率的高分辨处理方法。
背景技术
利用地质雷达(Ground-penetrating Radar,GPR)进行探测时,发射天线发射电磁波在介质中传播,经过存在电性差异的介质体或者界面时会产生反射信号,该反射信号被设置于检测表层的接收天线接收。通过分析该接收信号实现对介质内部分布规律的分析。地质雷达探测具有超浅层勘探的独特优势,是一种快速的高分辨率无损检测方法,在隧道以及铁路路基检测中具有广泛的应用。但是,受到介质耗散等因素的影响,电磁波在传播过程中会出现频率降低,幅度衰减等问题,使得接收信号的分辨率降低,从而影响后续的异常检测。目前用于提高地质雷达信号分辨率的方法主要有
现有技术1:谱白化类方法
(1)计算每个A-scan的一维傅里叶变换,得到信号的振幅谱和相位谱;
(2)在信号的有效频宽范围内,设定若干个带通滤波器,将幅度谱分为若干个子带;
(3)计算每个子带内的信号振幅谱加权系数,将每个子带内的信号谱乘以对应的子带加权系数,得到拓频后的频谱;
(4)基于拓频后的频谱,计算其一维傅里叶反变换,得到拓频后的A-scan信号;
该类方法还可以在时间域进行能量均衡。
现有技术1的缺点:
1)该类方法在提高地质雷达数据分辨率的同时,也放大了信号中的随机噪声;
2)人为设定各拓频子带宽度,没有统一的设置规则,易造成不同设置参数的拓频效果差异较大;
现有技术2:反褶积方法
(1)提取地质雷达信号子波表达式;
对每一个反射时间点进行如下操作:
(2)基于子波表达式,构造自相关矩阵;
(3)对自相关矩阵进行奇异值分解,并计算加权系数向量;
(4)计算加权系数向量与观测时间序列的内积得到该时间点的反射强度;
上述为时间域反褶积处理流程,还可以根据卷积定理,将时间域反褶积流程推广到频率域进行反褶积处理。
现有技术2的缺点:
1)地质雷达子波信号难以准确提取,因此影响后续的反射强度计算;
2)通常地质雷达数据中含有较强的随机噪声,这也使得反褶积效果受到影响。
发明内容
基于此,本公开揭示了一种用于提高地质雷达信号分辨率的信号处理方法,其特征在于,该方法包括以下步骤:
S100、对多道A-scan信号进行一维傅里叶变换,并在A-scan信号的有效频带宽度内确定A-scan信号的基础频段和各拓频频段;
S200、计算A-scan信号各拓频频段的加权系数;
S300、对A-scan信号进行连续小波变换,得到对应的时间-尺度域小波变换系数(时频分布);
S400、根据频率和尺度的换算关系,将不同时间-尺度范围内的小波变换系数乘以对应的加权系数;
S500、计算乘以对应的加权系数后小波变换系数对应的小波反变换,得到该A-scan信号对应的高分辨率信号;
S600、对每个A-scan信号重复步骤S300-S500。
本公开具有如下技术效果:
1)本公开的方法不需要估计地质雷达信号子波;
2)本公开参数设置简单,根据信号频谱分布自适应确定频谱划分和加权系数;
3)本公开的连续小波域信号自适应拓频方法可以有效的提高地质雷达剖面的时间分辨率以及增强深层弱反射的强度;
4)通过对频带的有效划分和加权处理,本公开所述方法对随机噪声有很好的抑制作用。
附图说明:
图1为本公开一个实施例中的GPR数据自适应频谱拓展流程;
图2为本公开一个实施例中的频谱划分示意;
图3(A)-图3(F)为本公开一个实施例的处理结果:图3(A)为隧道衬砌地质雷达数据B-scan剖面;图3(B)为连续小波域自适应频谱拓展结果;图3(C)为抽取的水平距离x=1.44m的A-scan信号及其拓频结果;图3(D)为x=1.44m处A-scan信号及其拓频结果的傅里叶频谱分析;图3(E)为抽取的水平距离x=4.83m的A-scan信号及其拓频结果;图3(F)为x=4.83m处A-scan信号及其拓频结果的傅里叶频谱分析;
图4(A)-图4(F)为本公开一个实施例中的处理结果:图4(A)铁路路基地质雷达数据B-scan剖面;图4(B)连续小波域自适应频谱拓展结果;图4(C)抽取的水平距离x=2m的A-scan信号及其拓频结果;图4(D)x=2m处A-scan信号及其拓频结果的傅里叶频谱分析;图4(E) 抽取的水平距离x=4m的A-scan信号及其拓频结果;图4(F)x=4m 处A-scan信号及其拓频结果的傅里叶频谱分析。
具体实施方式:
下面结合附图1至附图4(F)和具体实施方式对本发明做进一步详细的说明。
本发明提供一种用于隧道衬砌和铁路路基地质雷达数据的频谱拓展方法,包括对地质雷达数据的频谱分析、小波域自适应拓频处理等。
在一个实施例中,本公开揭示了一种用于提高地质雷达信号分辨率的信号处理方法,所述方法包括以下步骤:
S100、对多道A-scan信号进行一维傅里叶变换,并在A-scan信号的有效频带宽度内确定A-scan信号的基础频段和各拓频频段;
S200、计算A-scan信号各拓频频段的加权系数;
S300、对A-scan信号进行连续小波变换,得到对应的时间-尺度域小波变换系数(时频分布);
S400、根据频率和尺度的换算关系,对不同时间-尺度范围内的小波变换系数乘以对应的加权系数;
S500、计算乘以对应的加权系数后小波变换系数对应的小波反变换,得到A-scan信号对应的高分辨率信号;
S600、对每个A-scan信号重复步骤S300-S500。
在本实施例中,所述方法基于连续小波变换的多分辨率思想,在连续小波变换域实现对信号的分频拓展处理,以达到提高分辨率的目的。所述方法不需要估计子波信号,对频带能够实现有效划分,而且具有很好的抗噪性能,能够为后续的地质雷达数据处理步骤提供高保真的高分辨率数据。
本实施例的连续小波域信号自适应频率拓展方法可以有效的提高地质雷达数据的时间分辨率以及弱反射的幅度,有利于后期的异常检测。相比较于现有技术1,本实施例提供的方法能够在提高数据分辨率的同时,不增加噪声能量,具有较高的信噪比,且根据倍频程规则实现频带的稳定划分;相比较于现有技术2,本公开提供的方法不需要估计子波信号,仅需要基于傅里叶振幅谱对信号的小波变换结果进行有效的划分即可。
在一个实施例中,所述步骤S100具体包括如下步骤:
S101、对每一道A-scan信号进行一维傅里叶变换,然后计算对应傅里叶振幅谱的平均振幅谱;
S102、基于对应的平均振幅谱选定A-scan信号的基准频率点,作为标准,根据基准频率点确定基础频率段,然后在频率轴上确定各拓展频段;
所述拓展频段包括:高频段的一次、二次和高次谐波段,和低频段的一次次谐波、二次次谐波以及高次次谐波频率段。
在本实施例中,对若干个A-scan进行一维傅里叶变换,将选取的若干A-scan的傅里叶变换的模值对应点相加,然后除以A-scan个数,得到傅里叶振幅谱的平均振幅谱;基于对应的平均振幅谱选定信号的基准频率点,作为标准,根据基准频率点确定基础频率段,然后在频率轴上确定高频段的一次、二次和高次谐波段,和低频段的一次次谐波、二次次谐波以及高次次谐波频率段。
在一个实施例中,步骤S200所述的计算各拓频频段加权系数具体为:
S201、计算信号基础频段内频谱的平均能量密度P为
其中,F(ω)为A-scan的傅里叶变换,ω0表示基准频率点,ω0表示基础频段的下界;
S202、按照步骤S201中的公式,依次计算高频段的一次谐波段、二次谐波段、低频段的一次次谐波段和二次次谐波段的平均能量密度 P1,P2,P1,P2;
在一个实施例中,计算步骤S300所述的连续小波变换的关系式如下:
在一个实施例中,步骤S400所述的频率和尺度的换算关系如下:
其中,a为小波尺度因子,ω为信号频率,ωb为母小波函数的频率。
在一个实施例中,步骤S500所述的连续小波反变换的关系式如下:
在一个实施例中,所述A-scan信号的基准频率点选取为A-scan 信号平均振幅谱的最大值对应的频率值。
在一个实施例中,向A-scan信号的基准频率点的低频段方向幅度下降约30%时对应的频率定义为基础频率段。
在一个实施例中,在频率轴上小于基础频段的频率定义为低频,高于基础频段的频率定义为高频,然后向高频段随着频率升高依次确定一次谐波段(ω0,ω1),二次谐波段(ω1,ω2)以及高次谐波段,向低频段随着频率降低依次确定一次次谐波段(ω1,ω0),二次次谐波段(ω2,ω1)。
在一个实施例中,所述L个A-scan信号平均振幅谱的计算公式为:
其中:Fl(ω),l=1,2,…,L表示L个A-scan信号的傅里叶变换。
在一个实施例中,如图1所示,本发明所提供的地质雷达数据处理方法,具体通过如下步骤实施:
准备探地雷达数据,进行如下处理:
步骤1:对若干A-scan信号进行一维傅里叶变换,并确定各拓频频段,主要包括:
(1)对若干个A-scanf(t)进行一维傅里叶变换,然后计算对应傅里叶振幅谱的平均振幅谱,一维傅里叶变换由以下关系式决定:
选取L个A-scan信号,分别计算其傅里叶变换Fl(ω),l=1,2,…,L,平均振幅谱由以下关系式决定:
(2)参照图2,基于对应的平均振幅谱选定信号的基准频率点ω0,作为标准,ω0通常选取为平均振幅谱最大值对应的频率值,之后向低频段方向幅度下降约30%时对应的频率定义为ω0,确定基础频率段 (ω0,ω0),在频率轴上小于基础频段的频率定义为低频,高于基础频段的频率定义为高频,然后向高频段随着频率升高依次确定一次谐波段(ω0,ω1),二次谐波段(ω1,ω2)以及高次谐波段,向低频段随着频率降低依次确定一次次谐波段(ω1,ω0),二次次谐波段(ω2,ω1)等;其中,待扩展的高频段的一次和二次谐波频率段频宽分别为基础频段的一倍和二倍大小,更高次谐波以此类推;待扩展的低频段的一次次谐波和二次次谐波频宽分别为基础频段的一半和四分之一大小;
步骤2:计算各拓频频段的加权系数,其确定方法为:
(1)计算信号基础频段内频谱的平均能量密度P为
其中,F(ω)为A-scan的傅里叶变换;
(2)按照(1)中的公式,依次计算一次谐波段、二次谐波段、低频段一次次谐波段和二次次谐波段的平均能量密度P1,P2,P1,P2;
其余更高次的谐波和次谐波以此类推;
对B-scan剖面上的每一个A-scan重复以下步骤:
步骤3:对A-scan f(t)进行连续小波变换,得到对应的时间-尺度域分布剖面;
步骤4:根据频率和尺度的换算关系,对各个尺度范围内的小波变换系数乘以对应的加权系数,具体包括:
(2)计算对应加权后的时间-尺度域分布,对于一次谐波段内的连续小波变换分布,则有:
其余各次谐波和次谐波以此类推;
步骤5:基于时频域扩展后的连续小波域时间-尺度域分布,计算其对应的连续小波反变换,重构相应的高分辨率信号,连续小波反变换由以下公式决定;
本发明具有如下有益效果:
1)本发明的方法不需要估计地质雷达信号子波;
2)本发明参数设置简单,根据信号频谱分布自适应确定频谱划分和加权系数;
3)本发明的连续小波域信号自适应拓频方法可以有效的提高地质雷达剖面的时间分辨率以及增强深层弱反射的强度;
4)通过对频带的有效划分和加权处理,本发明的方法对随机噪声有很好的抑制作用。
下面利用本发明提供的数据处理方法对(1)实测的隧道衬砌地质雷达B-scan数据和(2)实测的铁路路基地质雷达B-scan数据进行了处理;处理流程包括步骤1-步骤5,并进行了单个A-scan的波形图及频谱对比。
图3(A)-图3(F)为实测的隧道衬砌地质雷达B-scan剖面及处理结果:
图3(A)为实测的B-scan数据,发射天线的主频为900MHz,隧道管片厚度为35厘米,其中分布两层钢筋,主筋半径为0.01米,配筋半径为0.005米。从图3(A)上可以较为明显地看到浅层的强能量层状杂波以及管片中钢筋反射的双曲线波形。而由于实际介质的频散作用,管片与后方围岩的界面反射在图3(A)上较弱,很难用于对隧道衬砌异常进行有效检测。图3B为高分辨拓频处理结果,可以看出剖面的时间方向分辨率得到了有效提高,而且深层的反射波能量得到了有效增强。图3(C)-图3(F)为抽取的A-scan信号,其中黑色虚线对应未做拓频处理的A-scan信号,黑色实线为拓频后的 A-scan信号。图3(C)为水平距离x=1.44m的A-scan信号,对比可以看出,经过连续小波域自适应拓频后,多层反射波的振幅得到了有效增强,突出了更多的反射界面。图3(D)为对应图3(C)的频谱分析,可以看到,经过拓频处理后,在有效频带内信号的低频、高频成分得到了增强,随机噪声占据的高频频段并没有得到显著提高,这验证了拓频结果的高可信度。图3(E)为水平距离x=4.83m处未拓频和已拓频的A-scan,图3(E)为对应的傅里叶振幅谱,可以看出,低频、高频成分得到了提升,因此增加了有效频宽。
图4(A)-图4(F)为实测的铁路路基地质雷达B-scan剖面及处理结果:
图4(A)为实测的用于铁路路基检测的B-scan数据,其中发射天线的主频为900MHz,地下分布有网状的钢筋结构。图中可以看到钢筋网反射的双曲线结构以及层状结构。但是,由于介质衰减、频散等因素的影响,双曲线结构模糊,且与水平界面反射的信号相互干扰。利用本发明专利的方法,对图4(A)B-scan剖面上的每个A-scan 进行拓频处理,提高其时间分辨率。图4(B)为对应的拓频处理剖面,可以看到钢筋反射的双曲线结构得到了有效增强,10ns左右的反射界面得以显现。图4(C)-图4(F)为抽取的A-scan信号,其中黑色虚线对应未做拓频处理的A-scan信号,黑色实线为拓频后的 A-scan信号。对比可以看出,基于小波域的自适应拓频技术有效的拓展了地质雷达信号的带宽,且将噪声限制在一定水平。这为后面的路基病害检测提供了有效的处理结果。
最后需要说明的是,以上模型和实际资料算例对本发明的目的,技术方案以及有益效果提供了进一步的验证,这仅属于本发明的具体实施算例,并不用于限定本发明的保护范围,在本发明的精神和原则之内,所做的任何修改,改进或等同替换等,均应在本发明的保护范围内。
Claims (9)
1.一种用于提高地质雷达信号分辨率的信号处理方法,其特征在于,所述方法包括以下步骤:
S100、对多道A-scan信号进行一维傅里叶变换,并在A-scan信号的有效频带宽度内确定A-scan信号的基础频段和各拓频频段;
S200、计算A-scan信号各拓频频段的加权系数;
S300、对A-scan信号进行连续小波变换,得到对应的时间-尺度域小波变换系数;
S400、根据频率和尺度的换算关系,对不同时间-尺度范围内的小波变换系数乘以对应的加权系数;
S500、计算所述乘以对应的加权系数后小波变换系数对应的连续小波反变换,得到A-scan信号对应的高分辨率信号;
S600、对每个A-scan信号重复步骤S300-S500;
其中,步骤S200所述计算A-scan信号各拓频频段加权系数包括以下步骤:
S201、计算A-scan信号基础频段内频谱的平均能量密度P,
其中,|F(ω)|为A-scan的平均振幅谱,ω0表示基准频率点,ω0表示基础频段的下界;
S202、按照步骤S201中的公式,依次计算A-scan信号高频段的一次谐波段、二次谐波段,低频段的一次次谐波段和二次次谐波段的平均能量密度P1,P2,P1,P2;
2.根据权利要求1所述的方法,其特征在于,所述步骤S100包括如下步骤:
S101、对多道A-scan信号进行一维傅里叶变换,然后计算对应傅里叶振幅谱的平均值,即为A-scan信号的平均振幅谱;
S102、基于对应的平均振幅谱选定A-scan信号的基准频率点,根据基准频率点确定基础频率段,然后在频率轴上确定各拓频频段;
所述拓频频段包括:高频段的一次、二次和高次谐波段,低频段的一次次谐波、二次次谐波以及高次次谐波频率段。
6.根据权利要求2所述的方法,其特征在于:所述A-scan信号的基准频率点ω0选取为A-scan信号平均振幅谱的最大值对应的频率值。
7.根据权利要求6所述的方法,其特征在于:在A-scan信号的平均振幅谱上,向A-scan信号的基准频率点的低频段方向幅度下降约30%时对应的频率定义为基础频率段的下界ω0。
8.根据权利要求7所述的方法,其特征在于:在频率轴上小于基础频段的频率定义为低频,高于基础频段的频率定义为高频,然后向高频段随着频率升高依次确定一次谐波段(ω0,ω1),二次谐波段(ω1,ω2)以及高次谐波段,向低频段随着频率降低依次确定一次次谐波段(ω1,ω0),二次次谐波段(ω2,ω1)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810868475.0A CN109283524B (zh) | 2018-08-01 | 2018-08-01 | 一种用于提高地质雷达信号分辨率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810868475.0A CN109283524B (zh) | 2018-08-01 | 2018-08-01 | 一种用于提高地质雷达信号分辨率的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109283524A CN109283524A (zh) | 2019-01-29 |
CN109283524B true CN109283524B (zh) | 2020-11-06 |
Family
ID=65183327
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810868475.0A Expired - Fee Related CN109283524B (zh) | 2018-08-01 | 2018-08-01 | 一种用于提高地质雷达信号分辨率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109283524B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113688692A (zh) * | 2021-07-29 | 2021-11-23 | 浙江大学 | 基于探地雷达时频尺度特征的供水管道漏损检测方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103485265A (zh) * | 2013-09-27 | 2014-01-01 | 华南理工大学 | 一种超宽带探地雷达公路质量探测方法及其探测装置 |
US20160097670A1 (en) * | 2014-10-01 | 2016-04-07 | Honeywell International Inc. | Resolution mode switching for pulsed radar |
CN104793203B (zh) * | 2015-04-17 | 2017-03-29 | 中南大学 | 一种用于多频多通道探地雷达的数据融合方法 |
CN105353373B (zh) * | 2015-12-16 | 2018-04-20 | 武汉大学 | 一种基于Hough变换探地雷达目标提取方法和装置 |
CN206020656U (zh) * | 2016-07-04 | 2017-03-15 | 杭州欧镭激光技术有限公司 | 一种用于激光扫描雷达的收发装置 |
CN106443771B (zh) * | 2016-08-05 | 2018-12-11 | 中国石油集团东方地球物理勘探有限责任公司 | 提高转换波地震数据分辨率的方法及其速度反演方法 |
CN107436450A (zh) * | 2017-07-26 | 2017-12-05 | 西安交通大学 | 一种基于连续小波变换的地震信号带宽拓展方法 |
-
2018
- 2018-08-01 CN CN201810868475.0A patent/CN109283524B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN109283524A (zh) | 2019-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Baili et al. | GPR signal de-noising by discrete wavelet transform | |
Wang et al. | Noise suppressing and direct wave arrivals removal in GPR data based on Shearlet transform | |
CN105334508B (zh) | 一种稀疏阵列宽带波束形成的栅瓣抑制方法 | |
CN104360251B (zh) | 一种变压器局部放电的超声波信号时延估计方法 | |
Rodriguez et al. | New analysis and extensions of split-spectrum processing algorithms | |
CN106842163B (zh) | 一种弹道目标回波信号时频特性估计方法 | |
CN106199532A (zh) | 基于混合傅立叶‑小波分析的探地雷达信号降噪方法 | |
Zhang et al. | Intensive interferences processing for GPR signal based on the wavelet transform and FK filtering | |
CN112285703B (zh) | 一种海杂波抑制及目标检测方法 | |
CN113887398A (zh) | 一种基于变分模态分解和奇异谱分析的gpr信号去噪方法 | |
CN109283524B (zh) | 一种用于提高地质雷达信号分辨率的方法 | |
Hao et al. | An automated GPR signal denoising scheme based on mode decomposition and principal component analysis | |
CN107121705B (zh) | 一种基于自动反相校正和峰度值比较的探地雷达回波信号去噪算法 | |
CN109283523B (zh) | 一种地质雷达B-scan数据处理方法 | |
Xin et al. | ISAR imaging of target with complex motion associated with the fractional Fourier transform | |
Karslı | Further improvement of temporal resolution of seismic data by autoregressive (AR) spectral extrapolation | |
Akbarpour et al. | Clutter removal of near-field UWB SAR imaging for pipeline penetrating radar | |
CN109975867B (zh) | 一种频率域信号混叠的地震数据拓频方法 | |
Perov et al. | Wavelet filtering of signals from ultrasonic flaw detector | |
CN110749923A (zh) | 一种基于范数方程提高分辨率的反褶积方法 | |
CN101930605A (zh) | 基于二维混合变换的sar图像目标提取方法及系统 | |
CN101639530B (zh) | 一种基于二维混合变换的sar回波信号去噪预处理方法 | |
CN104569945A (zh) | 一种用于合成孔径雷达图像的自适应滤波方法 | |
Muja et al. | Assessment of RF interference mitigation methods for automotive radars using real data | |
Yun et al. | Accurate Three-Dimensional Scattering Center Extraction for ISAR Image Using the Matched Filter-Based CLEAN Algorithm |
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: 20201106 |