CN105527617B - 一种基于鲁棒主元分析的探地雷达数据背景去除方法 - Google Patents
一种基于鲁棒主元分析的探地雷达数据背景去除方法 Download PDFInfo
- Publication number
- CN105527617B CN105527617B CN201610084826.XA CN201610084826A CN105527617B CN 105527617 B CN105527617 B CN 105527617B CN 201610084826 A CN201610084826 A CN 201610084826A CN 105527617 B CN105527617 B CN 105527617B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- data
- ground
- scan data
- 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
- 238000000513 principal component analysis Methods 0.000 title claims abstract description 31
- 230000001427 coherent effect Effects 0.000 title 1
- 238000000034 method Methods 0.000 claims abstract description 44
- 230000000149 penetrating effect Effects 0.000 claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims description 55
- 238000005457 optimization Methods 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 2
- 238000012847 principal component analysis method Methods 0.000 abstract description 15
- 238000001914 filtration Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation 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
- 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
- G01S7/414—Discriminating targets with respect to background clutter
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种基于鲁棒主元分析的探地雷达数据背景去除方法，涉及探地雷达技术领域，尤其涉及一种基于鲁棒主元分析的探地雷达数据背景去除方法。本发明要解决传统背景去除方法对探地雷达数据背景信号估计不准确、鲁棒性差的问题。本发明按以下步骤进行：一、探地雷达B‑SCAN数据标准化；二、探地雷达A‑SCAN数据去除直流偏移；三、利用滑动平均滤波法去除探地雷达B‑SCAN图像随机噪声；四、鲁棒主元分析去除背景。本发明采用鲁棒主元分析法取得良好效果，鲁棒性较强，本发明可应用于探地雷达技术领域。
Description
技术领域
本发明涉及探地雷达技术领域，尤其涉及一种基于鲁棒主元分析的探地雷达数据背景去除方法。
背景技术
对于收发天线共置的冲激脉冲探地雷达，随着雷达沿着水平测线移动，并不断向地下发射宽带电磁波，同时接收地下目标体的反射信号，在每一个位置可以得到一组采样数据，称为A-SCAN，由一系列A-SCAN数据组成二维探地雷达回波剖面图，称为B-SCAN图像。所以在进行后续的目标检测、识别以及解译工作之前，应采用有效的预处理方法，对探地雷达数据进行预处理，用以压制噪声和杂波等背景的干扰，尽可能突出目标信号。
经典的背景去除方法如主元分析方法，是将探地雷达B-SCAN数据矩阵进行奇异值分解，由于接收信号中，地面杂波分量较强，所以大奇异值对应的分量为主元，这些选定的主元所合成的信号即是地面杂波的估计值。这种方法的关键问题是如何选取主元。选择不当，则会产生较大误差，鲁棒性差。
考虑到探地雷达数据中背景杂波在水平方向通常是几乎不变的，比如当地表较为平整、介质较为均匀时，地表强反射波的时间延迟基本形同，故可以认为，探地雷达所采集的B-SCAN数据中，杂波分量可以用一个低秩矩阵来近似，而目标回波数据为少数的异常点，因此可以用一个稀疏矩阵来近似，因此应用鲁棒主元分析法能够有效的将背景杂波和目标回波分离，完成背景杂波的去除。
发明内容
本发明的基本原理：
首先对原始探地雷达B-SCAN数据进行标准化，由于原始数据可能存在直流偏移量，故将每一道A-SCAN数据减去自身均值；之后应用滑动平均滤波去除随机噪声，这样会使原始数据中的“毛刺”变得平滑；由于探地雷达数据中背景杂波在水平方向通常是变化缓慢的，比如当地表较为平整、介质较为均匀时，地表强反射波的时间延迟基本形同，故可以认为，杂波分量可以用一个低秩矩阵来近似，而目标回波数据为少数的异常点，因此可以用一个稀疏矩阵来近似，因此应用鲁棒主元分析法能够有效的将背景杂波和目标回波分离，完成背景杂波的去除。
本发明的目的是为了解决传统背景去除方法对探地雷达数据背景信号估计不准确、鲁棒性差的问题，提供一种基于鲁棒主元分析的探地雷达数据背景去除方法。本发明所述一种基于鲁棒主元分析的探地雷达数据背景去除方法，按以下步骤进行：
步骤一：探地雷达B-SCAN数据标准化
将探地雷达B-SCAN数据标准化处理，使二维数据矩阵中的值处于[-1,1]区间；
步骤二：探地雷达A-SCAN数据去除直流偏移；
探地雷达回波数据常常出现全部为正或者全部为负的情况，说明含有直流分量，需要进行压制，将每一道A-SCAN数据减去该道A-SCAN数据的均值；
步骤三：利用滑动平均滤波法去除探地雷达B-SCAN图像随机噪声；
利用一个滑动窗口，选择窗口宽度为5道，计算窗口内的5道A-SCAN数据的均值，用来代替窗口内上述5道中的中间道的A-SCAN数据；
步骤四：鲁棒主元分析去除背景
对于经过上述步骤处理之后的B-SCAN数据，进行鲁棒主元分析，将其分解成一个低秩矩阵L与一个稀疏矩阵S之和，其中低秩矩阵L作为背景杂波的估计矩阵，稀疏矩阵S作为探地雷达回波信号的估计矩阵。
本发明包括以下有益效果：
1、利用鲁棒主元分析法，相较于传统的主元分析法，避免因主元选取不当引起的较大误差；
2、传统的主元分析法，当杂波较强时，求解的结果会出现较大的误差，而鲁棒主元分析法仍能取得良好效果，鲁棒性较强；
3、本发明中，鲁棒主元分析法采用非精确拉格朗日乘子法，收敛速度快，运算精度较高，不需要占用过多的存储空间。
附图说明
图1为一种基于鲁棒主元分析的探地雷达数据背景去除方法总流程图；
图2为鲁棒主元分析算法流程图；
图3为原始B-SCAN数据灰度图；
图4为原始单道A-SCAN数据图；
图5为滑动平均滤波前后对比图；
图6为鲁棒主元分析处理后得到的稀疏矩阵灰度图；
图7为鲁棒主元分析处理后得到的低秩矩阵灰度图；
图8为稀疏矩阵中的单道A-SCAN数据图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合图1至8和具体实施方式对本发明作进一步详细的说明，其中图1为本发明所述的基于鲁棒主元分析的探地雷达数据背景去除方法总流程图。
具体实施方式一、本实施方式所述的一种基于鲁棒主元分析的探地雷达数据背景去除方法，按以下步骤进行：
步骤一：探地雷达B-SCAN数据标准化
将探地雷达B-SCAN数据标准化处理，使二维数据矩阵中的值处于[-1,1]区间；
步骤二：探地雷达A-SCAN数据去除直流偏移
探地雷达回波数据常常出现全部为正或者全部为负的情况，说明含有直流分量，需要进行压制，将每一道A-SCAN数据减去所有A-SCAN数据的均值；
步骤三：利用滑动平均滤波法去除探地雷达B-SCAN图像随机噪声
利用一个滑动窗口，选择窗口宽度为5道，计算窗口内的5道A-SCAN数据的均值，用来代替窗口内上述5道中的中间道的A-SCAN数据；
步骤四：鲁棒主元分析去除背景
对于经过上述步骤处理之后的B-SCAN数据，进行鲁棒主元分析，将其分解成一个低秩矩阵L与一个稀疏矩阵S之和，其中低秩矩阵L作为背景杂波的估计矩阵，稀疏矩阵S作为探地雷达回波信号的估计矩阵。
本实施方式包括以下有益效果：
2、利用鲁棒主元分析法，相较于传统的主元分析法，避免因主元选取不当引起的较大误差；
2、传统的主元分析法，当杂波较强时，求解的结果会出现较大的误差，而鲁棒主元分析法仍能取得良好效果，鲁棒性较强；
3、本发明中，鲁棒主元分析法采用非精确拉格朗日乘子法，收敛速度快，运算精度较高，不需要占用过多的存储空间。
具体实施方式二、本实施方式是对具体实施方式一所述的一种基于鲁棒主元分析的探地雷达数据背景去除方法的进一步说明，步骤一所述的将探地雷达B-SCAN数据标准化处理具体过程如下：
B-SCAN数据D为m×n的矩阵，且所有元素的值处于区间[MIN,MAX]，那么做以下标准化处理：
标准化后的B-SCAN数据取值范围变为[-1,1]。
具体实施方式三、本实施方式是对具体实施方式一或二所述的一种基于鲁棒主元分析的探地雷达数据背景去除方法的进一步说明，步骤二所述探地雷达A-SCAN数据去除直流偏移具体过程如下：
探地雷达B-SCAN数据中的每一列都称之为一道A-SCAN，将B-SCAN数据D写成A-SCAN的集合：
D＝[A1,A2,...,An]
其中Ai＝[ai1,ai2,...,aim]T表示i第道A-SCAN数据。对每一道A-SCAN数据做如下处理：
得到去除了直流分量的数据。
具体实施方式四、本实施方式是对具体实施方式一至三之一所述的一种基于鲁棒主元分析的探地雷达数据背景去除方法的进一步说明，步骤三所述利用滑动平均滤波法去除探地雷达B-SCAN图像随机噪声具体过程如下：
将窗口内的数据以列为单位取均值代替当前列，可以消除图像中的高频噪声干扰；算法设计如下：
对于m×n的数据矩阵，选择滑动窗口宽度为5道，随着窗口的滑动，对窗口内的5道A-SCAN求均值，代替当前窗口中间道的A-SCAN数据；公式如下：
对于边缘的A-SCAN数据，最前两道数据作如下处理：
A'1＝A1
最后两道数据，做相同处理：
A'n＝An
对于中间道数据，做如下处理：
得到的处理后的数据为D'＝[A'1,A'2,...,A'n]。
具体实施方式五、本实施方式是对具体实施方式一至四之一所述的一种基于鲁棒主元分析的探地雷达数据背景去除方法的进一步说明，步骤四所述鲁棒主元分析去除背景具体过程为：
(1)建立鲁棒主元分析模型
其中，D表示待处理的探地雷达B-SCAN数据矩阵，L和S分别表示待求解的低秩矩阵和稀疏矩阵，λ表示折中因子，rank(·)表示矩阵的秩，||·||0表示l0范数，表示使关于L和S的函数取得最小值；采用凸松弛法转化成如下公式：
其中，用核范数||·||*替代秩，l1范数||·||1代替l0范数；
(2)应用拉格朗日乘子法，构造拉格朗日函数：
其中，Y表示拉格朗日乘子，α>0表示惩罚因子，<·>是标准内积；将步骤(1)中的优化问题转化成当Y＝Yk,α＝αk时，对的优化问题进行交替式方法求解；Yk为第K步迭代的拉格朗日乘子，αk第K步迭代的惩罚因子；
(3)迭代更新求解
首先初始化参数Y和α，然后交替更新低秩矩阵L和稀疏矩阵S，迭代更新公式为
其中，运算符表示取使得关于L的函数取得最小值时的L值，表示第K步迭代系数，表示，第K步迭代，和λ有关的迭代系数，Lk+1和Sk+1分别表示第K步迭代所预期得到的L和S矩阵。
当矩阵L和S分别收敛于和时，对矩阵Y按照下式更新：
其中，和是指第K步迭代所求得的L和S的收敛解；Yk+1为第K步迭代得到的拉格朗日乘子；
对参数α按照下式更新，其中ρ>1为常数；ε>0为比较小的正数。
当满足精度要求时，停止迭代，即可得到矩阵L和S的近似解；
其中，αk+1为第K步迭代得到的惩罚因子，表示在求解之前的一个迭代步骤得到的收敛解；
(4)将所得的稀疏矩阵S作为目标回波信号的估计值，低秩矩阵L作为背景杂波信号的估计值，完成探地雷达数据背景去除。
为验证本发明的有益效果，作如下仿真实验：基于鲁棒主元分析的探地雷达数据背景去除
得到结果见图3至图8所示；从图3的原始B-SCAN数据图和图4单道A-SCAN数据图中看到，地表反射波较强，目标回波并不明显，被杂波淹没；图5是具体实施方式四的处理效果，可以看到，经过滑动平均滤波处理，单道A-SCAN数据波形变得平滑，噪声得到抑制；图6和图7是经过鲁棒主元分析处理后分别得到的稀疏矩阵和低秩矩阵灰度图，可以看到，图6中，目标回波得到突出和增强，证明了用稀疏矩阵来近似目标回波的合理性和有效性；图7中，背景杂波在水平方向变化缓慢，时间延迟基本相同，证明了用低秩矩阵来近似背景杂波的合理性和有效性；从图8经过鲁棒主元分析处理后的单道A-SCAN数据来看，背景杂波得到了抑制，目标信号得到突出和增强，也证明了鲁棒主元分析法，确实能有效的去除探地雷达数据的背景杂波。
Claims (4)
1.一种基于鲁棒主元分析的探地雷达数据背景去除方法，其特征在于它按以下步骤进行：
步骤一：探地雷达B-SCAN数据标准化
将探地雷达B-SCAN数据标准化处理，使其二维数据矩阵中的值处于[-1,1]区间；
步骤二：探地雷达A-SCAN数据去除直流偏移
探地雷达回波数据常常出现全部为正或者全部为负的情况，说明含有直流分量，需要进行压制，将每一道A-SCAN数据减去所有A-SCAN数据的均值；
步骤三：利用滑动平均滤波法去除探地雷达B-SCAN图像随机噪声
利用一个滑动窗口，选择窗口宽度为5道，计算窗口内的5道A-SCAN数据的均值，用来代替窗口内上述5道中的中间道的A-SCAN数据；
步骤四：鲁棒主元分析去除背景
对于经过上述步骤处理之后的B-SCAN数据，进行鲁棒主元分析，将其分解成一个低秩矩阵L与一个稀疏矩阵S之和，其中低秩矩阵L作为背景杂波的估计矩阵，稀疏矩阵S作为探地雷达回波信号的估计矩阵，鲁棒主元分析去除背景具体过程为：
(1)建立鲁棒主元分析模型
其中，D表示待处理的探地雷达B-SCAN数据矩阵，L和S分别表示待求解的低秩矩阵和稀疏矩阵，λ表示折中因子，rank(·)表示矩阵的秩，||·||0表示l0范数，表示使关于L和S的函数取得最小值；采用凸松弛法转化成如下公式：
其中，用核范数||·||*替代秩，l1范数||·||1代替l0范数；
(2)应用拉格朗日乘子法，构造拉格朗日函数：
<mrow>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mo>,</mo>
<mi>S</mi>
<mo>,</mo>
<mi>Y</mi>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>|</mo>
<mo>|</mo>
<mi>L</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mo>*</mo>
</msub>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>|</mo>
<mo>|</mo>
<mi>S</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
<mo>+</mo>
<mo><</mo>
<mi>Y</mi>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>-</mo>
<mi>L</mi>
<mo>-</mo>
<mi>S</mi>
<mo>)</mo>
</mrow>
<mo>></mo>
<mo>+</mo>
<mfrac>
<mi>&alpha;</mi>
<mn>2</mn>
</mfrac>
<mo>|</mo>
<mo>|</mo>
<mi>D</mi>
<mo>-</mo>
<mi>L</mi>
<mo>-</mo>
<mi>S</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mi>F</mi>
<mn>2</mn>
</msubsup>
</mrow>
其中，Y表示拉格朗日乘子，α>0表示惩罚因子，<·>是标准内积；将步骤(1)中的优化问题转化成当Y＝Yk,α＝αk时，对的优化问题进行交替式方法求解；Yk为第K步迭代的拉格朗日乘子，αk第K步迭代的惩罚因子；
(3)迭代更新求解
首先初始化参数Y和α，然后交替更新低秩矩阵L和稀疏矩阵S，迭代更新公式为
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>L</mi>
</munder>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mo>,</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>&alpha;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>D</mi>
<mrow>
<mn>1</mn>
<mo>/</mo>
<msub>
<mi>&alpha;</mi>
<mi>k</mi>
</msub>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<mo>/</mo>
<msub>
<mi>&alpha;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>S</mi>
</munder>
<mi>L</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<mi>S</mi>
<mo>,</mo>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>&alpha;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>&lambda;</mi>
<mo>/</mo>
<msub>
<mi>&alpha;</mi>
<mi>k</mi>
</msub>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>-</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<mo>/</mo>
<msub>
<mi>&alpha;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
1
其中，运算符表示取使得关于L的函数取得最小值时的L值，表示第K步迭代系数，表示第K步迭代和λ有关的迭代系数，Lk+1和Sk+1分别表示第K步迭代所预期得到的L和S矩阵；
当矩阵L和S分别收敛于和时，对矩阵Y按照下式更新：
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>-</mo>
<msubsup>
<mi>L</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>*</mo>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
其中，和是指第K步迭代所求得的L和S的收敛解；Yk+1为第K步迭代得到的拉格朗日乘子；
对参数α按照下式更新，其中ρ>1为常数；ε>0为小于10-4的正数；
当满足精度要求时，停止迭代，即可得到矩阵L和S的近似解；
其中，αk+1为第K步迭代得到的惩罚因子，表示在求解之前的一个迭代步骤得到的收敛解；
(4)将所得的稀疏矩阵S作为目标回波信号的估计值，低秩矩阵L作为背景杂波信号的估计值，完成探地雷达数据背景去除。
2.如权利要求1所述一种基于鲁棒主元分析的探地雷达数据背景去除方法，其特征在于步骤一所述的将探地雷达B-SCAN数据标准化处理具体过程如下：
B-SCAN数据D为m×n的矩阵，且所有元素的值处于区间[MIN,MAX]，那么做以下标准化处理：
<mrow>
<mfrac>
<mrow>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mi>M</mi>
<mi>I</mi>
<mi>N</mi>
<mo>+</mo>
<mi>M</mi>
<mi>A</mi>
<mi>X</mi>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
<mrow>
<mi>M</mi>
<mi>I</mi>
<mi>N</mi>
<mo>+</mo>
<mi>M</mi>
<mi>A</mi>
<mi>X</mi>
</mrow>
</mfrac>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>m</mi>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mi>n</mi>
</mrow>
标准化后的B-SCAN数据取值范围变为[-1,1]。
3.如权利要求1或2所述一种基于鲁棒主元分析的探地雷达数据背景去除方法，其特征在于步骤二所述探地雷达A-SCAN数据去除直流偏移具体过程如下：
探地雷达B-SCAN数据中的每一列都称之为一道A-SCAN，将B-SCAN数据D写成A-SCAN的集合：
D＝[A1,A2,...,An]
其中Ai＝[ai1,ai2,...,aim]T表示i第道A-SCAN数据，对每一道A-SCAN数据做如下处理：
<mrow>
<msubsup>
<mi>A</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>m</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
</mrow>
得到去除了直流分量的数据。
4.如权利要求3所述一种基于鲁棒主元分析的探地雷达数据背景去除方法，其特征在于步骤三所述利用滑动平均滤波法去除探地雷达B-SCAN图像随机噪声具体过程如下：
将窗口内的数据以列为单位取均值代替当前列，可以消除图像中的高频噪声干扰；算法设计如下：
对于m×n的数据矩阵，选择滑动窗口宽度为5道，随着窗口的滑动，对窗口内的5道A-SCAN求均值，代替当前窗口中间道的A-SCAN数据；公式如下：
对于边缘的A-SCAN数据，最前两道数据作如下处理：
A1'＝A1
<mrow>
<msubsup>
<mi>A</mi>
<mn>2</mn>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>3</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
最后两道数据，做相同处理：
A'n＝An
<mrow>
<msubsup>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>3</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
对于中间道数据，做如下处理：
<mrow>
<msubsup>
<mi>A</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>5</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mi>i</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</munderover>
<msub>
<mi>A</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<mn>3</mn>
<mo>&le;</mo>
<mi>k</mi>
<mo>&le;</mo>
<mi>n</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
得到的处理后的数据为D'＝[A1',A'2,...,A'n]。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610084826.XA CN105527617B (zh) | 2016-02-06 | 2016-02-06 | 一种基于鲁棒主元分析的探地雷达数据背景去除方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610084826.XA CN105527617B (zh) | 2016-02-06 | 2016-02-06 | 一种基于鲁棒主元分析的探地雷达数据背景去除方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105527617A CN105527617A (zh) | 2016-04-27 |
| CN105527617B true CN105527617B (zh) | 2017-11-07 |
Family
ID=55769954
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610084826.XA Expired - Fee Related CN105527617B (zh) | 2016-02-06 | 2016-02-06 | 一种基于鲁棒主元分析的探地雷达数据背景去除方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105527617B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106646383B (zh) * | 2016-09-28 | 2019-09-06 | 中国科学院电子学研究所 | 一种对月壤结构探测仪数据去除背景的方法 |
| CN108253961A (zh) * | 2016-12-29 | 2018-07-06 | 北京雷动云合智能技术有限公司 | 一种基于imu的轮式机器人定位方法 |
| CN107450054B (zh) * | 2017-07-14 | 2019-09-10 | 浙江省交通规划设计研究院 | 一种自适应探地雷达数据去噪方法 |
| CN109709544A (zh) * | 2019-01-04 | 2019-05-03 | 河北科技大学 | 一种探地雷达杂波去除方法 |
| CN111166322A (zh) * | 2020-01-17 | 2020-05-19 | 哈尔滨工业大学 | 心电信号特征波提取方法 |
| CN111239704B (zh) * | 2020-02-12 | 2021-09-21 | 中国科学院大气物理研究所 | 大气探测雷达信号识别方法、装置、设备及介质 |
| CN112666552B (zh) * | 2020-10-29 | 2023-04-25 | 西安交通大学 | 一种探地雷达数据背景杂波自适应去除方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102866429A (zh) * | 2012-04-28 | 2013-01-09 | 中国神华能源股份有限公司 | 一种地下水分布的确定方法 |
| CN104778298A (zh) * | 2015-01-26 | 2015-07-15 | 江南大学 | 基于egmm的高斯过程回归软测量建模方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4644816B2 (ja) * | 2006-04-14 | 2011-03-09 | 国立大学法人東北大学 | 地中レーダ装置及び画像信号処理方法 |
-
2016
- 2016-02-06 CN CN201610084826.XA patent/CN105527617B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102866429A (zh) * | 2012-04-28 | 2013-01-09 | 中国神华能源股份有限公司 | 一种地下水分布的确定方法 |
| CN104778298A (zh) * | 2015-01-26 | 2015-07-15 | 江南大学 | 基于egmm的高斯过程回归软测量建模方法 |
Non-Patent Citations (5)
| Title |
|---|
| GPR Anomaly Detection with Robust Principal Component Analysis;Matthew P. Masarik等;《Proc. of SPIE：Detection and Sensing of Mines, Explosive Objects, and Obscured Targets》;20150521;第9454卷;945414-1~945414-11页 * |
| Leveraging Robust Principal Component Analysis To Detect Buried Explosive Threats In Handheld Ground-Penetrating Radar Data;Dmitry Kalika等;《Proc. of SPIE：Detection and Sensing of Mines, Explosive Objects, and Obscured Targets》;20150521;第9454卷;94541D-1~94541D-10页 * |
| 探地雷达在混凝土施工质量无损检测中的应用;董延朋;《物探与化探》;20061031;第30卷(第5期);468-470页 * |
| 探地雷达属性技术及其在考古调查中的应用研究;赵文轲;《中国博士学位论文全文数据库 基础科学辑》;20130815(第08期);A011-197 页 * |
| 探地雷达成像算法研究及实现;刘永亮;《中国优秀硕士学位论文全文数据库 信息科技辑》;20061215(第12期);I136-985页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105527617A (zh) | 2016-04-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105527617B (zh) | 一种基于鲁棒主元分析的探地雷达数据背景去除方法 | |
| CN105005042B (zh) | 一种探地雷达地下目标定位方法 | |
| CN102788969B (zh) | 基于短时分数阶傅里叶变换的海面微动目标检测和特征提取方法 | |
| CN104360336B (zh) | 一种自适应提取雷达目标微动周期的新方法 | |
| CN108872962B (zh) | 基于分数阶傅里叶变换的激光雷达微弱信号提取和分解方法 | |
| CN107329137B (zh) | 基于box样条小波框架的探地雷达地下目标定位方法及系统 | |
| CN104122540A (zh) | 基于匹配字典和压缩感知的雷达一维距离像目标识别方法 | |
| CN109557429A (zh) | 基于改进小波阈值去噪的gis局部放电故障检测方法 | |
| CN106556820B (zh) | 一种直达波去除方法 | |
| CN105974468B (zh) | 一种能够同时进行五维地震数据重建和噪声压制的方法 | |
| CN102778680B (zh) | 基于参数化的匀加速运动刚体群目标成像方法 | |
| CN101872484B (zh) | 图像域中多个微弱目标航迹自适应生长检测方法 | |
| CN107390213B (zh) | 一种基于滑动时窗的探地雷达记录剖面的时延曲线提取方法 | |
| CN102721952A (zh) | 基于超分辨多普勒谱的地面运动目标稳健分类方法 | |
| CN101329402A (zh) | 基于改进Wedgelet的多尺度SAR图像边缘检测方法 | |
| CN108107429A (zh) | 基于最大似然估计的前视超分辨成像方法 | |
| CN110895331A (zh) | 基于结构化观测矩阵的脉冲多普勒雷达目标稀疏探测方法 | |
| CN103698765A (zh) | 一种isar成像方位定标方法 | |
| CN102156971A (zh) | 基于线状奇异性信息的sar图像相干斑抑制方法 | |
| CN107831473A (zh) | 基于高斯过程回归的距离‑瞬时多普勒图像序列降噪方法 | |
| CN105447318B (zh) | 一种微弱信号去噪方法及装置 | |
| CN109584256B (zh) | 一种基于霍夫直线检测的脉冲星色散值估计方法 | |
| Wang et al. | VAE-ResNet cascade network: An advanced algorithm for stochastic clutter suppression in ground penetrating radar data | |
| CN105389582A (zh) | 基于clean算法散射中心提取的isar目标轮廓提取方法 | |
| CN105717490B (zh) | 基于时频分析的lfm信号分离及参数估计方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into 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 |
Granted publication date: 20171107 Termination date: 20200206 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |