一种融合毫米波雷达和深度学习模型的无接触式心震信号检测方法及其装置与流程

1.本发明属于数字信号处理器应用技术领域，具体涉及一种融合毫米波雷达和深度学习模型的无接触式心震信号检测方法及其装置。


背景技术：

2.心震信号属于人体生物信号，是由人体发出的不稳定的低频微弱信号。可用心震图（seismocardiography，scg）表示心脏搏动对胸壁造成冲击而产生的局部震动。医学文献表明，scg可用于精确记录精准时间颗粒度下的心脏活动，包括心脏瓣膜的打开和关闭。这些测量值可用于检测和诊断多种心血管疾病，例如心肌梗塞（心脏病），冠心病，局部缺血和出血。
3.测量scg信号的标准方法依赖于可捕获胸壁微振动的加速度计。该过程通常要求用户脱下衬衫，仰卧，并使用胸带将加速度计固定在心脏的顶点附近，通常需要由医生在经过校准的医疗环境或受控环境中进行管理，这个过程十分麻烦费力，操作繁琐。
4.射频（radio frequency, rf）信号不仅可以实现无接触检测，而且rf信号内记录了大量有意义的特征，这意味着rf信号在感知监测人类生命体征领域有重大意义，但是传统的方法不能捕获隐藏在rf信号内的心脏运动细节。
5.有鉴于此，提出使用先进的数字信号处理手段，提取rf射频信号携带的有用信息得到心脏运动的scg信号，并实现心脏运动分类准则下的分类标注。


技术实现要素：

6.针对上述技术问题，本发明提供一种融合毫米波雷达和深度学习模型的无接触式心震信号检测方法及其装置，解决检测心震信号条件束缚多、成本高等问题；解决rf信号无法进行可视化查看等缺点，解决心脏运动难以检测分类的问题，为预防心脏疾病提供充足的判断数据和材料。
7.本发明的目的可以通过采取如图1所示的技术方案达到：步骤1：用多普勒毫米波雷达信号的发射与接收模块，获取雷达信号；步骤2：利用所述雷达信号，对心脏雷达信号进行定位与检测，获得心脏波束；步骤3：利用心震信号转换模块，由心脏波束信号获得与其对应的心震图信号；步骤4：由心震图信号分类标注模块，划分心脏在时序上的微运动状态，并判断心脏活动是否存在异常；步骤5：将心脏活动的分类结果在装置内进行存储，并且产生可视化结果推送至显示装置。
8.进一步地，所述的一种基于深度学习模型的无接触式心震信号检测毫米波雷达装置及其检测方法，其特征在于，所述雷达装置置于被检测者活动空间前方，通过发射天线向外发射固定频率的发射信号，经过所述被检测者会产生所述反射信号，所述装置无需接触
所述检测者。
9.进一步地，所述的心脏波束形成器，其特征在于，心脏波束形成器由如图2所示两个平行的信号处理流组成；步骤2-1：组合多通道采集数据，以十六个通道信号为例，组合公式：其中为发射射频信号的波长，、为第n个天线阵的坐标，x、y为投影坐标，为从第n个天线阵接收到的信号；步骤2-2：组合后的信号传入所述的第一个处理流提取心率、最大心率比；步骤2-3：组合后的信号传入所述的第二个处理流将十六个通道信号特征映射，在雷达检测区域内所有信号中确定心脏的射频信号；步骤2-4：合并两个处理流的输出：所述的波束形成器对三坐标参数化的时域信号进行快速傅里叶变换得到时域谱特性；由步骤2-2确定心率，由此得到四维特征数据（x坐标、y坐标、频率和心率）；步骤2-5：使用步骤2-4的四维输出结果，所述算法可以确定最佳的有最高的功率比的空间位置(即，心跳的周期最明显的地方)。
10.进一步地，所述的心脏波束形成器第一个提取心率的处理流，其特征在于，将原始数据d经过提取相位和微分变换后，送入一维卷积神经网络进行模板匹配；用随机梯度下降法求解所述过程，卷积神经网络训练时设定损失函数为：其中，是l2正则化函数，t是卷积滤波核，m是每个最大池化层输出元素的值，是调节超参数；相邻m的时间间隔为单周期心率；最后将卷积神经网络的最后一层最大池化层输出的数据为心率提取特征，使用频率分布直方图估计最大心率。
11.进一步地，所述的确定心脏射频信号的第二个处理流，其特征在于，由如下步骤组成：步骤2-3-1：使用波束形成器提升雷达精度和角度，并用傅里叶变换将原始数据进行特征空间映射；步骤2-3-2：提取时序数据的相位，并应用微分滤波器。
12.进一步地，所述的如图3所示的心脏波束与心震图信号转换模块，其特征在于，转换器由三模块神经网络组成，每个模块的神经网络由一维卷积层、激活层、批量标准化层和随机丢弃特征层组成；所述的三个模块神经网络具有较多数量的卷积核提取有效信息，卷积核越多越能得到更丰富的特征，此处以128个卷积核、32个卷积核和1个卷积核为例；在训练阶段，输入是rf信号，训练标签是同步采集的真实scg信号，损失函数是神经网络输出与真实scg信号之间的l2正则化函数；在测试阶段，仅需rf信号输入。
13.进一步地，所述的心震图信号分类标注模块，其特征在于，是一种适用于一维scg信号的u-net模型；u-net模型不仅能对输入的scg序列数据进行编码解码，还用层级融合，保留编码过程中的高维特征，以提高模型输出的准确度。
14.进一步地，一种装置，其特征在于，所述装置中存储有计算机程序，其中，所述计算机程序被设置为运行时执行上述任一项描述中的融合毫米波雷达和深度学习模型的无接触式心震信号检测方法。
15.与现有技术相比，本发明所提供的融合毫米波雷达和深度学习模型的无接触式心震信号检测方法及其装置具有以下优点效果：1)本发明提供的融合毫米波雷达和深度学习模型的无接触式心震信号检测方法及其装置能够减少检测心震信号的设施限制。
16.2)本发明提供的融合毫米波雷达和深度学习模型的无接触式心震信号检测方法及其装置能够降低检测心震信号的成本。
17.3)本发明提供的融合毫米波雷达和深度学习模型的无接触式心震信号检测方法及其装置能够解决rf信号无法进行可视化查看和利用率低的缺点。
18.4)本发明提供的融合毫米波雷达和深度学习模型的无接触式心震信号检测方法及其装置能够检测心脏运动并准确分类。
19.5)本发明提供的融合毫米波雷达和深度学习模型的无接触式心震信号检测方法及其装置能够为预防心脏疾病提供充足的判断数据和材料。
附图说明
20.为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
21.图1为本发明的结构示意图；图2为心脏波束形成器的两个信息处理流及其处理方法；图3为心脏波束信号与心震信号转化模块的网络结构图；图4为心震信号分类标注的两种可视化结果。
具体实施方式
22.为了使本发明实施例的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明实施例进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
23.请参阅图1-4，本发明还提供一种融合毫米波雷达和深度学习模型的无接触式心震信号检测方法，其包括步骤：步骤1：用多普勒毫米波雷达信号的发射与接收模块，获取雷达信号；具体的，所述雷达装置置于被检测者活动空间前方，通过发射天线向外发射固定频率的发射信号，经过所述被检测者会产生所述反射信号，所述装置无需接触所述检测者。
24.步骤2：利用所述雷达信号，对心脏雷达信号进行定位与检测，获得心脏波束；具体的，心脏波束形成器由两个平行的信号处理流组成；步骤2-1：组合多通道采集数据，以十六个通道信号为例，组合公式：
其中为发射射频信号的波长，、为第n个天线阵的坐标，x、y为投影坐标，为从第n个天线阵接收到的信号；步骤2-2：组合后的信号传入所述的第一个处理流提取心率、最大心率比；具体的，将原始数据d经过提取相位和微分变换后，送入一维卷积神经网络进行模板匹配；用随机梯度下降法求解所述过程，卷积神经网络训练时设定损失函数为：其中，是l2正则化函数，t是卷积滤波核，m是每个最大池化层输出元素的值，是调节超参数；相邻m的时间间隔为单周期心率；最后将卷积神经网络的最后一层最大池化层输出的数据为心率提取特征，使用频率分布直方图估计最大心率。
25.步骤2-3：组合后的信号传入所述的第二个处理流将十六个通道信号特征映射，在雷达检测区域内所有信号中确定心脏的射频信号；具体的，由如下步骤组成：步骤2-3-1：使用波束形成器提升雷达精度和角度，并用傅里叶变换将原始数据进行特征空间映射；步骤2-3-2：提取时序数据的相位，并应用微分滤波器。
26.步骤2-4：合并两个处理流的输出：所述的波束形成器对三坐标参数化的时域信号进行快速傅里叶变换得到时域谱特性；由步骤2-2确定心率，由此得到四维特征数据（x坐标、y坐标、频率和心率）；步骤2-5：使用步骤2-4的四维输出结果，所述算法可以确定最佳的有最高的功率比的空间位置(即，心跳的周期最明显的地方)。
27.步骤3：利用心震信号转换模块，由心脏波束信号获得与其对应的心震图信号；具体的，转换器由三模块神经网络组成，每个模块的神经网络由一维卷积层、激活层、批量标准化层和随机丢弃特征层组成；所述的三个模块神经网络具有较多数量的卷积核提取有效信息，卷积核越多越能得到更丰富的特征，此处以128个卷积核、32个卷积核和1个卷积核为例；在训练阶段，输入是rf信号，训练标签是同步采集的真实scg信号，损失函数是神经网络输出与真实scg信号之间的l2正则化函数；在测试阶段，仅需rf信号输入。
28.步骤4：由心震图信号分类标注模块，划分心脏在时序上的微运动状态，并判断心脏活动是否存在异常；具体的，是一种适用于一维scg信号的u-net模型；u-net模型不仅能对输入的scg序列数据进行编码解码，还用层级融合，保留编码过程中的高维特征，以提高模型输出的准确度。
29.步骤5：将心脏活动的分类结果在装置内进行存储，并且产生如图4所示的可视化结果推送至显示装置。
30.本说明书的描述中，参考术语“实施例”、“示例”等的描述意指结合该实施例或示
例描述的具体特征、结构、材料或特点包括于本发明的至少一个实施例或示例中。本发明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何一个或多个实施例或示例中以合适的方式结合。
31.以上内容仅仅是对本发明结构所作的举例和说明，所属本技术领域的技术人员对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，只要不偏离发明的结构或者超越本权利要求书所定义的范围，均应属于本发明的保护范围。