CN105259539B - 一种针对极地探冰雷达数据的分部式噪声抑制方法 - Google Patents
一种针对极地探冰雷达数据的分部式噪声抑制方法 Download PDFInfo
- Publication number
- CN105259539B CN105259539B CN201510716786.1A CN201510716786A CN105259539B CN 105259539 B CN105259539 B CN 105259539B CN 201510716786 A CN201510716786 A CN 201510716786A CN 105259539 B CN105259539 B CN 105259539B
- Authority
- CN
- China
- Prior art keywords
- mrow
- radar data
- branch
- ice radar
- msup
- 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
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/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
- Image Processing (AREA)
Abstract
本发明涉及极地探冰雷达数据的噪声抑制方法，具体涉及一种针对极地探冰雷达数据的分部式噪声抑制方法。包括步骤为：对待处理探冰雷达数据进行预滤波处理；基于滤波后各道雷达数据的幅度和相位，对待处理冰雷达数据进行分部划分；对各分部冰雷达数据，分别进行最优低秩滤波处理；对各分部滤波处理结果，进行累加求和，得到噪声抑制后的探冰雷达数据。本发明基于极地探冰雷达数据的幅度和相位信息，对冰雷达数据的不同介质层位进行分部式滤波处理，通过对不同介质层位的分部式噪声抑制处理，可有效反映冰下微弱信号的变化，去噪同时保留冰下微弱层位特征。
Description
技术领域
本发明涉及极地探冰雷达数据的噪声抑制方法，具体涉及一种针对极地探冰雷达数据的分部式噪声抑制方法。
背景技术
探冰雷达系统凭借探测深度深、分辨率高和数据信息丰富等优点，在各国的冰雪项目中发挥了重要作用。空气和冰盖交界面、冰盖内部层、冰下冰岩交界面、冰下湖是冰雷达探测的主要目标。由于探冰雷达在数据采集过程中会引入噪声干扰，通常噪声和干扰远强于深部目标的反射信号，难以提取深部层位特征。本发明主要涉及探冰雷达数据的随机噪声抑制方法。
近年来，针对随机噪声的抑制问题，人们提出了很多有效的噪声抑制方法，例如平均滤波、基于信号变换域特性的滤波(小波变换、曲波变换等)、基于图像局部平滑性的噪声抑制方法(总变分最小化、基于偏微分方程的方法等)、基于数据稀疏性的去噪方法(压缩感知、字典学习方法等)等。这些噪声抑制方法，多是对数据的全局性去噪处理，没有考虑含噪数据中不同介质的属性差异；虽然提高了图像的信噪比，但在某些情况下会丢失深部弱反射目标。
因此，如何在受噪声干扰的探冰雷达数据中，去噪同时有效保持层位边缘特征，尤其是深部微弱层位特征，对于探冰雷达成功用于极地冰川探测具有重要的意义。
发明内容
本发明的目的在于针对探冰雷达数据的不同介质属性，进行分部式滤波处理，实现噪声抑制的同时，保留层位边缘特征，即在提高信噪比的情况下，保证图像的空间分辨率。
本发明是通过以下技术方案实现的，包括以下步骤：
第一步，对待处理探冰雷达数据X(r,t)进行预滤波处理，得到预处理后的数据
所述的对待处理探冰雷达数据X(r,t)进行预滤波处理，具体是：对数据X(r,t)采用奇异值分解方法，进行低秩滤波处理。
第二步，对计算各道探冰雷达数据的幅度和相位，基于此对待处理探冰雷达数据X(r,t)进行分部划分，得到K个分部探冰雷达数据X(k)(r,t),(k＝1,2,...,K)。
所述的计算各道探冰雷达数据的幅度和相位，基于此对待处理探冰雷达数据X(r,t)进行分部划分，具体是：对滤波后各道探冰雷达数据，计算得到幅度序列和相位序列；绘制幅度序列和相位序列的包络，设定局部极值点判定阈值，得到各包络变化局部极大值；结合幅度序列和相位序列包络变化局部极大值分布，定义不同介质的分部参考位置，对待处理探冰雷达数据X(r,t)进行分部划分，得到K个分部探冰雷达数据X(k)(r,t),(k＝1,2,...,K)。
第三步，对各分部探冰雷达数据X(k)(r,t),(k＝1,2,...,K)，分别进行最优低秩滤波处理，得到各分部滤波处理结果
所述的低秩滤波处理，具体是：对X(k)(r,t),(k＝1,2,...,K)，分别构造其Casorati矩阵如下：
其中，代表在空间位置ri和采样时间tj的空间-时间采样值，N和M分别为第k个分部观测数据的空间和时间长度；对矩阵C(k)，采用奇异值分解方法进行最优低秩估计；假设C(k)的秩为L(k)，则其低秩估计结果如下：
其中，和分别是矩阵C(k)的第l个奇异值、第l个左奇异向量和第l个右奇异向量，所对应的矩阵表示形式分别为Σ(k)，U(k)和V(k)，H为取矩阵的共轭转置；L(k)为C(k)的秩，可由X(k)(r,t)中的噪声方差||ζ(k)||2通过下式估计获得
其中，为C(k)的第个特征值；由此，各分部数据的分部滤波处理结果为
第四步，将各分部滤波结果累加求和，得到噪声抑制后的处理结果为：
与现有技术相比，本发明的有益效果是：
(1)本发明一种针对极地探冰雷达数据的分部式噪声抑制方法，特别是采用预滤波处理技术，可有助于实现对不同介质层位的分部划分；
(2)基于探冰雷达数据的幅度和相位信息，可有效反映出冰下不同介质属性的差异，实现对冰下介质层位的分部划分；
(3)对不同介质层位进行不同程度的分部式滤波处理，可以对不同反射强度的层位特征给予不同的处理，有助于保留深部弱反射目标，去噪同时有效保留深部微弱层位结构。
附图说明
图1为本发明极地探冰雷达数据噪声抑制方法的流程图。
图2为本发明实施例中待处理的极地探冰雷达数据成像图。
图3为发明本实施例中待处理图像第900道的幅度和相位序列信号波形图。
图4为发明本实施例中待处理图像第900道的幅度和相位序列包络图。
图5为对本发明实施例进行噪声抑制后的探冰雷达数据成像图。
图5中，(a)为采用本发明方法的噪声抑制结果，(b)为采用小波变换的噪声抑制结果。
图6为对本发明实施例进行噪声抑制后第900道的信号波形图。
图6中，(a)为采用本发明方法的噪声抑制结果，(b)为采用小波变换的噪声抑制结果。
图7为对本发明实施例进行噪声抑制后第900道的灰度梯度图。
图7中，(a)为采用本发明方法的噪声抑制结果，(b)为采用小波变换的噪声抑制结果。
具体实施方式：
下面结合附图对本发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和过程，但本发明的保护范围不限于下述的实施例。
本发明中所使用的数据来源于中国科学院电子学研究所研制的探冰雷达，没有具体型号，没有批产。具体参数如下：工作带宽：100MHz，中心频率：150MHz，发射信号类型：线性调频信号，天线阵列：2个。
如图1所示，本发明极地探冰雷达数据噪声抑制方法实施例的流程如下：
第一步，对图2所示待处理探冰雷达数据X(r,t)(图像大小为1365×20480)进行预滤波处理，得到预处理后的数据
(1)对X(r,t)，采用奇异值分解方法，取前L＝20个奇异值，进行低秩滤波处理。
第二步，对计算各道探冰雷达数据的幅度和相位，基于此对待处理探冰雷达数据X(r,t)进行分部划分，得到K个分部探冰雷达数据X(k)(r,t),(k＝1,2,...,K)。
(1)对各道探冰雷达数据，计算得到幅度序列和相位序列，例如第900道数据的幅度序列和相位序列，如图3所示。
(2)绘制幅度序列和相位序列的包络，设定局部极值点判定阈值，得到各包络变化局部极大值，如图4所示。
(3)结合幅度序列和相位序列包络变化局部极大值分布，定义不同介质的分部参考位置，本实施例中，将待处理探冰雷达数据划分为地表、冰盖表层、冰盖内部层、冰盖底部、冰岩交界面和冰岩底部，可得6个分部冰雷达数据X(k)(r,t),(k＝1,2,...,6)。
第三步，对各分部探冰雷达数据X(k)(r,t),(k＝1,2,...,6)，分别进行最优低秩滤波处理，得到各分部滤波处理结果
(1)对X(k)(r,t),(k＝1,2,...,6），分别构造其Casorati矩阵如下：
其中，代表在空间位置ri和采样时间tj的空间-时间采样值，N和M分别为第k个分部观测数据的空间和时间长度。
(2)对矩阵C(k)，采用奇异值分解方法进行最优低秩估计：
其中，和分别是矩阵C(k)的第l个奇异值、第l个左奇异向量和第l个右奇异向量，所对应的矩阵表示形式分别为Σ(k)，U(k)和V(k)，H为取矩阵的共轭转置；L(k)为C(k)的秩，可由X(k)(r,t)中的噪声方差||ζ(k)||2通过下式估计获得，本实施例中L(1)＝2，L(2)＝10，L(3)＝10，L(4)＝2，L(5)＝10，L(6)＝2。
其中，为C(k)的第个特征值。
(3)基于奇异值分解结果，得到各分部数据的分部滤波处理结果为
第四步，将各分部滤波结果累加求和，得到噪声抑制后的处理结果为：
本实施例待处理探冰雷达数据第900道数据幅度序列和相位序列波形如图3所示，其对应包络如图4所示。本实施例噪声抑制结果、第900道数据噪声抑制后信号波形图、第900道数据噪声抑制后灰度梯度图分别如图5、图6和图7所示。与图5(b)、图6(b)和图7(b)小波滤波处理噪声抑制结果相比可见，本发明方法用于极地探冰雷达数据噪声抑制时，可以更有效的反映出不同介质层位的属性差异，通过对不同介质层位的分部式噪声抑制处理，反映出了冰下微弱信号变化，有效保留了冰下微弱层位特征，如图6(a)和图7(a)；在提高图像信噪比的同时，有效保证了图像的空间分辨率。
上述说明已经充分介绍了本发明的具体实施方式。需要指出的是，熟悉该领域的技术人员对本发明的具体实施方式所做的任何改动均不脱离本发明的权利要求书的范围。相应的，本发明的权利要求的范围也并不仅仅局限于前述具体实施方式。
Claims (3)
1.一种针对极地探冰雷达数据的分部式噪声抑制方法，其特征在于，包括以下步骤：
(1)对待处理探冰雷达数据X(r,t)进行预滤波处理，得到预处理后的数据
(2)对计算各道探冰雷达数据的幅度和相位，基于此对待处理探冰雷达数据X(r,t)进行分部划分，得到K个分部探冰雷达数据X(k)(r,t),(k＝1,2,...,K)；
(3)对各分部探冰雷达数据X(k)(r,t),(k＝1,2,...,K)，分别进行最优低秩滤波处理，得到各分部滤波处理结果
(4)将各分部滤波结果累加求和，得到噪声抑制后的处理结果为：
<mrow>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msup>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
所述步骤(3)的对各分部探冰雷达数据X(k)(r,t),(k＝1,2,...,K)，分别进行低秩滤波处理，具体是：对X(k)(r,t),(k＝1,2,...,K)，分别构造其Casorati矩阵如下：
其中，代表在空间位置ri和采样时间tj的空间-时间采样值，N和M分别为第k个分部观测数据的空间和时间长度；对矩阵C(k)，采用奇异值分解方法进行最优低秩估计；假设C(k)的秩为L(k)，则其低秩估计结果如下：
<mrow>
<msup>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>*</mo>
</mrow>
</msup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msup>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
</munderover>
<msubsup>
<mi>&sigma;</mi>
<mi>l</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msubsup>
<mi>u</mi>
<mi>l</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msubsup>
<mi>v</mi>
<mi>l</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mi>H</mi>
</mrow>
</msubsup>
<mo>,</mo>
</mrow>
其中，和分别是矩阵C(k)的第l个奇异值、第l个左奇异向量和第l个右奇异向量，所对应的矩阵表示形式分别为Σ(k)，U(k)和V(k)，H为取矩阵的共轭转置；L(k)为C(k)的秩，可由X(k)(r,t)中的噪声方差||ζ(k)||2通过下式估计获得
<mrow>
<msubsup>
<mi>&lambda;</mi>
<mrow>
<mover>
<mi>L</mi>
<mo>^</mo>
</mover>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msup>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mo>|</mo>
<mo>|</mo>
<msup>
<mi>&zeta;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
<mo>&le;</mo>
<msubsup>
<mi>&lambda;</mi>
<mover>
<mi>L</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msup>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中，为C(k)的第个特征值；由此，各分部数据的分部滤波处理结果为
2.根据权利要求1所述的针对极地探冰雷达数据的分部式噪声抑制方法，其特征在于，所述步骤(1)的对待处理探冰雷达数据X(r,t)进行预滤波处理，具体是：对数据X(r,t)采用奇异值分解方法，进行低秩滤波处理。
3.根据权利要求1所述的针对极地探冰雷达数据的分部式噪声抑制方法，其特征在于，所述步骤(2)的计算各道雷达数据的幅度和相位，基于此对待处理探冰雷达数据X(r,t)进行分部划分，具体是：对滤波后各道探冰雷达数据，计算得到幅度序列和相位序列；绘制幅度序列和相位序列的包络，设定局部极值点判定阈值，得到各包络变化局部极大值；结合幅度序列和相位序列包络变化局部极大值分布，定义不同介质的分部参考位置，对待处理探冰雷达数据X(r,t)进行分部划分。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510716786.1A CN105259539B (zh) | 2015-10-29 | 2015-10-29 | 一种针对极地探冰雷达数据的分部式噪声抑制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510716786.1A CN105259539B (zh) | 2015-10-29 | 2015-10-29 | 一种针对极地探冰雷达数据的分部式噪声抑制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105259539A CN105259539A (zh) | 2016-01-20 |
| CN105259539B true CN105259539B (zh) | 2017-09-29 |
Family
ID=55099300
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510716786.1A Active CN105259539B (zh) | 2015-10-29 | 2015-10-29 | 一种针对极地探冰雷达数据的分部式噪声抑制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105259539B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109492190B (zh) * | 2018-08-12 | 2022-11-15 | 中国科学院大学 | 一种基于分部式高斯混合模型的冰下层位检测方法 |
| CN109613490B (zh) * | 2018-12-21 | 2022-12-06 | 中国科学院国家天文台 | 一种测月雷达探测数据层位信号有效性的判定方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101561497B (zh) * | 2009-05-22 | 2011-09-21 | 西安电子科技大学 | 机载雷达杂波抑制方法 |
| CN102005033B (zh) * | 2010-11-16 | 2012-07-04 | 中国科学院遥感应用研究所 | 一种图像平滑抑制噪声方法 |
| FR2978273B1 (fr) * | 2011-07-22 | 2013-08-09 | Thales Sa | Procede de reduction du bruit dans une sequence d'images fluoroscopiques par filtrage temporel et spatial |
| CN103364770B (zh) * | 2013-07-14 | 2015-04-22 | 西安电子科技大学 | 基于矩阵填充的雷达目标检测系统及其检测方法 |
| CN103913733B (zh) * | 2014-04-14 | 2016-06-15 | 中国科学院电子学研究所 | 极地冰川厚度探测方法 |
| CN104463325A (zh) * | 2014-12-04 | 2015-03-25 | 中国科学院电子学研究所 | 一种极地探冰雷达原始数据噪声抑制方法 |
-
2015
- 2015-10-29 CN CN201510716786.1A patent/CN105259539B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105259539A (zh) | 2016-01-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105242245B (zh) | 基于极地探冰雷达数据低秩性和稀疏性的噪声抑制方法 | |
| CN105913393B (zh) | 一种自适应小波阈值图像去噪方法及装置 | |
| CN109658351B (zh) | 一种结合l0梯度约束和局部低秩矩阵恢复的高光谱图像去噪方法 | |
| CN101944230B (zh) | 基于多尺度的自然图像非局部均值去噪方法 | |
| CN106441288A (zh) | 一种加速度计的自适应小波降噪方法 | |
| CN108961181B (zh) | 一种基于shearlet变换的探地雷达图像去噪方法 | |
| CN107179550B (zh) | 一种数据驱动的地震信号零相位反褶积方法 | |
| CN102818629A (zh) | 基于平稳小波变换的微型光谱仪信号去噪方法 | |
| CN103020916A (zh) | 一种二维希尔伯特变换和bemd结合的图像去噪方法 | |
| CN113156513A (zh) | 一种基于注意力引导的卷积神经网络地震信号去噪方法 | |
| CN104820969A (zh) | 一种实时图像盲复原方法 | |
| CN103761719A (zh) | 一种基于邻域相关性的自适应小波阈值去噪方法 | |
| CN112087774A (zh) | 一种基于残差神经网络的通信辐射源个体识别方法 | |
| CN109472758B (zh) | 一种地震剖面图像纹理细节增强方法 | |
| CN104849757A (zh) | 消除地震信号中随机噪声系统及方法 | |
| CN112084845B (zh) | 基于多尺度小波系数自相关的低频1/f噪声消除方法 | |
| CN105259539B (zh) | 一种针对极地探冰雷达数据的分部式噪声抑制方法 | |
| CN103093432B (zh) | 基于极化分解和图像块相似性的极化sar图像降斑方法 | |
| CN112882115A (zh) | 基于gwo优化小波阈值的大地电磁信号去噪方法及系统 | |
| CN105550998A (zh) | 基于二代小波整数变换的图像增强方法及图像增强系统 | |
| CN105319593A (zh) | 基于曲波变换和奇异值分解的联合去噪方法 | |
| CN114779343A (zh) | 一种基于曲波变换-联合双边滤波的地震数据去噪方法 | |
| Wang et al. | An adaptive learning image denoising algorithm based on eigenvalue extraction and the GAN model | |
| CN117574062A (zh) | 基于vmd-dnn模型的小回线瞬变电磁信号去噪方法 | |
| CN108169795A (zh) | 基于随机采样的数据规则化方法 |
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 |