一种适用于穿墙雷达的目标检测方法与流程

本发明涉及穿墙雷达技术领域，尤其涉及一种适用于穿墙雷达的目标检测方法。

背景技术：

穿墙雷达能够对建筑物后的目标进行探测和区分，在生命迹象探测，动目标动作识别、多目标探测方面有着许多的具体应用。在军事上，穿墙雷达可以很好地穿透墙体，对建筑物后面的敌方士兵进行检测识别；在警备安防上，穿墙雷达可以获取建筑物内的犯罪分子的位置信息，在公安武警等执法部门进行侦查和抓捕行动中起到重要作用；在抗震救灾上，穿墙雷达可以灾后现场的建筑物进行探测，可以很好地辅助救灾人员的救援工作；目前国内外对穿墙雷达的已有较多研究，例如以色列的cameron公司开发的xaver系列雷达，针对不同用户需求，有1d手持穿墙雷达,2d成像雷达以及3d成像雷达，其分辨率可达到5cm，产品全面，性能较好，国内的有华诺星空的ce系列雷达，有1d和2d穿墙雷达，分辨率在30cm，性能较好。

专利号为201710632132公开了一种基于多功能超宽带穿墙雷达的动静目标检测方法及系统，包括确定穿墙雷达检测目标的模式，采用接收机的信号，共采集m*n组；m是有效脉冲长度，代表距离向数据；n是脉冲组数，代表速度向数据；该发明能够切换动静目标的检测方式，但不能区分显动目标与微动目标，检测出的结果偏差比较大。

技术实现要素：

本发明要解决的技术问题是提供一种适用于穿墙雷达的目标检测方法，该方法将显动目标和微动目标的检测流程分离，提供不同检测方法以提高目标的检测概率。

【微动目标】表示运动幅度不大于4cm的目标，例如静坐的人或者站立不动的人。

【显动目标】表示运动幅度大于4cm的目标，例如站立晃动或者行走的人。

为了解决上述技术问题，本发明采用以下技术方案：一种适用于穿墙雷达的目标检测方法，包括以下步骤：

s1，通过雷达向目标发射信号波，对接收到的回波信号与发射信号进行混频处理，得到原始回波数据，并对原始回波数据进行预处理；

s2，将预处理后的回波数据通过匹配滤波器进行脉冲压缩，得到回波数据的距离维信息；

s3，将脉冲压缩后的回波数据通过mti（动目标显示）滤波器对静止目标进行抑制，得到运动目标显示数据；

s4，将运动目标显示数据通过多普勒滤波器组，得到运动目标的多普勒频率-距离-角度维数据；

s5，在步骤s4的基础上对运动目标的多普勒频率-距离-角度维数据先进行三维恒虚警检测，然后对检测结果进行聚类，识别显动目标的数量和位置；

s6，在步骤s4的基础上对运动目标的多普勒频率-距离-角度维数据进行凝聚，确定目标的数量和位置，然后对凝聚之后的数据进行三维恒虚警检测，识别出微动目标。

进一步地，所述将原始回波数据预处理是将采集原始回波数据后通过一个幅度检测器，所述幅度检测器用于检测幅度是否饱和对原始回波数据增益进行调整。

进一步地，假设回波为x(t),匹配滤波器冲击响应为h(t),则所述脉冲压缩结果：,其中表示卷积运算。

进一步地，所述运动目标显示数据的计算公式为：

其中hmti(t)表示mti滤波器的单位冲激响应函数，y(t)为脉冲压缩结果。

进一步地，所述mti滤波器优选为iir滤波器。

进一步地，所述运动目标多普勒频率-距离-角度维数据的计算公式为：

其中y^kmtd(t)表示第k个多普勒滤波器的输出数据，表示第k个多普勒滤波器的单位冲激响应函数。

进一步地，所述多普勒滤波器组为fft滤波器组。

进一步地，所述三维恒虚警检测的步骤包括：

s51，求得三维数据内单个数据块近邻域的幅度均值；

s52，根据单个数据领域均值求得单个数据块远邻域的幅度均值；

s53，设定门限值，若单个数据块近邻域的幅度均值与远邻域的均值之间的比值大于门限值则表示单个数据块中有目标，若单个数据块近邻域的幅度均值与远邻域的均值之间的比值小于门限值则表示单个数据块中没有目标；

s54，计算出检测出来的目标个数和位置。

进一步地，所述聚类的步骤包括：

s511，在三维恒虚警检测结果的基础上，构建隶属度函数；

s512，通过隶属度函数构建代价函数；

s513，最小化代价函数即可得到聚类数目和位置。

通常的，会通过最大期望算法最小化代价函数。

进一步地，所述隶属度函数的公式为：

构建隶属度函数

其中，

n表示三维恒虚警检测结果的数量，m表示聚类后显动目标的个数。

进一步地，所述代价函数的公式为：

其中，是聚类表示矩阵，表示元素到聚类中心的距离，n表示三维恒虚警检测结果的数量，m表示聚类后显动目标的个数。

进一步地，通过动态规划算法构建值函数对所述多普勒频率-距离-角度维数据进行慢时间上的凝聚。

进一步地，所述值函数的公式为：

其中表示k时刻目标的值函数，表示k+1时刻的量测，且。

本发明的有益效果：

（1）本发明通过将显动目标和微动目标的进行分离检测，提高了检测精度和微动目标的检测概率。

（2）在将显动目标检测和微动目标检测分离后，针对微动目标增加了凝聚过程，增强了微动目标的信噪比，提高检测概率。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例的流程图；

图2为本发明实施例的数据预处理流程图；

图3为本发明实施例的fft滤波器幅频特性曲线图；

图4为本发明实施例的三维恒虚警立体结构图；

图5为本发明实施例的最大期望算法原理图。

具体实施方式

下面结合附图对发明进一步说明，但不用来限制本发明的范围。

如图1所示，本发明提供了一种适用于穿墙雷达的目标检测方法，包括以下步骤：

s1，通过雷达向目标发射信号波，对接收到的回波信号与发射信号进行混频处理，得到原始回波数据，并对原始回波数据进行预处理；

s2，将预处理后的回波数据通过匹配滤波器进行脉冲压缩，得到回波数据的距离维信息；

s3，将脉冲压缩后的回波数据通过mti滤波器对静止目标进行抑制，得到运动目标显示数据；

s4，将运动目标显示数据通过多普勒滤波器组，得到运动目标的多普勒频率-距离-角度维数据；

s5，在步骤s4的基础上对运动目标的多普勒频率-距离-角度维数据先进行三维恒虚警检测，然后对检测结果进行聚类，识别显动目标的数量和位置；

s6，在步骤s4的基础上对运动目标的多普勒频率-距离-角度维数据进行凝聚，确定目标的数量和位置，然后对凝聚之后的数据进行三维恒虚警检测，识别出微动目标。

如图2所示，若原始回波数据饱和，则会导致信号失真，检测结果不准确，将原始回波数据预处理是保证原始回波数据不饱和，具体的是将a/d采集器采集原始回波数据后通过一个幅度检测器，所述幅度检测器还能对原始回波数据增益进行调整。

通过设置有匹配滤波器对预处理后的回波数据进行脉冲压缩，假设回波为x(t)，匹配滤波器冲击响应为h(t)，则所述脉冲压缩结果：，其中表示卷积运算。

接下来将脉冲压缩后的结果通过动目标显示滤波器（mti），具体的，我们选用iir滤波器作为动目标显示滤波器使用，主要是为了对脉冲压缩后的结果中静目标进行抑制，得到运动目标显示数据。

运动目标显示数据的计算公式为：

其中hmti(t)表示iir滤波器的单位冲激响应函数，y(t)为脉冲压缩结果。

特别地，本实施例中多普勒滤波器组选用fft滤波器组，其幅频特性曲线如图3所示，对于8阶的fft滤波器将输出8组对应不同多普勒频率的距离维数据，其输出数据为：

k为多普勒通道数，y^kmtd(t)表示第k个多普勒滤波器的输出数据，表示第k个多普勒滤波器的单位冲激响应函数，对于本实施例选用的fft滤波器组而言：

其中n表示输入信号的慢时间维长度。

对于显动目标的处理：

对于上述通过fft滤波器组得到的多普勒频率-距离-角度维数，若m为角度维数据，n为距离维数据，k为多普勒通道数，则为n×m维矢量，表示第

k个多普勒滤波器输出数据，故为维数据块，在此基础上进行三维cfar（恒虚警检测），三维cfar立体窗结构如图4所示。

本实施例中三维恒虚警检测的原理是求得在维数据块内单个数据块的幅度值，若单个要判断某个数据块中是否有目标，主要的判断依据是该数据块中的幅度值是否比该数据块远邻域内的幅度值高，本实施例中为了检测的精准度，将求单个数据块的幅度值变成为求数据块近邻域内的幅度均值。

设a,b,c,分别代表检测单元的边长，d,e,f代表参考单元的边长，则对于维数据内单个数据块（i,j,k）近邻域的幅度均值为：

同样，我们可以求得远邻域内的幅度均值（噪声的基底能量）为：

若信杂比scr≥thre1，则表示检测单元内为目标，thre1表示第一预设门限。目标凝聚处理在三维cfar检测的基础上进行，通过对检测到的目标进行聚类得到目标的数量和位置，实例中选用fcm算法，具体实现步骤如下所示：

构建隶属度函数

满足

其中n表示三维恒虚警检测结果的数量，m表示聚类后目标的个数；

同时构建代价函数:

其中，是聚类表示矩阵，表示元素到聚类中心的距离；其中n表示三维恒虚警检测结果的数量，m表述聚类后目标的个数。

最小化代价函数即可得到目标的聚类数目和位置，如图5所示，此处可以使用em算法（最大期望算法）求解目标的聚类数和位置。

所述em算法流程如下所示：

（1）初始化目标点数m和聚类中心位置；

（2）根据初始化的目标点数或者上一步得到的目标点数m的和,求得最优的隶属度函数u:

（3）根据e步求出的求出聚类表示:

（4）求出当前的聚类值函数v:

（5）若，则将目标个数增加后（m=m+1）重复(2)-(4)步骤，直至，则停止迭代，并输出聚类数目m聚类表示矩阵θ。其中t0表示聚类收敛门限。

微动目标处理：

对于步骤s4中得到的多普勒频率-距离-角度维数据先进行慢时间上的凝聚，此处选用动态规划算法（不仅限于该算法）来完成该过程。

首先构建值函数：，其中表示k时刻目标的值函数，表示k+1时刻的量测，且,由此递推出任意时刻的值函数。

在得到任意时刻的值函数之后，需要做的就是确定目标的位置和数量，此时使用cfar（恒虚警检测）检测经过凝聚之后的值函数，若信杂比scr≥thre2，则表示检测单元内为目标，thre2表示第二预设门限。

本发明工作原理：通过雷达向目标发射信号波，对接收到的回波信号与发射信号进行混频处理，得到原始回波数据，并对原始回波数据进行预处理；将预处理后的回波数据通过匹配滤波器进行脉冲压缩，得到回波数据的距离维信息；将脉冲压缩后的回波数据通过mti滤波器对静止目标进行抑制，得到运动目标显示数据；将运动目标显示数据通过多普勒滤波器组，得到运动目标的多普勒频率-距离-角度维数据；

检测显动目标则是对运动目标的多普勒频率-距离-角度维数据先进行三维恒虚警检测，然后对检测结果进行聚类，识别显动目标的数量和位置；

检测微动目标则是对运动目标的多普勒频率-距离-角度维数据进行凝聚，确定目标的数量和位置，然后对凝聚之后的数据进行三维恒虚警检测，识别出微动目标。

本发明提供一种适于穿墙雷达的检测方法，通过将显动目标和微动目标的进行分离检测，提高了检测精度和微动目标的检测概率，在将显动目标检测和微动目标检测分离后，针对微动目标增加了凝聚过程，增强了微动目标的信噪比，提高检测概率。

以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等同物界定。