CN106680809B - 一种穿墙雷达自聚焦压缩感知成像方法 - Google Patents
一种穿墙雷达自聚焦压缩感知成像方法 Download PDFInfo
- Publication number
- CN106680809B CN106680809B CN201611222421.4A CN201611222421A CN106680809B CN 106680809 B CN106680809 B CN 106680809B CN 201611222421 A CN201611222421 A CN 201611222421A CN 106680809 B CN106680809 B CN 106680809B
- Authority
- CN
- China
- Prior art keywords
- wall
- scene
- iteration
- focusing
- matrix
- 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
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
- 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
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
技术领域
本发明涉及穿墙成像雷达技术领域,特别是一种穿墙雷达自聚焦压缩感知成像方法。
背景技术
穿墙成像雷达具有穿透常见的木材、砖石、混凝土等结构的墙体障碍物对其后隐藏的人员目标进行探测的能力,因此广泛应用于救援、监控、侦察等领域。传统雷达成像方法为了获得高分辨率与良好聚焦性的墙后场景图像,需要大带宽信号与密集采样等条件,使得穿墙探测的成本较高。由于室内同时出现大量目标的概率较低,墙后目标满足时间与空间分布上的稀疏性,故可以采用压缩感知算法,以少量的信号采样实现墙后场景的高质量重建,最终提高穿墙探测的效率。不同于自由空间探测,对墙后场景进行稀疏重构时,必须考虑墙体对雷达信号的迟滞效应才能建立准确的感知矩阵,即必须获知墙体的厚度与介电常数的信息。虽然已经存在一些关于墙体参数估计的研究,但目前其与墙后场景稀疏重构的研究是分离的,这意味着需要引入额外的操作进行墙体参数估计才能建立准确的感知矩阵,降低了墙后场景稀疏重构的效率。
穿墙成像探测时,未知的墙体参数会使得雷达成像结果中目标发生散焦与位置错误现象,传统的墙体参数估计算法需要反复调整天线位置,增加了操作复杂度,另一方面在对墙后场景进行成像时,为获得高分辨的目标像,传统穿墙成像雷达需要发射大带宽信号,因此会占用较多的频带资源,此外,接收端的密集采样也增加了系统成本。
发明内容
本发明所要解决的技术问题是,针对现有技术不足,提供一种穿墙雷达自聚焦压缩感知成像方法,将墙体参数估计与墙后场景成像统一到同一个处理框架中,在获得墙体参数估计的同时得到高分辨的目标像,减少系统发射端占用的频带资源与接收端的采样数,降低操作复杂度,提高探测效率。
为解决上述技术问题,本发明所采用的技术方案是:一种穿墙雷达自聚焦压缩感知成像方法,包括以下步骤:
1)记i为迭代次数索引,定义墙体厚度d(i)与介电常数为第i轮迭代中的超参数,迭代初始时刻,即i=0时,设定上述超参数的初始值为:d(0)=0,
2)若i=0,则转入步骤3);否则,利用上一次迭代获得的墙后场景估计值σ(i-1)更新本次迭代的超参数d(i)与σ(i-1)为第i-1次迭代中估计出的墙后场景;
3)利用超参数d(i)与更新本轮迭代中的墙后场景的探测基矩阵
4)建立如下的墙后场景稀疏重构模型:
其中,Θ为M′N′×MN的随机矩阵,代表Kronecker积;IM′为M′×M′的单位矩阵,为从N×N的单位矩阵中随机抽取N′行而形成的N′×N的测量矩阵,为从M×M的单位矩阵中随机抽取M′行而形成的M′×M的测量矩阵;s′为随机采样信号;λ为惩罚因子,σ(i)为本次迭代中估计出的墙后场景;M为频域采样点数目,N为方位向采样点数目;σ为向量化的探测平面网格的反射系数;
5)若i=0,则令i的值加1,并转入步骤2);否则记Δ为最近两次迭代获得的墙后场景的相对变化量:
6)如果Δ低于预设门限,则认为结果已经收敛,从而结束迭代操作,并输出最后的σ(i+1)作为墙后场景的重构结果;否则令i=i+1,并返回执行步骤2)。
所述预设门限设置为0.01。
利用超参数d(i)与更新本轮迭代中的墙后场景的探测基矩阵的表达式,其中为雷达在第n个方位向采样点上的探测基子矩阵,且有:
其中,θn,p=arctan[yp/(xn-xp)]为信号穿过墙体时折射角的近似值;rp=(xp,yp,zp)为第p个离散场景单元的空间位置矢量;rn=(xn,0,H)为天线在第n个方位向采样点时的空间位置矢量;km=2πfm/c为频率fm对应的波数。
与现有技术相比,本发明所具有的有益效果为:本发明在无需额外操作的情况下,仅利用少量采样信号即可完成墙体参数的估计与墙后目标的高分辨成像,减少了系统发射端占用的频带资源与接收端的采样数,降低了操作复杂度,提高了探测效率。
附图说明
图1为本发明原理框架图;
图2为本发明合成孔径模式下的穿墙探测信号传播模型;
图3为验证实验场景布局示意图,墙后两个三面角顶点坐标为(0.25,3.08)与(1.75,3.08);
图4(a)为不考虑墙体迟滞效应直接成像的结果,图4(b)为利用人工测量的墙体参数d=0.2m与εr=3进行补偿成像的结果。对比可以发现,考虑墙体补偿后,传统成像算法也能够修正墙后目标位置偏差,但所获得的墙后场景图像仍然存在目标像扩散以及栅旁瓣引起的虚假目标问题;
图5(a)为i=1时成像结果,获得的超参数估计值为d=0.05m,εr=1,目标位置估计为(0.29,3.29)与(1.71,3.29);
图5(b)为i=2时的成像结果,获得的超参数估计值为d=0.124m,εr=1.84,目标位置估计为(0.29,3.22)与(1.71,3.22);
图5(c)为i=3时的成像结果,获得的超参数估计值为d=0.204m,εr=2.75,目标位置估计为(0.29,3.08)与(1.71,3.08);
图5(d)为i=4时的成像结果,获得的超参数估计值为d=0.196m,εr=2.97,目标位置估计为(0.29,3.08)与(1.71,3.08)。
具体实施方式
本发明提出一种适用于穿墙雷达的自聚焦压缩感知成像技术,对墙体参数估计与墙后场景稀疏重构进行迭代优化,最终在获取墙体厚度与介电常数信息的同时得到高质量的墙后场景图像。本发明所用的信号模型与所提算法可阐述如下。
记fm为均匀采样模式下第m个频率采样点,则有:
fm=f0+MΔf,m=0,1,…,M-1 (1)
其中f0为初始发射频率,M为频域采样点数目,Δf为频域采样间隔。又将雷达测平面离散化为Nx×Ny的网格,其中Nx为方位向划分数,Ny为距离向划分数。记σ为向量化的探测平面网格的反射系数,则有
雷达成像的目的在于利用接收回波对σ进行估计。
假设雷达系统以合成孔径模式工作,如附图2所示,系统在第n个方位向采样位置时的第m个频域回波采样点s(m,n)可写作:
其中,rp=(xp,yp,zp)为第p个离散场景单元的空间位置矢量,σp为其反射系数,rn=(xn,0,H)为天线在第n个方位向采样点时的空间位置矢量,d为墙体厚度,εr为墙体的介电常数,为噪声,km=2πfm/c为频率fm对应的波数,c为光速。G(rn,rp,km,d,εr)为该探测场景的格林函数:
其中θn,p=arctan[yp/(xn-xp)]为信号穿过墙体时折射角的近似值,。当雷达系统在第n个方位向采样点位置完成全部的频率扫描时,获得的回波信号sn可写作:
其中
sn=[s(0,n) s(1,n) … s(M-1,n)]T (6)
综上,雷达系统完成全部方位向扫描后获得的回波s为:
其中
N为方位向采样点数目。
当雷达运动路线与成像场景确定时,Ψ中的各个空间位置矢量将成为已知常数,此时Ψ可视为仅与墙体的厚度d与墙体的介电常数εr相关的探测基矩阵,因此可以进一步将其(11)记作
常规雷达成像算法在对σ进行反演时,需要预先获知d与εr,否则获得的雷达图像中墙后目标位置将发生偏移,并且目标也会出现一定程度的散焦,如附图4(a)所示。为降低穿墙成像雷达成本,提高探测效率,增加穿墙成像雷达对环境的适应能力,本发明将墙体参数估计与墙后目标高分辨成像综合到一个统一的处理框架下,建立联合二者的迭代优化模型如下:
记i为迭代次数索引,定义墙体厚度d(i)与介电常数为第i轮迭代中的超参数,迭代初始时刻,即i=0时,设定上述超参数的初始值为:d(0)=0,
第一步,若i=0则转入第二步;否则,利用上一轮获得的墙后场景估计σ(i-1)更新超参数,可建立如下的优化模型:
该问题是一个无约束优化问题,可以采用经典的共轭梯度等算法进行求解。
第二步,根据式(4),利用超参数d(i)与更新本轮迭代中的墙后场景的探测基矩阵
第三步,由压缩感知理论可知,当信号具有稀疏性时,可以仅通过少量的采样即可实现高概率重构。大部分穿墙探测应用中,墙后目标满足时间与空间分布上的稀疏特性,因此可以建立如下的墙后场景稀疏重构模型:
其中λ为惩罚因子,σ(i)为本轮处理中估计出的墙后场景。Θ为M′N′×MN的随机矩阵:
代表Kronecker积,IM′为M′×M′的单位矩阵,为从N×N的单位矩阵中随机抽取N′行而形成的N′×N的测量矩阵,为从M×M的单位矩阵中随机抽取M′行而形成的M′×M的测量矩阵。而s′为随机采样信号
实际探测时,利用Θ定义的随机矩阵进行跳频发射获得随机采样信号s′,由于s′为M′×N′的矩阵,相对于均匀采样模式下M×N的接收信号,其采样规模已经降低,因而可以减少探测时占用的频带资源与所需的采样点数。根据压缩感知理论,M′与N′可根据雷达系统能够处理的最大目标数推算出。
式(13)中的问题是一个典型的稀疏重构问题,目前已经有大量公开的算法例如匹配追踪、基追踪、稀疏贝叶斯学习等可以用来求解该问题。
第四步,若i=0,则令并i=i+1转入第一步;否则记Δ为最近两次迭代获得的墙后场景的相对变化量
如果Δ低于某个预设门限,则可认为结果已经收敛,从而结束迭代操作,并输出最后的σ(i+1)作为墙后场景的重构结果,否则令i=i+1,并返回执行第一步。
Claims (3)
1.一种穿墙雷达自聚焦压缩感知成像方法,其特征在于,包括以下步骤:
1)记i为迭代次数索引,定义墙体厚度d(i)与介电常数为第i轮迭代中的超参数,迭代初始时刻,即i=0时,设定上述超参数的初始值为:d(0)=0,
2)若i=0,则转入步骤3);否则,利用上一次迭代获得的墙后场景估计值σ(i-1)更新本次迭代的超参数d(i)与σ(i-1)为第i-1次迭代中估计出的墙后场景估计值;
3)利用超参数d(i)与更新本轮迭代中的墙后场景的探测基矩阵
4)建立如下的墙后场景稀疏重构模型:
其中,Θ为M′N′×MN的随机矩阵, 代表Kronecker积;IM′为M′×M′的单位矩阵,为从N×N的单位矩阵中随机抽取N′行而形成的N′×N的测量矩阵,为从M×M的单位矩阵中随机抽取M′行而形成的M′×M的测量矩阵;s′为随机采样信号;λ为惩罚因子,σ(i)为本次迭代中估计出的墙后场景;M为频域采样点数目,N为方位向采样点数目;σ为向量化的探测平面网格的反射系数;
5)若i=0,则令i的值加1,并转入步骤2);否则记Δ为最近两次迭代获得的墙后场景的相对变化量:
6)如果Δ低于预设门限,则认为结果已经收敛,从而结束迭代操作,并输出最后的σ(i+1)作为墙后场景的重构结果;否则令i=i+1,并返回执行步骤2)。
2.根据权利要求1所述的穿墙雷达自聚焦压缩感知成像方法,其特征在于,所述预设门限设置为0.01。
3.根据权利要求1所述的穿墙雷达自聚焦压缩感知成像方法,其特征在于,利用超参数d(i)与更新本轮迭代中的墙后场景的探测基矩阵的表达式,
其中为雷达在第n个方位向采样点上的探测基子矩阵,且有:
其中,
θn,p=arctan[yp/(xn-xp)]为信号穿过墙体时折射角的近似值;
rp=(xp,yp,zp)为第p个离散场景单元的空间位置矢量;rn=(xn,0,H)为天线在第n个方位向采样点时的空间位置矢量;km=2πfm/c为频率fm对应的波数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611222421.4A CN106680809B (zh) | 2016-12-27 | 2016-12-27 | 一种穿墙雷达自聚焦压缩感知成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611222421.4A CN106680809B (zh) | 2016-12-27 | 2016-12-27 | 一种穿墙雷达自聚焦压缩感知成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106680809A CN106680809A (zh) | 2017-05-17 |
CN106680809B true CN106680809B (zh) | 2019-02-12 |
Family
ID=58870695
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611222421.4A Active CN106680809B (zh) | 2016-12-27 | 2016-12-27 | 一种穿墙雷达自聚焦压缩感知成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106680809B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107894591A (zh) * | 2017-09-30 | 2018-04-10 | 沈阳航空航天大学 | 基于压缩感知的穿墙雷达衍射层析成像方法 |
CN109581373A (zh) * | 2018-12-17 | 2019-04-05 | 中国人民解放军国防科技大学 | 一种全息穿透成像雷达的自聚焦成像方法 |
CN109683160B (zh) * | 2019-03-19 | 2019-06-18 | 长沙深之瞳信息科技有限公司 | 一种适用于穿墙雷达的目标检测方法 |
CN111766575B (zh) * | 2020-06-08 | 2023-04-21 | 桂林电子科技大学 | 一种穿墙雷达自聚焦稀疏成像方法及计算机设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102520405A (zh) * | 2011-12-16 | 2012-06-27 | 中国人民解放军国防科学技术大学 | 基于图像域墙体影响补偿的穿墙成像方法 |
CN103995256A (zh) * | 2014-05-29 | 2014-08-20 | 电子科技大学 | 一种适用于穿墙雷达成像的多面墙体快速补偿方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120235849A1 (en) * | 2010-11-19 | 2012-09-20 | Tatoian James Z | Method and apparatus for through the wall radar imaging |
-
2016
- 2016-12-27 CN CN201611222421.4A patent/CN106680809B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102520405A (zh) * | 2011-12-16 | 2012-06-27 | 中国人民解放军国防科学技术大学 | 基于图像域墙体影响补偿的穿墙成像方法 |
CN103995256A (zh) * | 2014-05-29 | 2014-08-20 | 电子科技大学 | 一种适用于穿墙雷达成像的多面墙体快速补偿方法 |
Non-Patent Citations (3)
Title |
---|
High-Resolution 3-D Imaging Algorithm without Derivative Operations for UWB Through-the-Wall Radars;Suouhei KIDERA et al.;《2008 IEEE Antennas and Propagation Society International Symposium》;20080909;全文 |
基于改进OMP的超宽带穿墙雷达稀疏成像方法;晋良念 等;《计算机技术与应用》;20151231;第41卷(第11期);全文 |
超宽带MIMO穿墙雷达信道建模与运动目标成像;孙鑫;《电子与信息学报》;20140831;第36卷(第8期);全文 |
Also Published As
Publication number | Publication date |
---|---|
CN106680809A (zh) | 2017-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106680809B (zh) | 一种穿墙雷达自聚焦压缩感知成像方法 | |
US9261592B2 (en) | Method and system for through-the-wall imaging using compressive sensing and MIMO antenna arrays | |
CN106646344B (zh) | 一种利用互质阵的波达方向估计方法 | |
US10042046B2 (en) | System and method for radar imaging using distributed arrays and compressive sensing | |
CN105911544B (zh) | 一种基于压缩感知技术的扩频穿墙雷达成像方法 | |
CN106772365A (zh) | 一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法 | |
CN107576961B (zh) | 一种互质降采样间歇合成孔径雷达稀疏成像方法 | |
CN106405548A (zh) | 基于多任务贝叶斯压缩感知的逆合成孔径雷达成像方法 | |
CN110780298B (zh) | 基于变分贝叶斯学习的多基isar融合成像方法 | |
Yang et al. | Multiple-measurement vector model and its application to through-the-wall radar imaging | |
JP5950835B2 (ja) | シーン内の反射物を再構成するシステム | |
CN104199029B (zh) | 一种提高压缩感知雷达目标成像性能的测量矩阵设计方法 | |
CN108845316A (zh) | 一种基于压缩感知理论的雷达稀疏探测方法 | |
CN106918810B (zh) | 一种存在阵元幅相误差时的微波关联成像方法 | |
CN108983168A (zh) | 基于孔径补全的压缩感知mimo雷达目标探测方法 | |
CN102879783A (zh) | 基于稀疏探测频率信号的isar成像方法 | |
CN107340517B (zh) | 一种多传感器多帧检测前跟踪方法 | |
CN110728706A (zh) | 基于深度学习的sar图像精细配准方法 | |
CN107656271B (zh) | 基于压缩感知重构的太赫兹雷达成像算法 | |
CN109507636A (zh) | 基于虚拟域信号重构的波达方向估计方法 | |
CN106707284A (zh) | 一种用于双基地逆合成孔径雷达的成像方法 | |
CN106990385B (zh) | 基于最小冗余线阵的非均匀面阵设计方法 | |
JP2011528784A (ja) | 広帯域rf探知 | |
Li et al. | Sparse MIMO architectures for through-the-wall imaging | |
Liu et al. | Sparsity-driven distributed array imaging |
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 |