CN112986983B - 一种全息穿透成像雷达杂波抑制方法 - Google Patents
一种全息穿透成像雷达杂波抑制方法 Download PDFInfo
- Publication number
- CN112986983B CN112986983B CN202110165788.1A CN202110165788A CN112986983B CN 112986983 B CN112986983 B CN 112986983B CN 202110165788 A CN202110165788 A CN 202110165788A CN 112986983 B CN112986983 B CN 112986983B
- Authority
- CN
- China
- Prior art keywords
- echo data
- matrix
- clutter
- target
- holographic
- 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 40
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000001629 suppression Effects 0.000 title claims abstract description 22
- 230000035515 penetration Effects 0.000 title claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims abstract description 43
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 19
- 230000000149 penetrating effect Effects 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000000203 mixture Substances 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/887—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种全息穿透成像雷达杂波抑制方法。技术方案是:首先,对接收到的全息穿透成像雷达回波数据进行奇异值分解;然后,以目标杂波比(TCR)最大作为准则求解重构系数矩阵;最后,利用求解出的重构系数矩阵,计算出杂波抑制后的回波数据。本发明能够有效抑制全息穿透成像雷达回波数据中的杂波分量,同时保留目标信息,实现介质中的埋藏弱目标增强成像。
Description
技术领域
本发明涉及全息穿透成像雷达技术领域,特别涉及一种全息穿透成像雷达的杂波抑制方法。
背景技术
全息穿透成像雷达是一种通过向被检测介质发射电磁波,探测介质中的电磁特性不连续的成像雷达设备,通过成像处理过程可对介质中的埋藏目标进行高分辨率二维成像,目前已被应用于建筑检查、无损探测和生物医学等领域。
然而,在全息穿透成像雷达工作过程中,由于探测天线与介质表面距离近且波束覆盖范围大,信号在首次穿透介质时大部分能量将直接返回,形成介质表面反射强杂波,其能量一般大于埋藏在介质内的目标信号,此外,收发天线之间的直耦波也会带来杂波干扰,最终导致成像结果中目标成分被遮蔽,因此,杂波抑制对于全息穿透成像雷达具有重要意义。
子空间分解技术将全息穿透成像雷达回波数据矩阵分解成杂波子空间与目标子空间,通过对目标子空间信息的重构,获得去除介质表面杂波与噪声的回波数据矩阵,达到杂波抑制的目的,实现埋藏弱目标增强成像。如何准确提取出目标子空间是通过子空间分解技术实现杂波抑制的难点,是全息穿透成像雷达埋藏弱目标探测的关键技术之一。在专利号为201918000350.1的文献中,提出的全息穿透成像雷达杂波抑制方法,取最大奇异值分量作为杂波分量,仅适用于特定应用场景,可能存在杂波抑制不彻底的情况。
发明内容
针对现有的全息穿透成像雷达回波数据成像方法中应用的子空间分解方法,无法准确区分目标子空间和杂波子空间这一难点,本发明提出了一种判断目标子空间成分的准则,利用这种准则可以准确提取出目标子空间,实现对全息穿透成像雷达回波数据的杂波抑制,提高了成像结果信杂比,改善了成像效果。
本发明的技术方案是:一种全息穿透成像雷达杂波抑制方法,具体包括如下步骤:
首先,对接收到的全息穿透成像雷达回波数据进行奇异值分解;然后,以目标杂波比(TCR)最大作为准则求解重构系数矩阵;最后,利用求解出的重构系数矩阵,计算出杂波抑制后的回波数据。
本发明还提供一种全息穿透成像雷达杂波抑制模块,其特征在于,具体包括:
对接收到的全息穿透成像雷达回波数据进行奇异值分解的子模块;以目标杂波比最大作为准则求解重构系数矩阵的子模块;利用求解出的重构系数矩阵,计算出杂波抑制后的回波数据的子模块。
所述模块采用现有的硬件,如FPGA或DSP数据处理平台即可实现。
本发明的有益效果是:本发明通过对雷达回波数据互相关特性、奇异值分布特性进行分析,归纳出左奇异向量方差与信号成分的组成规律,可以有效估计雷达回波数据奇异值分解后各分量中的目标和杂波成分,弥补了现有子空间分解方法难以准确区分目标子空间和杂波子空间的不足。实验结果表明本发明方法能够有效抑制全息穿透成像雷达回波数据中的杂波分量,同时保留目标信息,实现介质中的埋藏弱目标增强成像。
附图说明
图1是本发明提供的一种全息穿透成像雷达杂波抑制方法的流程示意图;
图2是成像实验场景中的介质与目标实物图;
图3是全息穿透成像雷达原始回波数据成像结果;
图4是利用现有子空间分解方法处理后的成像结果;
图5是利用本发明提供的方法处理后的成像结果。
具体实施方式
下面结合附图对本发明进行进一步的说明。
图1是本发明提供的一种全息穿透成像雷达杂波抑制方法的流程示意图,整个流程分为三步。第一步:雷达回波数据奇异值分解处理,第二步:求解重构系数矩阵;第三步:数据矩阵重构。详细描述如下:
第一步:雷达回波数据奇异值分解处理
根据奇异值分解原理,对雷达回波数据进行奇异值分解:
D=USVH
其中,D表示雷达回波数据,为M阶方阵,U=[u1,u2,...ui,...,uM]为M阶左奇异酉矩阵,V=[v1,v2,...vi,...,,vM]为M阶右奇异酉矩阵;S=diag(s1,s2,...si,...,sM)为M阶奇异值对角矩阵,其中,diag(·)表示对角矩阵,对角元素si表示雷达回波数据D的奇异值,满足s1≥s2≥...si...≥sM;H表示共轭转置。在实际工作时,如果雷达回波数据D,不是方阵,则可通过补零的方式获得方阵。
第二步:求解重构系数矩阵
通过最大化目标杂波比TCR,求解当目标杂波比TCR取得最大值时,对应的系数矩阵W为得到的重构系数矩阵即:
其中,目标杂波比TCR为:
上式中,σ2(ui)为ui的方差;Ii为第i个子矩阵Di'的所有像素值平方的和,
第三步:数据矩阵重构
利用重构系数矩阵对回波数据矩阵进行重构,得到杂波抑制后的回波数据/>
其中⊙表示矩阵的哈达玛积。
本发明的基本原理是:
对雷达回波数据D进行奇异值分解后,令B=DDH,则B为M阶矩阵,设dmn表示回波数据D的第m行、第n列元素值,bmn表示M阶矩阵B的第m行、第n列元素值,并且可计算出其各个元素为:
其中,上标*表示共轭,可见bmn表示回波数据D的互相关。
令λi表示B的特征值,与回波数据D的奇异值si的关系为根据奇异值分解原理,B×ui=λi×ui,进一步地,可以计算出:
令ui=[u1i,u2i,...,uMi]T,其中,T表示向量转置,可以得到:
进一步地,计算出ui的均值E(ui)和方差σ2(ui):
可以看出ui的方差σ2(ui)与B中列向量bl的均值相关,而bil是数据矩阵的互相关,表示回波数据中不同数据点之间的相关性。当回波数据中杂波很强时,回波数据主要由介质杂波组成,天线在不同位置信号相近,即dij≈d,d为常数,此时数据矩阵互相关bil也可认为是一个常值b,即bil=b。
代入E(ui),可以得到:
因此,可以近似认为Mb=λi,再将其代入σ2(ui),得到:
即当回波数据分量中只包含均匀介质回波时,σ2(ui)等于0,而当回波数据分量中包含目标时,天线在不同位置处信号存在较大差异,其相关性减弱,此时σ2(ui)>0。在实际情况中,当回波数据中存在强杂波掩盖目标信号时,回波数据分量主要由介质本身回波组成,此时σ2(ui)略大于0,当目标信号增强时,σ2(ui)随之也增大,因此,可将归一化后的σ2(ui)值作为回波数据分量中目标信号占比的衡量标准。
根据上述分析,本发明通过最大化目标杂波比TCR,求解当目标杂波比TCR取得最大值时,对应的系数矩阵W为得到的重构系数矩阵即:
其中,目标杂波比TCR为:
上式中,σ2(ui)为ui的方差;Ii为第i个子矩阵Di'的所有像素值平方的和,
利用重构系数矩阵对回波数据矩阵进行重构,即可得到杂波抑制后的回波数据/>
图2是成像实验中使用的介质与目标实物图,介质为混凝土,厚度为4cm;目标为“十”字形金属片,宽度为1cm,长度为10cm。实验时,目标紧贴在介质板的下表面,相当于目标位于介质中。
图3是全息穿透成像雷达原始回波数据成像结果,可以看出,混凝土介质存在很强的介质表面反射杂波,完全将目标信号掩盖。
图4是利用现有子空间分解方法处理后的成像结果,可以看出,使用现有奇异值方法处理后,图像左上角杂波得到较好的抑制,目标图像得到凸显,但图像右下角仍然存在较多杂波。
图5是利用本发明提供的方法处理后的成像结果。可以看出,通过改进子空间分解方法处理后,进一步抑制了图像右下角杂波,并保留了信号中的目标成分,验证了本发明所提方法的有效性。
以上所述的本发明实施方式,并不构成对本发明保护范围的限定,任何在本发明精神和原则之内所作的修改、等同替换和改进等,均应包含在本发明的权利要求保护范围之内。
Claims (1)
1.一种全息穿透成像雷达杂波抑制方法,其特征在于,具体包括如下步骤:
首先,对接收到的全息穿透成像雷达回波数据进行奇异值分解;然后,以目标杂波比最大作为准则求解重构系数矩阵;最后,利用求解出的重构系数矩阵,计算出杂波抑制后的回波数据;
其中,利用下式对雷达回波数据进行奇异值分解:
D=USVH
其中,D表示雷达回波数据,为M阶方阵,U=[u1,u2,...ui,...,uM]为M阶左奇异酉矩阵,V=[v1,v2,...vi,...,vM]为M阶右奇异酉矩阵;S=diag(s1,s2,...si,...,sM)为M阶奇异值对角矩阵,其中,diag(·)表示对角矩阵,对角元素si表示雷达回波数据D的奇异值,满足s1≥s2≥...si...≥sM;H表示共轭转置;
通过最大化目标杂波比,求解当目标杂波比取得最大值时,对应的系数矩阵W为得到的重构系数矩阵即:
其中,目标杂波比TCR为:
上式中,σ2(ui)为ui的方差;Ii为第i个子矩阵Di'的所有像素值平方的和,D′i=ui(wisi)vi H;
利用重构系数矩阵对回波数据矩阵进行重构,得到杂波抑制后的回波数据/>
其中,⊙表示矩阵的哈达玛积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110165788.1A CN112986983B (zh) | 2021-02-06 | 2021-02-06 | 一种全息穿透成像雷达杂波抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110165788.1A CN112986983B (zh) | 2021-02-06 | 2021-02-06 | 一种全息穿透成像雷达杂波抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112986983A CN112986983A (zh) | 2021-06-18 |
CN112986983B true CN112986983B (zh) | 2023-12-19 |
Family
ID=76348599
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110165788.1A Active CN112986983B (zh) | 2021-02-06 | 2021-02-06 | 一种全息穿透成像雷达杂波抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112986983B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114527462B (zh) * | 2022-04-22 | 2022-07-15 | 中国人民解放军国防科技大学 | 全息穿透成像雷达杂波抑制方法、装置、系统及介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998048260A1 (en) * | 1997-04-18 | 1998-10-29 | Raytheon Company | Method and system for remotely determining column density of trace gases |
EP1580573A1 (en) * | 2002-03-13 | 2005-09-28 | Raytheon Canada Limited | System and method for noise suppression in pre-processed radar data |
WO2013033394A2 (en) * | 2011-08-30 | 2013-03-07 | Watlow Electric Manufacturing Company | High definition heater system having a fluid medium |
CN104155632A (zh) * | 2014-07-18 | 2014-11-19 | 南京航空航天大学 | 一种基于局部相关性的改进子空间海杂波抑制方法 |
CN105929373A (zh) * | 2016-04-18 | 2016-09-07 | 北京理工大学 | 一种用于浅埋目标高分辨率透视成像雷达杂波抑制方法 |
CN107957574A (zh) * | 2017-12-28 | 2018-04-24 | 桂林电子科技大学 | 基于ifft和混合匹配追踪的时分地基mimo滑坡雷达成像方法 |
CN111580099A (zh) * | 2020-06-12 | 2020-08-25 | 南京信息工程大学 | 基于联合熵值的穿墙成像雷达墙体杂波抑制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090238426A1 (en) * | 2008-03-19 | 2009-09-24 | Uti Limited Partnership | System and Methods for Identifying an Object within a Complex Environment |
CN110009709B (zh) * | 2019-05-08 | 2023-07-07 | 上海联影医疗科技股份有限公司 | 医学图像成像方法和系统 |
-
2021
- 2021-02-06 CN CN202110165788.1A patent/CN112986983B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998048260A1 (en) * | 1997-04-18 | 1998-10-29 | Raytheon Company | Method and system for remotely determining column density of trace gases |
EP1580573A1 (en) * | 2002-03-13 | 2005-09-28 | Raytheon Canada Limited | System and method for noise suppression in pre-processed radar data |
WO2013033394A2 (en) * | 2011-08-30 | 2013-03-07 | Watlow Electric Manufacturing Company | High definition heater system having a fluid medium |
CN104155632A (zh) * | 2014-07-18 | 2014-11-19 | 南京航空航天大学 | 一种基于局部相关性的改进子空间海杂波抑制方法 |
CN105929373A (zh) * | 2016-04-18 | 2016-09-07 | 北京理工大学 | 一种用于浅埋目标高分辨率透视成像雷达杂波抑制方法 |
CN107957574A (zh) * | 2017-12-28 | 2018-04-24 | 桂林电子科技大学 | 基于ifft和混合匹配追踪的时分地基mimo滑坡雷达成像方法 |
CN111580099A (zh) * | 2020-06-12 | 2020-08-25 | 南京信息工程大学 | 基于联合熵值的穿墙成像雷达墙体杂波抑制方法 |
Non-Patent Citations (8)
Title |
---|
Front-Wall Clutter Removal in Through-the-Wall Radar Based on Weighted Nuclear Norm Minimization;Yi Zhou 等;《IEEE Geoscience and Remote Sensing Letters》;第1-5页 * |
Improvement of Surface Penetrating Radar Imaging by Suppressing Clutter Using Nonlinear Gain Control;Yawen Huang 等;《Twelfth International Conference on Digital Image Processing (ICDIP 2020)》;第1-5页 * |
Target-to-Clutter Ratio Enhancement of Images in Through-the-Wall Radar Using a Radiation Pattern-Based Delayed-Sum Algorithm;Youngjoon Lim 等;《JOURNAL OF ELECTROMAGNETIC ENGINEERING AND SCIENCE》;第405-410页 * |
基于K-均值聚类的SVD杂波抑制算法;黄凤青 等;《雷达科学与技术》;第611-617页 * |
基于组合滤波的压缩感知穿墙雷达杂波抑制;夏朝禹 等;《电讯技术》;第1186-1190页 * |
相控阵雷达资源优化管理的理论与方法;卢建斌;《中国博士学位论文全文数据库 信息科技辑》;第1-172页 * |
融合极化旋转域特征和超像素技术的极化SAR舰船检测;崔兴超 等;《雷达学报》;第1-15页 * |
超声血流成像中基于动态域的回归和奇异值分解的杂波抑制方法;肖磊 等;《计算机应用》;第265-269页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112986983A (zh) | 2021-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang | CFAR-based interference mitigation for FMCW automotive radar systems | |
US9250323B2 (en) | Target detection utilizing image array comparison | |
US9075129B2 (en) | Method and system for forming images by comparing subsets of image data | |
Sjögren et al. | Suppression of clutter in multichannel SAR GMTI | |
Ram et al. | Through-wall tracking of human movers using joint Doppler and array processing | |
Malanowski et al. | Detection of moving targets with continuous-wave noise radar: theory and measurements | |
DE60304692T2 (de) | System und verfahren zum erzeugen eines spektrums in radargeräten | |
CN105974405B (zh) | 基于幅度加权的探地雷达后向投影成像方法 | |
US6720910B2 (en) | Pri-staggered post-doppler adaptive monopulse processing for detection and location of a moving target in ground clutter | |
Wang et al. | Adaptive imaging for forward-looking ground penetrating radar | |
CN112986983B (zh) | 一种全息穿透成像雷达杂波抑制方法 | |
CN112098956B (zh) | 一种基于成像技术的交叉眼干扰对抗方法 | |
JP2009020015A (ja) | レーダ装置 | |
Mao et al. | An efficient anti-interference imaging technology for marine radar | |
CN106199538B (zh) | 用于提升扩频穿墙雷达跟踪动目标精度的杂波抑制方法 | |
Sharma et al. | Non-metallic pipe detection using SF-GPR: A new approach using neural network | |
Ram et al. | Human tracking using doppler processing and spatial beamforming | |
Agarwal et al. | Non-invasive concealed weapon detection and identification using V band millimeter wave imaging radar system | |
CN116540196A (zh) | 一种基于距离补偿和低秩稀疏分解的钢筋杂波抑制方法 | |
CN114152918A (zh) | 基于压缩感知的抗间歇式主瓣干扰方法 | |
Almutiry | Wideband tomographic super-resolution radar image | |
Jouadé et al. | High resolution radar focusing using spectral estimation methods in wide-band and near-field configurations: Application to millimeter-wave near-range imaging | |
Suwa et al. | Reconstruction of 3-D target geometry using radar movie | |
CN114527462B (zh) | 全息穿透成像雷达杂波抑制方法、装置、系统及介质 | |
Hu et al. | Grating lobe mitigation based on extended coherence factor in sparse MIMO UWB array |
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 |