穿墙雷达回波数据的处理方法技术

本发明专利技术公开了一种穿墙雷达回波数据的处理方法，包括获取雷达接收信号序列；进行逐列差分得到差分序列；在方向维计算能量得到能量序列；对差分序列逐行归一化得到归一化后的差分序列；在方向维累积求和得到序列；获取大于差分阈值的点，确定目标区间，保护区间和参考区间；判定目标区间是否出现目标并确定目标的初始位置；对序列进行滤波得到预处理后的数据；获得目标位置并显示目标轨迹，完成穿墙雷达回波数据的处理。本发明专利技术无需额外的硬件需求，适用性高，计算速度快，杂波抑制效果好，可应用于穿墙雷达运动目标的实时定位，无人机等动目标探测雷达系统的数据预处理中。

Processing method for echo data of through wall radar

The invention discloses a processing method for the echo data of a through-the-wall radar, which includes acquiring the sequence of radar received signals, obtaining the differential sequence by column difference, calculating the energy sequence in the direction dimension, normalizing the differential sequence row by row, and accumulating the sum sequence in the direction dimension. Target interval, protection interval and reference interval are determined by taking points larger than the difference threshold; target interval is determined and the initial position of the target is determined; pre-processed data is obtained by filtering the sequence; target position is obtained and target trajectory is displayed to complete the processing of radar echo data through the wall. The invention does not need additional hardware requirements, has high applicability, fast calculation speed and good clutter suppression effect, and can be applied to real-time positioning of moving targets of wall-penetrating radar and data pre-processing of moving targets detection radar systems such as unmanned aerial vehicles.

【技术实现步骤摘要】
穿墙雷达回波数据的处理方法
本专利技术具体涉及一种穿墙雷达回波数据的处理方法。
技术介绍
超宽带穿墙雷达是一种新型的对建筑物后的目标进行成像的雷达，由于其具有较高的距离分辨率和良好的穿透能力，在墙体后人质救援、生命探测、室内布局、反恐斗争等领域具有广泛的应用前景。为了抑制墙体，天线等产生的耦合信号，传统的实时处理方法有耦合对齐后再脉冲相消，该方法存在运动人体目标的特征信息无法完整保留、严重的多径虚假像干扰运动目标检测等问题。对于目标初始位置的确定，传统的恒虚警检测算法(CFAR)在一定程度上解决了这个问题，但其是在杂波服从高斯分布的假设下进行的。在实际的数据处理过程中，回波信号中高能量的位置依旧存在于雷达天线的耦合波附近。传统的脉冲对消法和指数平均法，虽然在一定程度上能去掉背景杂波的干扰。但是脉冲对消法存在抗噪声能力差，目标信息无法完整保留等缺陷。类似地，指数平均法存在权重系数无法自适应更新等问题。
技术实现思路
本专利技术的目的在于提供一种能够精准定位目标位置，且适用性好，计算速度快的穿墙雷达回波数据的处理方法。本专利技术提供的这种穿墙雷达回波数据的处理方法，包括如下步骤：S1.获取测量过程中雷达接收信号的所有A-scan序列x(k,n)，其中k＝1,2,...,K，n＝1,2,...,N，K和N分别为一道A-scan的采样点数和该次测量内接收到的A-scan道数；S2.对步骤S1得到的序列x(k,n)，逐列进行差分得到差分序列d(k,n1)，其中k＝1,2,...,K，n1＝1,2,...,N-1；S3.对步骤S2得到的差分序列d(k,n1)，在方向维计算能量得到能量序列E(k)，k＝1,2,...,K；S4.将步骤S2得到的差分序列进行逐行归一化，并得到归一化后的差分序列di(k,n1)，其中k＝1,2,...,K，n1＝1,2,...,N-1；S5.对步骤S4得到的归一化后的差分序列di(k,n1)，在方向维累积求和得到序列S(k)；S6.在步骤S5得到的序列S(k)中，获取大于差分阈值的点，从而确定目标区间，保护区间和参考区间；S7.在步骤S6确定的目标区间，保护区间和参考区间中，根据步骤S3得到的序列E(k)，判定目标区间是否出现目标，并将该目标区间作为目标的初始范围；S8.根据目标出现的范围确定滤波区域，并对滤波区域中序列x(k,n)进行滤波处理，得到预处理后的数据xd(k,n)；S9.根据步骤S8得到的预处理后的数据xd(k,n)获得目标位置并显示目标轨迹，从而完成穿墙雷达回波数据的处理。步骤S2所述的对步骤S1得到的序列x(k,n)逐列进行差分得到差分序列d(k,n1)，具体为采用如下算式计算d(k,n1)：d(k,i)＝x(k,i+1)-x(k,i)式中d(k,i)为差分序列d中第k行的第i列数据，且1≤i≤N-1。步骤S3所述的在方向维计算能量得到能量序列E(k)，具体为采用如下算式计算能量序列E(k)：式中j为能量序列E(k)中的第j个元素，且1≤j≤K。步骤S4所述的得到归一化后的差分序列di(k,n1)，具体为采用如下算式进行归一化：式中di(k,m)为归一化后的差分序列di中第k行的第m列数据，且1≤m≤N-1。步骤S5所述的在方向维累积求和得到序列S(k)，具体为采用如下算式计算序列S(k)：式中S(o)为序列S(k)中的第o个元素，且1≤o≤K。步骤S6所述的获取大于差分阈值的点并确定目标区间，保护区间和参考区间，具体为采用如下规则确定各个区间：R1.在序列S(k)中，找到大于事先设定的差分阈值的点temp；R2.采用如下规则划定目标区间，保护区间和参考区间：目标区间＝(temp-l1,temp+l1)保护区间＝(temp-l1-l2,temp-l1-1)+(temp+l1+1,temp+l1+l2)参考区间＝(temp-l1-l2-l3,temp-l1-l2-1)+(temp+l1+l2+1,temp+l1+l2+l3)其中l1、l2和l3均为事先设定的参数值。步骤S7所述的判定目标区间是否出现目标并将该目标区间作为目标的初始范围，具体为采用如下规则进行判定：若则认定目标出现在目标区间，并将该目标区间作为目标的初始范围；其中β为事先设定的能量系数，range_t表示目标区间，range_r表示参考区间。步骤S8所述的确定滤波区域并对滤波区域中进行滤波处理，具体为采用如下步骤确定滤波区域并进行滤波：A.采用如下公式确定滤波区域reg：式中pos为目标位置，r1、c1、r2和c2均为事先设定的区域值；B.在步骤A得到的滤波区域中，对序列x(k,n)进行二维离散傅里叶正变换：(k0,n0)＝reg(1,1)；u＝0,...,K1-1；v＝0,...,N1-1式中F(u,v)为x(k,n)在滤波区域中的傅里叶正变换，reg是一个二维矩阵，表示选定的滤波区域。K1为滤波区域reg的行宽，N1为滤波区域reg的列宽，reg(1,1)表示二维矩阵的第1行，第1列对应的坐标位置。C.按照如下公式设置二维滤波器：式中G(i,j)为K1×N1维矩阵，Fs为雷达接收信号的采样频率，Fl，Fh分别为距离维带通滤波器的上、下截止频率，其值与雷达发送信号的频谱相关，为向下取整，为向上取整；D.采用步骤C得到的二维滤波器，对步骤B得到的序列进行滤波，从而得到滤波后的矩阵Fd(u,v)＝G(i,j).×F(u,v)。该矩阵中，第u行第v列的值为矩阵G中第i行第j列的值与矩阵F中第u行第v列的值相乘。；E.对步骤D得到的滤波后的序列进行二维离散傅里叶反变换，从而得到最终的预处理后的数据：(k0,n0)＝reg(1,1)；k＝0,...,K1-1；n＝0,...,N1-1步骤S9所述的根据预处理后的数据xd(k,n)获得目标位置并显示目标轨迹，具体为根据预处理后的数据xd(k,n)，采用双站交叉定位算法获得目标位置，并显示目标轨迹，并当下一轮数据采集后，重复步骤S6～S9进行数据预处理和目标定位。本专利技术提供的这种穿墙雷达回波数据的处理方法，基于脉冲波穿墙雷达回波信号的时域表征形态，基于I-CFAR寻找目标的初始位置，并基于目标位置自适应选取二维区域，在区域内做二维FFT滤波，能有效地确定目标初始位置，并实时地滤除噪声和背景杂波的干扰；相对于传统的CFAR检测只适用于背景杂波服从高斯分布的不足，以及脉冲对消法不稳定、无法保留目标完整信息，指数平均无法自适应更新权重系数等缺点，本专利技术无需额外的硬件需求，适用性高，计算速度快，杂波抑制效果好，可应用于穿墙雷达运动目标的实时定位，无人机等动目标探测雷达系统的数据预处理中。附图说明图1为本专利技术方法的方法流程图。图2为本专利技术的实施例中目标真实的运动轨迹的示意图。图3为本专利技术的实施例中一次测量内所有A-scan波形的示意图。图4为本专利技术的实施例中一次测量内差分信号的波形的示意图。图5为本专利技术的实施例中一次测量内差分信号相加累积后的波形的示意图。图6为本专利技术的实施例中一次测量内方向维累积的差分信号能量分布的示意图。图7为本专利技术的实施例中一次测量内的行归一化后差分信号的波形的示意图。图8为本专利技术的实施例中一次测量内的行归一化后差分信号相加累积后的波形的示意图。图9为本专利技术本文档来自技高网...

【技术保护点】
1.一种穿墙雷达回波数据的处理方法，包括如下步骤：S1.获取测量过程中雷达接收信号的所有A‑scan序列x(k,n)，其中k＝1,2,...,K，n＝1,2,...,N，K和N分别为一道A‑scan的采样点数和该次测量内接收到的A‑scan道数；S2.对步骤S1得到的序列x(k,n)，逐列进行差分得到差分序列d(k,n1)，其中k＝1,2,...,K，n1＝1,2,...,N‑1；S3.对步骤S2得到的差分序列d(k,n1)，在方向维计算能量得到能量序列E(k)，k＝1,2,...,K；S4.将步骤S2得到的差分序列进行逐行归一化，并得到归一化后的差分序列di(k,n1)，其中k＝1,2,...,K，n1＝1,2,...,N‑1；S5.对步骤S4得到的归一化后的差分序列di(k,n1)，在方向维累积求和得到序列S(k)；S6.在步骤S5得到的序列S(k)中，获取大于差分阈值的点，从而确定目标区间，保护区间和参考区间；S7.在步骤S6确定的目标区间，保护区间和参考区间中，根据步骤S3得到的序列E(k)，判定目标区间是否出现目标，并将该目标区间作为目标的初始范围；S8.根据目标出现的范围确定滤波区域，并对滤波区域中序列x(k,n)进行滤波处理，得到预处理后的数据xd(k,n)；S9.根据步骤S8得到的预处理后的数据xd(k,n)获得目标位置并显示目标轨迹，从而完成穿墙雷达回波数据的处理。...

【技术特征摘要】
1.一种穿墙雷达回波数据的处理方法，包括如下步骤：S1.获取测量过程中雷达接收信号的所有A-scan序列x(k,n)，其中k＝1,2,...,K，n＝1,2,...,N，K和N分别为一道A-scan的采样点数和该次测量内接收到的A-scan道数；S2.对步骤S1得到的序列x(k,n)，逐列进行差分得到差分序列d(k,n1)，其中k＝1,2,...,K，n1＝1,2,...,N-1；S3.对步骤S2得到的差分序列d(k,n1)，在方向维计算能量得到能量序列E(k)，k＝1,2,...,K；S4.将步骤S2得到的差分序列进行逐行归一化，并得到归一化后的差分序列di(k,n1)，其中k＝1,2,...,K，n1＝1,2,...,N-1；S5.对步骤S4得到的归一化后的差分序列di(k,n1)，在方向维累积求和得到序列S(k)；S6.在步骤S5得到的序列S(k)中，获取大于差分阈值的点，从而确定目标区间，保护区间和参考区间；S7.在步骤S6确定的目标区间，保护区间和参考区间中，根据步骤S3得到的序列E(k)，判定目标区间是否出现目标，并将该目标区间作为目标的初始范围；S8.根据目标出现的范围确定滤波区域，并对滤波区域中序列x(k,n)进行滤波处理，得到预处理后的数据xd(k,n)；S9.根据步骤S8得到的预处理后的数据xd(k,n)获得目标位置并显示目标轨迹，从而完成穿墙雷达回波数据的处理。2.根据权利要求1所述的穿墙雷达回波数据的处理方法，其特征在于步骤S2所述的对步骤S1得到的序列x(k,n)逐列进行差分得到差分序列d(k,n1)，具体为采用如下算式计算d(k,n1)：d(k,i)＝x(k,i+1)-x(k,i)式中d(k,i)为差分序列d中第k行的第i列数据，且1≤i≤N-1。3.根据权利要求2所述的穿墙雷达回波数据的处理方法，其特征在于步骤S3所述的在方向维计算能量得到能量序列E(k)，具体为采用如下算式计算能量序列E(k)：式中j为能量序列E(k)中的第j个元素，且1≤j≤K。4.根据权利要求3所述的穿墙雷达回波数据的处理方法，其特征在于步骤S4所述的得到归一化后的差分序列di(k,n1)，具体为采用如下算式进行归一化：式中di(k,m)为归一化后的差分序列di中第k行的第m列数据，且1≤m≤N-1。5.根据权利要求4所述的穿墙雷达回波数据的处理方法，其特征在于步骤S5所述的在方向维累积求和得到序列S(k)，具体为采用如下算式计算序列S(k)：式中S(o)为序列S(k)中的第o个元素，且1≤o≤K。6.根据权利要求5所述的穿墙雷达回波数据的处理方法，其...

【专利技术属性】
技术研发人员：雷文太，席景春，候斐斐，李若楠，
申请(专利权)人：中南大学，
类型：发明
国别省市：湖南,43