CN109270526B - 一种针对墙体位置未知的角落目标定位方法 - Google Patents
一种针对墙体位置未知的角落目标定位方法 Download PDFInfo
- Publication number
- CN109270526B CN109270526B CN201811030816.3A CN201811030816A CN109270526B CN 109270526 B CN109270526 B CN 109270526B CN 201811030816 A CN201811030816 A CN 201811030816A CN 109270526 B CN109270526 B CN 109270526B
- Authority
- CN
- China
- Prior art keywords
- wall
- target
- radar
- corner
- array
- 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 22
- 238000003384 imaging method Methods 0.000 claims abstract description 25
- 230000005540 biological transmission Effects 0.000 claims abstract description 13
- 238000001514 detection method Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 2
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 230000001678 irradiating effect Effects 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 abstract description 4
- 238000002592 echocardiography Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 230000007123 defense Effects 0.000 description 2
- 241000588722 Escherichia Species 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
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/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
技术领域
本发明涉及多径利用雷达目标定位技术领域,特别是涉及目标位于雷达非直视区且墙体位置未知时的单目标定位方法。
背景技术
多径利用雷达目标定位技术领域主要利用电磁波在建筑物表面的镜面反射和衍射,对位于雷达非直视区的隐秘目标进行检测、定位和识别,在城市巷战、反恐作战、灾难救援等领域具有重大的应用价值。电磁波在拐角巷道内传播时,存在许多不同的传播路径,其中利用墙角的电磁波衍射可以探测到目标,同时,由于巷道光滑墙体的存在,电磁波从发射雷达发射出去之后会在墙体表面进行一次或者多次反射,以此也可以探测到目标。
国内外许多研究机构开展了非直视情况下利用多径回波对目标探测的研究。荷兰国防研究局利用多径回波产生的距离-多普勒图像区分不同的动目标(Deiana,D.,Kossen,A.S.,van Rossum,W.L.:'Multipath exploitation in an urban environment using aMIMO surveillance radar'.International Radar Symposium,Vilnius,Lithuania,2010,pp.1-4.),瑞典国防研究局利用多径回波成功获取了动目标的微多普勒特征(Gustafsson,M.,Andersson,Asa.,Johansson,T.,et al.:'Extraction of Human Micro-Doppler Signature in an Urban Environment Using a Sensing-Behind-the-CornerRadar',IEEE Geoscience and Remote Sensing Letters,2016,13,(2),pp.187-191).
上述研究只是利用多径探测到了目标,并未实现目标的定位。伊利诺大学芝加哥分校提出了利用点扩散函数和复合点扩散函数对隐秘目标进行定位的方法(Setlur,P.,Negishi,T.,Devroye,N.,et al.:'Multipath Exploitation in Non-LOS UrbanSynthetic Aperture Radar',IEEE Journal of Selected Topics in SignalProcessing,2014,8,(1),pp.137-152).伊尔默瑙工业大学提出了利用不同多径回波对动目标进行定位的方法(Zetik,R.,Eschrich,M.,Jovanoska,S.,et al.:'Looking behind acorner using multipath-exploiting UWB radar',IEEE Transactions on Aerospaceand Electronic Systems,2015,51,(3),pp.1916-1926),它从目标回波中提取出衍射和一次反射以及两者组合回波,然后利用场景几何结构对目标进行定位。但是上述方法都假设建筑物位置已知,在实际应用中,墙体位置一般未知。因此研究墙体位置未知环境下的非直视目标定位问题具有重要的价值。
发明内容
本发明为解决上述问题,提出了一种适用于墙体位置未知环境下的非直视目标定位方法,首先利用后向投影成像算法估计墙体位置,然后对目标进行成像,提取出雷达图像中的聚焦区域,然后根据聚焦区域中心到阵列中心的距离对聚焦区域进行排序,最后根据多径传输模型将聚焦区域与实际传输路径匹配,并按不同的路径结合墙体位置计算目标位置。
本发明技术方案如下:
一种针对墙体位置未知的角落目标定位方法,包括以下步骤:
步骤1:角落目标多径传输环境参数初始化
巷道的两面平行墙体Wall-1和Wall-2的位置记为D1、D2,与Wall-1垂直的墙体为Wall-3,两者组成的墙角位置记为(xc,yc),其中xc、yc分别为墙角的x轴坐标、y轴坐标。位于雷达非直视区的目标位置记为(xp,yp);雷达探测目标的两种路径记为Path-Ⅰ、Path-Ⅱ,Path-Ⅰ表示雷达发射的电磁波经过墙体Wall-2的一次镜面反射探测到目标的传输路径,Path-Ⅱ表示雷达发射的电磁波经过墙体Wall-2的一次镜面反射和墙体Wall-1的一次镜面反射探测到目标的传输路径;
步骤2:墙体位置获取
步骤2.1、将雷达阵列放置于与墙角共线的位置,波束方向对着墙体Wall-3方向;
步骤2.2、运用后向投影成像算法获取墙体Wall-3的图像,由于墙体具有较强的散射系数,图像中将有明显的Wall-3的图像;又由于后向投影成像算法基于距离的成像特性,在Wall-3关于阵列对称的位置,出现另一个虚拟墙体图像;
步骤2.3、利用霍夫变换直线提取方法,将上述图像中Wall-3所在直线提取出来,用符号Line1表示,同时将Wall-3的虚拟墙体图像所在直线提取出来,用符号Line2表示;由于Line1与Line2关于雷达阵列所在直线对称,且雷达阵列与墙角共线,可知两条直线的交点就是墙角所在位置;因此求解两条直线Line1与Line2的交点,即可获得墙角位置(xc,yc),由于墙角由两面相互垂直的墙体构成,所以墙体Wall-1的位置即可获得,为D1=xc;
步骤2.4、在雷达阵列中心不变的情况下旋转阵列,使其波束方向对着另一面平行墙体Wall-2;对该墙体用后向投影算法进行成像,对图像聚焦区域提取,即可获得墙体Wall-2的位置D2;
步骤3:采用雷达照射墙体,根据获得的回波信号进行目标位置获取。
进一步的,所述步骤3的具体方法为:
步骤3.1、对目标探测回波进行成像,获得成像结果,记为Io(·);
步骤3.2、提取出Io(·)的多个聚焦区域,并按照聚焦区域中心到阵列中心距离从小到大的顺序,对聚焦区域进行排序,得到聚焦区域序列,记为[A11,A12,A22],并将三个聚焦区域的中心分别记为(x1,y1)、(x2,y2)、(x3,y3);
步骤3.3、对于聚焦区域A11,根据几何关系,目标位置可以计算为:
步骤3.4、对于聚焦区域A22,根据几何关系,目标位置可以计算为:
步骤3.5、对于聚焦区域A12,根据几何关系,目标位置可由下面方程解出:
本发明的有益效果是:
本发明提出了一种适用于墙体位置未知环境下的非直视目标定位方法,首先通过雷达阵列旋转的方式,同时利用后向投影成像算法估计墙体位置,然后对目标进行成像,提取出雷达图像中的聚焦区域,然后根据聚焦区域中心到阵列中心的距离对聚焦区域进行排序,最后根据多径传输模型将聚焦区域与实际传输路径匹配,并按不同的路径结合墙体位置计算目标位置。本发明的优点是能准确估计非直视区目标的位置,且操作简单,不暴露操作人员位置。本发明可以应用于城市巷战、反恐作战等领域。
附图说明
图1为角落环境下多径传播模型。
图2为具体实施方式中仿真场景示意图。
图3为具体实施方式中对墙体Wall-3的成像结果
图4为具体实施方式中仿真场景墙体位置估计结果。
图5为具体实施方式中仿真场景目标位置估计结果。
具体实施方式
下面根据一个GPRMAX电磁仿真例子给出本发明的具体实施方式:
步骤1:角落目标多径传输环境参数初始化
对存在目标的待探测角落区域,其仿真场景如图2所示,坐标系零点位于场景左下角,墙体Wall-1与墙体Wall-2平行于y轴,且横坐标分别为D1=1.8m、D2=2.5m。墙体Wall-3垂直于Wall-1,组成的墙角坐标为(xc,yc)=(1.8m,1.5m)。单目标位于雷达非直视区,坐标为(xp,yp)=(2.0m,2.1m)。雷达为2发4收的时分MIMO雷达,其中两个发射天线位于阵列两端,与其邻近的接收天线相距3.75cm,四个接收天线位于阵列中间,相互之间间隔7.5cm。
步骤2:墙体位置获取
2.1、将雷达阵列放置于与墙角共线的位置,波束方向对着墙体Wall-3方向,如图2阵列A所示,其阵列中心坐标为(1.2343m,0.9343m),且与Wall-3夹角θ=45°。
2.2、运用后向投影成像算法获取墙体Wall-3的图像,如图3所示,墙体Wall-3在图像中有明显的聚焦区域,同时,在与阵列对称的位置,出现了Wall-3的虚拟图像。
2.3、利用霍夫变换直线提取方法,将上述图像中Wall-3所在直线提取出来,用符号Line1表示,同时将Wall-3的虚拟图像所在直线提取出来,用符号Line2表示。求解两条直线Line1与Line2的交点,即可获得墙角位置(xc,yc)=(1.8m,1.5m),由于墙角由两面相互垂直的墙体构成,所以墙体Wall-1的位置即可获得,为D1=xc=1.83m。
2.4、在雷达阵列中心不变的情况下旋转阵列,使其波束方向对着另一面平行墙体Wall-2。对该墙体用后向投影算法进行成像,对图像聚焦区域提取,即可获得墙体Wall-2的位置,最终仿真结果墙体位置估计如图4所示,黑色虚线即为三面墙体的实际位置。
步骤3:目标位置获取
3.1、对目标探测回波进行成像,获得成像结果如图5所示。
3.2、提取出图5中的多个聚焦区域,并按照聚焦区域中心到阵列中心距离从小到大的顺序,对聚焦区域进行排序,得到聚焦区域序列,记为[A11,A12,A22],且三个聚焦区域的中心分别为(x1,y1)=(2.9368m,2.0751m)、(x2,y2)=(3.2344m,2.0339m)、(x3,y3)=(3.4393m,2.0812m)。
3.3、对于聚焦区域A11,由式(1),目标位置可以计算为(xp1,yp1)=(2.0632m,2.0751m)
3.4、对于聚焦区域A22,由式(2),目标位置可以计算为(xp2,yp2)=(2.0781m,2.072m)
3.5、对于聚焦区域A12,由式(3),目标位置可由计算为(xp3,yp3)=(2.0393m,2.0812m)
3.6、对以上三个结果求平均,可得最终目标位置为(xp,yp)=(2.0602m,2.0761m),如图5中三角形所示,即为目标位置。
具体参数结果如表1所示,第一行为各个参数的真实值,第二行为本发明的估计值,第三行为两者的欧几里德距离。
由上表可以看出,最大的误差仅仅是0.065m。因此本发明提供的适用于墙体位置未知环境下的非直视目标定位方法不仅操作简单不会暴露目标,而且墙体位置及目标位置的估计误差都很小,验证了本发明的正确性和有效性。
Claims (2)
1.一种针对墙体位置未知的角落目标定位方法,该方法包括:
步骤1:角落目标多径传输环境参数初始化
巷道的两面平行墙体Wall-1和Wall-2的位置记为D1、D2,与Wall-1垂直的墙体为Wall-3,两者组成的墙角位置记为(xc,yc),其中xc、yc分别为墙角的x轴坐标、y轴坐标;位于雷达非直视区的目标位置记为(xp,yp);雷达探测目标的两种路径记为Path-Ⅰ、Path-Ⅱ,Path-Ⅰ表示雷达发射的电磁波经过墙体Wall-2的一次镜面反射探测到目标的传输路径,Path-Ⅱ表示雷达发射的电磁波经过墙体Wall-2的一次镜面反射和墙体Wall-1的一次镜面反射探测到目标的传输路径;所述雷达为MIMO雷达;
步骤2、获取墙体位置:
步骤2.1、将雷达阵列放置于与墙角共线的位置,波束方向对着墙体Wall-3方向;
步骤2.2、运用后向投影成像算法获取墙体Wall-3的图像;
步骤2.3、利用霍夫变换直线提取方法,将上述图像中Wall-3所在直线和Wall-3所产生的虚拟图像所在直线提取出来,求解两直线交点,得到墙体Wall-1位置;
步骤2.4、在雷达阵列中心不变的情况下旋转阵列,使其波束方向对着另一面平行墙体Wall-2;对该墙体用后向投影算法进行成像,对图像聚焦区域提取,获得墙体Wall-2的位置;
步骤3:采用雷达照射墙体,根据获得的回波信号进行目标位置获取;
步骤3.1、对目标探测回波进行成像,获得成像结果,记为Io(·);
步骤3.2、提取出Io(·)的多个聚焦区域,并按照聚焦区域中心到阵列中心距离从小到大的顺序,对聚焦区域进行排序,得到聚焦区域序列,记为[A11,A12,A22],并将三个聚焦区域的中心分别记为(x1,y1)、(x2,y2)、(x3,y3);
步骤3.3、对于聚焦区域A11,根据几何关系,目标位置计算为:
步骤3.4、对于聚焦区域A22,根据几何关系,目标位置计算为:
步骤3.5、对于聚焦区域A12,根据几何关系,目标位置由下面方程解出:
步骤3.6、对以上三个结果求平均,可得最终目标位置。
2.根据权利要求1所述的一种针对墙体位置未知的角落目标定位方法,其特征在于所述步骤2的具体方法为:
步骤2.1、将雷达阵列放置于与墙角共线的位置,波束方向对着墙体Wall-3方向;
步骤2.2、运用后向投影成像算法获取墙体Wall-3的图像,由于墙体具有较强的散射系数,图像中将有明显的Wall-3的图像;又由于后向投影成像算法基于距离的成像特性,在Wall-3关于阵列对称的位置,出现另一个虚拟墙体图像;
步骤2.3、利用霍夫变换直线提取方法,将上述图像中Wall-3所在直线提取出来,用符号Line1表示,同时将Wall-3的虚拟墙体图像所在直线提取出来,用符号Line2表示;由于Line1与Line2关于雷达阵列所在直线对称,且雷达阵列与墙角共线,可知两条直线的交点就是墙角所在位置;因此求解两条直线Line1与Line2的交点,即可获得墙角位置(xc,yc),由于墙角由两面相互垂直的墙体构成,所以墙体Wall-1的位置即可获得,为D1=xc;
步骤2.4、在雷达阵列中心不变的情况下旋转阵列,使其波束方向对着另一面平行墙体Wall-2;对该墙体用后向投影算法进行成像,对图像聚焦区域提取,即可获得墙体Wall-2的位置D2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811030816.3A CN109270526B (zh) | 2018-09-05 | 2018-09-05 | 一种针对墙体位置未知的角落目标定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811030816.3A CN109270526B (zh) | 2018-09-05 | 2018-09-05 | 一种针对墙体位置未知的角落目标定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109270526A CN109270526A (zh) | 2019-01-25 |
CN109270526B true CN109270526B (zh) | 2023-06-16 |
Family
ID=65188420
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811030816.3A Active CN109270526B (zh) | 2018-09-05 | 2018-09-05 | 一种针对墙体位置未知的角落目标定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109270526B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110865376A (zh) * | 2019-12-02 | 2020-03-06 | 电子科技大学 | 一种基于tbd的非直视拐角人体目标定位方法 |
CN111045006B (zh) * | 2019-12-23 | 2022-10-14 | 电子科技大学 | 一种基于多成像字典融合的拐角隐蔽目标成像方法 |
CN111551115A (zh) * | 2020-06-11 | 2020-08-18 | 中国二十二冶集团有限公司 | 测量视线受阻构件定位坐标的方法 |
EP3982160A1 (en) | 2020-10-08 | 2022-04-13 | Imec VZW | Method and system for indoor multipath ghosts recognition |
CN113960558B (zh) * | 2021-11-24 | 2023-10-24 | 成都理工大学 | 基于多输入多输出雷达的非视距目标定位方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105549011A (zh) * | 2015-12-16 | 2016-05-04 | 成都理工大学 | 一种基于mimo穿墙雷达的建筑物单边两点斜视成像方法 |
CN107589407A (zh) * | 2017-08-29 | 2018-01-16 | 电子科技大学 | 多发多收穿墙雷达成像后前后墙振荡一次多径抑制方法 |
-
2018
- 2018-09-05 CN CN201811030816.3A patent/CN109270526B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105549011A (zh) * | 2015-12-16 | 2016-05-04 | 成都理工大学 | 一种基于mimo穿墙雷达的建筑物单边两点斜视成像方法 |
CN107589407A (zh) * | 2017-08-29 | 2018-01-16 | 电子科技大学 | 多发多收穿墙雷达成像后前后墙振荡一次多径抑制方法 |
Non-Patent Citations (4)
Title |
---|
"Looking Behind a Corner Using Multipath-Exploiting UWB Radar";RUDOLF ZETIK et al.;《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》;20150731;第51卷(第3期);第1916-1926页 * |
"Single-Side Two-Location Spotlight Imaging for Building Based on MIMO Through-Wall-Radar";Yong Jia et al.;《sensors》;20160907;第1-17页 * |
"一种新的MIMO线阵穿墙成像模型及其环境参数估计";宋勇平 等;《电子与信息学报》;20141231;第36卷(第12期);第2980-2985页 * |
"基于改进Camshift的穿墙雷达运动人体目标成像跟踪算法";李松林 等;《计算机应用》;20180210;第38卷(第2期);第528-532、609页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109270526A (zh) | 2019-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109270526B (zh) | 一种针对墙体位置未知的角落目标定位方法 | |
Leigsnering et al. | Multipath exploitation and suppression for SAR imaging of building interiors: An overview of recent advances | |
Zhang et al. | Extending reliability of mmwave radar tracking and detection via fusion with camera | |
Setlur et al. | Multipath model and exploitation in through-the-wall and urban radar sensing | |
CN109581317B (zh) | 一种基于回波峰值匹配的角落目标定位方法 | |
Amin et al. | Through-the-wall radar imaging: theory and applications | |
Li et al. | Multiple targets localization behind L-shaped corner via UWB radar | |
CN105549011A (zh) | 一种基于mimo穿墙雷达的建筑物单边两点斜视成像方法 | |
CN107121665B (zh) | 一种基于稀疏阵的近场相干源的无源定位方法 | |
CN111045006B (zh) | 一种基于多成像字典融合的拐角隐蔽目标成像方法 | |
US20220114363A1 (en) | Method and System for Indoor Multipath Ghosts Recognition | |
Tahmoush et al. | Radar surveillance in urban environments | |
CN104166129A (zh) | 一种实波束雷达迭代最小均方误差角超分辨方法 | |
Li et al. | On the electromagnetic diffraction propagation model and applications | |
Guo et al. | MIMO radar localization of targets behind L-shaped corners | |
Smith et al. | Analysis and exploitation of multipath ghosts in radar target image classification | |
Zetik et al. | Cooperative localization and object recognition in autonomous UWB sensor networks | |
Yektakhah et al. | A method for detection of flat walls in through-the-wall SAR imaging | |
Long et al. | Ground moving target signal model and power calculation in forward scattering micro radar | |
CN113960558B (zh) | 基于多输入多输出雷达的非视距目标定位方法及系统 | |
Zhao et al. | Using sky-wave echoes information to extend HFSWR's maximum detection range | |
Gao et al. | Static background removal in vehicular radar: Filtering in azimuth-elevation-doppler domain | |
Skartados et al. | GPR Antenna localization based on A-Scans | |
Tran et al. | Toward passive RF tomography: Signal processing and experimental validation | |
Sakamoto et al. | A target tracking method with a single antenna using time-reversal UWB radar imaging in a multi-path environment |
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 |