CN103278818B - 基于相位统计参数加权的穿墙雷达成像方法 - Google Patents
基于相位统计参数加权的穿墙雷达成像方法 Download PDFInfo
- Publication number
- CN103278818B CN103278818B CN201310233525.5A CN201310233525A CN103278818B CN 103278818 B CN103278818 B CN 103278818B CN 201310233525 A CN201310233525 A CN 201310233525A CN 103278818 B CN103278818 B CN 103278818B
- Authority
- CN
- China
- Prior art keywords
- statistical parameter
- wall radar
- imaging
- echo data
- phase
- 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.)
- Expired - Fee Related
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种基于相位统计参数加权的穿墙雷达成像方法。技术方案是:第一步,对回波数据进行相位补偿,得到补偿后的成像目标区域的回波数据;第二步,求成像目标区域每个像素点对应的相位统计参数;第三步,利用标准后向投影成像算法进行成像处理,得到成像目标区域图像。第四步,将成像目标区域图像的每个像素点与该像素点对应的相位统计参数相乘,得到改善的成像目标区域图像。本发明通过合理利用穿墙雷达的回波数据的相位统计特性,对穿墙雷达图像进行加权处理,提高穿墙雷达图像质量。本发明计算简单,易于实现,可用于实际工程设计与信号处理。
Description
技术领域
本发明属于穿墙雷达成像技术领域,提出了一种基于相位统计参数加权的穿墙雷达成像方法,以提高穿墙雷达成像质量。
背景技术
穿墙雷达通常采用步进频率体制,发射具有良好穿透性的低频电磁波,能够穿透墙壁等障碍物,对隐蔽目标进行探测,在灾后救援、反恐维稳等方面具有广泛的应用前景,而穿墙雷达成像是穿墙雷达的重要发展方向。穿墙雷达成像是通过发射低频超宽带信号在保证穿透性的同时获得高的距离向分辨率,采用阵列天线在实现快速方位扫描的同时可获得高的方位向分辨率。
对墙后目标的高精度成像是提升穿墙雷达探测性能的关键。然而一方面电磁波在建筑物内传播不可避免存在衰减、多径、色散和折射等现象,另一方面由于采用超宽带信号,信号频率低端的波长往往是频率高端波长的几倍,由于结构的限制,阵列天线相邻阵元间的距离一般大于信号频率高端波长的一半,即阵列对高端频率信号是稀疏的。因此穿墙雷达图像质量受到图像散焦和假目标的困扰,从视觉效果上看,图像分辨率下降,对比度降低。
为了提高穿墙雷达成像质量,消除墙体对成像质量的影响,最有效的方法是测量墙体电磁参数,计算成像场景的电磁传播特性,在成像时对衰减、多径、色散和折射现象和阵列结构误差进行补偿。然而,不同墙体电磁参数随墙体材料不同和含水量不同而变化,建筑物内部电磁传播特性随建筑物结构不同而变化,实际应用时又很难对成像目标作现场测量,限制了该方法的实际应用效果。另外,该方法对抑制阵列稀疏性引起的栅瓣杂波也无能为力。
专利CN102520405A提出了一种基于图像域墙体影响补偿的穿墙成像方法,该方法先采用通用的雷达成像方法进行回波成像,再通过图像处理的方法消除墙体影响,提高图像质量。通常穿墙雷达回波非常微弱,信噪比低,先对回波进行处理,再利用处理后的回波进行成像,目前尚未查找到按照这种思路实现穿墙雷达成像的方法。
发明内容
本发明提出基于相位统计参数加权的穿墙雷达成像方法。该方法能提高穿墙雷达图像的分辨率和抑制杂波,同时具有方法简便、易于实现等优势。
本发明的技术方案是:
第一步,对穿墙雷达阵列天线接收到的步进频率波形回波数据进行基于自由空间传播历程的相位补偿,得到补偿后的成像目标区域的回波数据。
第二步,利用补偿后的成像目标区域的回波数据,求成像目标区域每个像素点对应的相位统计参数。该统计参数能表征墙体对穿透成像的影响和阵列稀疏性对步进频率信号雷达成像的影响。
第三步,利用标准后向投影成像算法(Back Projection,BP)对补偿后的成像目标区域的回波数据进行成像处理,得到成像目标区域图像。
第四步,将成像目标区域图像的每个像素点与该像素点对应的相位统计参数相乘,得到改善的成像目标区域图像。
本发明的有益效果:本发明通过合理利用穿墙雷达的回波数据的相位统计特性,设计了能有效表征墙体影响和阵列稀疏性影响的相位统计参数,对穿墙雷达图像进行加权处理,提高穿墙雷达图像质量。该方法既能改善穿墙雷达图像分辨率,又能抑制假目标,改善图像对比度。同时该方法思想简单,算法易于实现,可用于实际工程设计与信号处理。
附图说明
图1为本发明基于相位统计参数加权的穿墙雷达成像流程示意图;
图2为进行实验的场景图;
图3为传统成像算法获得的穿墙成像结果;
图4为利用本发明法获得的穿墙成像结果。
具体实施方式
下面结合附图对本发明进行详细介绍。
设穿墙雷达收发阵元数为M,发射信号步进数为K。第m个阵元接收到的第k个频率位于像素点(x,y)处的理想点目标的回波数据用Ek,m(x,y)表示;第k个步进对应的频率用fk表示,1≤k≤K;第m个阵元相位中心到位于像素点(x,y)处的点目标的斜距为rm(x,y),1≤m≤M。
图1为本发明流程图,具体实施步骤如下:
第一步,对穿墙雷达阵列天线接收到的步进频率波形回波数据进行基于自由空间传播历程的相位补偿,得到补偿后的成像目标区域的回波数据。
对接收到的每个频点、每个阵元的回波数据进行相位补偿,补偿回波数据自由空间传播路径相位延迟。具体方法为:将回波数据Ek,m(x,y)乘以相位补偿因子得到相位补偿后的回波数据:
(公式1)
第二步,利用补偿后的成像目标区域的回波数据,求成像目标区域每个像素点对应的相位统计参数。该统计参数能表征墙体对穿透成像的影响和阵列稀疏性对步进频率信号雷达成像的影响。
求第一步中得到的相位补偿后的回波数据ECk,m(x,y)的相位统计参数,作为最终图像的加权函数中的参数。具体实现步骤为:
步骤①:求像素点(x,y)所有阵元和频率点的回波数据瞬时相位的余弦函数标准差:
(公式2)
上式中,表示ECk,m(x,y)的相位。
步骤②:求像素点(x,y)所有阵元和频率点的回波数据瞬时相位的正弦函数标准差:
(公式3)
步骤③:由正弦函数标准差和余弦函数标准差定义像素点(x,y)回波数据相位的复数标准差:
(公式4)
步骤④:定义像素点(x,y)的相位统计参数:
(公式5)
第三步,利用标准后向投影成像算法(Back Projection,BP)对补偿后的成像目标区域的回波数据进行成像处理,得到成像目标区域图像。
由于在第一步预处理的过程中已经进行了回波数据相位延迟的补偿,利用BP算法时,对不同阵元不同频率的相位补偿后的回波数据ECk,m(x,y)直接相加就可得到目标区域图像。
其中wk为阵列加权系数,wm为频率加权系数。
第四步,将成像目标区域图像的每个像素点与该像素点对应的相位统计参数相乘,得到改善的成像目标区域图像。
对每个像素点(x,y),将第二步得到的相位统计参数W(x,y)与第三步得到的目标区域图像I(x,y)相乘得到改善的穿墙雷达图像。
Iw(x,y)=I(x,y)·W(x,y) (公式7)
图2至图4是进行实验的有关图像。
图2为进行实验的场景图。图中以雷达阵列方向为方位向,垂直于阵列方向为距离向,雷达阵列中点O为坐标原点,以包含成像目标点T和雷达阵列的平面为坐标平面建立直角坐标系。穿墙雷达天线阵列沿方位向排列,平行于一面厚度为0.3m的建筑物外墙。穿墙雷达天线阵列距离墙体23.7m。一人在建筑物内部,作为成像目标点T,位于(-0.4,28.3)m处。成像目标区域为[-5~5]m×[24~34]m。
图3为传统成像算法获得的穿墙成像结果。图中水平方向为方位向(单位:米),垂直方向为距离向(单位:米)。图像幅度以图像中最强目标为基准进行了归一化处理,单位为分贝(dB)。由图可见,得到的建筑物内部图像中出现了大量“虚假目标”,同时真实目标散焦严重,图像的整体质量不能令人满意。
图4为利用本发明法获得的穿墙成像结果。图中水平方向为方位向(单位:米),垂直方向为距离向(单位:米)。图像幅度以图像中最强目标为基准进行了归一化处理,单位为分贝(dB)。由图可见,采用本发明的方法后,墙体影响和阵列稀疏性影响得到校准,图像中“虚假目标”得到有效抑制,可以清楚看到建筑物内部的人体成像。同时图像整体质量良好。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (1)
1.一种基于相位统计参数加权的穿墙雷达成像方法,其特征在于,包括下述步骤:
第一步,对穿墙雷达阵列天线接收到的步进频率波形回波数据进行基于自由空间传播历程的相位补偿,得到补偿后的成像目标区域的回波数据;
第二步,利用补偿后的成像目标区域的回波数据,求成像目标区域每个像素点对应的相位统计参数;
第三步,利用标准后向投影成像算法对补偿后的成像目标区域的回波数据进行成像处理,得到成像目标区域图像;
第四步,将成像目标区域图像的每个像素点与该像素点对应的相位统计参数相乘,得到改善的成像目标区域图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310233525.5A CN103278818B (zh) | 2013-06-13 | 2013-06-13 | 基于相位统计参数加权的穿墙雷达成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310233525.5A CN103278818B (zh) | 2013-06-13 | 2013-06-13 | 基于相位统计参数加权的穿墙雷达成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103278818A CN103278818A (zh) | 2013-09-04 |
CN103278818B true CN103278818B (zh) | 2014-10-08 |
Family
ID=49061391
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310233525.5A Expired - Fee Related CN103278818B (zh) | 2013-06-13 | 2013-06-13 | 基于相位统计参数加权的穿墙雷达成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103278818B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108387896B (zh) * | 2018-01-03 | 2020-07-07 | 厦门大学 | 一种基于探地雷达回波数据的自动汇聚成像方法 |
CN109541585B (zh) * | 2018-10-09 | 2020-11-06 | 中国人民解放军第四军医大学 | 一种基于峰度评估的人体穿墙检测成像方法 |
CN109471097B (zh) * | 2018-11-14 | 2022-04-22 | 湖南华诺星空电子技术有限公司 | 一种穿墙雷达信号优化处理方法及装置 |
CN109298422B (zh) * | 2018-11-14 | 2022-06-28 | 湖南华诺星空电子技术有限公司 | 一种用于穿墙雷达的合成孔径成像优化处理方法、装置 |
CN109901161A (zh) * | 2019-03-06 | 2019-06-18 | 北京理工大学 | 一种基于多角度的透障雷达多径抑制方法 |
CN113238221B (zh) * | 2021-05-11 | 2023-09-22 | 沈阳航空航天大学 | 基于二维最小相位相干因子的mimo穿墙雷达成像方法 |
CN113341408B (zh) * | 2021-06-02 | 2022-07-01 | 中国人民解放军海军航空大学 | 一种基于穿墙雷达杂波抑制的成像方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102520405A (zh) * | 2011-12-16 | 2012-06-27 | 中国人民解放军国防科学技术大学 | 基于图像域墙体影响补偿的穿墙成像方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9057783B2 (en) * | 2011-01-18 | 2015-06-16 | The United States Of America As Represented By The Secretary Of The Army | Change detection method and system for use in detecting moving targets behind walls, barriers or otherwise visually obscured |
-
2013
- 2013-06-13 CN CN201310233525.5A patent/CN103278818B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102520405A (zh) * | 2011-12-16 | 2012-06-27 | 中国人民解放军国防科学技术大学 | 基于图像域墙体影响补偿的穿墙成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103278818A (zh) | 2013-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103278818B (zh) | 基于相位统计参数加权的穿墙雷达成像方法 | |
Wang | CFAR-based interference mitigation for FMCW automotive radar systems | |
CN110456348B (zh) | 多视向sar海浪谱数据融合的海浪截断波长补偿方法 | |
CN102333052B (zh) | 一种适用于浅海低频条件的水声信号盲解卷方法 | |
CN105137437A (zh) | 一种基于空域相位方差加权的目标检测方法 | |
CN113253213A (zh) | 高频地波雷达电离层回波空间分布特性测量方法及系统 | |
Martone et al. | Noncoherent approach for through-the-wall moving target indication | |
CN109471097B (zh) | 一种穿墙雷达信号优化处理方法及装置 | |
CN105301580A (zh) | 一种基于分裂阵互谱相位差方差加权的被动探测方法 | |
CN111580099A (zh) | 基于联合熵值的穿墙成像雷达墙体杂波抑制方法 | |
Mohsin Riaz et al. | Through‐wall image enhancement based on singular value decomposition | |
CN103558584A (zh) | 一种波达方向的检测方法及装置 | |
CN101592730B (zh) | 基于参数调节随机共振及后处理的传感器阵波束域微弱信号处理方法 | |
CN104020465B (zh) | 基于八单元小孔径圆阵天线的外辐射源雷达测角方法 | |
CN108872947A (zh) | 一种基于子空间技术的海杂波抑制方法 | |
CN110609085B (zh) | 一种基于矢量水听器的声学超材料声学性能测量方法 | |
Lei | A time-domain beamformer for UWB through-wall imaging | |
Lei et al. | Through-wall surveillance using ultra-wideband short pulse radar: Numerical simulation | |
Maaref et al. | Array‐Based Ultrawideband through‐Wall Radar: Prediction and Assessment of Real Radar Abilities | |
CN107728132B (zh) | 一种改善干扰阻塞算法输出信噪比的方法 | |
Anwar et al. | Clutter suppression in through-the-wall radar imaging using enhanced delay-and-sum beamformer | |
Chen et al. | A novel approach to mitigation of radar beam weighting effect on coherent radar imaging using VHF atmospheric radar | |
Shvets et al. | Lightning location and estimate of the lower ionosphere effective height using by dispersion properties of tweek atmospherics | |
Rahaman et al. | Performance Enhancement of Active Sonar System in Under Water Environment Using Spherical Hydrophone Array | |
Jia et al. | Improved cross-correlated back-projection algorithm for through-wall-radar imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141008 Termination date: 20160613 |