CN106093933B - 一种穿墙雷达成像后多目标侧墙一次多径抑制方法 - Google Patents
一种穿墙雷达成像后多目标侧墙一次多径抑制方法 Download PDFInfo
- Publication number
- CN106093933B CN106093933B CN201610698090.5A CN201610698090A CN106093933B CN 106093933 B CN106093933 B CN 106093933B CN 201610698090 A CN201610698090 A CN 201610698090A CN 106093933 B CN106093933 B CN 106093933B
- Authority
- CN
- China
- Prior art keywords
- vector
- multipath
- connected domain
- side wall
- radar
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000003384 imaging method Methods 0.000 title claims abstract description 25
- 239000013598 vector Substances 0.000 claims abstract description 144
- 239000011159 matrix material Substances 0.000 claims abstract description 47
- 230000005540 biological transmission Effects 0.000 claims abstract description 32
- 238000001514 detection method Methods 0.000 claims description 18
- 230000001629 suppression Effects 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000012163 sequencing technique Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 230000004927 fusion Effects 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 2
- 230000001174 ascending effect Effects 0.000 abstract description 2
- 230000000717 retained effect Effects 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 4
- 230000005764 inhibitory process Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002401 inhibitory effect Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000000452 restraining effect Effects 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
- G01S13/888—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons through wall detection
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开一种适用于穿墙雷达成像的多目标侧墙一次多径抑制方法,首先利用4连通域检测方法从二元原始雷达图像中提取出所有连通域;然后对所有连通域编号并依照它们的平均坐标离雷达阵列中心的欧式距离升序排序得到连通域向量;接着把连通域向量元素重组成两个子向量,分别按照到雷达的多径传输路径和直接回波传输路径对这两个子向量进行扩展;然后利用相似度函数计算两个子向量的相似度矩阵;最后使用相似度门限及判决规则的方法判定相似度矩阵的所有元素,得到真实目标所在的连通域并将多径假目标抑制掉。本发明能够对墙后隐蔽多目标进行检测,在保留目标的同时能够有效地将多径剔除。
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-wall radar imaging via pointspread functions[J].IEEE Transaction on Image Processing.2013.),它把多径假目标的值叠加到与它关联的真实目标上,在抑制多径的同时能够提高目标位置的信杂比。上述方法只能够应对单目标场景,或者多径与目标之间不存在重叠的场景,当处理存在目标与多径假目标重叠的多目标场景时,上述方法会造成多径假目标区域内的真实目标丢失,出现漏警。因此,研究密闭空间中存在重叠现象的多目标多径干扰抑制方法在穿墙雷达成像领域具有重要的价值。
发明内容
本发明为解决上述技术问题,提出了一种穿墙雷达成像后多目标侧墙一次多径抑制方法,首先通过门限检测将原始雷达图像进行二值化,使用连通域检测方法对该二值化图像进行连通域检测得到连通域向量,基于多径的产生原理使用相似度方法分析该连通域向量中两两元素间的关联性,判断出目标区域和多径区域,最后在保留目标的同时将多径假目标剔除,从而实现了多径抑制。
本发明采用的技术方案是:一种穿墙雷达成像后多目标侧墙一次多径抑制方法,包括:
S1、连通区域检测,根据设置的第一门限对原始雷达图像进行门限检测得到二元图像,对二元图像进行4-连通检测,并标记二元图像中的所有的N个连通域,按照平均坐标值到阵列中心的欧式距离由小到大对所有的N个连通域进行排序,得到一个连通域向量;
S2、构造相似度矩阵,将步骤S1得到的连通域向量分成两组,得到第一向量与第二向量;根据侧墙一次多径的形成原理对第一向量进行扩展,得到包含侧墙一次多径传输距离信息的第三向量;根据雷达的直接回波传输路径对第二向量进行扩展,得到包含直接回波传输距离信息的第四向量;根据相似度函数计算得到第三向量与第四向量的相似度矩阵;
S3、元素判决及多径抑制,根据第二门限判定相似度矩阵中的所有元素,得到目标区域。
进一步地,步骤S2中所述将步骤S1得到的所有连通域向量分成两组,具体为:所述第一向量为按照各个连通域的平均坐标到阵列中心的欧式距离由小到大的顺序连续地在连通域向量中取连通域,所述第二向量为按照平均坐标值到阵列中心的欧式距离由大到小的顺序连续地在连通域向量中取连通域,所述第一向量与第二向量所取连通域个数相等,记为N′,N′大于连通域向量中总的连通域个数的一半,并且小于总的连通域个数。
进一步地,步骤S2中所述的根据侧墙一次多径的形成原理对第一向量进行扩展,得到包含侧墙一次多径传输距离信息的第三向量;具体为:计算第一向量中每个连通域的平均坐标到K个雷达位置的侧墙一次多径传输距离,每个连通域得到的K个距离值组成一个向量ri;i=1,2,3…,N′,根据N′个向量ri得到包含侧墙一次多径传输距离信息的第三向量。
步骤S2中所述的根据雷达的直接回波传输路径对第二向量进行扩展,得到包含直接回波传输距离信息的第四向量;具体为:计算第二向量中每个连通域的平均坐标到K个雷达位置的直接回波距离,每个连通域得到的K个距离值组成一个向量ri′′;i′=1,2,3…,N′,根据N′个向量ri′′得到包含N′个直接回波传输距离信息的第四向量。
进一步地,步骤S2中所述的根据相似度函数计算得到第三向量与第四向量的相似度矩阵,具体为:根据相似度函数定义,第三向量与第四向量的相似度矩阵表达式为:
其中,sim(r1,r′2)表示计算向量r1和r′2之间相似度的相似度计算函数,其取值范围为[0,1]。
进一步地,步骤S3所述根据第二门限判定相似度矩阵中的所有元素,得到目标区域,具体为:
A1、对相似度矩阵S进行预处理,根据多径回波路径始终长于目标回波路径的特点,将相似度矩阵S中的上三角元素全部置零,得到一个下三角矩阵S′
A2、将矩阵S′中的所有非零元素按照取值从大到小依次放入行向量Sv:
A3、假如只考虑侧墙一次多径的影响,那么对于连通域向量D中的N个元素,最多存在N-1个目标,因此选取向量Sv中的前N-1个元素用来计算用于相似度检测的第二门限;
A4、将矩阵S′中大于第二门限的元素提取出来,并对提取出的元素按照以下步骤进行判决:
B1、如果sim(rk,r′p′)高于门限值,且矩阵S′中不存在第一个自变量为rp的元素,则将rk对应的连通域判定为目标区域,将r′p′对应的连通域判定为多径区域;
B2、如果同时存在sim(rk,r′p′)和sim(rp,r′q′)高于门限值,且r′p′与rp对应相同的连通域;根据sim(rk,r′p′)将rk对应的连通域判定为目标区域,将r′p′对应的连通域判定为多径区域;然后根据sim(rp,r′q′)将rp对应的连通域判定为目标区域,将r′q′对应的连通域判定为多径区域;最终得到rk对应的连通域与r′p′对应的连通域均为目标区域,r′q′对应的连通为多径区域;
B3、当矩阵S′中所有高于第二门限的元素判定完成之后,将判定为目标区域的连通域保留,同时将判定为多径区域的连通域置零,得到抑制了多径的后的二值化图像If(·);最后把原始图像Io(·)与If(·)进行乘法图像融合得到一幅保留所有目标并且消除了多径的雷达图像Ifinal(·)。
更进一步地,所述第二门限的计算式如下:
其中,j′=1,2,3…N-1。
本发明的有益效果:本发明的一种穿墙雷达成像后多目标侧墙一次多径抑制方法,首先利用4连通域检测方法从二元原始雷达图像中提取出所有连通域;然后对所有连通域编号并依照它们的平均坐标离雷达阵列中心的欧式距离升序排序得到连通域向量;接着把连通域向量元素重组成两个子向量,分别按照到雷达的多径传输路径和直接回波传输路径对这两个子向量进行扩展;然后利用相似度函数计算两个子向量的相似度矩阵;最后使用相似度门限及判决规则的方法判定相似度矩阵的所有元素,得到真实目标所在的连通域并将多径假目标抑制掉。本发明的方法能够有效识别雷达图像中的真实目标和它们的一阶侧墙多径假目标,在抑制多径假目标的同时不会造成真实目标的丢失。相比于现有技术中的多径抑制方法,本发明能够防止正常信噪比条件下目标的丢失,形成高质量的雷达图像,保证了穿墙雷达在密闭建筑物探测环境中的实用效果,为操作人员做出正确的决策提供了有力的保障。
附图说明
图1为侧墙一次多径抑制方法处理流程图。
图2为侧墙一次多径穿墙传播示意图。
图3为具体实施方式中仿真场景示意图。
图4为具体实施方式中原始雷达仿真图像。
图5为具体实施方式中4连通域检测及标记结果图。
图6为具体实施方式中多径抑制后二值化雷达图像。
图7为具体实施方式中最终结果图像。
具体实施方式
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。
如图1所示为本发明的处理流程图,本发明的技术方案为:一种穿墙雷达成像后多目标侧墙一次多径抑制方法,包括:
S1、连通区域检测,根据设置的第一门限对原始雷达图像进行门限检测得到二元图像,对二元图像进行4-连通检测,并标记二元图像中的所有连通域,按照平均坐标值到阵列中心的欧式距离由小到大对所有连通域进行排序,得到一个连通域向量。
设置第一门限THb,用THb对原始雷达图像Io(·)进行门限检测得到二元图像Ib(·),对Ib(·)进行4-连通(同一像素在上、下、左、右四个方向有连通)检测,并标记图像Ib(·)中的所有连通域。假设有N个连通域,计算每个连通域的平均坐标并按照平均坐标值到阵列中心的欧式距离由小到大对所有连通域进行排序,将排好序的连通域以符号D1-DN进行编号,得到一个连通域向量D=[D1,D2,…,DN]。
S2、构造相似度矩阵,将步骤S1得到的连通域向量分成两组,得到第一向量与第二向量;根据侧墙一次多径的形成原理对第一向量进行扩展,得到包含侧墙一次多径传输距离信息的第三向量;根据雷达的直接回波传输路径对第二向量进行扩展,得到包含直接回波传输距离信息的第四向量;根据相似度函数计算得到第三向量与第四向量的相似度矩阵。多径穿墙传播如图2所示,第三向量即由虚线表示的一阶多径,第四向量即由实现表示的直接路径。
将连通域向量D分成两组,具体为:所述第一向量为按照各个连通域的平均坐标值到阵列中心的欧式距离由小到大的顺序连续地在连通域向量中取连通域,所述第二向量为按照平均坐标值到阵列中心的欧式距离由大到小的顺序连续地在连通域向量中取连通域,所述第一向量与第二向量所取连通域个数相等,且第一向量或第二向量所取连通域的个数N′大于连通域向量中总的连通域个数的一半,并且小于总的连通域个数。假设第一组中包含元素D1-DN-1,用向量T=[D1,D2,…,DN-1]表示,即用T表示第一向量;第二组包含元素D2-DN,用向量M=[D2,D3,…,DN]表示,其中Dj代表第j个连通域,即用M表示第二向量;且第一向量T与第二向量M所取连通域个数相等,记为N′个,即在本实施例中N′=N-1。
接下来将向量T和M进行扩展;首先,根据侧墙一次多径的形成原理,计算向量T中每个元素所代表的连通域的平均坐标到各个雷达之间的侧墙一次多径传输距离。如:对于向量T中第i个连通域来说,它的平均坐标为假设雷达总数为K,第k个雷达的位置坐标为(xRk,0),则这个雷达到该连通域的多径传输距离为:
其中,d和ε分别为前墙厚度以及相对介电常数,xw代表的是左侧墙内表面的横坐标值,ΔxG为电磁波在从侧墙反射回接收雷达时在前墙中发生的水平偏移值,ΔxT为电磁波从目标反射回接收雷达时在前墙中发生的水平偏移值,表示向量T中第i个连通域的平均横坐标,表示向量T中第i个连通域的平均纵坐标,xRk表示第k个雷达的横坐标,k=1,2,3,···,K。
这样T中每个连通域对应K个雷达就有K个距离值,将第i个连通域计算得到的K个距离值放进向量ri中。对T中的每一个连通域都进行这样的距离计算后,得到一个包含侧墙一次多径传输距离信息的向量T′,T′的维度为1×(N-1)K,将T′表示如下:
其次,计算M中的每个连通域的平均坐标到K个雷达的直接回波距离,其中直接回波距离的计算公式为:
其中,表示向量M中第i′个连通域的平均横坐标,表示向量M中第i′个连通域的平均纵坐标。
同样地,对M中的每一个元素都进行这样的距离计算后得到一个包含直接回波传输距离信息的向量M′,M′的维度为1×(N-1)K,将M′表示如下:
定义任意两个向量x和y的相似度函数为:
其中,J代表向量中元素的个数。
基于这个相似度函数,T′和M′之间的相似度矩阵可以表示为:
由于多径回波路径始终长于目标回波路径,那么可将相似度矩阵S中的上三角元素全部置零,得到一个下三角矩阵:
S3、元素判决及多径抑制,根据第二门限判定相似度矩阵中的所有元素,得到目标区域。
为了便于后续分析,将矩阵S′的所有非零元素按照取值由大到小依次放入行向量Sv:
假如只考虑侧墙一次多径的影响,那么对于连通域向量D中的N个元素,最多存在N-1个目标,因此选取向量Sv中的前N-1个元素用来计算用于相似度检测的第二门限THsim,第二门限的计算公式为:
其中,j′=1,2,3…N-1。
将相似度矩阵S′中大于第二门限的元素提取出来,对这些提取出来的元素按照下面的判决规则进行判决。
B1、如果sim(rk,r′p′)高于门限值,且矩阵S′中不存在第一个自变量为rp的元素,则将rk对应的连通域判定为目标区域,将r′p′对应的连通域判定为多径区域;
B2、如果同时存在sim(rk,r′p′)和sim(rp,r′q′)高于门限值,且r′p′与rp对应相同的连通域,即r′p′与rp对应的是连通域向量D中的同一个连通域;根据sim(rk,r′p′)将rk对应的连通域判定为目标区域,将r′p′对应的连通域判定为多径区域;然后根据sim(rp,r′q′)将rp对应的连通域判定为目标区域,将r′q′对应的连通域判定为多径区域;最终得到rk、r′p′各自对应的连通域均为目标区域,r′q′对应的连通为多径区域;
B3、当矩阵S′中所有高于第二门限的元素判定完成之后,将判定为目标区域的连通域保留,同时将判定为多径区域的连通域置零,得到抑制了多径后的二值化图像If(·);最后把原始图像Io(·)与If(·)进行乘法图像融合得到一幅保留所有目标并且消除了多径的雷达图像Ifinal(·)。
下面根据一个MATLAB仿真例子给出本发明的具体实施方式。
对存在多个目标的待探测区域,其仿真场景如图3所示,坐标系零点位于前墙前表面的中心位置,31个收发共置雷达紧贴前墙放置,相邻两个雷达之间的水平间隔为0.1m,雷达的发射信号为中心频率1.5GHz、带宽1GHz的步进频信号,频率步进为2MHz。房间内表面的长度和宽度均为6m,四面墙的厚度与相对介电常数相等且为0.24m和7.6。三个目标分别位于(-0.5,3.0)m、(-1.5,4.0)m和(-2.4,3.9)m处,值得注意的是目标3与目标1的侧墙一次多径重合,原始雷达仿真图像如图4所示。
S1、使用门限THb对原始雷达仿真图像进行二值化检测,门限值取为雷达图像最大值像素点值的0.4倍,在这里它的取值为1.9134。门限检测后得到二值化的图像Ib(·)如图5所示。通过4-连通域检测的方法对图像Ib(·)进行检测并将所有连通域标记,可以得到五个连通域,其中,连通域1包含104个像素点,这些像素点的坐标分别为(-3.1,4.065)m,(-3.085,4.065)m…(-3.055,4.245)m,然后根据这些像素点坐标,计算求得连通域1的均值坐标为(-3.085,4.14)m,根据仿真场景示意图所建立的坐标系,雷达阵列中心位于(0,0)m,那么可以求得连通域1到雷达阵列中心点的距离为5.1630m。同理,可以得到其他连通域的坐标均值与到雷达阵列中心的距离,如表1所示:
表1其他连通域的坐标均值与到雷达阵列中心的距离
连通域编号 | x平均 | y平均 | 距离 |
连通域2 | -2.9350 | 4.1400 | 5.1630 |
连通域3 | -2.5150 | 4.1250 | 4.8312 |
连通域4 | -1.5550 | 4.2150 | 4.4927 |
连通域5 | -0.5350 | 3.2250 | 3.2691 |
将所有连通域按照其平均坐标到雷达阵列中心的欧式距离由小到大进行排序,得到五个连通域,依次编号为D1-D5,如图5所示。然后将这些连通域标号放进一个连通域向量D中。则有下面的对应关系:D1对应连通域5,D2对应连通域4,D3对应连通域3,D4对应连通域1,D5对应连通域2。
S2、将向量D分成向量T和向量M,其中向量T包含连通域D1-D4,向量M包含连通域D2-D5。对向量T进行扩展,计算T中的每一个连通域的中心坐标到31个雷达的侧墙一次多径传输距离。这样T的每一个元素就扩展成了一个1×31维的向量ri,i=1,2,3,4,扩展之后的向量记作T′=[r1r2r3r4]。同理,对向量M进行扩展,计算M中的每一个连通域的中心坐标到31个雷达的直接路径传输距离,则M的每一个元素也扩展成了一个1×31维的向量r′i,i=2,3,4,5,M向量扩展之后记作M′=[r′1r′2r′3r′4]。根据前面定义的相似度函数对向量T′和M′进行相似度计算,并且得到一个4×4维的相似度矩阵S,S的值如表2所示:
表2相似度矩阵S的值
D1 | D2 | D3 | D4 | |
D2 | 0.6177 | 0.3878 | 0.4665 | 0.4687 |
D3 | 0.9840 | 0.5106 | 0.6418 | 0.6449 |
D4 | 0.6110 | 0.7750 | 0.9150 | 0.5110 |
D5 | 0.4850 | 0.8990 | 0.6660 | 0.6630 |
因为目标的多径传播距离总是长于它的直接路径传播距离,则相似度矩阵S中对角线以上的元素全部近似为0后得到一个下三角矩阵S′,S′的值如表3所示:
表3矩阵S′的值
D1 | D2 | D3 | D4 | |
D2 | 0.6177 | 0 | 0 | 0 |
D3 | 0.9840 | 0.5106 | 0 | 0 |
D4 | 0.6114 | 0.7750 | 0.9150 | 0 |
D5 | 0.4853 | 0.8990 | 0.6660 | 0.6635 |
S3、元素判决及多径抑制
将下三角矩阵S′中的元素按值的大小进行降序排序得到一个1×10维的行向量Sv,
然后对该向量的前4个元素进行求均值计算后得到相似度检测门限THsim的值为0.8932,将矩阵S′中大于门限THsim的元素提取出来,它们依次为0.9840、0.9150和0.8990。所对应的连通域相似度函数为sim(r1,r′3)、sim(r2,r′5)和sim(r3,r′4),由判决准则可得连通域D1、D2和D3为真实目标,D4和D5分别为连通域D3和D2的侧墙一次多径假目标,连通域D1的侧墙一次多径假目标与连通域D3重合。经过多径抑制后的二值化雷达图像如图6所示,原始雷达仿真图像与多径抑制后的二值化雷达图像的非相干融合结果如图7所示。
由仿真结果可知,本发明提供的适用于穿墙雷达成像后多目标侧墙一次多径抑制方法不仅能够有效剔除多径,而且不会造成目标的丢失,验证了本发明的正确性和有效性。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (7)
1.一种穿墙雷达成像后多目标侧墙一次多径抑制方法,其特征在于,包括:
S1、连通区域检测,根据设置的第一门限对原始雷达图像进行门限检测得到二元图像,对二元图像进行4-连通检测,并标记二元图像中的所有的N个连通域,按照平均坐标值到阵列中心的欧式距离由小到大对所有的N个连通域进行排序,得到一个连通域向量;
S2、构造相似度矩阵,将步骤S1得到的连通域向量分成两组,得到第一向量与第二向量;根据侧墙一次多径的形成原理对第一向量进行扩展,得到包含侧墙一次多径传输距离信息的第三向量;根据雷达的直接回波传输路径对第二向量进行扩展,得到包含直接回波传输距离信息的第四向量;根据相似度函数计算得到第三向量与第四向量的相似度矩阵;
S3、元素判决及多径抑制,根据第二门限判定相似度矩阵中的所有元素,得到目标区域。
2.根据权利要求1所述的一种穿墙雷达成像后多目标侧墙一次多径抑制方法,其特征在于,步骤S2中所述将步骤S1得到的连通域向量分成两组,具体为:所述第一向量为按照平均坐标值到阵列中心的欧式距离由小到大的顺序连续地在连通域向量中取连通域,所述第二向量为按照平均坐标值到阵列中心的欧式距离由大到小的顺序连续地在连通域向量中取连通域,所述第一向量与第二向量所取连通域个数相等,记为N′,N′大于连通域向量中总的连通域个数的一半,并且小于总的连通域个数。
3.根据权利要求2所述的一种穿墙雷达成像后多目标侧墙一次多径抑制方法,其特征在于,步骤S2中所述的根据侧墙一次多径的形成原理对第一向量进行扩展,得到包含侧墙一次多径传输距离信息的第三向量;具体为:计算第一向量中每个连通域的平均坐标到K个雷达位置的侧墙一次多径传输距离,每个连通域得到的K个距离值组成一个向量ri;i=1,2,3…,N′,根据N′个向量ri得到包含侧墙一次多径传输距离信息的第三向量。
4.根据权利要求3所述的一种穿墙雷达成像后多目标侧墙一次多径抑制方法,其特征在于,步骤S2中所述的根据雷达的直接回波传输路径对第二向量进行扩展,得到包含直接回波传输距离信息的第四向量;具体为:计算第二向量中每个连通域的平均坐标到K个雷达位置的直接回波距离,每个连通域得到的K个距离值组成一个向量ri′′;i′=1,2,3…,N′,根据N′个向量ri′′得到包含N′个直接回波传输距离信息的第四向量。
5.根据权利要求4所述的一种穿墙雷达成像后多目标侧墙一次多径抑制方法,其特征在于,步骤S2中所述的根据相似度函数计算得到第三向量与第四向量的相似度矩阵,具体为:根据相似度函数定义,第三向量与第四向量的相似度矩阵表达式为:
其中,sim(r1,r2′)表示计算向量r1和r2′之间相似度的相似度计算函数,其取值范围为[0,1]。
6.根据权利要求5所述的一种穿墙雷达成像后多目标侧墙一次多径抑制方法,其特征在于,步骤S3所述根据第二门限判定相似度矩阵中的所有元素,得到目标区域,具体为:
A1、对相似度矩阵S进行预处理,根据多径回波路径始终长于目标回波路径的特点,将相似度矩阵S中的上三角元素全部置零,得到一个下三角矩阵S′
A2、将矩阵S′中的所有非零元素按照取值从大到小依次放入行向量Sv:
A3、假如只考虑侧墙一次多径的影响,那么对于连通域向量D中的N个元素,最多存在N-1个目标,因此选取向量Sv中的前N-1个元素用来计算用于相似度检测的第二门限;
A4、将矩阵S′中大于第二门限的元素提取出来,并对提取出的元素按照以下步骤进行判决:
B1、如果sim(rk,rp′′)高于门限值,且矩阵S′中不存在第一个自变量为rp的元素,则将rk对应的连通域判定为目标区域,将rp′′对应的连通域判定为多径区域;
B2、如果同时存在sim(rk,rp′′)和sim(rp,rq′′)高于门限值,且rp′′与rp对应相同的连通域;根据sim(rk,rp′′)将rk对应的连通域判定为目标区域,将rp′′对应的连通域判定为多径区域;然后根据sim(rp,rq′′)将rp对应的连通域判定为目标区域,将rq′′对应的连通域判定为多径区域;最终得到rk对应的连通域与rp′′对应的连通域均为目标区域,rq′′对应的连通域 为多径区域;
B3、当矩阵S′中所有高于第二门限的元素判定完成之后,将判定为目标区域的连通域保留,同时将判定为多径区域的连通域置零,得到抑制了多径后的二值化图像If(·);最后把原始图像Io(·)与If(·)进行乘法图像融合得到一幅保留所有目标并且消除了多径的雷达图像Ifinal(·)。
7.根据权利要求6所述的一种穿墙雷达成像后多目标侧墙一次多径抑制方法,其特征在于,所述第二门限的计算式如下:
其中,j′=1,2,3…N-1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610698090.5A CN106093933B (zh) | 2016-08-22 | 2016-08-22 | 一种穿墙雷达成像后多目标侧墙一次多径抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610698090.5A CN106093933B (zh) | 2016-08-22 | 2016-08-22 | 一种穿墙雷达成像后多目标侧墙一次多径抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106093933A CN106093933A (zh) | 2016-11-09 |
CN106093933B true CN106093933B (zh) | 2018-05-25 |
Family
ID=58070406
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610698090.5A Active CN106093933B (zh) | 2016-08-22 | 2016-08-22 | 一种穿墙雷达成像后多目标侧墙一次多径抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106093933B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6584369B2 (ja) * | 2015-08-20 | 2019-10-02 | 三菱電機株式会社 | 壁の背後に位置するシーンの雑音のない画像を求めるシステムおよび方法 |
CN107219523B (zh) * | 2017-05-22 | 2020-03-27 | 电子科技大学 | 一种穿墙雷达三维成像后多目标后墙一次多径抑制方法 |
CN109581317B (zh) * | 2018-12-24 | 2022-05-31 | 电子科技大学 | 一种基于回波峰值匹配的角落目标定位方法 |
CN111699408A (zh) * | 2019-05-31 | 2020-09-22 | 深圳市大疆创新科技有限公司 | 隧道场景检测方法及毫米波雷达 |
CN112213725B (zh) * | 2020-09-28 | 2022-10-25 | 森思泰克河北科技有限公司 | 车载雷达的多径虚警抑制方法、装置及终端设备 |
CN113109816B (zh) * | 2021-03-29 | 2023-10-03 | 广东工业大学 | 雷达回波图像的回波块追踪方法、装置及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103197289A (zh) * | 2013-04-03 | 2013-07-10 | 电子科技大学 | 一种适用于穿墙雷达成像的多径幻象抑制方法 |
CN104062638A (zh) * | 2014-06-18 | 2014-09-24 | 电子科技大学 | 一种用于穿墙雷达多目标检测的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9971019B2 (en) * | 2013-07-22 | 2018-05-15 | Mitsubishi Electric Research Laboratories, Inc. | System and method for through-the-wall-radar-imaging using total-variation denoising |
-
2016
- 2016-08-22 CN CN201610698090.5A patent/CN106093933B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103197289A (zh) * | 2013-04-03 | 2013-07-10 | 电子科技大学 | 一种适用于穿墙雷达成像的多径幻象抑制方法 |
CN104062638A (zh) * | 2014-06-18 | 2014-09-24 | 电子科技大学 | 一种用于穿墙雷达多目标检测的方法 |
Non-Patent Citations (1)
Title |
---|
一种适用于穿墙雷达成像的墙体补偿算法;姚雪 等;《雷达科学与技术》;20141231;第12卷(第6期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN106093933A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106093933B (zh) | 一种穿墙雷达成像后多目标侧墙一次多径抑制方法 | |
CN109188470B (zh) | 一种基于卷积神经网络的gnss欺骗干扰检测方法 | |
CN111123257B (zh) | 基于图时空网络的雷达动目标多帧联合检测方法 | |
CN107589407B (zh) | 多发多收穿墙雷达成像后前后墙振荡一次多径抑制方法 | |
CN111582089B (zh) | 基于卫星红外与可见光图像的海上目标信息融合方法 | |
CN109902715B (zh) | 一种基于上下文聚合网络的红外弱小目标检测方法 | |
CN104062638B (zh) | 一种用于穿墙雷达多目标检测的方法 | |
CN105549011B (zh) | 一种基于mimo穿墙雷达的建筑物单边两点斜视成像方法 | |
CN103197302B (zh) | 一种适用于穿墙雷达成像的目标位置提取方法 | |
CN106324595B (zh) | 一种基于多角度探测成像的多径假目标抑制方法 | |
US20070076917A1 (en) | Target detection improvements using temporal integrations and spatial fusion | |
CN109597065B (zh) | 一种用于穿墙雷达检测的虚警抑制方法、装置 | |
CN103745216B (zh) | 一种基于空域特性的雷达图像杂波抑制方法 | |
CN108983165B (zh) | 基于子站选择的多站雷达系统抗欺骗式干扰方法 | |
CN113376613B (zh) | 用于雷达探测的恒虚警检测方法、装置及电子设备 | |
CN111045006B (zh) | 一种基于多成像字典融合的拐角隐蔽目标成像方法 | |
Li et al. | Scale-adaptive human target tracking for through-wall imaging radar | |
CN115359258B (zh) | 一种成分不确定度测量的弱小目标检测方法及系统 | |
CN107219523B (zh) | 一种穿墙雷达三维成像后多目标后墙一次多径抑制方法 | |
Du et al. | NLOS target localization with an L-band UWB radar via grid matching | |
CN105866748B (zh) | 一种基于检测先验的固定窗长恒虚警检测方法 | |
CN113296089B (zh) | 用于多预警机目标跟踪系统的lmb密度融合方法及装置 | |
CN107219522B (zh) | 一种椭圆-双曲线联合的穿墙雷达目标定位方法 | |
CN113960558B (zh) | 基于多输入多输出雷达的非视距目标定位方法及系统 | |
Giacometti et al. | Localization of radar emitters from a single sensor using multipath and TDOA-AOA measurements in a naval context |
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 |