首页 / 专利库 / 电信 / 波达方向 / 基于生命体征SAR成像的人体目标身份识别方法及系统

基于生命体征SAR成像的人体目标身份识别方法及系统

阅读:471发布:2023-03-03

专利汇可以提供基于生命体征SAR成像的人体目标身份识别方法及系统专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种基于生命体征SAR成像的人体目标身份识别方法及系统,方法包括:采集待识别人员的呼吸 波形 ,提取呼吸波形的特征参数,并输入分类器中进行训练,获得身份识别模型;利用基于 相位 解调的生命体征SAR成像方法确定场景中人体目标的准确 位置 ;将雷达移动到确定好的人体目标位置正前方获取人体目标的呼吸波形;将待识别的呼吸波形输入到身份识别模型中,即可准确得到人体目标的身份。系统包括呼吸波形提取模 块 、特征提取模块、分类模型建立模块、成像模块、人体目标 定位 模块、身份识别模块。本发明方法的身份识别准确率高,可以准确获得场景中人体目标的位置及身份信息,鲁棒性好,且适用性更广。,下面是基于生命体征SAR成像的人体目标身份识别方法及系统专利的具体信息内容。

1.一种基于生命体征SAR成像的人体目标身份识别方法,其特征在于,包括以下步骤:
步骤1、将FMCW雷达静止置于待识别人员的正前方L米处,采集雷达回波信号,并对回波信号进行预处理获取呼吸波形Re(τa);
步骤2、提取步骤1得到的每一个呼吸波形Re(τa)的特征参数;
步骤3、根据步骤2中提取的呼吸波形的特征参数建立身份识别模型;
步骤4、待识别人员随机分布在某一场景中,FMCW雷达移动平台在所述场景内按直线轨迹运行采集回波信号,并对回波信号预处理后进行SAR成像,获得成像结果I;
步骤5、在步骤4的成像结果I中,定位人体目标位置Lx;
步骤6、确定人体目标位置后,将FMCW雷达移动到人体目标方位向位置Lx上,采集雷达回波信号,并对回波信号进行预处理获取人体目标呼吸波形Rex(τa);
步骤7、将步骤6中待识别的人体目标呼吸波形Rex(τa)输入至步骤3获得的身份识别模型,即可对待识别人体目标的身份进行识别。
2.根据权利要求1所述的基于生命体征SAR成像的人体目标身份识别方法,其特征在于,步骤1所述将FMCW雷达静止置于待识别人员的正前方L米处,采集雷达回波信号,并对回波信号进行预处理获取呼吸波形Re(τa),具体为:
步骤1-1、将FMCW雷达静止置于待识别人员的正前方L米处,采集雷达回波信号s(t);
步骤1-2、对回波信号s(t)进行预处理,具体为:
回波信号s(t)为:
式中,sT为雷达的发射信号, 为雷达的接收信号的复共轭,A1为信号sT的幅度,A2为信号 的幅度,K表示调频斜率,λ表示波长,R表示目标与雷达的距离,t表示波传播的“快时间”,c表示波传播的速度;
对回波信号s(t)进行傅里叶变换:
式中,A1A2为信号的幅度,T表示一个波形周期;f表示距离向频率
步骤1-3、根据傅里叶变换后的回波信号构建雷达回波矩阵,具体为:根据傅里叶变换后的回波信号构建Mp×Np雷达回波矩阵,该回波矩阵的行为距离向,每一行存储一个脉冲重复周期内得到的回波信号,列为方位向;其中Mp表示Mp个脉冲重复周期,Np表示一个脉冲重复周期内有Np个采样点Np;
步骤1-4、提取雷达回波矩阵中距离雷达L米处距离内的回波片段SH1,对该回波片段SH1进行相位解缠绕,获取该目标与雷达的距离R1:
式中,ang为相位α解缠绕后的相位角,
其中,相位角 a为回波片段SH1中复数元素的实部,b为回波片段SH1中复
数元素的虚部
步骤1-5、求取人体目标的呼吸波形Re(τa),具体为:将目标与雷达的距离R1对应到每一个方位向采样的慢时间τa获得R1(τa),R1(τa)包含雷达与目标的正对距离R0和人体目标的呼吸波形Re(τa),因此:
Re(τa)=R1(τa)-R0。
3.根据权利要求1所述的基于生命体征SAR成像的人体目标身份识别方法,其特征在于,步骤2所述提取步骤1得到的每一个呼吸波形Re(τa)的特征参数,具体为:
步骤2-1、求取呼吸频率Bf,具体为:对呼吸波形Re(τa)进行FFT变换,即可获得呼吸频率Bf;
步骤2-2、求取平均呼气开始时间周期Tex,具体为:
步骤2-2-1、获取呼吸波形Re(τa)的位移峰值集合:
Pex=[px1 px2 px3 px4 … … pxN]
式中,pxi为第i个呼吸波形的峰值,N表示有N个呼吸波形;
步骤2-2-2、根据位移峰值对应的时间索引求取呼气开始时间周期Texi为:
Texi=tpx(i+1)-tpx(i)
式中,tpxi为第i个呼吸波形的峰值pxi对应的时间索引;
步骤2-2-3、求取平均呼气开始时间周期Tex为:
步骤2-3、求取平均呼吸深度depth,具体为:
步骤2-3-1、获取呼吸波形Re(τa)的最低点集合:
Nin=[nx1 nx2 nx3 nx4 … nxN]
式中,nxi为第i个呼吸波形的最低点,N表示有N个呼吸波形;
步骤2-3-2、求取呼吸深度dxi为:
dxi=pxi-nxi
步骤2-3-3、求取平均呼吸深度depth为:
步骤2-4、求取平均呼气开始时间周期Tex的标准差,具体为:
步骤2-5、求取平均吸气开始时间周期Tin,具体为:
步骤2-5-1、根据呼吸波形最低点对应的时间索引求取吸气开始时间周期Tini为:
Tini=tnx(i+1)-tnxi
式中,tnxi为第i个呼吸波形的最低点nxi对应的时间索引;
步骤2-5-2、求取平均吸气开始时间周期Tin为:
步骤2-6、求取平均吸气开始时间周期Tin的标准差 具体为:
步骤2-7、求取平均呼气速度vexi,具体为:
步骤2-8、求取平均吸气速度vini,具体为:
步骤2-9、求取平均面积比率,具体为:
步骤2-9-1、提取呼吸段中位移值为峰值q%的点为:
p0.qi=0.q*[px1 px2 px3 px4 … … pxN]
提取呼吸段中位移值比最低点高s%的点:
n1.si=1.s*[nx1 nx2 nx3 nx4 … nxN]
步骤2-9-2、根据步骤2-9-1获得呼气-吸气面积A_exi的四个边界点为:[p0.qi n1.si p0.qi n1.si+1];吸气-呼气面积A_ini的四个边界点为:[p0.qi n1.si+1 p0.qi+1 n1.si+1];
步骤2-9-3、计算面积比,并求均值r1为:
步骤2-10、求取每一个呼吸波形峰值前后1秒的欧几里得距离比的均值r2为:
式中,da表示顶点出现一秒之后的数据点到顶点的欧几里德距离,dp表示顶点出现一秒之前的数据点到顶点的欧几里德距离。
4.根据权利要求3所述的基于生命体征SAR成像的人体目标身份识别方法,其特征在于,步骤2-9-1所述q%=70%;s%=30%。
5.根据权利要求2或3所述的基于生命体征SAR成像的人体目标身份识别方法,其特征在于,步骤3所述根据步骤2中提取的呼吸波形的特征参数建立身份识别模型,具体为:将呼吸波形Re(τa)作为样本,将每个样本的特征及其对应的身份标签输入至分类器进行训练,获得身份识别模型。
6.根据权利要求1或2所述的基于生命体征SAR成像的人体目标身份识别方法,其特征在于,步骤4所述FMCW雷达移动平台在所述场景内按直线轨迹运行采集回波信号,并对回波信号预处理后进行SAR成像,获得成像结果I,具体为:
步骤4-1、FMCW雷达沿直线轨迹运行采集原始回波信号s'(t),运行的方向为方位向,FMCW雷达照射方向垂直于运行轨迹,照射方向为距离向;
步骤4-2、对原始回波信号s'(t)进行傅里叶变换后构建雷达回波矩阵,过程与步骤1-
2、1-3相同;
步骤4-3、利用后向投影BP算法对步骤4-2的雷达回波矩阵进行SAR成像,获得成像结果I。
7.根据权利要求6所述的基于生命体征SAR成像的人体目标身份识别方法,其特征在于,步骤5所述在步骤4的成像结果I中,定位人体目标位置Lx,具体为:
步骤5-1、在成像图I中,获取所有局部能量最大的点即所有目标;
步骤5-2、定位各个目标的方位向位置Ln,具体为:针对每一个目标,将其对应的局部能量最大点的行坐标值与方位向分辨率ρa的乘积,作为该目标的方位向位置Ln;其中方位向分辨率ρa为:
式中,La为天线的真实孔径;
步骤5-3、定位各个目标的距离向位置Hn即目标所在距离门Hn,具体为:针对每一个目标,将其对应的局部能量最大点的列坐标值与距离向分辨率ρr的乘积,作为该目标的距离向位置Hn;其中,距离向分辨率ρr为:
式中,c为光速,Bp为FMCW雷达信号带宽;
步骤5-4、提取目标所在距离门Hn内预处理后的回波片段SHn,对该回波片段SHn进行相位解缠绕,获取该目标与雷达的距离R:
式中,ang为相位角α解缠绕后的相位角,
其中,相位角 a为回波片段SHn中复数元素的实部,b为回波片段SHn中复
数元素的虚部;
步骤5-5、对步骤5-4中目标与雷达的距离R进行去拟合处理,获取目标的微多普勒信息,具体为:
步骤5-5-1、将目标与雷达的距离R对应到每一个方位向采样的慢时间τa获得R(τa),R(τa)包含由雷达运动产生的距离Rx(τa)和目标的微多普勒信息Re'(τa);其中:
式中,R0表示雷达在目标正前方时,两者的距离向长度,τa表示方位向运动的慢时间,v表示雷达运动的速度;
步骤5-5-2、根据目标与雷达的距离R对应的所有R(τa),利用偏差平方和最小的原则获得线性拟合曲线函数fit(τa),并由fit(τa)代替Rx(τa);
步骤5-5-3、根据R(τa)和fit(τa)求取目标的微多普勒信息Re'(τa)为:
Re'(τa)=R(τa)-fit(τa)
步骤5-6、对目标的微多普勒信息Re'(τa)进行平滑处理,所用平滑公式为:
式中,Res(τa)m为平滑后的目标微多普勒波形,m为目标的微多普勒信息Re'(τa)中方位向慢时间的第m个采样点,m∈[1,Mn];Mn为采样点数;Re'(τa)m为目标微多普勒信息Re'(τa)的第m个采样点;Re'(τa)m+k为目标微多普勒信息Re'(τa)的第m+k个采样点;K为每个采样点的加权点数,可以通过设定K的值来设置平滑的点数;
步骤5-7、假设雷达运动到目标所在方位向位置Ln的时刻为ta,根据ta时刻Res(τa)m波形判断定位人体目标:若在ta时刻附近Res(τa)m波形中存在有规律的波形,则表明此目标为人类目标,并记录该目标对应的方位向位置Lx;若在ta时刻附近Res(τa)m波形与其他时刻的波形没有太大差异,则表明此目标是非生命体目标。
8.根据权利要求7所述的基于生命体征SAR成像的人体目标身份识别方法,其特征在于,步骤6所述确定人体目标位置后,将FMCW雷达移动到人体目标方位向位置Lx上,采集雷达回波信号,对回波信号进行预处理获取呼吸波形Rex(τa),具体为:
步骤6-1、将FMCW雷达移动到目标的方位向位置Lx上,采集回波数据;
步骤6-2、对步骤6-1采集的回波数据进行傅里叶变换后,构建雷达回波矩阵,具体过程与步骤1-3相同;
步骤6-3、定位人体目标的距离向位置Hx即人体目标所在距离门Hx,具体为:将步骤6-2中雷达回波矩阵中能量最大点的列坐标值与距离向分辨率ρr的乘积,作为该人体目标的距离向位置Hx;
步骤6-4、提取人体目标所在距离门Hx内预处理后的回波片段SHx,对该回波片段SHx进行相位解缠绕,获取人体目标与雷达的距离Rx:
步骤6-5、求取人体目标的呼吸波形Rex(τa),具体为:将目标与雷达的距离Rx对应到每一个方位向采样的慢时间τa获得Rx(τa),Rx(τa)包含雷达与目标的正对距离R0和人体目标的呼吸波形Rex(τa),则:
Rex(τa)=Rx(τa)-R0。
9.一种实现权利要求1至8任意一项所述的基于生命体征SAR成像的人体目标身份识别方法的系统,其特征在于,包括:
呼吸波形提取模,用于采集所有待识别人员的目标呼吸波形;
特征提取模块,用于提取人体目标呼吸波形的特征参数;
分类模型建立模块,用于根据呼吸波形的特征参数建立身份识别模型;
成像模块,用于对场景进行SAR成像;
人体目标定位模块,用于获取SAR成像结果中人体目标的位置;
身份识别模块,用于将场景中采集到的人体目标呼吸波形输入至身份识别模型,获取该人体目标的身份。

说明书全文

基于生命体征SAR成像的人体目标身份识别方法及系统

技术领域

[0001] 本发明属于雷达技术领域,特别涉及一种基于生命体征SAR成像的人体目标身份识别方法及系统。

背景技术

[0002] 移动雷达平台通过天线的方位移动,能够获得更大的合成孔径和更好的空间分辨率,同时雷达平台使用的天线比传统波束扫描雷达所使用的天线要小很多,便于部署和装载。如今调频连续波雷达技术的发展日益成熟,现代调频连续波雷达具有轻量化、低成本和高分辨率等优势,其在地球科学、安防和生命体征监测等领域的应用中具有极大的潜
[0003] 随着雷达的不断发展,对场景中人体目标的定位及身份识别的需求日益增加。因为其在家庭安全、用户认证和健康监测等方面有很大的应用前景。现有的FMCW雷达系统,受限于精度,无法直接从成像中得到人体目标的具体位置,需要结合CW雷达系统,测量目标的微多普勒信息,才能得到生命体目标位置,这种方法增加了系统的复杂度以及成本;目前,根据呼吸模式来确定人体目标身份这类研究尚不成熟,主要存在识别准确率低,数据集不足等问题。因此,快速准确的获取人体目标的位置并进一步识别人体目标身份有着很大的应用前景。

发明内容

[0004] 本发明的目的在于提供一种基于生命体征SAR成像的人体目标身份识别方法及系统。
[0005] 实现本发明目的的技术解决方案为:一种基于生命体征SAR成像的人体目标身份识别方法,包括以下步骤:
[0006] 步骤1、将FMCW雷达静止置于待识别人员的正前方L米处,采集雷达回波信号,并对回波信号进行预处理获取呼吸波形Re(τa);
[0007] 步骤2、提取步骤1得到的每一个呼吸波形Re(τa)的特征参数;
[0008] 步骤3、根据步骤2中提取的呼吸波形的特征参数建立身份识别模型;
[0009] 步骤4、待识别人员随机分布在某一场景中,FMCW雷达移动平台在所述场景内按直线轨迹运行采集回波信号,并对回波信号预处理后进行SAR成像,获得成像结果I;
[0010] 步骤5、在步骤4的成像结果I中,定位人体目标位置Lx;
[0011] 步骤6、确定人体目标位置后,将FMCW雷达移动到人体目标方位向位置Lx上,采集雷达回波信号,并对回波信号进行预处理获取人体目标呼吸波形Rex(τa);
[0012] 步骤7、将步骤6中待识别的人体目标呼吸波形Rex(τa)输入至步骤3获得的身份识别模型,即可对待识别人体目标的身份进行识别。
[0013] 实现上述基于生命体征SAR成像的人体目标身份识别方法的系统,包括:
[0014] 呼吸波形提取模,用于采集所有待识别人员的目标呼吸波形;
[0015] 特征提取模块,用于提取人体目标呼吸波形的特征参数;
[0016] 分类模型建立模块,用于根据呼吸波形的特征参数建立身份识别模型;
[0017] 成像模块,用于对场景进行SAR成像;
[0018] 人体目标定位模块,用于获取SAR成像结果中人体目标的位置;
[0019] 身份识别模块,用于将场景中采集到的人体目标呼吸波形输入至身份识别模型,获取该人体目标的身份。
[0020] 本发明与现有技术相比,其显著优点为:1)本发明通过FMCW移动平台,对场景区域进行生命体征SAR成像,从而确定人体目标位置,不需要事先知道人体目标的位置,克服局限性,使适用场景更加广泛;2)本发明利用雷达传感器获取人体目标的呼吸信号,与传统的接触式监测相比,设备简单易操作,能减少人体的不适感,从而不在实验方面影响实验者的呼吸,减小测量误差,并且可以克服很多的局限性,如被要求接触处受到烧伤或其他因素导致接触片无法直接接触等;3)本发明选择了合适的特征参数来表征不同的个人身份,通过机器学习方法训练出可以进行人体目标身份识别的模型,最终的识别结果较好,准确率较高;4)本发明的方法简单有效,设备简单易实现,低成本,易操作,性能可靠。
[0021] 下面结合附图对本发明做进一步描述。

附图说明

[0022] 图1为本发明基于生命体征SAR成像的人体目标身份识别方法流程图
[0023] 图2为本发明实施例中呼吸波形的动态分割面积比特征示意图。
[0024] 图3为本发明实施例中部吸满时呼吸特征示意图。
[0025] 图4为本发明实施例中SVM混淆矩阵示意图。
[0026] 图5为本发明实施例中实验场景设计图。
[0027] 图6为本发明实施例中BP算法成像图。
[0028] 图7为本发明实施例中解缠绕后的相位信号示意图。
[0029] 图8为本发明实施例中10阶线性拟合后的解缠绕相位信号示意图。
[0030] 图9为本发明实施例中原始呼吸信号示意图。
[0031] 图10为本发明实施例中平滑后的呼吸信号示意图。
[0032] 图11为本发明实施例中最终定位及识别结果图。

具体实施方式

[0033] 结合图1,本发明基于生命体征SAR成像的人体目标身份识别方法,包括以下步骤:
[0034] 步骤1、将FMCW雷达静止置于待识别人员的正前方L米处,采集雷达回波信号,并对回波信号进行预处理获取呼吸波形Re(τa);
[0035] 步骤2、提取步骤1得到的每一个呼吸波形Re(τa)的特征参数;
[0036] 步骤3、根据步骤2中提取的呼吸波形的特征参数建立身份识别模型;
[0037] 步骤4、待识别人员随机分布在某一场景中,FMCW雷达移动平台在所述场景内按直线轨迹运行采集回波信号,并对回波信号预处理后进行SAR成像,获得成像结果I;
[0038] 步骤5、在步骤4的成像结果I中,定位人体目标位置Lx;
[0039] 步骤6、确定人体目标位置后,将FMCW雷达移动到人体目标方位向位置Lx上,采集雷达回波信号,并对回波信号进行预处理获取人体目标呼吸波形Rex(τa);
[0040] 步骤7、将步骤6中待识别的人体目标呼吸波形Rex(τa)输入至步骤3获得的身份识别模型,即可对待识别人体目标的身份进行识别。
[0041] 进一步地,步骤1将FMCW雷达静止置于待识别人员的正前方L米处,采集雷达回波信号,并对回波信号进行预处理获取呼吸波形Re(τa),具体为:
[0042] 步骤1-1、将FMCW雷达静止置于待识别人员的正前方L米处,采集雷达回波信号s(t);
[0043] 步骤1-2、对回波信号s(t)进行预处理,具体为:
[0044] 回波信号s(t)为:
[0045]
[0046] 式中,sT为雷达的发射信号, 为雷达的接收信号的复共轭,A1为信号sT的幅度,A2为信号 的幅度,K表示调频斜率,λ表示波长,R表示目标与雷达的距离,t表示波传播的“快时间”,c表示波传播的速度;
[0047] 对回波信号s(t)进行傅里叶变换:
[0048]
[0049] 式中,A1A2为信号的幅度,T表示一个波形周期;f表示距离向频率
[0050] 步骤1-3、根据傅里叶变换后的回波信号构建雷达回波矩阵,具体为:根据傅里叶变换后的回波信号构建Mp×Np雷达回波矩阵,该回波矩阵的行为距离向,每一行存储一个脉冲重复周期内得到的回波信号,列为方位向;其中Mp表示Mp个脉冲重复周期,Np表示一个脉冲重复周期内有Np个采样点Np;
[0051] 步骤1-4、提取雷达回波矩阵中距离雷达L米处距离内的回波片段SH1,对该回波片段SH1进行相位解缠绕,获取该目标与雷达的距离R1:
[0052]
[0053] 式中,ang为相位α解缠绕后的相位角,
[0054]
[0055] 其中,相位角 a为回波片段SH1中复数元素的实部,b为回波片段SH1中复数元素的虚部
[0056] 步骤1-5、求取人体目标的呼吸波形Re(τa),具体为:将目标与雷达的距离R1对应到每一个方位向采样的慢时间τa获得R1(τa),R1(τa)包含雷达与目标的正对距离R0和人体目标的呼吸波形Re(τa),因此:
[0057] Re(τa)=R1(τa)-R0。
[0058] 进一步地,步骤2提取步骤1得到的每一个呼吸波形Re(τa)的特征参数,具体为:
[0059] 步骤2-1、求取呼吸频率Bf,具体为:对呼吸波形Re(τa)进行FFT变换,即可获得呼吸频率Bf;
[0060] 步骤2-2、求取平均呼气开始时间周期Tex,具体为:
[0061] 步骤2-2-1、获取呼吸波形Re(τa)的位移峰值集合:
[0062] Pex=[px1 px2 px3 px4 … … pxN]
[0063] 式中,pxi为第i个呼吸波形的峰值,N表示有N个呼吸波形;
[0064] 步骤2-2-2、根据位移峰值对应的时间索引求取呼气开始时间周期Texi为:
[0065] Texi=tpx(i+1)-tpx(i)
[0066] 式中,tpxi为第i个呼吸波形的峰值pxi对应的时间索引;
[0067] 步骤2-2-3、求取平均呼气开始时间周期Tex为:
[0068]
[0069] 步骤2-3、求取平均呼吸深度depth,具体为:
[0070] 步骤2-3-1、获取呼吸波形Re(τa)的最低点集合:
[0071] Nin=[nx1 nx2 nx3 nx4 ... nxN]
[0072] 式中,nxi为第i个呼吸波形的最低点,N表示有N个呼吸波形;
[0073] 步骤2-3-2、求取呼吸深度dxi为:
[0074] dxi=pxi-nxi
[0075] 步骤2-3-3、求取平均呼吸深度depth为:
[0076]
[0077] 步骤2-4、求取平均呼气开始时间周期Tex的标准差,具体为:
[0078]
[0079] 步骤2-5、求取平均吸气开始时间周期Tin,具体为:
[0080] 步骤2-5-1、根据呼吸波形最低点对应的时间索引求取吸气开始时间周期Tini为:
[0081] Tini=tnx(i+1)-tnxi
[0082] 式中,tnxi为第i个呼吸波形的最低点nxi对应的时间索引;
[0083] 步骤2-5-2、求取平均吸气开始时间周期Tin为:
[0084]
[0085] 步骤2-6、求取平均吸气开始时间周期Tin的标准差 具体为:
[0086]
[0087] 步骤2-7、求取平均呼气速度vexi,具体为:
[0088]
[0089] 步骤2-8、求取平均吸气速度vini,具体为:
[0090]
[0091] 步骤2-9、求取平均面积比率,具体为:
[0092] 步骤2-9-1、提取呼吸段中位移值为峰值q%的点为:
[0093] p0.qi=0.q*[px1 px2 px3 px4 ... ... pxN]
[0094] 提取呼吸段中位移值比最低点高s%的点:
[0095] n1.si=1.s*[nx1 nx2 nx3 nx4 … nxN]
[0096] 步骤2-9-2、根据步骤2-9-1获得呼气-吸气面积A-exi的四个边界点为:[p0.qi n1.si p0.qi n1.si+1];吸气-呼气面积A-ini的四个边界点为:[p0.qi n1.si+1 p0.qi+1 n1.si+1];
[0097] 步骤2-9-3、计算面积比,并求均值r1为:
[0098]
[0099] 步骤2-10、求取每一个呼吸波形峰值前后1秒的欧几里得距离比的均值r2为:
[0100]
[0101] 式中,da表示顶点出现一秒之后的数据点到顶点的欧几里德距离,dp表示顶点出现一秒之前的数据点到顶点的欧几里德距离。
[0102] 示例性优选地,步骤2-9-1中q%=70%;s%=30%。
[0103] 进一步地,步骤3根据步骤2中提取的呼吸波形的特征参数建立身份识别模型,具体为:将呼吸波形Re(τa)作为样本,将每个样本的特征及其对应的身份标签输入至分类器进行训练,获得身份识别模型。
[0104] 进一步地,步骤4中FMCW雷达移动平台在所述场景内按直线轨迹运行采集回波信号,并对回波信号预处理后进行SAR成像,获得成像结果I,具体为:
[0105] 步骤4-1、FMCW雷达沿直线轨迹运行采集原始回波信号s'(t),运行的方向为方位向,FMCW雷达照射方向垂直于运行轨迹,照射方向为距离向;
[0106] 步骤4-2、对原始回波信号s'(t)进行傅里叶变换后构建雷达回波矩阵,过程与步骤1-2、1-3相同;
[0107] 步骤4-3、利用后向投影BP算法对步骤4-2的雷达回波矩阵进行SAR成像,获得成像结果I。
[0108] 进一步地,步骤5在步骤4的成像结果I中,定位人体目标位置Lx,具体为:
[0109] 步骤5-1、在成像图I中,获取所有局部能量最大的点即所有目标;
[0110] 步骤5-2、定位各个目标的方位向位置Ln,具体为:针对每一个目标,将其对应的局部能量最大点的行坐标值与方位向分辨率ρa的乘积,作为该目标的方位向位置Ln;其中方位向分辨率ρa为:
[0111]
[0112] 式中,La为天线的真实孔径;
[0113] 步骤5-3、定位各个目标的距离向位置Hn即目标所在距离门Hn,具体为:针对每一个目标,将其对应的局部能量最大点的列坐标值与距离向分辨率ρr的乘积,作为该目标的距离向位置Hn;其中,距离向分辨率ρr为:
[0114]
[0115] 式中,c为光速,Bp为FMCW雷达信号带宽;
[0116] 步骤5-4、提取目标所在距离门Hn内预处理后的回波片段SHn,对该回波片段SHn进行相位解缠绕,获取该目标与雷达的距离R:
[0117]
[0118] 式中,ang为相位角α解缠绕后的相位角,
[0119]
[0120] 其中,相位角 a为回波片段SHn中复数元素的实部,b为回波片段SHn中复数元素的虚部;
[0121] 步骤5-5、对步骤5-4中目标与雷达的距离R进行去拟合处理,获取目标的微多普勒信息,具体为:
[0122] 步骤5-5-1、将目标与雷达的距离R对应到每一个方位向采样的慢时间τa获得R(τa),R(τa)包含由雷达运动产生的距离Rx(τa)和目标的微多普勒信息Re'(τa);其中:
[0123]
[0124] 式中,R0表示雷达在目标正前方时,两者的距离向长度,τa表示方位向运动的慢时间,v表示雷达运动的速度;
[0125] 步骤5-5-2、根据目标与雷达的距离R对应的所有R(τa),利用偏差平方和最小的原则获得线性拟合曲线函数fit(τa),并由fit(τa)代替Rx(τa);
[0126] 步骤5-5-3、根据R(τa)和fit(τa)求取目标的微多普勒信息Re'(τa)为:
[0127] Re'(τa)=R(τa)-fit(τa)
[0128] 步骤5-6、对目标的微多普勒信息Re'(τa)进行平滑处理,所用平滑公式为:
[0129]
[0130] 式中,Res(τa)m为平滑后的目标微多普勒波形,m为目标的微多普勒信息Re'(τa)中方位向慢时间的第m个采样点,m∈[1,Mn];Mn为采样点数;Re'(τa)m为目标微多普勒信息Re'(τa)的第m个采样点;Re'(τa)m+k为目标微多普勒信息Re'(τa)的第m+k个采样点;K为每个采样点的加权点数,可以通过设定K的值来设置平滑的点数;
[0131] 步骤5-7、假设雷达运动到目标所在方位向位置Ln的时刻为ta,根据ta时刻Res(τa)m波形判断定位人体目标:若在ta时刻附近Res(τa)m波形中存在有规律的波形,则表明此目标为人类目标,并记录该目标对应的方位向位置Lx;若在ta时刻附近Res(τa)m波形与其他时刻的波形没有太大差异,则表明此目标是非生命体目标。
[0132] 进一步地,步骤6确定人体目标位置后,将FMCW雷达移动到人体目标方位向位置Lx上,采集雷达回波信号,对回波信号进行预处理获取呼吸波形Rex(τa),具体为:
[0133] 步骤6-1、将FMCW雷达移动到目标的方位向位置Lx上,采集回波数据;
[0134] 步骤6-2、对步骤6-1采集的回波数据进行傅里叶变换后,构建雷达回波矩阵,具体过程与步骤1-3相同;
[0135] 步骤6-3、定位人体目标的距离向位置Hx即人体目标所在距离门Hx,具体为:将步骤6-2中雷达回波矩阵中能量最大点的列坐标值与距离向分辨率ρr的乘积,作为该人体目标的距离向位置Hx;
[0136] 步骤6-4、提取人体目标所在距离门Hx内预处理后的回波片段SHx,对该回波片段SHx进行相位解缠绕,获取人体目标与雷达的距离Rx:
[0137]
[0138] 步骤6-5、求取人体目标的呼吸波形Rex(τa),具体为:将目标与雷达的距离Rx对应到每一个方位向采样的慢时间τa获得Rx(τa),Rx(τa)包含雷达与目标的正对距离R0和人体目标的呼吸波形Rex(τa),则:
[0139] Rex(τa)=Rx(τa)-R0。
[0140] 实现上述基于生命体征SAR成像的人体目标身份识别方法的系统,包括:
[0141] 呼吸波形提取模块,用于采集所有待识别人员的目标呼吸波形;
[0142] 特征提取模块,用于提取人体目标呼吸波形的特征参数;
[0143] 分类模型建立模块,用于根据呼吸波形的特征参数建立身份识别模型;
[0144] 成像模块,用于对场景进行SAR成像;
[0145] 人体目标定位模块,用于获取SAR成像结果中人体目标的位置;
[0146] 身份识别模块,用于将场景中采集到的人体目标呼吸波形输入至身份识别模型,获取该人体目标的身份。
[0147] 下面结合实施例对本发明作进一步详细的描述。
[0148] 实施例
[0149] 结合图1,本发明基于生命体征SAR成像的人体目标身份识别方法,包括以下内容:
[0150] 1、在本实施例中采用一个基于移动平台的FMCW雷达系统,系统雷达的载频为5.8GHz,发射信号的带宽为320MHz,天线的增益为11.3dB,半功率角为46°,采样频率
192kHz。
[0151] 2、将FMCW雷达静止放置在实验人员正前方1米处,采集10个实验人员的呼吸波形数据,每包呼吸数据为30秒,每个实验人员采集50包,共500包数据。
[0152] 3、提取呼吸波形的特征参数,其中呼吸波形的动态分割面积比特征如图2所示,肺部吸满时前后一秒的呼吸特征如图3所示。
[0153] 4、将10个实验人员分别贴上标签1-10,然后将提取的特征参数和对应的标签送入分类器进行训练,获得身份识别模型,分类结果混淆矩阵如图4所示。
[0154] 5、设置场景为方位向长3.6m,距离向长5m的矩形区域,在场景区域内共放置了两个目标,其中目标L1为反射体,目标L2为生命体目标,其身份是标签为1的实验人员,目标距离场景左边延分别为1.2m和2.4m,距离场景下边沿的距离都是1m,具体的实验场景设计如图5所示。
[0155] 6、让雷达延3.6m长的场景下边沿(方位向)由左至右以速v=0.2m/s进行运动,对回波数据进行预处理后使用BP算法对回波矩阵进行SAR成像,成像结果I如图6所示,图中显示了两个目标L1和L2。
[0156] 7、首先由于两个目标的距离向距离相同,故提取任一目标所在距离门H1的回波片段SH1,处理后得雷达与目标的距离R(τa)如图7所示;其次,为了消除由雷达运动产生的距离信息Rx(τa),计算一个10阶的拟合函数fit(τa),结果如图8所示;最后,用R(τa)减去拟合函数fit(τa)即可得到目标的微多普勒信息Re'(τa),如图9所示;对Re'(τa)做250点的平滑处理,得到平滑后的目标微多普勒波形Res(τa)如图10所示;从图10可以看出,虚线框中的12到27秒,可以清楚地看到规则的振动信号,振幅比前12秒大10倍,而12到27秒的时间,正是雷达运动到目标L2附近的时间,也就是说,目标L2是图6中的人体目标。
[0157] 8、将FMCW移动到人体目标L2的位置,静止的测量目标的呼吸波形。
[0158] 9、将测量得到的人体目标L2的呼吸波形输入到身份识别模型,即可获得到人体目标L2的身份。本实施例最后的识别结果如图11所示,人体目标L2是标签为1的实验人员。
[0159] 本发明通过雷达传感器对生理信号进行采集,提取人体目标呼吸波形的特征参数,训练出身份识别模型,接着利用基于相位解调的生命体征SAR成像方法对人体目标进行定位,最后提取人体目标呼吸波形的特征参数送入身份识别模型中,识别出人体目标的身份。本发明的人体目标定位效果好,身份识别准确率高,鲁棒性好,且适用性更广。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈