CN107479043A - 基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法 - Google Patents
基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法 Download PDFInfo
- Publication number
- CN107479043A CN107479043A CN201710694893.8A CN201710694893A CN107479043A CN 107479043 A CN107479043 A CN 107479043A CN 201710694893 A CN201710694893 A CN 201710694893A CN 107479043 A CN107479043 A CN 107479043A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msqrt
- msup
- 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.)
- Granted
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
- 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
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9027—Pattern recognition for feature extraction
-
- 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
-
- 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
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9004—SAR image acquisition techniques
- G01S13/9017—SAR image acquisition techniques with time domain processing of the SAR signals in azimuth
-
- 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
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
-
- 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/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
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)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,应用于合成孔径穿墙雷达一次多径抑制技术领域,首先,根据建筑物结构特征、穿墙雷达离前墙位置以及各反射路径回波,建立目标直接回波字典与多径假目标字典;其次,对回波数据,分别使用目标直接回波字典与多径假目标字典进行后向投影成像处理,得到目标空间图像和多径假目标空间图像;最后,将目标空间图像和多径假目标空间图像进行非相干融合,实现建筑物内一次多径假目标的抑制;相比现有技术复杂度更低,多径抑制效果较好。
Description
技术领域
本发明属于合成孔径穿墙雷达一次多径抑制技术领域,特别涉及合成孔径穿墙雷达一次多径抑制假目标剔除技术。
背景技术
利用穿墙雷达系统对封闭空间内的隐蔽目标进行探测时,电磁波除了直接通过目标反射回接收天线外,还将在建筑物墙体内表面、地面发生反射,最后传输到接收天线,在雷达图像中则表现为多径假目标,导致虚警,严重影响了穿墙雷达的探测性能。
合成孔径穿墙雷达领域的多径可按电磁波在建筑物内表面的反射情况进行分类,主要包括:前墙振荡多径、侧墙一次多径、侧墙二次多径、后墙一次多径、后墙二次多径、前后墙振荡多径。二次多径及前后墙振荡多径分布在成像区域之外,不影响房间内隐蔽目标的检测;前墙振荡多径附着于目标图像区域后方,可通过图像域滤波的方法解决。本发明主要解决侧墙一次与后墙一次多径假目标的抑制问题,也适用于其他多径假目标的剔除。
国内外许多研究机构开展了基于合成孔径穿墙雷达的多径建模及假目标剔除方法的研究。其中,美国维纳诺瓦大学大学提出了一种基于压缩感知的多径抑制方法(M.Leigsnering,F.Ahmad,M.Amin,and A.Zoubir,“Multipath exploitation inthrough-the-wall radar imaging using sparse reconstruction,”IEEE Transactionson Geoscience and Remote Sensing,vol.50,no.2,April 2014,pp.920–939.),它将原始数据稀疏化后,利用最优化求解方法进行成像场景恢复,可实现多径的抑制。该方法需要求解复杂的最优化问题,会耗费大量的计算资源及运行时间,不能够实时处理。法赫德国王石油矿产大学提出了一种将子孔径成像与压缩感知结合的多径抑制方法(A.Abdalla,A.Muqaibel,and A.Suhail,“Aspect dependent multipath ghost suppression in TWRIunder compressive sensing framework,”in Proceedings of the InternationalConference on Communication,Signal Processing,and their Application,Sharjah,February 2015,pp.1–6),该方法虽然能够减少运算数据量及抑制多径假目标,但是需进行雷达选择和优化求解,会导致操作复杂并且不能满足实时性的要求。
发明内容
为解决上述技术问题,本发明提出了一种基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,通过建立目标直接回波字典矩阵与多径字典矩阵;并将两类字典用于后向投影成像,得到目标空间图像和多径假目标空间图像;通过将这些图像进行非相干融合,实现建筑物内一次多径假目标的抑制。
本发明采用的技术方案为:基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,首先,根据建筑物结构特征、穿墙雷达离前墙位置以及各反射路径回波,建立目标直接回波字典矩阵与多径假目标字典矩阵;其次,对回波数据,分别使用目标直接回波字典矩阵与多径假目标字典矩阵进行后向投影成像处理,得到目标空间图像和多径假目标空间图像;最后,将目标空间图像和多径假目标空间图像进行非相干融合,得到抑制多径后的图像。
进一步地,包括以下步骤:
S1、根据穿墙雷达离前墙位置、建筑物布局信息以及目标直接反射路径回波,建立目标直接回波字典;
根据穿墙雷达离前墙位置、建筑物布局信息以及左侧墙一次多径回波、后墙一次多径回波、右侧墙一次多径回波,建立多径假目标字典;
S2、利用穿墙雷达对建筑物内隐蔽目标进行探测,得到回波数据;
S3、对步骤S2得到的回波数据,根据目标直接回波字典和多径假目标字典,分别进行后向投影成像处理;得到目标直接回波字典以及多径假目标字典对应的图像;
S4、分别对目标直接回波字典以及多径假目标字典对应的图像进行归一化处理,得到归一化后的目标直接回波字典以及多径假目标字典对应的图像;
S5、将步骤S4得到的归一化后的目标直接回波字典以及多径假目标字典对应的图像,进行非相干融合;再对得到的非相干融合图像进行归一化处理,得到抑制多径后的图像。
更进一步地,所述多径假目标字典包括:左侧墙一次多径假目标字典、后墙一次多径假目标字典以及右侧墙一次多径假目标字典。
更进一步地,步骤S5所述的非相干融合公式为:
其中,|.|为取绝对值操作,Ifusion(x,y)为非相干融合图像中位于(x,y)的像素点的取值,为目标直接回波字典对应的图像位于(x,y)的像素点的归一化取值,为左侧墙一次多径假目标字典对应的图像位于(x,y)的像素点的归一化取值,为后墙一次多径假目标字典对应的图像位于(x,y)的像素点的归一化取值,为右侧墙一次多径假目标字典对应的图像位于(x,y)的像素点的归一化取值。
本发明的有益效果:本发明的基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,通过直接反射路径回波、左侧墙一次多径回波、后墙一次多径回波、右侧墙一次多径回波、以及穿墙雷达离前墙位置和建筑物布局信息建立目标直接回波字典、左侧墙一次多径字典、后墙一次多径字典以及右侧墙一次多径字典;然后利用多幅图像连乘实现多径抑制;本申请的方法相比现有技术复杂度更低,多径抑制效果较好。
附图说明
图1为封闭空间一次多径剔除方法处理流程图。
图2为具体实施方式中仿真场景示意图。
图3为具体实施方式中四个成像空间的归一化图像。
图4为具体实施方式中多径抑制后的最终图像。
具体实施方式
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。
本实施例设置如图2所示的仿真场景,其中,坐标系零点位于前墙前表面中心位置,四面墙体的厚度均为0.24m,介电常数为8.6,房间的内部范围大小为6m×5.52m,前墙前表面位于y=0m处,左侧墙内表面位于xl=-3m处,右侧墙内表面位于xr=3m处,后墙内表面位于yB=5.76m处。
仿真中,发射信号为步进变频连续波,信号频率范围为1.1GHz-2.1GHz,频率步进为2MHz。单发单收穿墙雷达从左往右依次放置在31个位置上,相邻位置的横坐标间隔为0.1m,所有位置的纵坐标为0m。最左边位置为(-1.5m,0m),最右边位置为(1.5m,0m),31个位置的横轴中心与坐标原点重合。
如图1所示为本申请的方案流程图,本发明的技术方案为:基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,首先,根据建筑物结构特征、穿墙雷达离前墙位置以及各反射路径回波,建立目标直接回波字典矩阵与多径字典矩阵;其次,对回波数据,分别使用目标直接回波字典矩阵与多径字典矩阵进行后向投影成像处理,得到目标空间图像和多径假目标空间图像;最后,将目标空间图像和多径假目标空间图像进行非相干融合,实现建筑物内一次多径假目标的抑制。具体包括以下步骤:
S1、根据穿墙雷达离前墙位置、建筑物布局信息以及目标直接反射路径回波,建立目标直接回波字典AT;
根据穿墙雷达离前墙位置、建筑物布局信息以及左侧墙一次多径回波、后墙一次多径回波、右侧墙一次多径回波,建立多径假目标字典AJ,i,i=1,2,3;即,AJ,1表示左侧墙一次多径字典,AJ,2表示后墙一次多径字典,AJ,3表示右侧墙一次多径字典;
S2、利用穿墙雷达对建筑物内隐蔽目标进行探测成像,得到回波数据;
假设得到的回波数据表示为{s(n,t),n=1,2,...,N},每幅图像包含201×201个像素点,具体是:方位向201个,距离向201个,因此共有201×201个像素点;N为穿墙雷达系统放置的位置数。
S3、对步骤S2得到的回波数据,根据步骤S1得到的目标直接回波字典AT和多径假目标字典AJ,i,i=1,2,3,分别进行后向成影成像处理(即BP成像);得到目标直接回波字典以及多径假目标字典对应的图像;即对应的四幅图像IT、IJ,1、IJ,2、IJ,3;
对应的四幅图像IT、IJ,1、IJ,2、IJ,3,位于(x,y)的像素点的取值计算公式为:
其中,s(n,t)为利用穿墙雷达对建筑物内隐蔽目标进行探测成像得到的回波数据,t为时间变量,N为穿墙雷达系统放置的位置数,n=1,2,...,N,τT,n为目标直接反射路径回波聚焦延时、τJ,1,n为左侧墙一次多径回波聚焦延时、τJ,2,n为后墙一次多径回波聚焦延时和τJ,3,n为右侧墙一次多径回波聚焦延时。
τT,n的计算公式如下所示:
其中,d和ε分别为前墙厚度以及相对介电常数,ΔxT为电磁波从目标反射回接收雷达时在前墙中发生的水平偏移值,N为穿墙雷达系统放置的位置数,(Xn,Yn)为穿墙雷达系统放置的第n个位置处对应的坐标,c为光速。
τJ,1,n的计算公式如下所示:
其中,d和ε分别为前墙厚度以及相对介电常数,ΔxJ,1为电磁波在从左侧墙反射回接收雷达时在前墙中发生的水平偏移值,xl为左侧墙内表面的x轴坐标,N为穿墙雷达系统放置的位置数,(Xn,Yn)为穿墙雷达系统放置的第n个位置处对应的坐标,c为光速。
τJ,2,n的计算公式如下所示:
其中,d和ε分别为前墙厚度以及相对介电常数,ΔxJ,2为电磁波在从后墙反射回接收雷达时在前墙中发生的水平偏移值,yB为后墙内表面的y轴坐标,N为穿墙雷达系统放置的位置数,(Xn,Yn)为穿墙雷达系统放置的第n个位置处对应的坐标,c为光速。
τJ,3,n的计算公式如下所示:
其中,d和ε分别为前墙厚度以及相对介电常数,ΔxJ,3为电磁波在从右侧墙反射回接收雷达时在前墙中发生的水平偏移值,xr为右侧墙内表面的x轴坐标,N为穿墙雷达系统放置的位置数,(Xn,Yn)为穿墙雷达系统放置的第n个位置处对应的坐标,c为光速。
S4、分别对目标直接回波字典以及多径假目标字典对应的图像进行归一化处理,得到归一化后的目标直接回波字典以及多径假目标字典对应的图像;即如图3所示,四个图中,真实目标位置所在区域像素点的强度值均较大。
图像归一化时,位于(x,y)的像素点的归一化取值为:
其中,IT为目标直接回波字典对应的图像,IJ,1为左侧墙一次多径假目标字典对应的图像,IJ,2为后墙一次多径假目标字典对应的图像,IJ,3为右侧墙一次多径假目标字典对应的图像,IT(x,y)为目标直接回波字典对应的图像位于(x,y)的像素点取值,IJ,1(x,y)为左侧墙一次多径假目标字典对应的图像位于(x,y)的像素点取值,IJ,2(x,y)为后墙一次多径假目标字典对应的图像位于(x,y)的像素点取值,IJ,3(x,y)为右侧墙一次多径假目标字典对应的图像位于(x,y)的像素点取值,max(IT)为IT中的最大像素值,max(IJ,1)为IJ,1中的最大像素值,max(IJ,2)为IJ,2中的最大像素值,max(IJ,3)为图像IJ,3中的最大像素值。
S5、将步骤S4得到的归一化后的目标直接回波字典以及多径假目标字典对应的图像,进行非相干融合;再对得到的非相干融合图像Ifusion进行归一化处理,得到抑制多径后的图像Ifinal,如图4所示,从图中可以看到,三个多径假目标均被抑制了。
非相干图像融合的计算公式为:
其中,*为乘法运算符号,|.|为取绝对值操作,Ifusion(x,y)为非相干融合图像Ifusion中位于(x,y)的像素点的取值,对图像Ifusion的归一化处理表示为:
其中,Ifinal(x,y)为归一化后的非相干融合后的图像Ifinal中位于(x,y)的像素点的取值,max(Ifusion)为Ifusion的最大像素值。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (10)
1.基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,其特征在于,包括:首先,根据建筑物结构特征、穿墙雷达离前墙位置以及各反射路径回波,建立目标直接回波字典与多径假目标字典;其次,对回波数据,分别使用目标直接回波字典与多径假目标字典进行后向投影成像处理,得到目标空间图像和多径假目标空间图像;最后,将目标空间图像和多径假目标空间图像进行非相干融合,得到抑制多径后的图像。
2.根据权利要求1所述的基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,其特征在于,包括以下步骤:
S1、根据穿墙雷达离前墙位置、建筑物布局信息以及目标直接反射路径回波,建立目标直接回波字典;
根据穿墙雷达离前墙位置、建筑物布局信息以及左侧墙一次多径回波、后墙一次多径回波、右侧墙一次多径回波,建立多径假目标字典;
S2、利用穿墙雷达对建筑物内隐蔽目标进行探测,得到回波数据;
S3、对步骤S2得到的回波数据,根据目标直接回波字典和多径假目标字典,分别进行后向投影成像处理;得到目标直接回波字典以及多径假目标字典对应的图像;
S4、分别对目标直接回波字典以及多径假目标字典对应的图像进行归一化处理,得到归一化后的目标直接回波字典以及多径假目标字典对应的图像;
S5、将步骤S4得到的归一化后的目标直接回波字典以及多径假目标字典对应的图像,进行非相干融合;再对得到的非相干融合图像进行归一化处理,得到抑制多径后的图像。
3.根据权利要求2所述的基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,其特征在于,所述多径假目标字典包括:左侧墙一次多径假目标字典、后墙一次多径假目标字典以及右侧墙一次多径假目标字典。
4.根据权利要求3所述的基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,其特征在于,步骤S5所述的非相干融合公式为:
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
<mi>s</mi>
<mi>i</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>|</mo>
<msub>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mi>T</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>*</mo>
<mo>|</mo>
<msub>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>*</mo>
<mo>|</mo>
<msub>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>*</mo>
<mo>|</mo>
<msub>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
其中,|.|为取绝对值操作,Ifusion(x,y)为非相干融合图像中位于(x,y)的像素点的取值,为目标直接回波字典对应的图像位于(x,y)的像素点的归一化取值,为左侧墙一次多径假目标字典对应的图像位于(x,y)的像素点的归一化取值,为后墙一次多径假目标字典对应的图像位于(x,y)的像素点的归一化取值,为右侧墙一次多径假目标字典对应的图像位于(x,y)的像素点的归一化取值。
5.根据权利要求4所述的基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,其特征在于,所述以及各自的计算式为:
<mrow>
<msub>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mi>T</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>I</mi>
<mi>T</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>max</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>max</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
1
<mrow>
<msub>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>max</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>max</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中,IT为目标直接回波字典对应的图像,IJ,1为左侧墙一次多径假目标字典对应的图像,IJ,2为后墙一次多径假目标字典对应的图像,IJ,3为右侧墙一次多径假目标字典对应的图像,IT(x,y)为目标直接回波字典对应的图像位于(x,y)的像素点取值,IJ,1(x,y)为左侧墙一次多径假目标字典对应的图像位于(x,y)的像素点取值,IJ,2(x,y)为后墙一次多径假目标字典对应的图像位于(x,y)的像素点取值,IJ,3(x,y)为右侧墙一次多径假目标字典对应的图像位于(x,y)的像素点取值,max(IT)为IT中的最大像素值,max(IJ,1)为IJ,1中的最大像素值,max(IJ,2)为IJ,2中的最大像素值,max(IJ,3)为图像IJ,3中的最大像素值。
6.根据权利要求5所述的基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,其特征在于,所述IT(x,y)、IJ,1(x,y)、IJ,2(x,y)以及IJ,3(x,y)各自的计算式为:
<mrow>
<msub>
<mi>I</mi>
<mi>T</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>3</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
其中,s(n,t)为利用穿墙雷达对建筑物内隐蔽目标进行探测成像得到的回波数据,t为时间变量,N为穿墙雷达系统放置的位置数,n=1,2,...,N,τT,n为目标直接反射路径回波聚焦延时,τJ,1,n为左侧墙一次多径回波聚焦延时,τJ,2,n为后墙一次多径回波聚焦延时,τJ,3,n为右侧墙一次多径回波聚焦延时。
7.根据权利要求6所述的基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,其特征在于,所述τT,n的计算公式如下所示:
<mrow>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mi>&epsiv;</mi>
</msqrt>
<msqrt>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>&Delta;x</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>+</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;x</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
</mrow>
其中,d和ε分别为前墙厚度以及相对介电常数,ΔxT为电磁波从目标反射回接收雷达时在前墙中发生的水平偏移值,N为穿墙雷达系统放置的位置数,(Xn,Yn)为穿墙雷达系统放置的第n个位置处对应的坐标,c为光速。
8.根据权利要求6所述的基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,其特征在于,所述τJ,1,n的计算公式如下所示:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mi>&epsiv;</mi>
</msqrt>
<msqrt>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>&Delta;x</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>+</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;x</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mi>&epsiv;</mi>
</msqrt>
<msqrt>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>&Delta;x</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>+</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>x</mi>
<mi>l</mi>
</msub>
<mo>-</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;x</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,d和ε分别为前墙厚度以及相对介电常数,ΔxJ,1为电磁波在从左侧墙反射回接收雷达时在前墙中发生的水平偏移值,xl为左侧墙内表面的x轴坐标,N为穿墙雷达系统放置的位置数,(Xn,Yn)为穿墙雷达系统放置的第n个位置处对应的坐标,c为光速。
9.根据权利要求6所述的基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,其特征在于,所述τJ,2,n的计算公式如下所示:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mi>&epsiv;</mi>
</msqrt>
<msqrt>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>&Delta;x</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>+</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;x</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mi>&epsiv;</mi>
</msqrt>
<msqrt>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>&Delta;x</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>+</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;x</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>y</mi>
<mi>B</mi>
</msub>
<mo>-</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,d和ε分别为前墙厚度以及相对介电常数,ΔxJ,2为电磁波在从后墙反射回接收雷达时在前墙中发生的水平偏移值,yB为后墙内表面的y轴坐标,N为穿墙雷达系统放置的位置数,(Xn,Yn)为穿墙雷达系统放置的第n个位置处对应的坐标,c为光速。
10.根据权利要求6所述的基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法,其特征在于,所述τJ,3,n的计算公式如下所示:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>3</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mi>&epsiv;</mi>
</msqrt>
<msqrt>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>&Delta;x</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>+</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;x</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mi>&epsiv;</mi>
</msqrt>
<msqrt>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>&Delta;x</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>+</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>-</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;x</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,d和ε分别为前墙厚度以及相对介电常数,ΔxJ,3为电磁波在从右侧墙反射回接收雷达时在前墙中发生的水平偏移值,xr为右侧墙内表面的x轴坐标,N为穿墙雷达系统放置的位置数,(Xn,Yn)为穿墙雷达系统放置的第n个位置处对应的坐标,c为光速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710694893.8A CN107479043B (zh) | 2017-08-15 | 2017-08-15 | 基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710694893.8A CN107479043B (zh) | 2017-08-15 | 2017-08-15 | 基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107479043A true CN107479043A (zh) | 2017-12-15 |
CN107479043B CN107479043B (zh) | 2020-10-30 |
Family
ID=60600471
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710694893.8A Active CN107479043B (zh) | 2017-08-15 | 2017-08-15 | 基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107479043B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108828592A (zh) * | 2018-04-25 | 2018-11-16 | 中国科学院电子学研究所 | 基于mimo矩形平面阵列的方位向成像方法 |
CN109061622A (zh) * | 2018-06-22 | 2018-12-21 | 电子科技大学 | 一种基于毫米波雷达的隐蔽目标多径探测方法 |
CN109298418A (zh) * | 2018-09-30 | 2019-02-01 | 湖南华诺星空电子技术有限公司 | 基于建筑物内部结构特征的雷达检测虚警抑制方法及装置 |
CN109581317A (zh) * | 2018-12-24 | 2019-04-05 | 电子科技大学 | 一种基于回波峰值匹配的角落目标定位方法 |
CN109597065A (zh) * | 2018-12-11 | 2019-04-09 | 湖南华诺星空电子技术有限公司 | 一种用于穿墙雷达检测的虚警抑制方法、装置 |
CN109696671A (zh) * | 2018-12-29 | 2019-04-30 | 南京汇君半导体科技有限公司 | 基于组结构与相关性学习的高分辨穿墙雷达成像方法 |
CN111045006A (zh) * | 2019-12-23 | 2020-04-21 | 电子科技大学 | 一种基于多成像字典融合的拐角隐蔽目标成像方法 |
CN112698326A (zh) * | 2020-11-30 | 2021-04-23 | 浙江华消科技有限公司 | 雷达穿墙探测对象的方法、雷达探测设备及电子装置 |
CN114066792A (zh) * | 2021-11-16 | 2022-02-18 | 沈阳航空航天大学 | 一种基于多分辨率融合卷积神经网络的穿墙雷达成像方法 |
CN114066792B (zh) * | 2021-11-16 | 2024-10-25 | 沈阳航空航天大学 | 一种基于多分辨率融合卷积神经网络的穿墙雷达成像方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130222172A1 (en) * | 2012-02-28 | 2013-08-29 | L-3 Communications Cyterra Corporation | Determining penetrability of a barrier |
CN103995256A (zh) * | 2014-05-29 | 2014-08-20 | 电子科技大学 | 一种适用于穿墙雷达成像的多面墙体快速补偿方法 |
CN106226765A (zh) * | 2016-09-12 | 2016-12-14 | 桂林电子科技大学 | 一种建筑物布局成像方法及系统 |
CN106324595A (zh) * | 2016-09-14 | 2017-01-11 | 电子科技大学 | 一种基于多角度探测成像的多径假目标抑制方法 |
CN106772365A (zh) * | 2016-11-25 | 2017-05-31 | 南京理工大学 | 一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法 |
-
2017
- 2017-08-15 CN CN201710694893.8A patent/CN107479043B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130222172A1 (en) * | 2012-02-28 | 2013-08-29 | L-3 Communications Cyterra Corporation | Determining penetrability of a barrier |
CN103995256A (zh) * | 2014-05-29 | 2014-08-20 | 电子科技大学 | 一种适用于穿墙雷达成像的多面墙体快速补偿方法 |
CN106226765A (zh) * | 2016-09-12 | 2016-12-14 | 桂林电子科技大学 | 一种建筑物布局成像方法及系统 |
CN106324595A (zh) * | 2016-09-14 | 2017-01-11 | 电子科技大学 | 一种基于多角度探测成像的多径假目标抑制方法 |
CN106772365A (zh) * | 2016-11-25 | 2017-05-31 | 南京理工大学 | 一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法 |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108828592B (zh) * | 2018-04-25 | 2020-12-22 | 中国科学院电子学研究所 | 基于mimo矩形平面阵列的方位向成像方法 |
CN108828592A (zh) * | 2018-04-25 | 2018-11-16 | 中国科学院电子学研究所 | 基于mimo矩形平面阵列的方位向成像方法 |
CN109061622A (zh) * | 2018-06-22 | 2018-12-21 | 电子科技大学 | 一种基于毫米波雷达的隐蔽目标多径探测方法 |
CN109061622B (zh) * | 2018-06-22 | 2022-11-08 | 电子科技大学 | 一种基于毫米波雷达的隐蔽目标多径探测方法 |
CN109298418A (zh) * | 2018-09-30 | 2019-02-01 | 湖南华诺星空电子技术有限公司 | 基于建筑物内部结构特征的雷达检测虚警抑制方法及装置 |
CN109298418B (zh) * | 2018-09-30 | 2022-08-05 | 湖南华诺星空电子技术有限公司 | 基于建筑物内部结构特征的雷达检测虚警抑制方法及装置 |
CN109597065B (zh) * | 2018-12-11 | 2022-09-09 | 湖南华诺星空电子技术有限公司 | 一种用于穿墙雷达检测的虚警抑制方法、装置 |
CN109597065A (zh) * | 2018-12-11 | 2019-04-09 | 湖南华诺星空电子技术有限公司 | 一种用于穿墙雷达检测的虚警抑制方法、装置 |
CN109581317B (zh) * | 2018-12-24 | 2022-05-31 | 电子科技大学 | 一种基于回波峰值匹配的角落目标定位方法 |
CN109581317A (zh) * | 2018-12-24 | 2019-04-05 | 电子科技大学 | 一种基于回波峰值匹配的角落目标定位方法 |
CN109696671A (zh) * | 2018-12-29 | 2019-04-30 | 南京汇君半导体科技有限公司 | 基于组结构与相关性学习的高分辨穿墙雷达成像方法 |
CN111045006A (zh) * | 2019-12-23 | 2020-04-21 | 电子科技大学 | 一种基于多成像字典融合的拐角隐蔽目标成像方法 |
CN111045006B (zh) * | 2019-12-23 | 2022-10-14 | 电子科技大学 | 一种基于多成像字典融合的拐角隐蔽目标成像方法 |
CN112698326A (zh) * | 2020-11-30 | 2021-04-23 | 浙江华消科技有限公司 | 雷达穿墙探测对象的方法、雷达探测设备及电子装置 |
CN114066792A (zh) * | 2021-11-16 | 2022-02-18 | 沈阳航空航天大学 | 一种基于多分辨率融合卷积神经网络的穿墙雷达成像方法 |
CN114066792B (zh) * | 2021-11-16 | 2024-10-25 | 沈阳航空航天大学 | 一种基于多分辨率融合卷积神经网络的穿墙雷达成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107479043B (zh) | 2020-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107479043A (zh) | 基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法 | |
CN105954745B (zh) | 一种适用于穿墙雷达多径幻象抑制的成像方法 | |
CN106324595B (zh) | 一种基于多角度探测成像的多径假目标抑制方法 | |
CN105549011B (zh) | 一种基于mimo穿墙雷达的建筑物单边两点斜视成像方法 | |
CN103197289B (zh) | 一种适用于穿墙雷达成像的多径幻象抑制方法 | |
CN111123257A (zh) | 基于图时空网络的雷达动目标多帧联合检测方法 | |
CN113192646B (zh) | 目标检测模型构建方法及不同目标间距离监控方法、装置 | |
CN106093934B (zh) | 基于改进动态规划的穿墙雷达成像后多目标位置估计方法 | |
CN109471097B (zh) | 一种穿墙雷达信号优化处理方法及装置 | |
CN107589407A (zh) | 多发多收穿墙雷达成像后前后墙振荡一次多径抑制方法 | |
CN111340882B (zh) | 基于图像的室内定位方法及装置 | |
CN110794471A (zh) | 一种毫米波稀疏阵列远程监视成像方法及系统 | |
CN108872977A (zh) | 基于单通道超宽带雷达的生命体双站协同探测方法 | |
CN111045006A (zh) | 一种基于多成像字典融合的拐角隐蔽目标成像方法 | |
CN112560619A (zh) | 一种基于多聚焦图像融合的多距离鸟类精准识别方法 | |
Jia et al. | Multipath ghost and side/grating lobe suppression based on stacked generative adversarial nets in MIMO through-wall radar imaging | |
CN108919259A (zh) | 一种基于多径利用的穿墙雷达建筑墙体位置获取方法 | |
JP7142400B1 (ja) | 多入力多出力レーダーに基づく見通し外ターゲット測位の方法及びシステム | |
CN113050083A (zh) | 一种基于点云的超宽带雷达人体姿态重构方法 | |
CN110927681A (zh) | 一种信号处理方法、装置、终端及存储介质 | |
He et al. | 3d radio imaging under low-rank constraint | |
CN107479052A (zh) | 基于广义高斯分布模型的地面隐蔽目标检测方法 | |
CN106291480A (zh) | 一种基于节点耦合的多径抑制算法 | |
CN116453160A (zh) | 基于毫米波成像的可疑物品无停留式检测方法及系统 | |
Yan et al. | Multipath ghosts location and sub-aperture based suppression algorithm for TWIR |
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 |