CN107219523B - 一种穿墙雷达三维成像后多目标后墙一次多径抑制方法 - Google Patents
一种穿墙雷达三维成像后多目标后墙一次多径抑制方法 Download PDFInfo
- Publication number
- CN107219523B CN107219523B CN201710361374.XA CN201710361374A CN107219523B CN 107219523 B CN107219523 B CN 107219523B CN 201710361374 A CN201710361374 A CN 201710361374A CN 107219523 B CN107219523 B CN 107219523B
- Authority
- CN
- China
- Prior art keywords
- multipath
- point
- objects
- radar
- wall
- 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
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/887—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
- G01S13/888—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons through wall detection
-
- 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
- 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/021—Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals
-
- 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/411—Identification of targets based on measurements of radar reflectivity
- G01S7/412—Identification of targets based on measurements of radar reflectivity based on a comparison between measured values and known or stored values
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种穿墙雷达三维成像后多目标后墙一次多径抑制方法,首先从二值化后的原始雷达三维图像中提取出所有图像中物体;然后提取出所有物体重心坐标,并依照它们的重心坐标离雷达阵列中心的欧式距离进行升序排序;其次计算出所有物体到每个天线阵元的直接路径传播距离和多径路径传播距离;然后定义相关匹配系数,并两两计算物体直接的直接路径传播距离和多径路径传播距离的相关匹配系数;最后定义判决门限及规则,根据两个物体的相关匹配系数判定它们是否具有目标‑多径关系,找出图像中所有多径后将他们剔除。本发明能够有效地将穿墙雷达三维图像中的后墙一次多径剔除。
Description
技术领域
本发明涉及穿墙雷达三维成像技术领域,特别是涉及后墙一次多径目标的多径假目标的抑制方法。
背景技术
穿墙雷达三维成像技术主要是利用电磁波穿透建筑物墙体等障碍物,对密闭的建筑物内的静止或移动目标进行三维检测、定位和识别,同时对建筑物布局进行成像。电磁波穿透建筑物墙体后,在密闭建筑物空间内并不沿着一条路线传播,而是存在许多不同的传播路径,其中从发射雷达打到目标然后直接反射回接收雷达的路径称为直接路径,由于建筑物内存在墙体、天花板、地板等反射面,从雷达发射出的电磁波会在墙体等表面进行一次或者多次反射后再经过目标反射回接收雷达,或者先打到目标,然后再反射到内墙面,最后反射回接收雷达,这就会造成雷达收到多次反射电磁波,在雷达成像中表现为多径假目标。由于电磁波在每一次墙面反射时强度都会衰减,所以二次及高次多径的强度都很弱,因此通常只考虑一次多径。多径假目标会导致虚警,严重影响穿墙雷达的成像质量。因此,密闭建筑物中穿墙成像多径幻象抑制十分重要。
当前,关于多径干扰抑制方法中,文献“Z.X.Li,Y.Jia,et al,“A novel approachof multi-path suppression based on sub-aperture imaging in through-wall-radarimaging”,IEEE Radar Conference,2013.”提出了一种基于子孔径成像的多径抑制方法,该方法基于多径位置随着天线阵列移动而移动,但目标位置始终不变的多径特性,通过将不同天线阵列的雷达成像进行相乘融合来达到多径抑制的目的。在基于点扩散函数的多径抑制方法(文献Setlur P,Alli G,Nuzzo L.Multipath exploitation in through-wallradar imaging via point spread functions[J].IEEE Transaction on ImageProcessing.2013.)中,其把多径假目标反推并与其对应的真实目标关联,将多径假目标的幅度值转而叠加到该真实目标上,在抑制多径幅度的同时能够提高目标位置的信杂比。而在一种基于压缩感知的多径抑制方法(文献J.Wang,P.Wang,Y.Li,Q.Song and Z.Zhou,“AMultipath Suppression Technique for Through-the-wall Radar”,IEEEInternational Conference on Ultra-Wideband,2013.)中,其将雷达回波看作是目标空间对环境传输函数的映射,因此在已知建筑物环境的情况下,可以重构传输函数,对雷达回波进行最小范数反转,从而重构真实雷达图像,避免多径假目标。上述已有多径干扰抑制方法只能够应对穿墙雷达二维成像对应的场景,当处理三维图像时,上述方法并不适用。因此,研究密闭空间中三维成像的多目标多径干扰抑制方法在穿墙雷达成像领域具有重要的价值。
发明内容
本发明提供一种适用于密闭建筑物空间三维成像多目标后墙一次多径抑制方法,本发明首先将原始的雷达三维图像进行归一化,然后通过门限检测将其进行二值化,再对该图像进行三维目标提取,得到图像中所有物体(包括真实目标和多径幻想)的三维重心坐标;其次基于多径的产生原理使用相关匹配方法分析所有物体中两两元素间的关联性,判断出目标区域和多径区域;最后在保留目标的同时将多径假目标剔除,从而实现了多径抑制。
本发明的一种穿墙雷达三维成像后多目标后墙一次多径抑制方法,包括以下步骤:
步骤1:三维图像中物体提取。
对原始雷达三维图像进行归一化处理,得到归一化后的原始雷达三维图像Iorg(·),然后对图像Iorg(·)进行二值化,即设置门限Tr,用Tr对图像Iorg(·)进行门限检测得到二元图像Ib(·),对图像Ib(·)进行三维成像物体提取,提取出所有物体的重心坐标,并计算各重心坐标到阵列中心的欧式距离,对所有物体按照欧式距离进行升序排序并按序编号,得到图像Ib(·)的所有物体的距离矩阵P=[P1P2…PK],其中K表示图像Ib(·)中的物体个数,Pk表示第k个物体。
步骤2:相关匹配系数的计算。
计算全部K个物体到所有天线阵元的直接传播距离:
其中
其中,M,N分别为天线收发阵元个数,第m个发射天线位于点Tm,第n个接收天线位于点Rn,εr为前墙相对介电常数,Ak,t和Bk,t为电磁波沿直接路径传播时在前墙后表面的折射点,为各点间对应的长度值,点Pk则对应第k个物体的重心。
其次计算全部K个物体到所有天线阵元的多径传播距离:
其中
可知的取值范围为(0,1],且当取到一个较大的接近1的值时,即说明第j个物体的直接传播距离和第i个物体的多径传播距离匹配程度很高,相关性很强,则可以判定第j个物体与第i个物体有多径-目标关系,即第j个物体是第i个物体的一个后墙一次多径。
由于物体的多径传播距离始终大于其直接传播距离,所以只需要计算当j>i时的相关匹配系数即可。计算得到相关匹配系数矩阵μ:
步骤3:判决准则及多径抑制。
由于K个物体中,最多有K-1个多径,可知对于真正的多径而言,其μ值应该大于最大的K个μ值的平均值,故可以设定判定门限δ为:相关匹配系数矩阵μ中前K个最大值的均值。
最后,将系数矩阵μ中的非零值与门限δ进行对比,若有则可判定第j个物体(Pj)为第i个物体(Pi)的多径。当全部比较完毕后,除了被判定为多径的物体之外,剩下的物体全部为真实目标。最后在二值三维图中去除多径,并将这个多径抑制后的二值三维图与原三维雷达图像进行非相干融合,得到多径抑制后的雷达图像If(·),即If(·)=Ib(·)⊙Iorg(·)。
附图说明
图1为具体实施方式中仿真场景示意图。
图2为具体实施方式中仿真阵列阵元分布示意图。
图3为三维后墙一次多径穿墙传播示意图。
图4为具体实施方式中原始雷达仿真图像。
图5为三维后墙一次多径抑制方法处理流程图。
图6为具体实施方式中最终结果图像。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合实施方式和附图,对本发明作进一步地详细描述。
对存在多个目标的待探测区域,其仿真场景如图1所示,坐标系零点位于前墙前表面的中心位置,4个发射天线和16个接收天线紧贴前墙放置,天线分布图如图2所示,雷达的发射信号为中心频率1.8GHz、带宽0.8GHz的步进频信号,频率步进为2MHz。前墙的厚度为0.24m,后墙前表面位于z=10m处,后墙高度为2m,两墙的相对介电常数相等且均为8.6。两个目标P1、P2分别位于(-1.0,1.3,6.0)m和(1.0,1.3,8.0)m处,其后墙一次多径穿墙传播示意图如图3所示,其中点Tm、Tn分别为发射、接收天线的坐标,点A、B为电磁波沿直接路径传播时在前墙后表面的折射点,点C为电磁波沿多径路径传播时在前墙后表面的折射点,点F为电磁波沿多径路径传播时在后墙前表面的反射点。且目标P1、P2的原始的雷达三维成像如图4所示。
参见图5,根据本发明的后墙一次多径抑制方法,对图1所示的仿真场景的具体处理步骤如下:
步骤1:三维图像中物体提取。
101:使用门限Tr对雷达原始图像(三维图像)进行二值化检测,得到二值化雷达图像,其中门限值取为雷达原始图像最大值像素点值的0.5倍;
102:提取出所有物体的重心坐标,得到四个物体重心坐标,根据仿真场景示意图所建立的坐标系,雷达阵列中心位于(0.0,1.3,0.0)m,则按照物体到阵列中心的距离进行升序排序得到P=[P1 P2 P3 P4],物体的重心坐标和到阵列中心的距离如下表所示:
步骤2:相关匹配系数的计算。
由于多径的传播延时总是大于直接回波传播延时,故无需计算所有的仅计算当j>i的相关匹配系数值。本实施例中,三维图像中共有四个物体,故计算 的值即可。按照相关匹配系数的计算公式,得到所有相关系数值如下表所示:
步骤3:多径判决及抑制。
301:求取门限δ。
303:将删除掉多径幻象后的二值化雷达图像与原始雷达图像进行非相干融合,即可得到多径抑制后的雷达图像,结果如图6所示。
由图6可知,本发明提供的适用于穿墙雷达三维成像后多目标后墙一次多径抑制方法可以有效剔除多径,得到正确的探测结果,验证了本发明的正确性和有效性。
以上所述,仅为本发明的具体实施方式,本说明书中所公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换;所公开的所有特征、或所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以任何方式组合。
Claims (2)
1.一种穿墙雷达三维成像后多目标后墙一次多径抑制方法,其特征在于,包括以下步骤:
步骤1:三维图像中物体提取:
对原始雷达三维图像进行图像预处理,得到二值化的三维图像Ib(·);
对图像Ib(·)进行三维成像物体提取,提取出所有物体的重心坐标,并计算各重心坐标到阵列中心的欧式距离;
对所有物体按照欧式距离进行升序排序并按序编号,得到图像Ib(·)的所有物体的距离矩阵P=[P1 P2…PK],其中K表示图像Ib(·)中的物体个数;
步骤2:相关匹配系数的计算:
εr为前墙相对介电常数,表示点Tm与Ak,t的长度值,表示点Ak,t与物体pk的重心的长度值,表示物体pk的重心与点Bk,t的长度值,表示点Bk,t与Rn之间的长度值,点Tm表示第m个发射天线的位置,点Rn表示第n个接收天线的位置,点Ak,t和Bk,t表示电磁波沿直接路径传播时在前墙后表面的折射点,t表示电磁波发射时间;
表示点Tm与点Cm,t的长度值,表示点Cm,t与F的的长度值,表示点F与物体pk的重心的长度值,点Cm,t为电磁波沿多径路径传播时在前墙后表面的折射点,点F为电磁波沿多径路径传播时在后墙前表面的反射点;
步骤3:判决准则及多径抑制:
基于相关匹配系数矩阵μ设置门限δ:将相关匹配系数矩阵μ中前K个最大值的均值作为门限δ;
从图像Ib(·)中删除判定为多径的物体,再与原始雷达三维图像进行非相干融合,得到多径抑制后的雷达图像。
2.如权利要求1所述的方法,其特征在于,步骤3中,设置门限δ具体为:对相关匹配系数矩阵μ中的非零值进行降序排序并取前K个值的均值作为门限δ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710361374.XA CN107219523B (zh) | 2017-05-22 | 2017-05-22 | 一种穿墙雷达三维成像后多目标后墙一次多径抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710361374.XA CN107219523B (zh) | 2017-05-22 | 2017-05-22 | 一种穿墙雷达三维成像后多目标后墙一次多径抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107219523A CN107219523A (zh) | 2017-09-29 |
CN107219523B true CN107219523B (zh) | 2020-03-27 |
Family
ID=59945481
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710361374.XA Active CN107219523B (zh) | 2017-05-22 | 2017-05-22 | 一种穿墙雷达三维成像后多目标后墙一次多径抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107219523B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108761418B (zh) * | 2018-06-11 | 2022-06-03 | 成都纳雷科技有限公司 | 一种基于微多普勒特征的雷达多径目标剔除方法 |
CN109541585B (zh) * | 2018-10-09 | 2020-11-06 | 中国人民解放军第四军医大学 | 一种基于峰度评估的人体穿墙检测成像方法 |
CN114152943B (zh) * | 2021-12-15 | 2023-04-21 | 电子科技大学 | 一种基于超宽带穿墙雷达的两阶段墙体参数估计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104062638A (zh) * | 2014-06-18 | 2014-09-24 | 电子科技大学 | 一种用于穿墙雷达多目标检测的方法 |
CN105954745A (zh) * | 2016-04-29 | 2016-09-21 | 电子科技大学 | 一种适用于穿墙雷达多径幻象抑制的成像方法 |
CN106093933A (zh) * | 2016-08-22 | 2016-11-09 | 电子科技大学 | 一种穿墙雷达成像后多目标侧墙一次多径抑制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010095946A1 (en) * | 2009-02-20 | 2010-08-26 | Nederlandse Organisatie Voor Toegepast-Natuurwetenschappelijk Onderzoek Tno | A method of detecting a scatterer in a structure, a radar system and a computer program product |
US8912950B2 (en) * | 2009-08-03 | 2014-12-16 | Raytheon Company | Interference mitigation in through the wall radar |
-
2017
- 2017-05-22 CN CN201710361374.XA patent/CN107219523B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104062638A (zh) * | 2014-06-18 | 2014-09-24 | 电子科技大学 | 一种用于穿墙雷达多目标检测的方法 |
CN105954745A (zh) * | 2016-04-29 | 2016-09-21 | 电子科技大学 | 一种适用于穿墙雷达多径幻象抑制的成像方法 |
CN106093933A (zh) * | 2016-08-22 | 2016-11-09 | 电子科技大学 | 一种穿墙雷达成像后多目标侧墙一次多径抑制方法 |
Non-Patent Citations (3)
Title |
---|
Multipath ghosts location and sub-aperture based suppression algorithm for TWIR;Dong Yan et al.;《2016 IEEE Radar Conference (RadarConf)》;20160609 * |
Three-dimensional human imaging for through-the-wall radar;Lingjiang Kong et al.;《2009 IEEE Radar Conference》;20090526 * |
穿墙雷达室内多径机理分析与抑制方法;申文婷 等;《雷达科学与技术》;20161215;605-613 * |
Also Published As
Publication number | Publication date |
---|---|
CN107219523A (zh) | 2017-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107589407B (zh) | 多发多收穿墙雷达成像后前后墙振荡一次多径抑制方法 | |
CN107861123B (zh) | 一种穿墙雷达在复杂环境下对多运动目标实时跟踪的方法 | |
CN106093933B (zh) | 一种穿墙雷达成像后多目标侧墙一次多径抑制方法 | |
CN107219523B (zh) | 一种穿墙雷达三维成像后多目标后墙一次多径抑制方法 | |
US9971019B2 (en) | System and method for through-the-wall-radar-imaging using total-variation denoising | |
CN109031277B (zh) | 一种穿墙雷达多目标图像域稳健跟踪方法 | |
US12096237B2 (en) | Method for predicting structure of indoor space using radio propagation channel analysis through deep learning | |
CN113376613B (zh) | 用于雷达探测的恒虚警检测方法、装置及电子设备 | |
Li et al. | Scale-adaptive human target tracking for through-wall imaging radar | |
KR101897763B1 (ko) | 레이더를 이용한 거리 측정 방법 및 장치 | |
Gingras et al. | Electromagnetic matched-field processing: Basic concepts and tropospheric simulations | |
Jana et al. | Compensation of wall effect for through wall tracking of moving targets | |
Liu et al. | Clutter reduction and target tracking in through-the-wall radar | |
CN113960558B (zh) | 基于多输入多输出雷达的非视距目标定位方法及系统 | |
Wang et al. | Fast detection of GPR objects with cross correlation and hough transform | |
CN108919259A (zh) | 一种基于多径利用的穿墙雷达建筑墙体位置获取方法 | |
CN106680809B (zh) | 一种穿墙雷达自聚焦压缩感知成像方法 | |
CN113064160B (zh) | 一种u型建筑布局下目标定位方法 | |
Yektakhah et al. | A method for detection of flat walls in through-the-wall SAR imaging | |
Kocur et al. | Through wall tracking of moving targets by M-sequence UWB radar | |
Luo et al. | An effective multipath ghost recognition method for sparse MIMO radar | |
CN110133641B (zh) | 一种尺度自适应的穿墙成像雷达目标跟踪方法 | |
CN114720944A (zh) | 一种基于svd迭代熵值的墙体杂波抑制方法 | |
Singh et al. | A novel technique for contrast target detection in through-the-wall radar images | |
Kidera et al. | Ellipse based image extrapolation method with RPM imaging for through-the-wall UWB radar |
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 |