CN108828590A - 一种低复杂度熵值扩展穿墙雷达成像方法 - Google Patents
一种低复杂度熵值扩展穿墙雷达成像方法 Download PDFInfo
- Publication number
- CN108828590A CN108828590A CN201810718283.1A CN201810718283A CN108828590A CN 108828590 A CN108828590 A CN 108828590A CN 201810718283 A CN201810718283 A CN 201810718283A CN 108828590 A CN108828590 A CN 108828590A
- Authority
- CN
- China
- Prior art keywords
- entropy
- matrix
- echo
- wall
- signal
- 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.)
- Pending
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/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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种穿墙雷达成像方法,属于穿墙雷达成像算法领域,该方法对天线阵元所接收的回波信号矩阵先进行降维处理,得到低维度的子矩阵,然后根据离散无记忆信源熵值原理,分别对子矩阵进行扩展,利用熵值对每个矩阵设置不同的门限,消除墙体杂波的影响,实现目标成像。该方法能够大幅降低计算复杂度,减少算法中数据占据的内存。
Description
技术领域
本发明涉及一种低复杂度熵值扩展穿墙雷达成像方法,属于穿墙雷达成像算法领域。
背景技术
探测建筑物内或障碍物后的目标问题广泛存在于城市作战、反恐维稳和灾后救援等现代社会的诸多方面。当今的社会整体形式稳定,但是局部地区依然频发冲突乃至爆发局部战争。从伊拉克战争之后,世界各国已经意识到城市作战将在未来的战争中占据越来越重要的地位。对于作战的双方,对周边复杂环境(建筑物,建筑物,掩体)的探测能力技高一筹,便能在战争中获得先机。所以说无论是军用方面还是民用方面,都迫切地需要穿墙成像技术。
穿墙成像过程中,接收天线阵元接收到的回波信号必然会包含大量强杂波信号(如前墙体散射与后墙体散射),这些杂波信号相对于目标信号幅度较大,以至于淹没所需的目标信号,从而严重影响成像效果。因此,为了实现穿墙雷达的目标最优成像,我们必须在成像之前先对回波进行一系列处理,尽最大努力消除不需要的其他信号,使得目标能够清晰地展现出来。
对上面所述问题,国内外的科研人员进行了深入的研究,获取不少的成果,并且提出了不少优良成像算法。其中,背景对消法能够很好地消除墙体散射等杂波信号,对墙后目标进行成像。但是在实际应用中,这种方法不易实现,因为我们不可能将墙后物体移除后对背景进行测量。基于时间门限的算法,是通过设置时间门限将接收回波信号的前部滤除,以此消除墙体散射信号,达到穿墙成像的目的。然而时间门限选择不当,往往会导致部分目标信号也被消除。基于熵值扩展的算法,能够有效抑制墙体强杂波,提高成像精度,但该方法对硬件要求很高、成本较大,穿墙雷达设计将较复杂,尤其该算法计算复杂度大,占用系统内存较多,难于满足实时性需求。因此,如何在不减少天线阵元的数量情况下,又能降低计算复杂度,而通过算法处理消除杂波信号的影响,以保证和提高成像精度,将是我们必须面对的关键问题。
发明内容
本发明为了解决现有技术中存在的问题,提供一种能够大幅降低计算复杂度,减少算法中数据占据内存的熵值扩展穿墙雷达成像方法。
为了达到上述目的,本发明提出的技术方案为:一种低复杂度熵值扩展穿墙雷达成像方法,其特征在于,包括如下步骤:
步骤一、利用均匀分布于测线的自发自收天线阵元,依次发射并接收信号;
步骤二、对天线阵元所接收的回波信号矩阵进行分块处理,获得一系列的子矩阵;
步骤三、然后利用熵值扩展的算法对每个子矩阵进行扩展,并计算扩展后的熵值;
步骤四、对熵值设定相应的门限,消除墙体散射等杂波信号,实现穿墙目标成像。
对上述技术方案的进一步设计为:所述天线阵元数目为N,天线阵元沿平行于墙体侧边的侧线均匀分布。
对回波信号矩阵进行分块处理时,按照矩阵的列进行分块。
矩阵分块的具体方法为:设回波信号矩阵为e,
(aij,i=1,2,…,M,j=1,2,…,N),为各接收天线在阵元j处的采样点i时刻接收到回波信号的强度,N为天线阵元回波信号的数目,M为采样次数;
对回波矩阵e按照列进行分块,得到的子矩阵b1和b2,其大小分别为为M*N1,M*N2。其中各矩阵列数满足:
N1+N2=N
那么回波数据矩阵就转化为b1,b2,如下表示:
由复杂度公式可知,对子矩阵进行扩展时,可以大幅减少计算复杂度,并且能减少数据所占内存,同时算法的处理时间获得提高。另外,根据离散无记忆扩展信源熵值的理论可知,对子矩阵的扩展次数即为熵值扩大倍数,因此能够增大门限可调范围,提高成像的精度,所需天线阵元数目也会大大减少。与其他成像方法相比,本方法在保障计算量低的同时,达到成像精度高,阵元数目少的特点。因此在完成穿墙成像任务中,对硬件设备的要求较低,能够大大简化穿墙雷达硬件设备的设计。
附图说明
图1是本实施例穿墙成像模型结构示意图;
图2是在采样点数与扩展次数一定,算法复杂度与阵元之间的关系图;
图3是在采样点数与阵元数据一定,算法复杂度与扩展次数之间的关系图;
图4是阵元数目与扩展次数同时改变,计算复杂度的变化图;
图5是原始数据成像图;
图6是采用本实施例方法进行杂波滤除后的成像图。
具体实施方式
下面结合附图以及具体实施例对本发明进行详细说明。
实施例
下面结合附图对本发明的实现步骤做出详细的说明:
本实施例的低复杂度熵值扩展穿墙雷达成像方法,具体步骤为:
步骤1,建立穿墙成像模型
参照图1,设置收发天线阵元数目N,沿平行于墙体的侧边均匀分布,其与墙体之间距离h。发射信号是Ricker子波,中心频率为f。前墙体和后墙体的材料为混凝土,是均匀无耗介质,厚度d,相对介电常数εr。目标是半径为a的圆球,为理想电导体,球心距离墙壁δ。
步骤2,获得原始信号数据
原始回波数据由仿真软件GprMax来获得。对N组回波信号分别进行采样,记采样次数为M。回波信号数据可组成M*N维的矩阵:
(aij,i=1,2,…,M,j=1,2,…,N),为各接收天线在阵元j处的采样点i时刻接收到回波信号的强度。
步骤3,回波矩阵按列分块并分别进行扩展
回波矩阵e按照列进行分块,即对回波矩阵进行分解,得到低维度的子矩阵b1,b2。其大小分别为为M*N1,M*N2。其中各矩阵列数满足:
N1+N2=N
那么回波数据矩阵就转化为b1,b2,如下表示:
利用熵值扩展的算法对矩阵b1进行处理得到矩阵b1z,利用熵值扩展的算法对矩阵b2进行处理得到矩阵b2z。
将b1z与b2z进行合并得到
bz=(b1z,b2z)
由图5与图6可以看出,本算法可以达到滤除杂波的效果,结合表1与图2、图3、图4来看,算法复杂度低,避免了不必要的计算步骤。同时,从占据内存上来看也是非常明显的。理论上,扩展次数越大,门限可调范围越大,成像精度也越高,同时也能够大大减少所需阵元数目,然而在一定程度上会增大计算量,因此,实际工程中根据具体情况而定,但是对于本算法可以更进一步的进行扩展。因为将总阵元数目分为两个比较小的阵元数目,复杂度度大幅下降,另外图像得到进一步优化,杂波进一步被滤除。
表1为算法的计算复杂度以及在程序运行中关键数据所占据的内存大小于扩展次数的关系。
时间算法复杂度的计算公式为o(M*NL),其中M为采样数,N为阵元数,L为扩展次数,由图2看出,在扩展次数一定情况下,随着L的增加,计算复杂度会急剧上升,并且随着算法的进行,数据所需的内存也是非常巨大的。由图3看出,在阵元数目一定情况下,随着N的增加,计算复杂度上升巨大,并且显示扩展次数比阵元数目更具有影响。图4更能体现出三者之间的关系。
经过本发明算法的改进后,计算复杂度公式变为o(M*(N1L+N2L)),由于N1+N2=N,那么随着扩展次数的增加,那么
N1L+N2L<<NL
因此本发明的方法可大幅降低计算复杂度。
本发明的不局限于上述各实施例,凡采用等同替换方式得到的技术方案均落在本发明要求保护的范围内。
Claims (4)
1.一种低复杂度熵值扩展穿墙雷达成像方法,其特征在于,包括如下步骤:
步骤一、利用均匀分布于测线的自发自收天线阵元,依次发射并接收信号;
步骤二、对天线阵元所接收的回波信号矩阵进行分块处理,获得一系列的子矩阵;
步骤三、然后利用熵值扩展的算法对每个子矩阵进行扩展,并计算扩展后的熵值;
步骤四、对熵值设定相应的门限,消除墙体散射等杂波信号,实现穿墙目标成像。
2.根据权利要求1所述低复杂度熵值扩展穿墙雷达成像方法,其特征在于,所述天线阵元数目为N,天线阵元沿平行于墙体侧边的侧线均匀分布。
3.根据权利要求2所述低复杂度熵值扩展穿墙雷达成像方法,其特征在于,对回波信号矩阵进行分块处理时,按照矩阵的列进行分块。
4.根据权利要求3所述低复杂度熵值扩展穿墙雷达成像方法,其特征在于,设回波信号矩阵为e,
(aij,i=1,2,…,M,j=1,2,…,N),为各接收天线在阵元j处的采样点i时刻接收到回波信号的强度,N为天线阵元回波信号的数目,M为采样次数;
对回波矩阵e按照列进行分块,得到的子矩阵b1和b2,其大小分别为为M*N1,M*N2。其中各矩阵列数满足:
N1+N2=N
那么回波数据矩阵就转化为b1,b2,如下表示:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810718283.1A CN108828590A (zh) | 2018-07-03 | 2018-07-03 | 一种低复杂度熵值扩展穿墙雷达成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810718283.1A CN108828590A (zh) | 2018-07-03 | 2018-07-03 | 一种低复杂度熵值扩展穿墙雷达成像方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108828590A true CN108828590A (zh) | 2018-11-16 |
Family
ID=64133655
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810718283.1A Pending CN108828590A (zh) | 2018-07-03 | 2018-07-03 | 一种低复杂度熵值扩展穿墙雷达成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108828590A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110187340A (zh) * | 2019-06-17 | 2019-08-30 | 中国电子科技集团公司信息科学研究院 | 一种基于熵的探测目标的信息表征方法及系统 |
CN111580099A (zh) * | 2020-06-12 | 2020-08-25 | 南京信息工程大学 | 基于联合熵值的穿墙成像雷达墙体杂波抑制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101915943A (zh) * | 2010-08-10 | 2010-12-15 | 中南大学 | 均匀背景介质的介电常数和隐蔽目标参数的联合反演方法 |
CN104297749A (zh) * | 2014-10-31 | 2015-01-21 | 南京信息工程大学 | 一种基于回波扩展熵的穿墙雷达成像方法 |
CN104965203A (zh) * | 2015-05-25 | 2015-10-07 | 北京理工雷科电子信息技术有限公司 | 一种目标透视雷达成像方法 |
CN106407363A (zh) * | 2016-09-08 | 2017-02-15 | 电子科技大学 | 一种基于信息熵的超高维数据降维算法 |
-
2018
- 2018-07-03 CN CN201810718283.1A patent/CN108828590A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101915943A (zh) * | 2010-08-10 | 2010-12-15 | 中南大学 | 均匀背景介质的介电常数和隐蔽目标参数的联合反演方法 |
CN104297749A (zh) * | 2014-10-31 | 2015-01-21 | 南京信息工程大学 | 一种基于回波扩展熵的穿墙雷达成像方法 |
CN104965203A (zh) * | 2015-05-25 | 2015-10-07 | 北京理工雷科电子信息技术有限公司 | 一种目标透视雷达成像方法 |
CN106407363A (zh) * | 2016-09-08 | 2017-02-15 | 电子科技大学 | 一种基于信息熵的超高维数据降维算法 |
Non-Patent Citations (3)
Title |
---|
会飞的牛牛: "矩阵运算中选择分块矩阵策略的研究", 《矩阵运算中选择分块矩阵策略的研究》 * |
李家强等: "基于扩展信源熵值的穿墙成像雷达墙体强杂波抑制", 《物理学报》 * |
管涛: "《数据分析基础及模型》", 31 July 2015, 合肥工业大学出版社 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110187340A (zh) * | 2019-06-17 | 2019-08-30 | 中国电子科技集团公司信息科学研究院 | 一种基于熵的探测目标的信息表征方法及系统 |
CN111580099A (zh) * | 2020-06-12 | 2020-08-25 | 南京信息工程大学 | 基于联合熵值的穿墙成像雷达墙体杂波抑制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102243303B (zh) | 一种基于呼吸特征的静止人体穿墙定位方法 | |
US8077072B2 (en) | Static RF imaging for inside walls of a premises | |
CN109375204B (zh) | 基于雷达的目标检测方法、系统、设备及介质 | |
CN109597065B (zh) | 一种用于穿墙雷达检测的虚警抑制方法、装置 | |
CN109919870B (zh) | 一种基于bm3d的sar图像相干斑抑制方法 | |
CN108120976B (zh) | 一种基于多普勒通道特性的地杂波谱泄露抑制方法 | |
CN103197302B (zh) | 一种适用于穿墙雷达成像的目标位置提取方法 | |
CN108828590A (zh) | 一种低复杂度熵值扩展穿墙雷达成像方法 | |
CN103197289A (zh) | 一种适用于穿墙雷达成像的多径幻象抑制方法 | |
CN105549011A (zh) | 一种基于mimo穿墙雷达的建筑物单边两点斜视成像方法 | |
CN111580099A (zh) | 基于联合熵值的穿墙成像雷达墙体杂波抑制方法 | |
CN107942329A (zh) | 机动平台单通道sar对海面舰船目标检测方法 | |
CN109541585A (zh) | 一种基于峰度评估的人体穿墙检测成像方法 | |
CN108802725A (zh) | 一种浅层穿透雷达合成孔径成像方法 | |
CN105093187A (zh) | 一种基于目标轨迹提取的穿墙雷达的成像预处理方法 | |
CN105513076A (zh) | 基于目标坐标分布特性的弱小目标恒虚警检测方法 | |
Sun et al. | Moving target imaging and ghost mitigation in through-the-wall sensing application | |
CN104732190B (zh) | 一种基于正交纹理相关分析的合成孔径声纳目标检测方法 | |
Luo et al. | Improved SMC-PHD filter for multi-target track-before-detect | |
CN110221289A (zh) | 用于三坐标有源相控阵雷达的目标检测方法 | |
CN110632592B (zh) | 一种手持式穿墙雷达的虚警消除方法 | |
Tiwari et al. | Development of a model for detection and estimation of depth of shallow buried non-metallic landmine at microwave X-band frequency | |
CN112904334B (zh) | 一种基于互相关的探地雷达后向投影快速成像方法 | |
CN103675811B (zh) | 基于阵列雷达的扩展动目标穿墙成像方法 | |
Smitha et al. | Clutter reduction using background subtraction of ground penetrating radar for landmine detection |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181116 |