CN109471097B - 一种穿墙雷达信号优化处理方法及装置 - Google Patents

一种穿墙雷达信号优化处理方法及装置 Download PDF

Info

Publication number
CN109471097B
CN109471097B CN201811353797.8A CN201811353797A CN109471097B CN 109471097 B CN109471097 B CN 109471097B CN 201811353797 A CN201811353797 A CN 201811353797A CN 109471097 B CN109471097 B CN 109471097B
Authority
CN
China
Prior art keywords
aperture
data
imaging
imaging result
compressed
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
Application number
CN201811353797.8A
Other languages
English (en)
Other versions
CN109471097A (zh
Inventor
唐良勇
龚赟
胡俊
王生水
韩乃军
韩明华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huanuo Xingkong Technology Co.,Ltd.
Hunan Huanuo Xingkong Electronic Technology Co ltd
Original Assignee
HUNAN NOVASKY ELECTRONIC TECHNOLOGY CO LTD
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by HUNAN NOVASKY ELECTRONIC TECHNOLOGY CO LTD filed Critical HUNAN NOVASKY ELECTRONIC TECHNOLOGY CO LTD
Priority to CN201811353797.8A priority Critical patent/CN109471097B/zh
Publication of CN109471097A publication Critical patent/CN109471097A/zh
Application granted granted Critical
Publication of CN109471097B publication Critical patent/CN109471097B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开一种穿墙雷达信号优化处理方法及装置，该方法包括步骤包括：S1.使用穿墙雷达对待测区域进行目标检测，接收对应多个孔径的多通道回波信号；S2.从多通道回波信号中分别获取不同孔径的孔径域数据构成多个子孔径数据组，将多通道回波信号对应的全孔径数据以及各子孔径数据组分别进行加权成像后比较加权成像结果，由各加权成像结果中相似部分得到优化处理后的成像数据信号；该装置包括雷达探测模块、信号优化处理模块；本发明具有实现方法简单、成本低、能够实现实时信号处理同时抑制杂波，且优化性能好、效率高等优点。

Description

一种穿墙雷达信号优化处理方法及装置
技术领域
本发明涉及雷达信号处理技术领域，尤其涉及一种穿墙雷达信号优化处理方法及装置。
背景技术
穿墙雷达利用电磁波的穿透性及传输特性，可穿过墙体等非金属介质对墙后区域进行探测，通过对墙后区域的回波信号进行一系列数据处理，可探测跟踪多个隐蔽人体目标，因而被广泛应用于城市巷战、反恐斗争、和灾后救援等军事民事领域。在穿墙雷达探测过程中，由于电磁波在封闭的建筑物内部传播时会发生各种反射、散射和折射，导致回波信号中存在较大旁瓣、栅瓣及多径干扰，从而引起虚警进而影响建筑物内目标的精确检测和判断，因此针对穿墙雷达的探测信号需要优化信号质量，其中关键即为抑制杂波。
现有技术中对雷达信号通常是直接采用后向投影(BP，Back Projection)算法进行成像，但是BP算法对于由环境、目标产生的旁瓣、栅瓣及多径干扰抑制效果差，尤其是在处于较为封闭的室内时易造成虚景，而针对雷达信号的杂波抑制目前通常是基于仿真阶段的处理方法，并不是针对实际探测环境，所取的环境都较为理想，而实际应用时雷达所处探测环境较为复杂，该类雷达信号处理方式的处理性能不佳，并不能有效的滤除实际探测环境中的杂波信号。
有从业者提出使用多径机理分析来实现穿墙雷达室内多径抑制，即首先分析室内多径信号传播模型，得到室内常见的两种多径回波信号分量：目标和墙体之间的多径分量、目标之间的多径分量；然后基于双圆解析表达式求得任意两圆交点的位置，并详细分析交点位置与孔径大小以及孔径中心位置之间的关联性，再通过求解交点聚焦后散射区域的中心位置得到虚假像的位置随子孔径移动而改变的结论；最后利用虚假像的中心位置以及能量在不同位置或不同大小的子孔径内变化大的特点，采用子孔径双层融合法抑制多径虚假像。但是上述雷达处理方法基于多径机理来实现多径抑制，需要依赖模型分析及大量的计算过程，实现复杂、成本较高且效率低，不适用于实时性要求高的场合。
发明内容
本发明要解决的技术问题就在于：针对现有技术存在的技术问题，本发明提供一种实现方法简单、成本低、能够实现实时信号处理同时抑制杂波，且优化性能好、效率高的穿墙雷达信号优化处理方法及装置。
为解决上述技术问题，本发明提出的技术方案为：
一种穿墙雷达信号优化处理方法，步骤包括：
S1.雷达探测：使用穿墙雷达对待测区域进行目标检测，接收对应多个孔径的多通道回波信号；
S2.信号优化处理：从所述多通道回波信号中分别获取不同孔径的孔径域数据构成多个子孔径数据组，将所述多通道回波信号对应的全孔径数据以及各所述子孔径数据组分别进行加权成像后比较各加权成像结果，由各所述加权成像结果中相似部分得到优化处理后的成像数据信号。
作为本发明方法的进一步改进，所述步骤S2中进行加权成像包括：将所述全孔径数据进行成像后得到全孔径成像结果，以及将各所述子孔径数据组进行成像后，对应得到多个压缩孔径成像结果，所述全孔径成像结果、各所述压缩孔径成像结果分别与对应的加权因子进行加权后，得到全孔径加权成像结果以及各压缩孔径加权成像结果。
作为本发明方法的进一步改进，还包括配置所述加权因子的步骤，所述配置所述加权因子的步骤包括：根据孔径数据的成像结果以及对应的孔径数量确定得到初始加权因子，取所述初始加权因子的n次方得到最终的加权因子，其中n>1。
作为本发明方法的进一步改进，所述全孔径成像结果的初始加权因子按照式
Figure BDA0001865559870000021
计算得到；
第l个所述压缩孔径成像结果的初始加权因子按照式
Figure BDA0001865559870000022
计算得到；
其中，zk(x,y)为第k个孔径的孔径成像结果，K为雷达孔径数量，L为每个子孔径数据组中子孔径数量。
作为本发明方法的进一步改进：取所述初始加权因子的3次方得到最终的加权因子。
作为本发明方法的进一步改进，所述步骤S2中采用BP算法进行成像，步骤包括：
S21.将成像区域划分为多个像素点；
S22.计算每一个像素点与所有收发天线组合之间的双程时延，以及根据各像素点到发射天线、接收天线的距离计算每一个像素点的相位补偿；
S23.根据计算得到的各像素点与所有收发天线组合之间的所述双程时延、各像素点的所述相位补偿，计算得到每个孔径数据的后向投影成像结果。
作为本发明方法的进一步改进，所述步骤S2中比较各加权成像结果步骤包括：分别将全孔径加权成像结果中每个像素点与各压缩孔径加权成像结果中对应像素点进行比较，所述全孔径加权成像结果由所述全孔径数据进行加权成像得到，各所述压缩孔径加权成像结果由各所述子孔径数据组进行加权得到，每次比较时取值最小的像素点输出，得到所述优化处理后的成像数据信号。
作为本发明方法的进一步改进，所述比较各加权成像结果具体步骤包括：将所述全孔径加权成像结果中每个像素点与第一组压缩孔径加权成像结果中对应像素点进行比较，每次比较时取值最小的像素点输出，得到第一次比较后的成像结果；将所述第一次比较后的成像结果中每个像素点与第二组压缩孔径加权成像结果中对应像素点进行比较，每次比较时取值最小的像素点输出，得到第二次比较后的成像结果，以此类推，直至完成所有压缩孔径加权成像结果的比较，得到所述优化处理后的成像数据信号。
作为本发明方法的进一步改进，所述步骤S2中选取相同数量的不同子孔径数据构成各个子孔径数据组，具体包括：从雷达的K个孔径中分别选取L个不同孔径数据形成L个子孔径数据组，每个子孔径数据组对应构造得到一个压缩孔径，第l个压缩孔径为：
Al＝<al1,al2,...,alL>
其中，alp为第p个孔径位置信息的索引矢量，1≤alp≤K，且当p≠q时，alp≠alq。
一种穿墙雷达信号处理装置，包括：
雷达探测模块，用于使用穿墙雷达对待测区域进行目标检测，接收对应多个孔径的多通道回波信号；
信号优化处理模块，用于从所述多通道回波信号中分别获取不同孔径的孔径域数据构成多组子孔径数据，将所述多通道回波信号对应的全孔径数据以及各所述子孔径数据分别进行加权成像后比较各加权成像结果，由各所述加权成像结果中相似部分得到优化处理后的成像数据信号。
与现有技术相比，本发明的优点在于：
1、本发明穿墙雷达信号优化处理方法及装置，通过对穿墙雷达实测信号进行实时处理，由全孔径和不同压缩孔径形成的图像数据进行比较，来获取各图像数据中相似部分作为最终的成像数据，可以充分利用不同压缩孔径的信息实现高效的穿墙雷达信号优化处理，同时有效抑制目标杂波。
2、本发明穿墙雷达信号优化处理方法及装置，利用不同压缩孔径形成数据图像进行最小值提取，通过使得每个像素点取全孔径加权成像结果与各压缩孔径加权成像结果中最小值，充分利用孔径信息的同时实现最小旁瓣抑制处理，可以充分利用孔径信息的同时有效抑制杂波。
3、本发明穿墙雷达信号优化处理方法及装置，通过取相干因子CF的n次方作为各孔径成像结果的加权因子，可以使得目标成像结果更为聚焦，进一步提高优化处理性能；进一步基于相干因子CF的3次方对各孔径成像结果进行加权，可以获得目标成像结果最佳性能的聚焦，进而可以实现最佳性能的优化处理。
附图说明
图1是本实施例穿墙雷达信号优化处理方法的实现流程示意图。
图2是本发明具体应用实施例中实现信号优化处理的实现流程示意图。
图3是本发明具体应用实施例中得到的各孔径成像结果示意图。
图4是本发明具体应用实施例中得到的各孔径加权成像结果示意图。
图5是本发明具体应用实施例中得到的各次迭代比较处理后的结果示意图。
具体实施方式
以下结合说明书附图和具体优选的实施例对本发明作进一步描述，但并不因此而限制本发明的保护范围。
如图1所示，本实施例穿墙雷达信号优化处理方法步骤包括：
S1.雷达探测：使用穿墙雷达对待测区域进行目标检测，接收对应多个孔径的多通道回波信号；
S2.信号优化处理：从多通道回波信号中分别获取不同孔径的孔径域数据构成多个子孔径数据组，即构成多个不同的压缩孔径，将多通道回波信号对应的全孔径数据以及各子孔径数据组分别进行加权成像后比较各加权成像结果，由各加权成像结果中相似部分得到优化处理后的成像数据信号。
由于在图像形成过程中，目标的主要响应是连贯的且被标准化，由不同稀疏数据组生成的图像中目标的位置和峰值在是相同的，同时在不同稀疏数组中缺失的数据是不同的，因此旁瓣的位置和峰值与谷值的幅度也不一样。本实施例通过对穿墙雷达实测信号进行实时处理，由全孔径和不同压缩孔径形成的图像数据进行比较来获取最终的成像数据，可以充分利用不同压缩孔径的信息实现穿墙雷达信号优化处理，有效抑制目标杂波。
本实施例具体将穿墙雷达安装在待测区域外部的墙体附近，雷达等效孔径总数为K，由穿墙雷达对待测区域进行目标检测，通过雷达系统获得经过信号预处理后的多通道回波信号echos(t)，从多通道回波信号echos(t)的实测数据中选取不同孔径的孔径域数据来构建压缩孔径。
本实施例步骤S2中选取相同数量的不同子孔径数据构成各个子孔径数据组，具体包括：从雷达的K个孔径中分别选取L个不同孔径数据形成L个子孔径数据组，每个子孔径数据组对应构造得到一个压缩孔径，即图像形成时选取部分孔径位置信息，对实测数据选取多组相同数量、不同孔径的孔径域的数据构造新的随机压缩孔径，其中所选取的各子孔径不重复，每个压缩孔径即为从成像的k个孔径中选取L个孔径形成，L的值取为：
L＝p·K (1)
其中0＜p＜1，构造压缩孔径的个数具体为H。
选取部分孔径构成压缩孔径进行成像后，丢弃其它的K-L个孔径位置信息，其中第l个压缩孔径表示为：
Al＝<al1,al2,...,alL> (2)
其中，alp为第p个孔径位置信息的索引矢量，1≤alp≤K，且当p≠q时，alp≠alq。
本实施例步骤S2中进行加权成像包括：将全孔径数据进行成像后得到全孔径成像结果，以及将各子孔径数据组进行成像后，对应得到多个压缩孔径成像结果，全孔径成像结果、各压缩孔径成像结果分别与对应的加权因子进行加权后，得到全孔径加权成像结果以及各压缩孔径加权成像结果。通过将全孔径、各压缩孔径数据进行加权成像后比较，可以利用各压缩孔径数据的加权成像结果去除雷达图像信号中的杂波，实现杂波抑制。
本实施例步骤S2中对全孔径以及压缩孔径数据分别采用后向投影BP算法进行成像，采用BP算法进行成像步骤包括：
S21.将成像区域划分为多个像素点；
S22.计算每一个像素点与所有收发天线组合之间的双程时延，以及根据各像素点到发射天线、接收天线的距离计算每一个像素点的相位补偿；
S23.根据计算得到的各像素点与所有收发天线组合之间的双程时延、各像素点的相位补偿，计算得到每个孔径数据的后向投影成像结果。
在具体应用实施例中，对全孔径以及压缩孔径数据采用后向投影BP算法进行成像的详细流程为：
①首先将成像区域划分为M×N个像素点，其中ym(m＝1,2,...,M)和xn(n＝1,2,...,N)分别表示距离向和方位向上像素点的坐标值；
②计算每一个像素点(xn,ym)与所有收发天线组合之间的双程时延：
Figure BDA0001865559870000061
其中，k＝1,2,...,K表示第k个孔径，(xT(k),yT(k))表示发射天线坐标，(xR(k),yR(k))表示接收天线坐标，，c为光速且具体c＝3.0×108m/s；
Figure BDA0001865559870000062
表示像素点到发射天线的距离；
Figure BDA0001865559870000063
表示像素点到接收天线的距离；
③计算每一个像素点的相位补偿：
phasek(xn,ym)＝exp(j2πfc·(dT+dR)/c) (4)
其中，fc为载频；
④计算每个孔径后向投影成像结果：
zk(x,y)＝echos(t-τk(x,y))×phasek(x,y) (5)
通过上述得到全孔径数据、各压缩孔径数据的后向投影成像结果后，分别与对应的加权因子进行加权，得到对应的加权成像结果。
本实施例中，还包括配置加权因子的步骤，配置加权因子的步骤包括：根据孔径数据的成像结果以及对应的孔径数量确定得到初始加权因子，取初始加权因子的n次方得到最终的加权因子，其中n>1。
本实施例全孔径成像结果的初始加权因子具体按照式(6)计算得到，第l个压缩孔径成像结果的初始加权因子具体按照式(7)计算得到。
Figure BDA0001865559870000064
Figure BDA0001865559870000065
取式(6)的n次方作为对应全孔径成像结果的加权因子、取式(7)的n次方作为对应压缩孔径成像结果的加权因子，相比于传统的直接使用式(6)、(7)作为加权因子，可以使得目标成像结果更为聚焦，进一步提高优化处理性能。
本实施例具体取初始加权因子的3次方得到最终的加权因子，即对应全孔径成像结果的加权因子为：
Figure BDA0001865559870000071
第l个构造的压缩孔径成像结果的加权因子为：
Figure BDA0001865559870000072
则计算全孔径加权成像结果为：
Figure BDA0001865559870000073
计算子孔径加权成像结果为：
Figure BDA0001865559870000074
l表示第第l个构造压缩孔径。
基于相干因子CF的3次方对各孔径成像结果进行加权，可以获得目标成像结果最佳性能的聚焦，进而可以实现最佳性能的优化处理。
本实施例步骤S2中比较各加权成像结果步骤包括：分别将全孔径加权成像结果中每个像素点与各压缩孔径加权成像结果中对应像素点进行比较，全孔径加权成像结果由全孔径数据进行加权成像得到，各压缩孔径加权成像结果由各子孔径数据组进行加权得到，每次比较时取值最小的像素点输出，得到优化处理后的成像数据信号。本实施例利用不同压缩孔径形成数据图像进行最小值提取，通过使得每个像素点取全孔径加权成像结果与各压缩孔径加权成像结果中最小值，可以充分利用孔径信息的同时实现最小旁瓣抑制处理，有效抑制杂波。
本实施例中，比较各加权成像结果具体步骤包括：将全孔径加权成像结果中每个像素点与第一组压缩孔径加权成像结果中对应像素点进行比较，每次比较时取值最小的像素点输出，得到第一次比较后的成像结果；将第一次比较后的成像结果中每个像素点与第二组压缩孔径加权成像结果中对应像素点进行比较，每次比较时取值最小的像素点输出，得到第二次比较后的成像结果，以此类推，直至完成所有压缩孔径加权成像结果的比较，得到优化处理后的成像数据信号。
本实施例具体将得到的全孔径成像数据I(x,y)与第l个构造的压缩子孔径成像数据Il(x,y)对应像素按照式(12)进行比较，找出每个像素点的最小值构成新的成像数据
I(xn,ym)＝min<I(xn,ym)，Il(xn,ym)>,1≤l≤H
按照上述迭代H次后就可得到经过优化后可抑制杂波的成像数据信号。
如图2所示，本发明在具体应用实施例中首先通过穿墙雷达对待测区域进行探测，接收多通道回波数据；接收到回波数据后进行信号预处理，将经过信号预处理的雷达全孔径回波数据分成全孔径和H个由L个不同孔径构成的压缩孔径；对全孔径数据、H个孔径构成的压缩孔径数据采用BP算法进行成像，得到的全孔径BP成像结果、压缩孔径BP成像结果如图3所示，其中图3(a)对应为全孔径BP成像结果，图3(b)为第l个压缩孔径BP成像结果；计算加权因子CF3，并使用加权因子CF3对图3所示的全孔径原始BP成像结果和压缩孔径原始BP成像结果进行加权，得到对应的全孔径加权成像结果、压缩孔径加权成像结果，如图4所示，其中图4(a)对应为全孔径加权成像结果，图4(b)对应为第l个压缩孔径加权成像结果；然后依次对全孔径加权成像结果、各压缩孔径加权成像结果进行最小旁瓣抑制处理，先对如图4(a)所示的全孔径加权成像结果和第一个子孔径加权成像结果采用最小像素比较法得到第1次迭代结果，第二个子孔径加权成像结果再与第1次迭代结果采用最小像素比较法进行比较得到新的结果，以此类推，经过H次比较后得到最终的优化后处理结果，即得到抑制杂波后的最佳目标成像结果，如图5所示，其中图5(a)对应为第1次迭代结果，图5(b)对应为第2次迭代结果，图(c)对应为最终优化后处理结果。
本实施例穿墙雷达信号处理装置，包括：
雷达探测模块，用于使用穿墙雷达对待测区域进行目标检测，接收对应多个孔径的多通道回波信号；
信号优化处理模块，用于从多通道回波信号中分别获取不同孔径的孔径域数据构成多组子孔径数据，将多通道回波信号对应的全孔径数据以及各子孔径数据分别进行加权成像后比较各加权成像结果，由各加权成像结果中相似部分得到优化处理后的成像数据信号。
本实施例信号优化处理模块中进行加权成像包括：将全孔径数据进行成像后得到全孔径成像结果，以及将各子孔径数据组进行成像后，对应得到多个压缩孔径成像结果，全孔径成像结果、各压缩孔径成像结果分别与对应的加权因子进行加权后，得到全孔径加权成像结果以及各压缩孔径加权成像结果。
本实施例信号优化处理模块还包括配置加权因子，包括：根据孔径数据的成像结果以及对应的孔径数量确定得到初始加权因子，取初始加权因子的n次方得到最终的加权因子，其中n>1。
本实施例中，信号优化处理模块中采用BP算法进行成像，包括：
第一单元，用于将成像区域划分为多个像素点；
第二单元，用于计算每一个像素点与所有收发天线组合之间的双程时延，以及根据各像素点到发射天线、接收天线的距离计算每一个像素点的相位补偿；
第三单元，用于根据计算得到的各像素点与所有收发天线组合之间的双程时延、各像素点的相位补偿，计算得到每个孔径数据的后向投影成像结果。
本实施例信号优化处理模块中比较各加权成像结果步骤包括：分别将全孔径加权成像结果中每个像素点与各压缩孔径加权成像结果中对应像素点进行比较，全孔径加权成像结果由全孔径数据进行加权成像得到，各压缩孔径加权成像结果由各子孔径数据组进行加权得到，每次比较时取值最小的像素点输出，得到优化处理后的成像数据信号。
本实施例信号优化处理模块中比较各加权成像结果具体步骤包括：将全孔径加权成像结果中每个像素点与第一组压缩孔径加权成像结果中对应像素点进行比较，每次比较时取值最小的像素点输出，得到第一次比较后的成像结果；将第一次比较后的成像结果中每个像素点与第二组压缩孔径加权成像结果中对应像素点进行比较，每次比较时取值最小的像素点输出，得到第二次比较后的成像结果，以此类推，直至完成所有压缩孔径加权成像结果的比较，得到优化处理后的成像数据信号。
本实施例穿墙雷达信号处理装置与上述穿墙雷达信号处理方法为一一对应的装置，具体实现原理、所能达到的技术效果均相同，在此不再一一赘述。
上述只是本发明的较佳实施例，并非对本发明作任何形式上的限制。虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明。因此，凡是未脱离本发明技术方案的内容，依据本发明技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均应落在本发明技术方案保护的范围内。

Claims (9)

1.一种穿墙雷达信号优化处理方法，其特征在于，步骤包括：
S1.雷达探测：使用穿墙雷达对待测区域进行目标检测，接收对应多个孔径的多通道回波信号；
S2.信号优化处理：从所述多通道回波信号中分别获取不同孔径的孔径域数据构成多个子孔径数据组，将所述多通道回波信号对应的全孔径数据以及各所述子孔径数据组分别进行加权成像后比较各加权成像结果，由各所述加权成像结果中相似部分得到优化处理后的成像数据信号；
所述步骤S2中进行加权成像包括：将所述全孔径数据进行成像后得到全孔径成像结果，以及将各所述子孔径数据组进行成像后，对应得到多个压缩孔径成像结果，所述全孔径成像结果、各所述压缩孔径成像结果分别与对应的加权因子进行加权后，得到全孔径加权成像结果以及各压缩孔径加权成像结果。
2.根据权利要求1所述的穿墙雷达信号优化处理方法，其特征在于，还包括配置所述加权因子的步骤，所述配置所述加权因子的步骤包括：根据孔径数据的成像结果以及对应的孔径数量确定得到初始加权因子，取所述初始加权因子的n次方得到最终的加权因子，其中n>1。
3.根据权利要求2所述的穿墙雷达信号优化处理方法，其特征在于，所述全孔径成像结果的初始加权因子按照式
Figure FDA0003500430710000011
计算得到；
第l个所述压缩孔径成像结果的初始加权因子按照式
Figure FDA0003500430710000012
计算得到；
其中，zk(x,y)为第k个孔径的孔径成像结果，K为雷达孔径数量，L为每个子孔径数据组中子孔径数量；al1为第1个孔径位置信息的索引矢量，alL为第L个孔径位置信息的索引矢量。
4.根据权利要求2或3所述的穿墙雷达信号优化处理方法，其特征在于：取所述初始加权因子的3次方得到最终的加权因子。
5.根据权利要求1～3中任意一项所述的穿墙雷达信号优化处理方法，其特征在于，所述步骤S2中采用BP算法进行成像，步骤包括：
S21.将成像区域划分为多个像素点；
S22.计算每一个像素点与所有收发天线组合之间的双程时延，以及根据各像素点到发射天线、接收天线的距离计算每一个像素点的相位补偿；
S23.根据计算得到的各像素点与所有收发天线组合之间的所述双程时延、各像素点的所述相位补偿，计算得到每个孔径数据的后向投影成像结果。
6.根据权利要求1～3中任意一项所述的穿墙雷达信号优化处理方法，其特征在于，所述步骤S2中比较各加权成像结果步骤包括：分别将全孔径加权成像结果中每个像素点与各压缩孔径加权成像结果中对应像素点进行比较，所述全孔径加权成像结果由所述全孔径数据进行加权成像得到，各所述压缩孔径加权成像结果由各所述子孔径数据组进行加权得到，每次比较时取值最小的像素点输出，得到所述优化处理后的成像数据信号。
7.根据权利要求6所述的穿墙雷达信号优化处理方法，其特征在于，所述比较各加权成像结果具体步骤包括：将所述全孔径加权成像结果中每个像素点与第一组压缩孔径加权成像结果中对应像素点进行比较，每次比较时取值最小的像素点输出，得到第一次比较后的成像结果；将所述第一次比较后的成像结果中每个像素点与第二组压缩孔径加权成像结果中对应像素点进行比较，每次比较时取值最小的像素点输出，得到第二次比较后的成像结果，以此类推，直至完成所有压缩孔径加权成像结果的比较，得到所述优化处理后的成像数据信号。
8.根据权利要求1～3中任意一项所述的穿墙雷达信号优化处理方法，其特征在于，所述步骤S2中选取相同数量的不同子孔径数据构成各个子孔径数据组，具体包括：从雷达的K个孔径中分别选取L个不同孔径数据形成L个子孔径数据组，每个子孔径数据组对应构造得到一个压缩孔径，第l个压缩孔径为：
Al＝<al1,al2,...,alL>
其中，alp为第p个孔径位置信息的索引矢量，1≤alp≤K，且当p≠q时，alp≠alq。
9.一种穿墙雷达信号优化处理装置，其特征在于，包括：
雷达探测模块，用于使用穿墙雷达对待测区域进行目标检测，接收对应多个孔径的多通道回波信号；
信号优化处理模块，用于从所述多通道回波信号中分别获取不同孔径的孔径域数据构成多组子孔径数据，将所述多通道回波信号对应的全孔径数据以及各所述子孔径数据分别进行加权成像后比较各加权成像结果，由各所述加权成像结果中相似部分得到优化处理后的成像数据信号；所述信号优化处理模块中进行加权成像包括：将所述全孔径数据进行成像后得到全孔径成像结果，以及将各所述子孔径数据组进行成像后，对应得到多个压缩孔径成像结果，所述全孔径成像结果、各所述压缩孔径成像结果分别与对应的加权因子进行加权后，得到全孔径加权成像结果以及各压缩孔径加权成像结果。
CN201811353797.8A 2018-11-14 2018-11-14 一种穿墙雷达信号优化处理方法及装置 Active CN109471097B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811353797.8A CN109471097B (zh) 2018-11-14 2018-11-14 一种穿墙雷达信号优化处理方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811353797.8A CN109471097B (zh) 2018-11-14 2018-11-14 一种穿墙雷达信号优化处理方法及装置

Publications (2)

Publication Number Publication Date
CN109471097A CN109471097A (zh) 2019-03-15
CN109471097B true CN109471097B (zh) 2022-04-22

Family

ID=65672891

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811353797.8A Active CN109471097B (zh) 2018-11-14 2018-11-14 一种穿墙雷达信号优化处理方法及装置

Country Status (1)

Country Link
CN (1) CN109471097B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111025256A (zh) * 2019-12-26 2020-04-17 湖南华诺星空电子技术有限公司 一种机载雷达的微弱生命体征信号的检测方法及系统
CN113238221B (zh) * 2021-05-11 2023-09-22 沈阳航空航天大学 基于二维最小相位相干因子的mimo穿墙雷达成像方法
CN113341408B (zh) * 2021-06-02 2022-07-01 中国人民解放军海军航空大学 一种基于穿墙雷达杂波抑制的成像方法及系统
CN114296074B (zh) * 2021-12-29 2022-10-14 湖南华诺星空电子技术有限公司 基于飞行平台的穿透式侦察雷达成像方法及雷达系统

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2317334A2 (en) * 2009-11-02 2011-05-04 Furuno Electric Co., Ltd. Method and program for displaying radar image, and target object detection signal processing device
CN102721965A (zh) * 2012-06-08 2012-10-10 北京航空航天大学 用于合成孔径雷达成像处理的全向旁瓣抑制方法
CN102721964A (zh) * 2012-04-28 2012-10-10 南京航空航天大学 一种多子孔径自聚焦方法
CN103278818A (zh) * 2013-06-13 2013-09-04 中国人民解放军国防科学技术大学 基于相位统计参数加权的穿墙雷达成像方法
CN103885061A (zh) * 2014-03-27 2014-06-25 西安电子科技大学 基于改进相位梯度自聚焦的双基sar运动补偿方法
CN104502911A (zh) * 2014-12-25 2015-04-08 湖南华诺星空电子技术有限公司 一种穿墙成像雷达的墙壁参数估计方法
KR20160038464A (ko) * 2014-09-30 2016-04-07 국방과학연구소 고해상도 거리 측면도와 역합성 개구면 레이더 영상의 융합을 통한 표적 식별 방법 및 그 장치
CN105954745A (zh) * 2016-04-29 2016-09-21 电子科技大学 一种适用于穿墙雷达多径幻象抑制的成像方法
CN106772365A (zh) * 2016-11-25 2017-05-31 南京理工大学 一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法
CN108120980A (zh) * 2017-12-13 2018-06-05 南京航空航天大学 一种星载sar多模式成像信号处理算法的fpga的实现方法
CN108802725A (zh) * 2017-04-27 2018-11-13 南京理工大学 一种浅层穿透雷达合成孔径成像方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7914454B2 (en) * 2004-06-25 2011-03-29 Wilk Ultrasound Of Canada, Inc. Real-time 3D ultrasonic imaging apparatus and method
US20120235849A1 (en) * 2010-11-19 2012-09-20 Tatoian James Z Method and apparatus for through the wall radar imaging

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2317334A2 (en) * 2009-11-02 2011-05-04 Furuno Electric Co., Ltd. Method and program for displaying radar image, and target object detection signal processing device
CN102721964A (zh) * 2012-04-28 2012-10-10 南京航空航天大学 一种多子孔径自聚焦方法
CN102721965A (zh) * 2012-06-08 2012-10-10 北京航空航天大学 用于合成孔径雷达成像处理的全向旁瓣抑制方法
CN103278818A (zh) * 2013-06-13 2013-09-04 中国人民解放军国防科学技术大学 基于相位统计参数加权的穿墙雷达成像方法
CN103885061A (zh) * 2014-03-27 2014-06-25 西安电子科技大学 基于改进相位梯度自聚焦的双基sar运动补偿方法
KR20160038464A (ko) * 2014-09-30 2016-04-07 국방과학연구소 고해상도 거리 측면도와 역합성 개구면 레이더 영상의 융합을 통한 표적 식별 방법 및 그 장치
CN104502911A (zh) * 2014-12-25 2015-04-08 湖南华诺星空电子技术有限公司 一种穿墙成像雷达的墙壁参数估计方法
CN105954745A (zh) * 2016-04-29 2016-09-21 电子科技大学 一种适用于穿墙雷达多径幻象抑制的成像方法
CN106772365A (zh) * 2016-11-25 2017-05-31 南京理工大学 一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法
CN108802725A (zh) * 2017-04-27 2018-11-13 南京理工大学 一种浅层穿透雷达合成孔径成像方法
CN108120980A (zh) * 2017-12-13 2018-06-05 南京航空航天大学 一种星载sar多模式成像信号处理算法的fpga的实现方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Grating lobe suppression in sparse array-based ultrawideband through-wall imaging radar;Xin Tu 等;《IEEE Antennas and Wireless Propagation Letters》;20151012;第1020-1023页 *
基于WA-EMD算法的脉冲式超宽带雷达多目标生命体征检测;唐良勇 等;《南京理工大学学报》;20170430;第41卷(第2期);第198-206页 *
基于子孔径方位特性的穿墙雷达虚拟多视角成像方法;任楠 等;《电讯技术》;20180331;第58卷(第3期);第276-283页 *
穿墙雷达室内多径机理分析与抑制方法;申文婷 等;《雷达科学与技术》;20161231;第14卷(第6期);第605-613页 *

Also Published As

Publication number Publication date
CN109471097A (zh) 2019-03-15

Similar Documents

Publication Publication Date Title
CN109471097B (zh) 一种穿墙雷达信号优化处理方法及装置
Leigsnering et al. Multipath exploitation in through-the-wall radar imaging using sparse reconstruction
CN105549011B (zh) 一种基于mimo穿墙雷达的建筑物单边两点斜视成像方法
CN107678028B (zh) 低信噪比条件下的微波凝视关联成像方法
CN110596646B (zh) 一种基于mimo体制的提高雷达角度分辨率的布局及方法
CN109298418B (zh) 基于建筑物内部结构特征的雷达检测虚警抑制方法及装置
CN106324595B (zh) 一种基于多角度探测成像的多径假目标抑制方法
CN107561507B (zh) 一种外辐射源雷达的杂波对消方法
CN108776336A (zh) 一种基于emd的自适应穿墙雷达静止人体目标定位方法
CN103018727A (zh) 一种基于样本训练的机载雷达非平稳杂波抑制方法
CN111580099A (zh) 基于联合熵值的穿墙成像雷达墙体杂波抑制方法
CN109765529B (zh) 一种基于数字波束形成的毫米波雷达抗干扰方法及系统
CN109298422B (zh) 一种用于穿墙雷达的合成孔径成像优化处理方法、装置
CN109298417B (zh) 一种基于雷达信号处理的建筑物内部结构探测方法及装置
CN103543451A (zh) 一种基于压缩感知的多径虚像抑制sar后处理方法
CN105954745A (zh) 一种适用于穿墙雷达多径幻象抑制的成像方法
CN113238221B (zh) 基于二维最小相位相干因子的mimo穿墙雷达成像方法
CN108845318B (zh) 基于Relax算法的星载高分宽幅成像方法
Ma et al. Spatiotemporal two-dimensional deconvolution beam imaging technology
Jagadesh et al. Modeling Target Detection and Performance Analysis of Electronic Countermeasures for Phased Radar.
CN116930963A (zh) 一种基于无线通信系统的穿墙成像方法
CN103969647A (zh) 基于图像域滤波的稀疏阵列建筑布局成像方法
Leigsnering et al. CS based wall ringing and reverberation mitigation for through-the-wall radar imaging
US7876256B2 (en) Antenna back-lobe rejection
Boudamouz et al. Through the wall radar imaging with mimo beamforming processing-simulation and experimental results

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
CP01 Change in the name or title of a patent holder

Address after: 410205 B7 Lugu Enterprise Plaza, 27 Wenxuan Road, Changsha High-tech Development Zone, Hunan Province

Patentee after: Huanuo Xingkong Technology Co.,Ltd.

Address before: 410205 B7 Lugu Enterprise Plaza, 27 Wenxuan Road, Changsha High-tech Development Zone, Hunan Province

Patentee before: Hunan Huanuo Xingkong Electronic Technology Co.,Ltd.

Address after: 410205 B7 Lugu Enterprise Plaza, 27 Wenxuan Road, Changsha High-tech Development Zone, Hunan Province

Patentee after: Hunan Huanuo Xingkong Electronic Technology Co.,Ltd.

Address before: 410205 B7 Lugu Enterprise Plaza, 27 Wenxuan Road, Changsha High-tech Development Zone, Hunan Province

Patentee before: HUNAN NOVASKY ELECTRONIC TECHNOLOGY Co.,Ltd.

CP01 Change in the name or title of a patent holder